Simulation of laser propagation in a plasma with a frequency wave equation

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, we use an iterative Krylov method preconditioned by a separable matrix. The corresponding linear system is solved with a block cyclic reduction method. Some enlightenments 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.