Commit c8ccb16d authored by MIJIEUX Thomas's avatar MIJIEUX Thomas

Add test case in RESULTS.org for influence of restart parameter

parent 4108a3aa
...@@ -143,3 +143,28 @@ ggplot(df, aes(x=nb_mvp)) + ...@@ -143,3 +143,28 @@ ggplot(df, aes(x=nb_mvp)) +
scale_y_log10() + geom_hline(aes(yintercept=1e-6, color="threshold")) + 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)")) ggtitle("young1c (nrhs=6, m=90 k=5)") + ylab(TeX("$\\eta_b$ (min, max)"))
#+end_src #+end_src
** Influence of restart parameter (nrhs=6, m=90)
*** run test case
#+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_core/
mkdir -p ../data/res
#nbRHS=6; maxSpace=90
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A QRDR -r DEFLATED -p 5 -u -o "BGMRES-DR(05)"
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A QRDR -r DEFLATED -p 10 -u -o "BGMRES-DR(10)"
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A QRDR -r DEFLATED -p 15 -u -o "BGMRES-DR(15)"
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -A QRDR -r DEFLATED -p 20 -u -o "BGMRES-DR(20)"
#+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*
library(ggplot2)
library(latex2exp)
df <- read.table("./build/src/data/res/BGMRES-DR(05).res", header=T)
df <- rbind(df, read.table("./build/src/data/res/BGMRES-DR(10).res", header=T))
df <- rbind(df, read.table("./build/src/data/res/BGMRES-DR(15).res", header=T))
df <- rbind(df, read.table("./build/src/data/res/BGMRES-DR(20).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,10,15,20))") + ylab(TeX("$\\eta_b$ (min, max)"))
#+end_src
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