From a45b1d7df54f09b78efd3333be8d8cc4a89dccdf Mon Sep 17 00:00:00 2001 From: ziane Date: Tue, 22 Sep 2015 18:36:07 +0200 Subject: [PATCH] figure --- IJHPCN/paper.tex | 71 +++++-------------------------------------- IJHPCN/snes_ex14.pdf | Bin 0 -> 7298 bytes 2 files changed, 7 insertions(+), 64 deletions(-) create mode 100644 IJHPCN/snes_ex14.pdf diff --git a/IJHPCN/paper.tex b/IJHPCN/paper.tex index db5f791..2d0bd55 100644 --- a/IJHPCN/paper.tex +++ b/IJHPCN/paper.tex @@ -1035,80 +1035,23 @@ cores to more than 16 with 8,192 cores. \subsection{Influence of parameters for TSIRM} In this section we present some experimental results in order to study the influence of some parameters on the TSIRM algorithm. We conducted experiments on $16$ cores to solve 3D problems of size $200,000$ components per core. We solved nonlinear problems token from examples of PETSc. We fixed some parameters of the TSIRM algorithm as follows: the nonlinear systems are solved with a precision of $10^{-8}$, block Jacobi preconditioner is used, the tolerance threshold $\epsilon_{tsirm}$ is $10^{-8}$ , the maximum number of iterations $max\_iter_{tsirm}$ is set to $10,000$ iterations, the FGMRES method is used as the inner solver with a tolerance threshold $\epsilon_{kryl}=10^{-10}$ and the least-squares problem is solved with a precision $\epsilon_{ls}=10^{-40}$ in the minimization process. -%time mpirun ../ex48 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -snes_rtol 1.e-8 -snes_monitor -ksp_type tsirm -ksp_pc_type bjacobi -pc_type ksp -ksp_tsirm_tol 1e-8 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 30 -ksp_tsirm_inner_restarts 30 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 15 -ksp_tsirm_size_ls 10 +%time mpirun ../ex34 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -ksp_type tsirm -ksp_pc_type asm -pc_type ksp -ksp_tsirm_tol 1e-10 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 30 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 20 -ksp_tsirm_size_ls 10 \begin{figure}[htbp] \centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_cgls_iter_total} -\caption{Number of total iterations using two different methods for the minimization: LSQR and CGLS.} -\label{fig:cgls-iter} + \includegraphics[width=0.5\textwidth]{ksp_tsirm_cgls} +\caption{Number of total iterations using two different methods for the minimization: CGLS and LSQR.} +\label{fig:cgls} \end{figure} \begin{figure}[htbp] \centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_cgls_time} -\caption{Execution time in seconds using two different methods for the minimization: LSQR and CGLS.} -\label{fig:cgls-time} + \includegraphics[width=0.5\textwidth]{snes_ex14} +\caption{Total number of iterations in example {\it snes ex14} of PETSc by varyin the number of inner iterations and the size of the least-squares problem.} +\label{fig:snes_ex14} \end{figure} -%time mpirun ../ex35 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -snes_rtol 1.e-8 -snes_monitor -ksp_type tsirm -ksp_pc_type bjacobi -pc_type ksp -ksp_tsirm_tol 1e-8 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 30 -ksp_tsirm_inner_restarts 38 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 15 -ksp_tsirm_size_ls 10 -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_inner_restarts_iter_total} -\caption{Number of total iterations with variation of restarts in the inner solver FGMRES.} -\label{fig:inner_restarts_iter_total} -\end{figure} -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_inner_restarts_time} -\caption{Execution time in seconds with variation of restarts in the inner solver FGMRES.} -\label{fig:inner_restarts_time} -\end{figure} -%time mpirun ../ex14 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -snes_rtol 1.e-8 -snes_monitor -ksp_type tsirm -ksp_pc_type bjacobi -pc_type ksp -ksp_tsirm_tol 1e-8 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 1000 -ksp_tsirm_inner_restarts 30 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 15 -ksp_tsirm_size_ls 10 -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_max_inner_iter} -\caption{Number of total iterations with variation of number of inner iterations.} -\label{fig:max_inner_iter} -\end{figure} - -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_max_inner_time} -\caption{Execution time in seconds with variation of number of inner iterations.} -\label{fig:max_inner_time} -\end{figure} - -%time mpirun ../ex14 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -snes_rtol 1.e-8 -snes_monitor -ksp_type tsirm -ksp_pc_type bjacobi -pc_type ksp -ksp_tsirm_tol 1e-8 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 30 -ksp_tsirm_inner_restarts 30 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 5 -ksp_tsirm_size_ls 10 -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_maxiter_ls_iter} -\caption{Number of total iterations with variation of number of iterations in the minimization process.} -\label{fig:maxiter_ls_iter} -\end{figure} - -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_maxiter_ls_time} -\caption{Execution time in seconds with variation of number of iterations in the minimization process.} -\label{fig:maxiter_ls_time} -\end{figure} - -%time mpirun ../ex14 -da_grid_x 147 -da_grid_y 147 -da_grid_z 147 -snes_rtol 1.e-8 -snes_monitor -ksp_type tsirm -ksp_pc_type bjacobi -pc_type ksp -ksp_tsirm_tol 1e-8 -ksp_tsirm_maxiter 10000 -ksp_ksp_type fgmres -ksp_tsirm_max_inner_iter 30 -ksp_tsirm_inner_restarts 30 -ksp_tsirm_inner_tol 1e-10 -ksp_tsirm_cgls 0 -ksp_tsirm_tol_ls 1.e-40 -ksp_tsirm_maxiter_ls 15 -ksp_tsirm_size_ls 2 -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_size_ls_iter} -\caption{Number of total iterations with variation of the size of the least-squares problem in the minimization process.} -\label{fig:size_ls_iter} -\end{figure} - -\begin{figure}[htbp] -\centering - \includegraphics[angle=-90,width=0.5\textwidth]{ksp_tsirm_size_ls_time} -\caption{Execution time in seconds with variation of the size of the least-squares problem in the minimization process.} -\label{fig:size_ls_time} -\end{figure} %%ENDNEW diff --git a/IJHPCN/snes_ex14.pdf b/IJHPCN/snes_ex14.pdf new file mode 100644 index 0000000000000000000000000000000000000000..4ec9de40f287177f367f6cc9c9422fe21b65ee18 GIT binary patch literal 7298 zcmb_B2|SeD*Of$LsaHvndXR=Rp4l0U>`RPnMUfga55_XXEFp!`W=o|+N_H*wB5O)Y zDrAYul2D>hQi{s|nTeLS_xry8_xrswKfk&6Z1>!A&b{}XL)qHM1cgz@Bb003#H1n! z02*MqdLndmKud_}&h-F@FvJ!#VR&;PHUOG2{vsCi^({~J^)oS_vTtQI1Ft*IpU#HV)M&-iEV2GcGsKT%W}VI zSS>$aF;&%*!BoA$P|FWHEp#?q98FVG%^FgkGd}6@yy(-S;U&`3$@$K6Wu8RH2IM+* z;c68)!i_rNh88}L@>(TPP~@sqm2mRhlDNPZC}B-jVB^Y|2N>o1It#Bp@LD5!@som( zVn=y)<$A{x$~PT)Dih)ZitY_1Y0H-?U09I3cff751d4g=v`i=WY#KhHLI#)RB|Rrq z;>|~Uji+jQhQekldUg@`~yCeZsTakSRtKOqxN`>9BEhA+gIcAiI@A#nTw?=6j#V=i{+5+%5#n`%=_%x zy*0X<;(G6fUoFFnji*Sa82Jr%MH$EPglpbP=TTHo->qln7ieu@L5&IiortTK@U|D# zLoTx!AKcPp4w+oOtK2ZJ-?h#Az{^*KbxB1lFP(Vk5rsafH007ml2A@RsBE$+zK>FQ zz}qKQJb?l5$j5Gmdiq2xwElv>T>Mw@)(Hbqf~TRx$f-TYTZV~=j`O#4-Lh;|nwS1Y zvCTd*#r#YoCZj*&*gAckrs*?pWsBuc`z5~ByPS4B&t&h+hNi#xhZvfh|rPEJ7NrLL>9}PYsK`A$pJT9X|bOc9AB;QR8`pM%DZ+2OYWHD z&6vgW_Z$XF?)dqOa$TzobyjKgv&Zd2-lNb@iY|AO$CNBDP83JnXYEr_ulGM9bF>0w zDQ}c3J0gFl#C*bXvY?7+TK%BxI?;q!H8DOlHdOwhZ^vlW$~KSpKKWhsA8(-7oRv*D zWZJX%;xX?Au7kDYZZY%NlZ#F@J*H|6dgUmIIUEX;@^2x7p;TeYrIq9ax2s*X-sR6i z`xYFzNo{UisXFddGMIppaM2%L8gmw(EcZM`Z}pXuSQGJLF>1_`4#hj4nGs^V?M>Dv zn<)iPT7ly78^zjYgUQRq&Joj$_3PKkD=8&x`=D{~%O0=!x7HL0ebl>>m;TN?Tlk>i zww%XXR^?UZ_|+Xv8TBC?0liHMp)K#p==ZA9+pLStYkT`JFU+o-Uf+?Bc;w^m$)h_{ zu{lP{&S!IEo@R%IR4-RnX?75ekn(t@a`MuPl{>N%d|F3rIE!&hXP zd7l2WlWbQ0_0^s-9H*-W55%r!Nhpjz=00-{eWbKT^zp8R!@?RttPhK$L`y$V;rOzd z6S-h_x*t{e&GCsx$Dig^hFt}5DHesMTPx)HEa;%Dn+Vp4>Ol%{LG^?NV_DPOX=oS5*sb014P!RY}E1%WqGK zQa#!8<*THh?!Mls4VDS@Eevgv3PO&%MGhdp1-uG$w-J9H#TnV&9XBVDut;ql)}hVn zMrGmcCCcnR&zm)O%TJAei~Q8t`qzTp0qa9Hj(=&slHHTr5w3e@VUqh$XFJC|T7=xb z_rdN6SJ&1KrJVW)*iBP`FRiw3(0G1P?Imu7Cw}p_;RmHkFY58!)Ufw@NtrZed+F9g zX^H8tJ|>GSk&6*0y=@h;i&6?oO6ABbk}neLAJ)&XlJ(2~Xq2QB6ulJN_#rU&(6?3N zlOKVzu^!QiSoUY5%q3YAOI*w(+NV6Qg#7ALc1l?4xeTc1u1alMale$(ZVfpuX{@V>PD+ni*-{-iCt#E_AC;Em;#g zin|u4bv;_x?#goc{JJ)zJFlw3Ys=Hj+aJZ%zTE+}1~u+iZ8V)yf|&G~M;0FOA7=h@ z%>Tt`YT^(xF#OL)+dj}20chB77RKEa6<6awnc0QUsY*U*E?25?x=k~%ZSV6mVBgUkLf4}e} zei{EJ_s2>z6htE5F8}!bMaB;<0E@-`pcX&{R5%gDz^4L*#Q{Vt2GGD` zw(}EPcMcp!7z`SZ0QEUE{vm-eaC}pJ%^-%m$IO)POFj6(5SznrXM!wv>I<$hFfxE%o z5SyhcumYU8LpT7#_s6es&7aTM?|W+hYi18`4$0)+CtRjwAc&Hz@Y6&po?tulnI-4w~Z&mRhA7 z;lBnCzJr`Nkz<8LTNh+Bo;k2)N62*Ow4y$hO1&3tJpDP%Ik&E{?WKyl=xSlvt~Te& zwYB-yrIwmGEs(oPQCVtd_8w@rTjtspeGAKWMzKm-A@50Psgu}A-TQk^-Mr#1Qol6r zp?~2Axoa=f9u1N{-*bfao0L;dCryL^#x}8>u`!8)4#Jb6ho$!d-`}bbau*S9T$pkp zU?U|bv&|ZJvB-X@wBbSPSXGUP%+8CdCq@8F-CL1z2-CXBRn5-qh$uj~q!RmZGh3wb@fgYvrdgnChYRD;_W z$}FoBU!7pAx#0e%UYSdBHNpi(#VT!Jw(gR(L(-WOtT)-~_0FiOswOJ6dtKMu+EeuU z)Ahddt%aO;!UHBHyk$PmsU|t}0U##1@eMmXX(@^Fk z?q1b_$~Dd&Y;j4kj6>Cl@9PrpZ=QIwaczOgGaDC|$-xzQ<;`7F{?Z|#j>AbxMF=3w zAj82XkFrgDbOoZO&ft80Yo?5>#Fwm9$O&xcoC2JI#s#taW&Mvdt4AWqT_sj4!T$Up zlWODI;~l%)>!emT8C%-sE(yM&E~Sv2B9|JyxEWOy`e9dBvfSZ!e?7jPmbABCd#F%E z#^s&#)p$(okvgfu{yJ6<%69}e&n_jtc6n+WBh5x}+m525g14_LcgXCZT10<&emrsb z>ycp}hl8v1L?RZc(2pUbE~G9lmTtZ_W!!^o#Xoats(n1(8`v_nBEcsvMJ6Vpzdlad z5bbFrQR>$fAiBCQ@@B`inDtT!(RT}#Ms8Au+yfDg1zW1g!Dn_*8)1{&I{h|p?m1CJ_;>w03TaRWBb?AmDR(bgi zlnvPXkT_j^!KU|B^5w(lD!WMPTka&@T#t*E*=DTS9`~T*T)fc~y>8P34OEC~`Ju

