To solve the neutron diffusion equation many linear systems has to be solved using preconditioned Krylov methods.

Traditional preconditioners based on incomplete factorizations are expensive in terms of memory.

Domain decomposition preconditioners are studied including substructuring preconditioners and additive Schwarz preconditioners.

2D and 3D benchmarks have been studied obtaining better performance results than usual preconditioners.