]> AND Private Git Repository - GMRES2stage.git/blobdiff - paper.tex
Logo AND Algorithmique Numérique Distribuée

Private GIT Repository
sdjfo
[GMRES2stage.git] / paper.tex
index 98f208f53ee4693c769022865d26765750903215..16b2de588eaeb4977240ffb3598fae1e9dd79cf3 100644 (file)
--- a/paper.tex
+++ b/paper.tex
@@ -884,21 +884,22 @@ torso3             & 2D/3D problem & 259,156 & 4,429,042 \\
 \label{tab:01}
 \end{center}
 \end{table}
-Chosen parameters  are detailed below.   As by default  the restart of  GMRES is
-performed  every 30  iterations,  we have  chosen  to stop  the  GMRES every  30
-iterations (\emph{i.e.} $max\_iter_{kryl}=30$).  $s$ is set to 8. CGLS is chosen
-to   minimize  the   least-squares  problem   with  the   following  parameters:
+Chosen parameters  are detailed below.   
+We have  stopped  the  GMRES every  30
+iterations (\emph{i.e.}, $max\_iter_{kryl}=30$), which is the default 
+setting of GMRES.  $s$, for its part, has been set to 8. CGLS 
+ minimizes  the   least-squares  problem   with  parameters
 $\epsilon_{ls}=1e-40$ and $max\_iter_{ls}=20$.  The external precision is set to
-$\epsilon_{tsirm}=1e-10$.  Those  experiments have been performed  on a Intel(R)
-Core(TM) i7-3630QM CPU @ 2.40GHz with the version 3.5.1 of PETSc.
+$\epsilon_{tsirm}=1e-10$.  These  experiments have been performed  on an Intel(R)
+Core(TM) i7-3630QM CPU @ 2.40GHz with the 3.5.1 version  of PETSc.
 
 
-In  Table~\ref{tab:02}, some  experiments comparing  the solving  of  the linear
-systems obtained with the previous matrices  with a GMRES variant and with TSIRM
-are given. In the  second column, it can be noticed that  either GMRES or FGMRES
-(Flexible GMRES)~\cite{Saad:1993} is used to solve the linear system.  According
-to  the matrices,  different preconditioners  are  used.  With  TSIRM, the  same
-solver and the  same preconditionner are used.  This Table  shows that TSIRM can
+Experiments comparing 
+a GMRES variant with TSIRM in the resolution of linear systems are given in  Table~\ref{tab:02}. 
+The  second column describes whether GMRES or FGMRES
+(Flexible GMRES~\cite{Saad:1993}) has been used for linear systems solving.  
+Different preconditioners  have been used according to the matrices.  With  TSIRM, the  same
+solver and the  same preconditionner are used.  This table  shows that TSIRM can
 drastically reduce  the number of iterations  to reach the  convergence when the
 number of iterations for  the normal GMRES is more or less  greater than 500. In
 fact this also depends on two parameters: the number of iterations to stop GMRES
@@ -923,7 +924,7 @@ torso3             & fgmres / sor  & 37.70 & 565 & 34.97 & 510 \\
 \hline
 
 \end{tabular}
-\caption{Comparison of (F)GMRES and TSIRM with (F)GMRES in sequential with some matrices, time is expressed in seconds.}
+\caption{Comparison between sequential standalone (F)GMRES and TSIRM with (F)GMRES (time in seconds).}
 \label{tab:02}
 \end{center}
 \end{table}
@@ -933,10 +934,10 @@ torso3             & fgmres / sor  & 37.70 & 565 & 34.97 & 510 \\
 
 
 In order to perform larger experiments, we have tested some example applications
-of  PETSc. Those  applications are  available in  the \emph{ksp}  part  which is
+of  PETSc. Those  applications are  available in  the \emph{ksp}  part,  which is
 suited for scalable linear equations solvers:
 \begin{itemize}
-\item ex15  is an example  which solves in  parallel an operator using  a finite
+\item ex15  is an example  that solves in  parallel an operator using  a finite
   difference  scheme.   The  diagonal  is  equal to  4  and  4  extra-diagonals
   representing the neighbors in each directions  are equal to -1. This example is
   used  in many  physical phenomena, for  example, heat  and fluid  flow, wave