  The aim of this work is to perform numerical simulations of the propagation
of a laser in a plasma. At each time step, one has to solve a Helmholtz
equation in a domain which consists in some hundreds of millions of cells. To
solve this huge linear system, one uses a iterative Krylov method with a
preconditioning by a separable matrix. The corresponding linear system is
solved with a block cyclic reduction method. Some enlightments on the parallel
implementation are also given. Lastly, numerical results are presented
including some features concerning the scalability of the numerical method on a
parallel architecture.
