From c58e487c34213347636db95b9de6e3db660b6c4b Mon Sep 17 00:00:00 2001
From: couturie <raphael.couturier@univ-fcomte.Fr>
Date: Tue, 8 Sep 2015 19:56:48 +0200
Subject: [PATCH] new

---
 IJHPCN/nb_iter_sec_ex45.pdf       | Bin 7345 -> 0 bytes
 IJHPCN/nb_iter_sec_ex45_curie.txt |   2 +-
 IJHPCN/paper.tex                  |  20 ++++++++++++++------
 3 files changed, 15 insertions(+), 7 deletions(-)
 delete mode 100644 IJHPCN/nb_iter_sec_ex45.pdf

diff --git a/IJHPCN/nb_iter_sec_ex45.pdf b/IJHPCN/nb_iter_sec_ex45.pdf
deleted file mode 100644
index 5c0750f0043e4feb6f5d896ad7cbdbcbfb6785cb..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 7345
zcmb_h2|SeB`@e0fNhPGT=p~F|=A9X{V<cO$RF<+t#(0NeHZx<rmedW2x}`+v%9^C4
za<8qjw1{rfDxp$Z6-voce&?NNx%dA5%ja`vK4#AQoagyI&v}+}oM&{_T3F$Tcrr%k
za`&!m3<V-UBA)<^u`z6i2>oRK5EW1yVJjYAhDadTiqDoI76?b=i?CQ2sZ4^f1(>jm
zLw6d5DnwO<hi5Zu8&6f2?&`5!k+)K72lbwd0c4}U_DJ4z<>Fntn)=#u;~e%{e*T<t
zl@@ILp=#Cfq+1_52LxEWC?}GR$?4j}cG;t!)m7TN%;sQ5w?%J$zV#zzSG=d<f&0OO
zBU{62R37el@cW^WvSrU4ZftIAo6*^jp3!g^ch0E%#f5E(L+wxZm$fDwOtk7v^6}Z*
zkmy+Hr5KlM9LvKjSgf^!Lb?)6DUUJ2*u?7RaL5I#3Y!h@ytw6>l$>|K^<|T{K+jaW
zguuyZOupFSTDU#$vm!NMUwW{Unv8n&PQAh1eO=Lv8Lj*8*>&?tv8S(plUrZi&0lSn
z^n1ztg|?KETbVUh#|%mc2OaOEh4kENdE;96A~i3X5g5!Ln0+;WZ!n9j#^D^*O&v0H
z^>g+Pq2~E();?L(aVW?3X!BEKQupyhhU$XW26C9y@Kx_*ZAV+}t~~CY<$1o=TsK{-
zqWIq~*7F7F$HSKe@-0Fy)TUiDO*_8ZHDG&QRhm75`|V^%ZAHuJ`RT`TrztQ}b3px-
zE<Yr{WQ@S&tdFQ2B6@H`IDy)!?mo+41ow;jLS6o9u%N4VQtus@cN%dcsZe*i`9;ZD
zys!aY4b1C?RQr}@l^Ist%Z*Nd?zpEk*RG}L*^Z!fi(Pl=W9#eAv9OO`^^V5tXm{iw
ze78>xUWHq77QcBBzsGd6d-m$j9ZQ}BosWuc)ik>LzR1ok;r`9J(VqvsgqgF4IiE7{
zmxJ76>yAfY@vKe7!s0voAM9#>etF_XQ$c&t%By>qSS}At4${0Ed}PG^61}aYX1|F}
z{E)$;jD4(EpG!`TuKv_C&>7wT=ed__uG2n+$6b1CzSPTl)7fx(yvHPGYS|jAj+L}m
z@E`hvFNR(Xy{(J5J-gb|O|xN--uZpSufyooigmH+aXCp#`|Ea}5!JCtg5i%bY4wvN
zm$8eYN_K^flt$BsoteSDWr;$G<&m9=8}9#kV&uZ^p~~Ke{e%6xHwbqv?3=AmieV!{
z-?6d)Beayz0z$uObRuSqME+J~&f#JNUd#4FV0*-u$6hH4gFFaennEB$hE#$l224s2
zp$sxa>ktDwB2rPXgo8*yg|3x|I8KPn0|qq=!OlpS46s{<$*i4Zph_VkiUUNfL_nNK
zls|x?h0Vb;K$Jpgu{&Y_Y0NBC#h2m-XuvK=ZkQB;ZF#;@$OD6#?TGP&075Vd`JE_=
zKDKXsku+z^*nE-SIK-Is_K1w_3uYkV7s5B)v>yS!Edv9uS!)jwS*$T5<*OP~<w(Pz
za*p3%<jWkNrDNzn&Z3WFlPlR$1a%B-<FwM<#Yuk^!Vf`YJPup#E(=7;k?_PakpyC(
z5OR0OXMj6-5~<AGpDlq5Nw6JzY)m9l;SD@rnZMKnL!zOI{@WX!1%DU{g$n&&{=UP(
z5FjcQT``SB0-xxZ3RX;$>l@dKV$uy!{}KTI7v6Es@37GMFK@v0)ege9@i-nyt_$Tw
zQO7vt))OE@lmcV`FPeHna4ztKe!vkzGoh6C_1P)dM~3D(ni?RI<N1cV_PeC{i|6H!
z4r~GP&mmd~O2q?5hj1h!9imZ55RFXqM9(KbDP%~-z-CeoS~h?RGM+77h4B3R#~#t>
zqI0My2DSl-%Hx;`{rCv*80;iN1nWWV00?`YR0@_O534+6u<IBM#gG9|K@NcO(}T@m
zbJ!BLf~{d2*bcUb9bqTf8D_&iFbDR95m>Go=E7WF2m*6K`NBL{01IIuPl&)GSOh-B
zuo#i>M82>Dmcmjb1QEhgUKlLpv!(v9RLtfeund;@OAzoKDuN>riD;SJLEvCU3?U+#
z4&Rl5&Pv3OCmi1^m=G6mL3-#mf#y6agF1c_k+cjW7afbjms9@xLmv3!``DR}|5rd_
zFz77S*Yn>cb%wK(^@7VH#qx5kY+bK)&DyL)V^t+H%gIZOCz2DkJlvwE1s%3ou|jG2
z3AIc`|G;C88>ZC!Auz*rPMQ)K=(1wBGI8_vl655&FW-bpbu-$h3;%e{U<|aCr>FNE
z936Q!Fmu<bm)3n@#;+dwA`X%ecb_eJl-;MF1?3C{<qU3f6|}M&RFd*lhc_u&P$PXg
zF2sg>%gLg+6@!qgQCjD_ut3`ngGc<HAo+Q@+{Q&YT562)GP{w-jlLq0=-%6^2cx^E
zX#3^m6(F;o=wRK_f6H3E-fy^m#72}y!mUfqzS_9z&pkTYI$$o&604W8;mVd9dwJ^W
zFBl~27;Jic<?s`mCyywweQBh9ssYLDkd9Nm)*h!jdMc#Rx4X0T;u$Yxqe<0Q4;zGs
ztg4k*`nWp2O-wjbzis)z6%8ljv-u7A_5&foIRz!obV}xeHMqlz3RVS7Ed5oIo2};;
zfbUvkpFL4;!}4D-*-vh%F5mk$x2hogmeOjA?MJgzZD;0m>uFFj7MLaHDy-Xjbh^cc
z#C)#fptDJjVMm15%$_LI;b@Pc2{Fqz)zI9hxaPTk@_p~9w7fvMDsTUdSKIEaw0@+U
zWtEVgXJ<WX_G!b1p63sEg21zN6``A7|1_1py@+owELh}2$@p1)my<d3e%#F|OZ-QU
z+vrqN$@teU)1!0lYkL0lz~8lMr^3SLYAuqPU5vKG_3LaaH_hvQenw@XN10@On-Jf#
zWx=EK#ihoTn^FpXeK;zUkh1Xo2ew-9`z?uuI#Ttszb0K1KGXgEU7_sy;m39C6Nwj9
zj4yTmwCM1X84iS+-Z@<JglAK`wQwIeo4nrD)%6~pa!z<+&TB}aIXHP;K%o2X<6RMF
zs%)*!T}d}DT43BxZc$-mr8gAxDvj<9n&z>T?__e0@vwD^wR*AA<g7o6XPc=$UODd)
z-_GO6`fc%fgIWqPPZbn@Sy`s>3*qV0+g1H4H*(U4_x_wc-?nhUHMgY!XAJriuQ{}J
zB;4lDTa*=;^YYXIrAZ%}VyCi{9!~jm(s?;m;c$YB4g0E!spb>?mQA8RtIc~iha0x#
zb$(bjgXyiLSkaud-(7aMb{0<6-%=q*e}fPHQA78IMUk`G4QN(@8k44ERHmw#()RZk
zmXY3GN`H3a_uiyw=U2f2x`Qp-D-YrdX&%m|1{Q*;udZpBN_u<m(h}2!Y_lu*oih{W
zgj#2H9a9@76(_t|@nPuM=w*YCCK^>YJ43a9vkSyOcTwkg2=5=6Sy8YR4vxg{J&o+>
zXcRq+wJ^A!j`8a`pgq56Pr#iHy@GnPD@%Pi8?JdC*chQ_eJc9F(<qO^^kw}H#wU3b
zl5{7|?5P-VP>B(?v5z*H4VvDHU;946_L&3Kq`(n2*xhnq7B`T6P;F;`WvP{3{HqRk
zqmp~7btMMp-poudCu=+<P2ZKK64~ELy3!d@I!k^2ucInOHFG~Mo|Dp`UHts^!-V4%
zC9~hI_<V0p(9MbSI1_3^6tBfKw=DLF{@L-<ldRlYH*tzy7^!lfRHz&^_25Z3=#5c}
zrA7aizMFj?e;PDDowW0WW4pe-{@}~BT+1@{{h@4FDdkj9XsK$ip{DxI2Kxy})gSPn
zlEASfCBg;f?XY*B=)o8WS(iL(^_8M$<xauf`|tmejZ0V;Q__>5Dx8u0wlVEco@>h9
z<yRV73R=%4yvPpy6ggC;Kon435N}Xc?mciCpR43i+<3nt;8T?yvS^QCXQ@h0{OGG_
zbIZr}+doG3-?cuq;PpVxjXgP)^SghNsO8Oj>c>qHHk|3&bfGm~H<MrH9@QOw)cO2C
z^)*gjra{%o3@88UV!K-{)Jm0?+jmJ9*t|WvJd1w*^#mM?*2CXgaa%)ML8*FU@4KM{
zozF3iX_?yF9<LNvcNO)e?5^XMjo!Psang+gzaZJOe2&cTv`#bU&ZSi<H(pX|4B@<6
z!e!cIoWXRIcog$LvIfFdAL{(L(dP(J*n(-gcuqk*X`AMnR`2J+7vu`>)>FUybfFB#
z-&*5+w;;UPu`=?FP!gaK<9Ob_>Gh3w<+s?sR>xK=Mx<a)ugqvQnx=2Rb#hvfUDJe_
z+qp_3E}NG!3wko%wK!xn&2^U?Ke$u>#jUt{PONtNI*qjWDDmN0@4S<oGe)O*vLZff
zsPwxgyepr4<DR%ZzJMT(ZBqRsd`*oSUt}7k`TDW?WLm6W<#U50B>l|YTQ9zjdvEZv
zC6d>ke#VdW-oZB{-p2Rn>AvoH9Pf!OGU^0Rhtva$Rrf0=WjEhw=x=RPIq0Bz+fc{8
z`^0YtLo1bX&gu`e?)-h8VE+g8T`9EKsAJE6Dmk&~NLT+G%{%XJ-7rmEOg|eqbvaXy
zr*@*fJd3`B7;W3>XKOOC3-`3l0q=m-uF{Sl-J%u}U}?29Ewkx)iAu~b;k}!;*-yqF
zy1hLm&oayZ1+6cya1K;(*h!zZ_+E$8WNP3s?&s<R(}1Svm^`U=0+;$YNO`1fBF2F`
z#fY7JQsa<bQN2lqQ-$(HL$5xYoys2d+J~1nZKIS|4SPn4b+<dv_J`ik-FN#;i6{Bd
zl9~ad<mY=>7p95&B}>ip{hoX{Wl6YpRHaEfBAglV_>soErJMX#Yx<U~7#X=+bYSkB
zrRy9WvL(sP&?|9+{<?|GH20LCmp!Ixnb_;gX5o`fe+zT#*f%jU*L-#(e*Nu3jPTj*
zmT~pNJFWYk1=mgAroUlvCG~CJ*$P^9+m7E{4PHFX-hSupL!+fjJ_}CyWvvYqeu%R_
z`(bii{^gj$VD;3@LuE&b+e}I~Z_pp<digrv$W(W#t+=P>&-*SmnyI=p_m8B)j%5T*
zzvv>1a9oz7LqeNMLTZs$p|{3>OYM4xoJWTYSCq!izHCS2Ps3~|PViJtH7VKUy=}3p
zqDqaMEoRe%Rd+@f44n48*fn^0%bdr7?_Xxm*ja6xYA2B-|2(N#DXQ=Gv2*6U#78w=
zCQCj@bc)I?M2T*Yvyaysrt$L3b*{Bqi#Oc7_PA<hq%eLbW#Gv9vq6RKp<RN<y$==?
zpHY4-<a*lE4~Zk1W4+1hQ;Jt(n(R81&bMba-!#HUMMa7%FP6`+UZ8X{NZ)7ulZK0l
zsz>i@b<+iIV$F~xHOqz*3Stwfu6OaBxpB~>a;@rs>nDGH{?f;+sq5oWjTt%Z8$0o5
z9~kId)bbff{B*4F3GS##t~b5vR?RB&mt99M=XHh-?23#Te6;H1q2|QbePwRvT9AVl
zYvyI%cPoFk_Fj5x!?A?sjT=cux@-3JZzBihNCI7VPYJ*9%3%%NQhhcvL2>cJ%RZ(i
zAEP2ksUB*}F9@Dai#?|FGV>%oCC;vFc)Rl9o0|{rygaDnMymHp8*XlX%HDFMe@(Dy
zH+HvIJ7)3ly%EB^^JRvM8QMK(XM`NtU(8zge$_p*w;BT<CegQEeNWMtn^NXgx-ml0
zXNvemjC<Q1S*q#H-hj+H25DHWGc(egOW|M!v#Xy>wv1eR0$XcPJ89O{!gU!t%W`hH
z3g?-K%rJG>g(Y+``%<<qk9f8G!_eZ2t=Dv-uzu5$s#@N}4XT*bn)`*tQ_fz0P;vb9
z$?P=~-z4u0jdglg{wcMzz=T0ej_P{~b<a<pk+#2Iy*_Hzj%t+~ESKr@_w!GYvjfZv
zCM4Ex{g`q3+0g0ZrTbi-PnZ<BW~SOBO{YkrX)dX&*>U^FPf58K2G1?+J4s71Q#eIa
z@Yd<n=-2McA5d2?JQ**jR9&K~ywKFof&Y4IensnL$H;rDN4#S6B|~R=>#uk{x@%-4
z@_cnnTRZv}g*Rv3IV?0#rc^tJm~x}1PsH%9*WEO2H+irz-@ICD{;vbJwe7jD*jU_t
z+Rx9XyC%jj+{k{uwW8tJyj;sLhuj^LFy@WdRLMlbUXLpQJ@xIz&$flzC!ZT^5WniU
z$vP1I@CB<^^X()fVcXQk(rQ)nwmsGX6FXWe*NOMVCFdtkyx)yEQKeYhpgQ-=rHY?7
zKh<A&@Qp~I;CL%$Lcr1>Tl5z8W2ZuXDgXXLHhx|Ib_;N|w;R6#5XlVk*DHXrX&AVq
z1R^rX2fcZiV0$XhV<Dce33da;o?tIFNBnuKA|!}Ygo86DB9O!I#hS7(#-?o0B0+Fr
zY%va8Q9y@L3W2+siI%yeos|~K5XKA>hy|c)0EG$oLMb!M1S=O|f-y>lv5=fa7HEQX
zwO<R#L?V75PX>|j3_KA>G&F#?5^#kK6-fdiBA$u`kP=@m(-GW+z5s$J6Rf{XCT7BL
zXlN)t)DSO{_`yU5g8>srFo}c%95`vXP{t0!38hPcKDi-0*RNtiDPC?9o+A>#JXe^A
zC&1&1#v2Jdkq|-^vwcLtG81fYFwd7sqH}$?BnksZM<`Srk;o<B2;hYy(P=~?hXFph
zq;YlP@67x`AD93N4S^rPhzUAEzQG(sVuH0^YX@2Ti=?tIjSPr^Cs48Q*z*4^t3Yi2
z*$W_`V{R}8JWPAg7T|_kfTk3SL?BXd1O|>kcP3Jq6bh3_)+dmd1cEXA7fuj-AZiXe
z%%J0piaG})5Y7p~i6Li*&V-0eAWvp8DbN?u?>N5`6-k^yO0fKd!D7Bh29fbpNX(XT
z{Q2OH&4;MpiH>m^!)Tno^U%NTGq%{?{(ttGK!Bz-%Ez=2ae@V)S7BrE|4EvEJI=nG
zuLg^QC46~>_;O%`kD#_mL8ytqVqXrEE0PG<AivpSG3cqVWgzWgke67NF_@I`WPF4r
zuRk-s?8|`9L+gpjXAAvIuwgh~gv$=*%djlHWnbdJ=dr#F<Pz{V6fokTsUebl0dhiu
z#t~T@P-DRo9`Y9mISUH$RS=AQ$2xt$(!T?ZXAPP<{|)R1fPcsxISK4HAdTR!yBP}u
z{Y$TZ59r$#4h{wrY$&=x|LurC`%Gge#Q0vL_fK2RGSD7F!voxk)~+5}yH1E0g6I%h
z_GD1`5V{orKCn5+SuB_+fxsZpiBt+jzAa!01nfAHNTN3$j=vT$U~xusOAib{yD9+T
z$iDKTp#x;IV3|Ka7|)DlU)#xH5`+sHphPkT-Fv^hAfVrn0&$@+8j-*Nt#b4Q3CC$v
z8UeIef20u!hM?L1BaJ|(|4o(v*7<in>fdBRPZ_jbf5an_sG!&W7aH}S_y~XJBN_gU
zMk0_w8}>&$3egb6d7LJbutC320)7&NU?*M#_+5bv!Hy!ZnW9a9M+mkNaz&850STI^
z7LbQA!-7h&B3RHY$#ion-HKvJwlHIu(TIkY6b6k<C$Ro^4o~#A1*wcJk;wx@We`ai
L9UaRxR+#?)2fbmt

diff --git a/IJHPCN/nb_iter_sec_ex45_curie.txt b/IJHPCN/nb_iter_sec_ex45_curie.txt
index 114b677..02a00b9 100644
--- a/IJHPCN/nb_iter_sec_ex45_curie.txt
+++ b/IJHPCN/nb_iter_sec_ex45_curie.txt
@@ -1,5 +1,5 @@
 #			fgmres/asm	tsirm/asm		fgmres/hypre
 512		123.55			227.454			0.0698
-2048	104.34			143.29			0.0268
+2048	104.34			173.68			0.0268
 4096	94.26				153.09			
 8192	72.08				142.40
diff --git a/IJHPCN/paper.tex b/IJHPCN/paper.tex
index 3aa76de..7b78160 100644
--- a/IJHPCN/paper.tex
+++ b/IJHPCN/paper.tex
@@ -638,7 +638,7 @@ preconditioners in PETSc, readers are referred to~\cite{petsc-web-page}.
 \hline
 
 \end{tabular}
-\caption{Comparison of FGMRES and TSIRM with FGMRES for example ex15 of PETSc with two preconditioners (mg and sor) having 25,000 components per core on Juqueen ($\epsilon_{tsirm}=1e-3$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
+\caption{Comparison of FGMRES and TSIRM with FGMRES for example ex15 of PETSc/KSP with two preconditioners (mg and sor) having 25,000 components per core on Juqueen ($\epsilon_{tsirm}=1e-3$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
 \label{tab:03}
 \end{center}
 \end{table*}
@@ -710,7 +710,7 @@ interesting.
 \hline
 
 \end{tabular}
-\caption{Comparison of FGMRES  and TSIRM with FGMRES algorithms for ex54 of Petsc (both with the MG preconditioner) with 25,000 components per core on Curie ($max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
+\caption{Comparison of FGMRES  and TSIRM with FGMRES algorithms for ex54 of PETSc/KSP (both with the MG preconditioner) with 25,000 components per core on Curie ($max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
 \label{tab:04}
 \end{center}
 \end{table*}
@@ -769,7 +769,7 @@ taken into account with TSIRM.
 \hline
 
 \end{tabular}
-\caption{Comparison of FGMRES  and TSIRM for ex54 of PETSc (both with the MG preconditioner) with 204,919,225 components on Curie with different number of cores ($\epsilon_{tsirm}=5e-5$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
+\caption{Comparison of FGMRES  and TSIRM for ex54 of PETSc/KSP (both with the MG preconditioner) with 204,919,225 components on Curie with different number of cores ($\epsilon_{tsirm}=5e-5$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
 \label{tab:05}
 \end{center}
 \end{table*}
@@ -784,7 +784,7 @@ taken into account with TSIRM.
 
 Concerning the  experiments some  other remarks are  interesting.
 \begin{itemize}
-\item We have tested other examples  of PETSc (ex29, ex45, ex49).  For all these
+\item We have tested other examples  of PETSc/KSP (ex29, ex45, ex49).  For all these
   examples,  we have also  obtained similar  gains between  GMRES and  TSIRM but
   those  examples are  not scalable  with many  cores. In  general, we  had some
   problems with more than $4,096$ cores.
@@ -815,18 +815,26 @@ Concerning the  experiments some  other remarks are  interesting.
 \cline{2-7}
                     & Time  & \# Iter.  & Time  & \# Iter. & Time  & \# Iter.   \\\hline \hline
    512              & 5.54      & 685    & 2.5 &       570  & 128.9 & 9     \\
-   2048             & 14.95     & 1,560  &  5.2 &     746   & 335.7 & 9 \\
+   2048             & 14.95     & 1,560  &  4.32 &     746   & 335.7 & 9 \\
    4096             & 25.13    & 2,369   & 5.61 &   859  &   >1000  & -- \\
    8192             & 44.35   & 3,197   &  7.6  &  1083 & >1000 &  --   \\
 
 \hline
 
 \end{tabular}
-\caption{Comparison of FGMRES  and TSIRM for ex45 of PETSc with two preconditioner (ASM and HYPRE)  having 25,000 components per core on Curie ($\epsilon_{tsirm}=1e-10$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
+\caption{Comparison of FGMRES  and TSIRM for ex45 of PETSc/KSP with two preconditioner (ASM and HYPRE)  having 25,000 components per core on Curie ($\epsilon_{tsirm}=1e-10$, $max\_iter_{kryl}=30$, $s=12$, $max\_iter_{ls}=15$, $\epsilon_{ls}=1e-40$),  time is expressed in seconds.}
 \label{tab:06}
 \end{center}
 \end{table*}
 
+
+\begin{figure}[htbp]
+\centering
+  \includegraphics[width=0.5\textwidth]{nb_iter_sec_ex45_curie}
+\caption{Number of iterations per second with ex45 and the same parameters as in Table~\ref{tab:06} (weak scaling)}
+\label{fig:03}
+\end{figure}
+
 %%ENDNEW
 
 %%%*********************************************************
-- 
2.39.5