Conman: vectorizing a finite element code for incompressible two-dimensional convection in the Earth's mantle

https://doi.org/10.1016/0031-9201(90)90225-MGet rights and content

Abstract

We discuss some simple concepts for vectorizing scientific codes, then apply these concepts to ConMan, a finite element code for simulations of mantle convection. We demonstrate that large speed-ups, close to the theoretical limit of the machine, are possible for entire codes, not just specially constructed routines. Although our specific code uses the finite element method, the vectorizing concepts discussed are widely applicable.

References (12)

  • D.S. Malkus et al.

    Mixed finite element methods-reduced and selective integration techniques: a unification of concepts

    Comput. Meth. Appl. Mech. Eng.

    (1978)
  • C.C. Ashcraft et al.

    Progress in sparse matrix methods for large linear systems on vector supercomputers

    Int. J. Supercomput. Appl.

    (1987)
  • A. Brooks

    A Petrov-Galerkin Finite-element Formulation for Convection Dominated Flows

  • J.J. Dongarra

    Performance of various computers using standard linear equation software in a FORTRAN environment

    Argonne National Laboratory Tech. Memo

    (1987)
  • J.J. Dongarra et al.

    Squeezing the most out of an algorithm in Cray FORTRAN

    ACM Trans. Math. Software

    (1984)
  • T.J.R. Hughes

    The Finite Element Method

There are more references available in the full text version of this article.

Cited by (154)

  • Tectono-Convective Modes on Earth and Other Terrestrial Bodies

    2023, Dynamics of Plate Tectonics and Mantle Convection
  • Implementation of a Volume-of-Fluid method in a finite element code with applications to thermochemical convection in a density stratified fluid in the Earth's mantle

    2019, Computers and Fluids
    Citation Excerpt :

    In particular, there have been a sequence of codes developed over this period of time that are now freely available to any individual who wishes to study mantle dynamics. They include HC [25,26,71], ConMan [38], CitCom S [48,75,85], Citcom CU [51,84] and ASPECT [29,41]. These codes, as well as others, can be downloaded from the Computational Infrastructure for Geodynamics (CIG) at U.C. Davis.1

  • Effects of temporal plume-slab interaction on the partial melting of the subducted oceanic crust

    2015, Journal of Asian Earth Sciences
    Citation Excerpt :

    All of the model domains consist of four-node quadrilateral elements. The finite element code ConMan (King et␣al., 1990) was used to calculate the governing equations and rheology, described above. To minimize the effect of the initial conditions on the thermal and flow structures of the subduction zone, we allowed subduction for 30 My using constant subduction parameters (near-steady state), and injection of the plume blob into the mantle wedge was followed.

View all citing articles on Scopus

Present address: Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A.

View full text