MLD2P4: A package of parallel algebraic multilevel domain decomposition preconditioners in Fortran 95

Domain decomposition ideas have long been an essential tool for the solution of PDEs on parallel computers. In recent years many research efforts have been focused on employing recursively domain decomposition methods to obtain multilevel preconditioners to be used with Krylov solvers. In this context, we developed MLD2P4 (MultiLevel Domain Decomposition Parallel Preconditioners Package based on PSBLAS), a package of parallel multilevel preconditioners that combines Additive Schwarz domain decomposition methods with a smoothed aggregation technique to build a hierarchy of coarse-level corrections in an algebraic way. The design of MLD2P4 was guided by objectives such as extensibility, flexibility, performance, portability and ease of use. They were achieved by following an object-oriented approach while using the Fortran 95 language, as well as by employing the PSBLAS library as basic framework. In this paper we present MLD2P4 focusing on its design principles, software architecture and use.