Commit 502d1121 authored by MIJIEUX Thomas's avatar MIJIEUX Thomas

update results.org and TODO.org

parent cc7783e8
......@@ -85,8 +85,8 @@ ggplot(df, aes(x=nb_mvp, y=maxRes, color=name)) +
#+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_basic/
mkdir -p ../data/res
./testMatrixMarketChoice -t BLOCK -s CGS -m 900 -A STDDR -u -o "Full_GELS"
./testMatrixMarketChoice -t BLOCK -s CGS -m 900 -A QR -u -o "GELS_with_Incremental_QR_factorization"
./testMatrixMarketChoice -t BLOCK -s CGS -m 700 -A STDDR -u -o "Full_GELS"
./testMatrixMarketChoice -t BLOCK -s CGS -m 700 -A QR -u -o "GELS_with_Incremental_QR_factorization"
#+end_src
*** plot the graphic
#+begin_src R :results output graphics :file (org-babel-temp-file "figure" ".png") :exports both :width 600 :height 400 :session *R*
......@@ -96,6 +96,7 @@ df <- rbind(df, read.table("./build/src/data/res/GELS_with_Incremental_QR_factor
ggplot(df, aes(x=global_iteration, y=least_square_time, color=name)) +
geom_line() + ggtitle("Least Square duration (young1c)")
#+end_src
** Influence of the Algorithm
*** run test case
#+begin_src sh :session *TEST* :results silent
......@@ -202,3 +203,68 @@ ggplot(df, aes(x=nb_mvp)) +
scale_y_log10() + geom_hline(aes(yintercept=1e-6, color="threshold")) +
ggtitle("young1c (nrhs=6, m=90 k=(5,10,15,20))") + ylab(TeX("$\\eta_b$ (min, max)"))
#+end_src
** csc lightInTissue 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 -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
*** 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/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
** 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)"))
#+end_src
** 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
......@@ -10,7 +10,7 @@ export WORKDIR=/home/tmijieux/fabulous # for example
#+end_src
* SANDBOX
** csc lightInTissue
** csc lightInTissue timing
*** run test
#+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_basic/
......@@ -56,3 +56,30 @@ ggplot(df, aes(x=nb_mvp)) +
ggtitle("young1c (nrhs=6, m=90 k=5)") + ylab(TeX("$\\eta_b$ (min, max)"))
#+end_src
#+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]]
......@@ -21,16 +21,16 @@
*** DONE Add timer for orthogonalization, least square and matrix vector product
*** DONE Add global timer
*** TODO Improve timer semantics and logs
*** TODO Fix problem of least square time measure in QR+DR version
*** TODO Fix problem of least square time measure in QR+DR version; IBQRDR
(actual factorization if performed in notify_factorization_end() )
*** TODO Add flops/s counter
*** TODO print flops
** DONE Improve RESULTS.org
Eventually, anyone must be able to gather all interesting results into RESULTS.org
just by evaluating code blocks from RESULTS.org and/or tangling RESULTS.org
* DONE Implement IB+DR
* TODO Implement IB+DR+QR
** tpqrt and tpmqrt kernels
* TODO iterated orthogonalization stop criterion
* TODO find out how to link fabulous with parallel mkl with spack
* TODO parallel(distributed) test case
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment