]> AND Private Git Repository - Krylov_multi.git/commitdiff
Logo AND Algorithmique Numérique Distribuée

Private GIT Repository
new
authorcouturie <couturie@extinction>
Tue, 17 Dec 2013 13:28:45 +0000 (14:28 +0100)
committercouturie <couturie@extinction>
Tue, 17 Dec 2013 13:28:45 +0000 (14:28 +0100)
biblio.bib
krylov_multi.tex

index 8d0b469e7ed28153e064c1713afe76e2fb474f43..d1f3b0090530934908f05724e1bfe8da92d4de38 100644 (file)
@@ -6,4 +6,70 @@
   pages={9--29},
   year={1993},
   publisher={Elsevier}
+}
+
+@article{o1985multi,
+  title={Multi-splittings of matrices and parallel solution of linear systems},
+  author={O'Leary, Dianne P. and White, Robert E.},
+  journal={SIAM Journal on Algebraic Discrete Methods},
+  volume={6},
+  number={4},
+  pages={630--640},
+  year={1985},
+  publisher={SIAM}
+}
+
+
+@article{frommer1992h,
+  title={H-splittings and two-stage iterative methods},
+  author={Frommer, Andreas and Szyld, Daniel B.},
+  journal={Numerische Mathematik},
+  volume={63},
+  number={1},
+  pages={345--356},
+  year={1992},
+  publisher={Springer-Verlag}
+}
+
+@article{bru1995parallel,
+  title={Parallel, synchronous and asynchronous two-stage multisplitting methods},
+  author={Bru, Rafael and Migall{\'o}n, Violeta and Penad{\'e}s, Jos{\'e} and Szyld, Daniel B},
+  journal={Electronic Transactions on Numerical Analysis},
+  volume={3},
+  pages={24--38},
+  year={1995}
+}
+
+@article{bai1999block,
+  title={Block and asynchronous two-stage methods for mildly nonlinear systems},
+  author={Bai, Zhong-Zhi and Migall{\'o}n, Violeta and Penad{\'e}s, Jos{\'e} and Szyld, Daniel B.},
+  journal={Numerische Mathematik},
+  volume={82},
+  number={1},
+  pages={1--20},
+  year={1999},
+  publisher={Springer}
+}
+
+@article{bahi2000asynchronous,
+  title={Asynchronous iterative algorithms for nonexpansive linear systems},
+  author={Bahi, Jacques M.},
+  journal={Journal of Parallel and Distributed Computing},
+  volume={60},
+  number={1},
+  pages={92--112},
+  year={2000},
+  publisher={Elsevier}
+}
+
+
+@article{couturier2008gremlins,
+  title={GREMLINS: a large sparse linear solver for grid environment},
+  author={Couturier, Rapha{\"e}l and Denis, Christophe and J{\'e}z{\'e}quel, Fabienne},
+  journal={Parallel Computing},
+  volume={34},
+  number={6},
+  pages={380--391},
+  year={2008},
+  publisher={Elsevier}
 }
\ No newline at end of file
index 8a64840d9ddb1bad5ccfe26e06b39383b362c9de..0fd3b7919cddd12d22e8c9750d7b800db22081e0 100644 (file)
 \begin{abstract}
 In  this  paper we  revist  the  krylov  multisplitting algorithm  presented  in
 \cite{huang1993krylov}  which  uses  a  scalar  method to  minimize  the  krylov
-iterations computed by a multisplitting algorithm. Our new algorithm is simply a
-parallel multisplitting algorithm  with few blocks of large  size and a parallel
-krylov  minimization  is used  to  improve  the  convergence. Some  large  scale
-experiments  with a 3D  Poisson problem  are presented.  They show  the obtained
-improvements compared to a classical GMRES both in terms of number of iterations
-and execution times.
+iterations computed by a multisplitting algorithm. Our new algorithm is based on
+a  parallel multisplitting  algorithm  with few  blocks  of large  size using  a
+parallel GMRES method inside each block and on a parallel krylov minimization in
+order to improve the convergence. Some large scale experiments with a 3D Poisson
+problem  are presented.   They  show  the obtained  improvements  compared to  a
+classical GMRES both in terms of number of iterations and execution times.
 \end{abstract}
 
 \section{Introduction}
 
 Iterative methods are used to solve  large sparse linear systems of equations of
 the form  $Ax=b$ because they are  easier to parallelize than  direct ones. Many
-iterative  methods have  been proposed  and adpated  by many  researchers.  When
+iterative  methods have  been proposed  and  adapted by  many researchers.  When
 solving large  linear systems  with many cores,  iterative methods  often suffer
 from  scalability  problems.    This  is  due  to  their   need  for  collective
-communications to perform matrix-vector products and reduction operations.
+communications  to  perform  matrix-vector  products and  reduction  operations.
+Preconditionners can be  used in order to increase  the convergence of iterative
+solvers.   However, most  of the  good preconditionners  are not  sclalable when
+thousands of cores are used.
+
+
+A completer...
+On ne peut pas parler de tout...
+
+\section{Related works}
+
+
+A general framework  for studying parallel multisplitting has  been presented in
+\cite{o1985multi} by O'Leary and White. Convergence conditions are given for the
+most general case.  Many authors improved multisplitting algorithms by proposing
+for  example  a  asynchronous  version  \cite{bru1995parallel}  and  convergence
+condition  \cite{bai1999block,bahi2000asynchronous}   in  this  case   or  other
+two-stage algorithms~\cite{frommer1992h,bru1995parallel}
+
+In  \cite{huang1993krylov},  the  authors  proposed  a  parallel  multisplitting
+algorithm in which all the tasks except  one are devoted to solve a sub-block of
+the splitting  and to send their  local solution to  the first task which  is in
+charge to  combine the vectors at  each iteration.  These vectors  form a Krylov
+basis for  which the first tasks minimize  the error function over  the basis to
+increase the convergence, then the other tasks receive the update solution until
+convergence of the global system. 
+
+
+
+In \cite{couturier2008gremlins}, the  authors proposed practical implementations
+of multisplitting algorithms that take benefit from multisplitting algorithms to
+solve large scale linear systems. Inner  solvers could be based on scalar direct
+method with the LU method or scalar iterative one with GMRES.
 
 \bibliographystyle{plain}
 \bibliography{biblio}