SANDBOX.org 3.47 KB
Newer Older
1 2 3 4 5 6 7
# -*- mode: org -*-
# -*- coding: utf-8 -*-

#+TITLE: FABuLOuS sandbox

* Setup workstation
  Set this variable to your work directory (where you downloaded fabulous):
8
  #+begin_src sh :session *TEST* :results silent
9
export WORKDIR=/home/tmijieux/fabulous # for example
10
  #+end_src
11 12

* SANDBOX
13
** csc lightInTissue timing
14
*** run test
15 16 17 18 19 20 21
    #+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_basic/
mkdir -p ../data/res
#nbRHS=6; maxSpace=90
./testMatrixMarketChoice -f ../data/New/lightINtissue.mtx -k CSC -n 6 -m 90 -A IBDR \
                         -r DEFLATED -p 5 -e 1e-4 -u -o "IB-BGMRES-DR"
    #+end_src
22
*** plot
23 24 25 26 27 28 29 30 31 32 33 34 35 36
    #+begin_src R :results output graphics :file (org-babel-temp-file "figure" ".png") :exports both :width 600 :height 400 :session *R*
library(ggplot2)
library(latex2exp)
df <- read.table("./build/src/data/res/IB-BGMRES-DR.res", header=T)
#df <- rbind(df, read.table("./build/src/data/res/IB-BGMRES.res", header=T))
#df <- rbind(df, read.table("./build/src/data/res/BGMRES-DR.res", header=T))
#df <- rbind(df, read.table("./build/src/data/res/IB-BGMRES-DR.res", header=T))
ggplot(df, aes(x=nb_mvp)) +
    geom_line(aes(y=time, color="total")) +
    geom_line(aes(y=least_square_time, color="gels")) +
    geom_line(aes(y=mvp_spent, color="mvp")) +
    geom_line(aes(y=ortho_spent, color="ortho")) +
    ggtitle("lightInTissue (nrhs=6, m=90 k=5)") + ylab(TeX("$\\eta_b$ (min, max)"))
    #+end_src
37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
** qr ib dr
*** run test
    #+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_basic/
mkdir -p ../data/res
#nbRHS=6; maxSpace=90
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A IBDR -r DEFLATED -p 5 -u -o "IB-BGMRES-DR"
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A QRIBDR -r DEFLATED -p 5 -u -o "QR-IB-BGMRES-DR"
    #+end_src
*** plot
    #+begin_src R :results output graphics :file (org-babel-temp-file "figure" ".png") :exports both :width 600 :height 400 :session *R*
library(ggplot2)
library(latex2exp)
df <- read.table("./build/src/data/res/IB-BGMRES-DR.res", header=T)
df <- rbind(df, read.table("./build/src/data/res/QR-IB-BGMRES-DR.res", header=T))
ggplot(df, aes(x=nb_mvp)) +
    geom_line(aes(y=maxRes, color=name)) +
    geom_line(aes(y=minRes, color=name)) +
    scale_y_log10() + geom_hline(aes(yintercept=1e-6, color="threshold")) +
    ggtitle("young1c (nrhs=6, m=90 k=5)") + ylab(TeX("$\\eta_b$ (min, max)"))
57

58
    #+end_src
59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85

    #+results:
    [[file:/tmp/babel-4496LoC/figure44968FD.png]]
** qr ib dr timing
*** run test
    #+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_basic/
mkdir -p ../data/res
#nbRHS=6; maxSpace=90
./testMatrixMarketChoice -n 100 -t RUHE -s MGS -m 700 -e 1e-4 -A IBDR -r DEFLATED -p 5 -u -o "IB-BGMRES-DR"
./testMatrixMarketChoice -n 100 -t RUHE -s MGS -m 700 -e 1e-4 -A QRIBDR -r DEFLATED -p 5 -u -o "QR-IB-BGMRES-DR"
    #+end_src
*** plot
    #+begin_src R :results output graphics :file (org-babel-temp-file "figure" ".png") :exports both :width 600 :height 400 :session *R*
library(ggplot2)
library(latex2exp)
df <- read.table("./build/src/data/res/IB-BGMRES-DR.res", header=T)
df <- rbind(df, read.table("./build/src/data/res/QR-IB-BGMRES-DR.res", header=T))
ggplot(df, aes(x=nb_mvp)) +
    geom_line(aes(y=least_square_time, color=name)) +
    geom_line(aes(y=least_square_time, color=name)) +
    ggtitle("young1c (nrhs=100, m=500 k=5, max_mvp=10000)") + ylab(TeX("$\\eta_b$ (min, max)"))

    #+end_src

    #+results:
    [[file:/tmp/babel-4496LoC/figure4496wnn.png]]