Parallel GeoFEST for regional faulted deformation
The GeoFEST application, which can accurately simulate deformation and stress in a faulted region of the earth's crust and mantle in the presence of heterogeneous rheology, faulting, tectonic loading, and viscoelastic relaxation, is a scalable parallel application well suited for the largest supercomputers. In order to achieve such scaling, the authors implemented a carefully choreographed algorithm using the Message Passing interface (MPI) calls, made use of parallel partitioning and data management functions from the PYRAMID and ParMetis libraries, and designed a technique for using the strain energy density from an initial scratch elastic solution to direct PYRAMID in performing parallel adaptive mesh refinement. An example problem follows that illustrates how an initial geometry and mesh designed on a laptop computer is automatically processed on hundreds to thousands of processors using millions of refined finite elements. Copyright (c) 2009 John Wiley & Sons, Ltd.
adaptive mesh refinement finite element earthquake