Parallel multistage preconditioners by Hierarchical Interface Decomposition on “T2K Open Super Computer (Todai Combined Cluster)” with Hybrid parallel programming models

In this work, parallel preconditioning methods based on ldquoHierarchical Interface Decomposition (HID)rdquo and hybrid parallel programming models were applied to finite-element based simulations of linear elasticity problems in media with heterogeneous material properties. Reverse Cuthill-McKee reordering with cyclic multicoloring (CM-RCM) was applied for parallelism through OpenMP. The developed code has been tested on the ldquoT2K Open Super Computer (Todai Combined Cluster)rdquo using up to 512 cores. Preconditioners based on HID provide a scalable performance and robustness in comparison to conventional localized block Jacobi preconditioners. Performance of Hybrid 4x4 parallel programming model is competitive with that of Flat MPI.

[1]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[2]  Rolf Rabenseifner,et al.  Communication Bandwidth of Parallel Programming Models on Hybrid Architectures , 2009, ISHPC.

[3]  Kengo Nakajima.  Parallel Multistage Preconditioners Based on a Hierarchical Graph Decomposition for SMP Cluster Architectures with a Hybrid Parallel Programming Model , 2007, HPCC.

[4]  Shao-Liang Zhang,et al.  GPBi-CG: Generalized Product-type Methods Based on Bi-CG for Solving Nonsymmetric Linear Systems , 1997, SIAM J. Sci. Comput..

[5]  Kengo Nakajima,et al.  The Impact of Parallel Programming Models on the Linear Algebra Performance for Finite Element Simulations , 2006 .

[6]  K. Nakajima.  Parallel Iterative Solvers of GeoFEM with Selective Blocking Preconditioning for Nonlinear Contact Problems on the Earth Simulator , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[7]  Kengo Nakajima.  Strategies for Preconditioning Methods of Parallel Iterative Solvers for Finite-Element Applications in Geophysics , 2009 .

[8]  Kengo Nakajima.  Parallel Iterative Solvers of GeoFEM with Selective Blocking Preconditioning for Nonlinear Contact Problems on the Earth Simulator , 2003, SC.

[9]  Yousef Saad,et al.  A Parallel Multistage ILU Factorization Based on a Hierarchical Graph Decomposition , 2006, SIAM J. Sci. Comput..

[10]  Kengo Nakajima,et al.  The Impact of Parallel Programming Models on the Performance of Iterative Linear Solvers for Finite Element Applications , 2006, VECPAR.