Commit 9236607e authored by MIJIEUX Thomas's avatar MIJIEUX Thomas

Add test case for young1c (paper figure2 example8)

parent cbdc117a
......@@ -83,14 +83,13 @@ cd ${WORKDIR}/build/src/test_core/
./testMatrixMarketChoice -t BLOCK -s CGS -m 900 -A STDDR -u -o "Basic_GELS"
./testMatrixMarketChoice -t BLOCK -s CGS -m 900 -A QR -u -o "QR_factorization"
#+end_src
**** plot the graphic
*** 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)
df <- read.table("./build/src/data/res/Basic_GELS.res", header=T)
df <- rbind(df, read.table("./build/src/data/res/QR_factorization.res", header=T))
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
......@@ -116,3 +115,26 @@ ggplot(df, aes(x=nb_mvp)) +
scale_y_log10() +
ggtitle("coneSphere problem with CGS RUHE (STD, STD+DR and IB algorithm)")
#+end_src
** young1c (nrhs=6, m=90, k=5)
*** run test case
#+begin_src sh :session *TEST* :results silent
cd ${WORKDIR}/build/src/test_core/
#nbRHS=6; maxSpace=90
./testMatrixMarketChoice -n 6 -t BLOCK -s MGS -m 90 -e 1e-6 -u -o "BGMRES"
./testMatrixMarketChoice -n 6 -t BLOCK -s MGS -m 90 -e 1e-6 -A IB -u -o "IB-BGMRES"
./testMatrixMarketChoice -n 6 -t RUHE -s MGS -m 90 -e 1e-6 -r DEFLATED -p 5 -u -o "BGMRES-DR"
#+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.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))
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
......@@ -11,7 +11,7 @@
* HOLD Reproducible installation
** DONE Instruction in INSTALL.org
** DONE Update spack package
** HOLD merge request to make Chameleon headers compatible with C++ complex types
** DONE merge request to make Chameleon headers compatible with C++ complex types
* TODO reproducible results and visualization
** TODO Improve Logging
*** DONE Add timer for iterations (global)
......
......@@ -56,7 +56,7 @@ void run_test_algo(const struct gengetopt_args_info &info, Algo algo, Restart re
Block XExact{dim, nbRHS};
Block RHS{dim, nbRHS};
std::vector<P> epsilon{1e-4};
std::vector<P> epsilon{static_cast<P>(info.epsilon_arg)};
RandomMatrixLoader rnld;
rnld.LoadMatrix(XExact);
......
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