Efficient Parallel Resolution of The Simplified Transport Equations in Mixed-Dual Formulation
RAMET, Pierre
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
Parallel tools for Numerical Algorithms and Resolution of essentially Hyperbolic problems [BACCHUS]
Voir plus >
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
Parallel tools for Numerical Algorithms and Resolution of essentially Hyperbolic problems [BACCHUS]
RAMET, Pierre
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
Parallel tools for Numerical Algorithms and Resolution of essentially Hyperbolic problems [BACCHUS]
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
Parallel tools for Numerical Algorithms and Resolution of essentially Hyperbolic problems [BACCHUS]
ROMAN, Jean
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
High-End Parallel Algorithms for Challenging Numerical Simulations [HiePACS]
< Réduire
Laboratoire Bordelais de Recherche en Informatique [LaBRI]
High-End Parallel Algorithms for Challenging Numerical Simulations [HiePACS]
Langue
en
Article de revue
Ce document a été publié dans
Journal of Computational Physics. 2011, vol. 230, n° 5, p. 2004-2020
Elsevier
Résumé en anglais
A reactivity computation consists of computing the highest eigenvalue of a generalized eigenvalue problem, for which an inverse power algorithm is commonly used. Very fine modelizations are difficult to treat for our ...Lire la suite >
A reactivity computation consists of computing the highest eigenvalue of a generalized eigenvalue problem, for which an inverse power algorithm is commonly used. Very fine modelizations are difficult to treat for our sequential solver, based on the simplified transport equations, in terms of memory consumption and computational time. A first implementation of a Lagrangian based domain decomposition method brings to a poor parallel efficiency because of an increase in the power iterations. In order to obtain a high parallel efficiency, we improve the parallelization scheme by changing the location of the loop over the subdomains in the overall algorithm and by benefiting from the characteristics of the Raviart-Thomas finite element. The new parallel algorithm still allows us to locally adapt the numerical scheme (mesh, finite element order). However, it can be significantly optimized for the matching grid case. The good behavior of the new parallelization scheme is demonstrated for the matching grid case on several hundreds of nodes for computations based on a pin-by-pin discretization.< Réduire
Origine
Importé de halUnités de recherche