G6%J&lYcmpdK8SK%$)vPU{sjL5OzG7Qk*{fu$j;`31Mx{p zjrwy(^bHRL1@V|COl;;^3q^7Oj6k`va!Jwx|Gdp^YrQ~mBBb=-vAWvkn?7e-YsAT`85VZ3qsyISN7c+O#|Kdvc8%v{@+1+_)46$h!)3-UVJ9%Uc%?==`$+r6^5nD?dMOvF>k zNOrUJ+N!7x+&)jg~~BstQq_QB8|; z`cr_EE%TO4`{Nok&s{uZR`*Oe6Jsr#-MhoaIe{}z#d1MqdOc{I99$O^RJG?tw7NI8 zNn*=|qKKX-^hIeEU3a^Uo@W*2{G}1Mr_)eVbWW8^&rW23Pf*#2NFF6etbI9`94PM$Y(QwMGTSrHZp%MH(mBpX6is9SLONL%t935Qu$V;{jc4>crRGwwdHl^ z6$JCm$H5KPy5G-vvR}Q|yh7&ilIgE5%QDlSTR&=weZh;~TwRGLWT=~~H@gQX-+h)g zwbvOlb?;lxlxS^S%R{BR4~72H3D4K7dA5pDE9aWtzmKdhrE*X$Cg4P^-m0q>^g|n+ zZ6y6(6r74;e7P~dC99!9=)88hb!*k!x!XHR`I*q~rGX#`{xvC@O^<)R=Z+N1ACEZ! zjs5YM>ud;svj{JU3%Ek=3?>;lT3&?&7<4kyfj~i1d<`KFhS_d5WVhSOp0?YIMxrA( z=puABP~n;cLIqHLQEq|;AqRl7BeJ5Qt)+<~pCCXhz{l4I&eDMZA8#f{D}am?uxPcQf0zZVVw zNgYi3eg(54;)t^Q7fXWxA7*0Vq1;M&)^Pk-8iv#MuD_VDRdYFXy{UycxRR2LfjBr$Bf_ z<3|aL{avOVW0J8x4sW1=GeK*hO)ywP-T#|n ZJHNZf;ZoULfrE&eL;^xt*?5x);y;PgUjYCB literal 0 HcmV?d00001 -- 2.39.5