A parallel implementation of a two-level overlapping Schwarz method with energy-minimizing coarse space based on trilinos. (English) Zbl 1355.65168
Summary: We describe a new implementation of a two-level overlapping Schwarz preconditioner with energy-minimizing coarse space (GDSW: generalized Dryja-Smith-Widlund) and show numerical results for an additive and a hybrid additive-multiplicative version. Our parallel implementation makes use of the Trilinos software library and provides a framework for parallel two-level Schwarz methods. We show parallel scalability for two- and three-dimensional scalar second-order elliptic and linear elasticity problems for several thousands of cores. We also discuss techniques for the parallel construction of coarse spaces which are also of interest for other parallel preconditioners and discretization methods using energy minimizing coarse functions. We finally show an application in monolithic fluid-structure interaction, where significant improvements are achieved compared to a standard algebraic, one-level overlapping Schwarz method.
MSC:
| 65N55 | Multigrid methods; domain decomposition for boundary value problems involving PDEs |
| 65Y05 | Parallel numerical computation |
| 35J25 | Boundary value problems for second-order elliptic equations |
| 74B05 | Classical linear elasticity |
| 74S30 | Other numerical methods in solid mechanics (MSC2010) |
| 65F08 | Preconditioners for iterative methods |
Keywords:
domain decomposition; scientific software; parallel; overlapping Schwarz; fluid-structure interaction; GDSW: generalized Dryja-Smith-Widlund; numerical results; second-order elliptic; linear elasticity; preconditionersReferences:
| [1] | J. Aarnes and T. Y. Hou, {\it Multiscale domain decomposition methods for elliptic problems with high aspect ratios}, Acta Math. Appl. Sin. Engl. Ser., 18 (2002), pp. 63-76, . · Zbl 1003.65142 |
| [2] | P. R. Amestoy, I. S. Duff, J.-Y. L’Excellent, and J. Koster, {\it A fully asynchronous multifrontal solver using distributed dynamic scheduling}, SIAM J. Matrix Anal. Appl., 23 (2001), pp. 15-41, . · Zbl 0992.65018 |
| [3] | E. Anderson, Z. Bai, C. Bischof, L. S. Blackford, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, and D. Sorensen, {\it LAPACK Users’ Guide}, 3rd ed., SIAM, Philadelphia, 1999, . · Zbl 0934.65030 |
| [4] | S. Badia, A. F. Martín, and J. Principe, {\it Multilevel balancing domain decomposition at extreme scales}, SIAM J. Sci. Comput., 38 (2016), pp. C22-C52, . · Zbl 1334.65217 |
| [5] | S. Badia, A. Quaini, and A. Quarteroni, {\it Splitting methods based on algebraic factorization for fluid-structure interaction}, SIAM J. Sci. Comput., 30 (2008), pp. 1778-1805, . · Zbl 1368.74021 |
| [6] | A. H. Baker, A. Klawonn, T. Kolev, M. Lanser, O. Rheinbach, and U. M. Yang, {\it Scalability of classical algebraic multigrid for elasticity to half a million parallel tasks}, in Software for Exascale Computing - SPPEXA 2013-2015, Lect. Notes Comput. Sci. Eng. 113, H.-J. Bungartz, P. Neumann, and E. W. Nagel, eds., Springer, Cham, Switzerland, 2016, pp. 113-140, ; also available as TUBAF Preprint 2015-14 online from . |
| [7] | G. Ballard, C. Siefert, and J. Hu, {\it Reducing communication costs for sparse matrix multiplication within algebraic multigrid}, SIAM J. Sci. Comput., 38 (2016), pp. C203-C231, . · Zbl 1339.65058 |
| [8] | D. Balzani, D. Böse, D. Brands, R. Erbel, A. Klawonn, O. Rheinbach, and J. Schröder, {\it Parallel simulation of patient-specific atherosclerotic arteries for the enhancement of intravascular ultrasound diagnostics}, Eng. Computation, 29 (2012), pp. 888-906, . |
| [9] | D. Balzani, S. Deparis, S. Fausten, D. Forti, A. Heinlein, A. Klawonn, A. Quarteroni, O. Rheinbach, and J. Schröder, {\it Aspects of arterial wall simulations: Nonlinear anisotropic material models and fluid structure interaction}, in Proceedings of the WCCM XI, 2014. |
| [10] | D. Balzani, S. Deparis, S. Fausten, D. Forti, A. Heinlein, A. Klawonn, A. Quarteroni, O. Rheinbach, and J. Schröder, {\it Numerical modeling of fluid-structure interaction in arteries with anisotropic polyconvex hyperelastic and anisotropic viscoelastic material models at finite strains}, Int. J. Numer. Methods Biomed. Eng., 32 (2016), e02756, . |
| [11] | D. Balzani, P. Neff, J. Schröder, and G. A. Holzapfel, {\it A polyconvex framework for soft biological tissues. Adjustment to experimental data}, Internat. J. Solids Struct., 43 (2006), pp. 6052-6070, . · Zbl 1120.74632 |
| [12] | D. Brands, A. Klawonn, O. Rheinbach, and J. Schröder, {\it Modelling and convergence in arterial wall simulations using a parallel FETI solution strategy}, Comput. Methods Biomech. Biomed. Engin., 11 (2008), pp. 569-583, . |
| [13] | M. Buck, O. Iliev, and H. Andrä, {\it Multiscale finite elements for linear elasticity: Oscillatory boundary conditions}, in Domain Decomposition Methods in Science and Engineering XXI, Lect. Notes Comput. Sci. Eng. 98, J. Erhal, M. J. Gander, L. Halpern, G. Pichot, T. Sassi, and O. Widlund, eds., Springer, Cham, Switzerland, 2014, pp. 237-245, . · Zbl 1455.74012 |
| [14] | P. Crosetto, S. Deparis, G. Fourestey, and A. Quarteroni, {\it Parallel algorithms for fluid-structure interaction problems in haemodynamics}, SIAM J. Sci. Comput., 33 (2011), pp. 1598-1622, . · Zbl 1417.92008 |
| [15] | T. A. Davis and I. S. Duff, {\it An unsymmetric-pattern multifrontal method for sparse LU factorization}, SIAM J. Matrix Anal. Appl., 18 (1997), pp. 140-158, . · Zbl 0884.65021 |
| [16] | S. Deparis, D. Forti, G. Grandperrin, and A. Quarteroni, {\it FaCSI: A block parallel preconditioner for fluid-structure interaction in hemodynamics}, J. Comput. Phys., 327 (2016), pp. 700-718, . · Zbl 1373.74036 |
| [17] | C. R. Dohrmann, A. Klawonn, and O. B. Widlund, {\it Domain decomposition for less regular subdomains: Overlapping Schwarz in two dimensions}, SIAM J. Numer. Anal., 46 (2008), pp. 2153-2168, . · Zbl 1183.65160 |
| [18] | C. R. Dohrmann, A. Klawonn, and O. B. Widlund, {\it A family of energy minimizing coarse spaces for overlapping Schwarz preconditioners}, in Domain Decomposition Methods in Science and Engineering XVII, Lect. Notes Comput. Sci. Eng. 60, Springer, Berlin, 2008, pp. 247-254, . · Zbl 1359.65289 |
| [19] | C. R. Dohrmann and O. B. Widlund, {\it Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity}, Internat. J. Numer. Methods Engrg., 82 (2010), pp. 157-183, . · Zbl 1188.74053 |
| [20] | C. R. Dohrmann and O. B. Widlund, {\it Lower dimensional coarse spaces for domain decomposition}, in Domain Decomposition Methods in Science and Engineering XXI, Lect. Notes Comput. Sci. Eng. 98, J. Erhel, J. M. Gander, L. Halpern, G. Pichot, T. Sassi, and O. Widlund, eds., Springer, Cham, Switzerland, 2014, pp. 527-535, . · Zbl 1382.65437 |
| [21] | M. Dryja, B. F. Smith, and O. B. Widlund, {\it Schwarz analysis of iterative substructuring algorithms for elliptic problems in three dimensions}, SIAM J. Numer. Anal., 31 (1994), pp. 1662-1694, . · Zbl 0818.65114 |
| [22] | M. Á. Fernández and M. Moubachir, {\it A Newton method using exact Jacobians for solving fluid-structure coupling}, Comput. & Structures, 83 (2005), pp. 127-142, . |
| [23] | L. Formaggia, M. Fernandez, A. Gauthier, J. F. Gerbeau, C. Prud’homme, and A. Veneziani, {\it The LifeV Project}, , 2015. |
| [24] | P. C. Franzone, L. F. Pavarino, and S. Scacchi, {\it Parallel multilevel solvers for the cardiac electro-mechanical coupling}, Appl. Numer. Math., 95 (2015), pp. 140-153, . · Zbl 1320.65138 |
| [25] | I. Guèye, S. E. Arem, F. Feyel, F.-X. Roux, and G. Cailletaud, {\it A new parallel sparse direct solver: Presentation and numerical experiments in large-scale structural mechanics parallel computing}, Internat. J. Numer. Methods Engrg., 88 (2011), pp. 370-384, . · Zbl 1242.74122 |
| [26] | A. Heinlein, {\it Parallel Overlapping Schwarz Preconditioners and Multiscale Discretizations with Applications to Fluid-Structure Interaction and Highly Heterogeneous Problems}, Ph.D. thesis, Universität zu Köln, Köln, Germany, 2016. |
| [27] | A. Heinlein, U. Hetmaniuk, A. Klawonn, and O. Rheinbach, {\it The approximate component mode synthesis special finite element method in two dimensions: Parallel implementation and numerical results}, J. Comput. Appl. Math., 289 (2015), pp. 116-133, . · Zbl 1326.65160 |
| [28] | A. Heinlein, A. Klawonn, J. Knepper, and O. Rheinbach, {\it Multiscale Coarse Spaces for Overlapping Schwarz Methods Based on the ACMS Space in \textup2D}, Tech. Rep., Fakultät für Mathematik und Informatik, Technische Universität Bergakademie Freiberg, Freiberg, Germany, 2016, submitted; available as TUBAF Preprint 2016-09 online from . · Zbl 1448.65263 |
| [29] | A. Heinlein, A. Klawonn, and O. Rheinbach, {\it Parallel overlapping Schwarz with an energy-minimizing coarse space}, in Proceedings of the 23rd International Conference on Domain Decomposition Methods, Lect. Notes Comput. Sci. Eng., Springer, Cham, Switzerland; to appear; available as TUBAF Preprint 2015-17 online from . · Zbl 1367.65043 |
| [30] | A. Heinlein, A. Klawonn, and O. Rheinbach, {\it Parallel two-level overlapping Schwarz methods in fluid-structure interaction}, in Proceedings of ENUMATH 2015, Lect. Notes Comput. Sci. Eng., Springer, Cham, Switzerland, to appear; available as TUBAF Preprint 2015-15 online from . · Zbl 1387.74108 |
| [31] | M. A. Heroux, R. A. Bartlett, V. E. Howle, R. J. Hoekstra, J. J. Hu, T. G. Kolda, R. B. Lehoucq, K. R. Long, R. P. Pawlowski, E. T. Phipps, A. G. Salinger, H. K. Thornquist, R. S. Tuminaro, J. M. Willenbring, A. Williams, and K. S. Stanley, {\it An overview of the Trilinos project}, ACM Trans. Math. Software, 31 (2005), pp. 397-423, . · Zbl 1136.65354 |
| [32] | P. Hüser, {\it GDSW - ein überlappendes Zwei-Level-Schwarz-Verfahren mit energieminimierendem groben Gitter}, Master’s thesis, Fakultät für Mathematik, Universität Duisburg-Essen, Duisburg, Germany, 2012. |
| [33] | P. Jolivet, {\it Méthodes de décomposition de domaine. Application au calcul haute performance}, Ph.D. thesis, Université de Grenoble, Grenoble, France, 2014. |
| [34] | P. Jolivet, F. Hecht, F. Nataf, and C. Prud’homme, {\it Scalable domain decomposition preconditioners for heterogeneous elliptic problems}, in Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis (SC ’13), ACM, New York, 2013, pp. 80:1-80:11, . |
| [35] | G. Karypis, K. Schloegel, and V. Kumar, {\it ParMETIS - Parallel Graph Partitioning and Sparse Matrix Ordering. Version 3.2}, Tech. Rep., Department of Computer Science and Engineering, University of Minnesota, Minneapolis, MN, 2011. |
| [36] | A. Klawonn, M. Lanser, and O. Rheinbach, {\it Toward extremely scalable nonlinear domain decomposition methods for elliptic partial differential equations}, SIAM J. Sci. Comput., 37 (2015), pp. C667-C696, . · Zbl 1329.65294 |
| [37] | A. Klawonn, M. Lanser, and O. Rheinbach, {\it \(FE^2\) TI: Computational scale bridging for dual-phase steels}, in Parallel Computing: On the Road to Exascale, Advances in Parallel Computing 27, IOS Press, Amsterdam, 2016, pp. 797-806; also available as TUBAF Preprint 2015-12 online from . |
| [38] | A. Klawonn, M. Lanser, O. Rheinbach, H. Stengel, and G. Wellein, {\it Hybrid MPI/OpenMP parallelization in FETI-DP methods}, in Recent Trends in Computational Engineering - CE2014, Lect. Notes Comput. Sci. Eng. 105, M. Mehl, M. Bischoff, and M. Schäfer, eds., Springer, Cham, Switzerland, 2015, pp. 67-84, . |
| [39] | A. Klawonn and O. Rheinbach, {\it Deflation, projector preconditioning, and balancing in iterative substructuring methods: Connections and new results}, SIAM J. Sci. Comput., 34 (2012), pp. A459-A484, . · Zbl 1248.65129 |
| [40] | F. Kong and X.-C. Cai, {\it A scalable Schwarz method for \textup3D linear elasticity problems on domains with complex geometry}, poster presented at The International Conference for High Performance Computing, Networking, Storage and Analysis (SC13), Denver, CO, 2013, . |
| [41] | J. Mandel and C. R. Dohrmann, {\it Convergence of a balancing domain decomposition by constraints and energy minimization}, Numer. Linear Algebra Appl., 10 (2003), pp. 639-659, . · Zbl 1071.65558 |
| [42] | M. Munteanu, L. F. Pavarino, and S. Scacchi, {\it A scalable Newton-Krylov-Schwarz method for the bidomain reaction-diffusion system}, SIAM J. Sci. Comput., 31 (2009), pp. 3861-3883, . · Zbl 1205.65261 |
| [43] | S. V. Patankar and D. B. Spalding, {\it A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows}, Int. J. Heat Mass Transfer, 15 (1972), pp. 1787-1806, . . · Zbl 0246.76080 |
| [44] | S. Scacchi, {\it A hybrid multilevel Schwarz method for the bidomain model}, Comput. Methods Appl. Mech. Engrg., 197 (2008), pp. 4051-4061, . · Zbl 1194.78048 |
| [45] | S. Scacchi, {\it A multilevel hybrid Newton-Krylov-Schwarz method for the bidomain model of electrocardiology}, Comput. Methods Appl. Mech. Engrg., 200 (2011), pp. 717-725, . · Zbl 1225.92011 |
| [46] | C. Siefert, developer of the Trilinos packages \tt ML, \tt MueLu, \tt Ifpack, Ifpack2, and EpetraExt, {\it private communication}, 2015. |
| [47] | M. Stephan and J. Docter, {\it JUQUEEN: IBM Blue Gene/Q Supercomputer System at the Jülich Supercomputing Centre}, J. Large-Scale Res. Facil., 1 (2015), A1, . |
| [48] | A. Toselli and O. Widlund, {\it Domain Decomposition Methods–Algorithms and Theory}, Springer Ser. Comput. Math. 34, Springer-Verlag, Berlin, 2005, . · Zbl 1069.65138 |
| [49] | Y. Wu and X.-C. Cai, {\it A fully implicit domain decomposition based ALE framework for three-dimensional fluid-structure interaction with application in blood flow computation}, J. Comput. Phys., 258 (2014), pp. 524-537, . · Zbl 1349.76552 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
