Skip to content

Latest commit

 

History

History
90 lines (54 loc) · 3.39 KB

README.md

File metadata and controls

90 lines (54 loc) · 3.39 KB
Raw

libCEED + PETSc Examples

PETSc v3.17 or a development version of PETSc at commit 0e95d842 or later is required.

CEED bakeoff problems with raw mesh management - bpsraw

This code solves the CEED bakeoff problems on a structured grid generated and referenced using only low-level communication primitives.

To build, run make bpsraw

To run, ./bpsraw -ceed [ceed-resource] -problem bp[1-6] -degree [degree]

In addition to the common arguments, the following arguments may be set:

  • -local - Target number of locally owned DoFs per process

CEED bakeoff problems with DMPlex - bps

This code solves the CEED bakeoff problems on a unstructured grid using DMPlex. This example requires a PETSc version later than 3.11.3.

To build, run make bps

To run, ./bps -ceed [ceed-resource] -problem bp[1-6] -degree [degree]

In addition to the common arguments, the following arguments may be set:

  • -mesh - Read mesh from file
  • -cells - Number of cells per dimension

Running a suite

Some run-time arguments can be passed lists, which allows a single mpiexec invocation to run many experiments. For example

mpiexec -n 64 ./bps -problem bp1,bp2,bp3,bp4 -degree 2,3,5,7                \
    -ceed /cpu/self/opt/serial,/cpu/self/xsmm/serial,/cpu/self/xsmm/blocked \
    -local_nodes 600,20000 | tee bps.log

which will sample from the 4*4*3=48 specified combinations, each of which will run a problem-size sweep of 600, 1200, 2400, 4800, 9600, 192000 FEM nodes per MPI rank. The resulting log file can be read by the Python plotting scripts in benchmarks/.

CEED bakeoff problems with DMPlex and PCMG - multigrid

This code solves the CEED bakeoff problems on a unstructured grid using DMPlex with p-multigrid implemented in PCMG. This example requires a PETSc version later than 3.11.3.

To build, run make multigrid

To run, ./multigrid -ceed [ceed-resource] -problem bp[1-6] -degree [degree]

In addition to the common arguments, the following arguments may be set:

  • -mesh - Read mesh from file
  • -cells - Number of cells per dimension

Command line arguments

The following arguments can be specified for all of the above examples:

  • -ceed - CEED resource specifier
  • -problem - CEED benchmark problem to solve
  • -degree - Polynomial degree of tensor product basis
  • -q_extra - Number of extra quadrature points
  • -test - Testing mode (do not print unless error is large)
  • -benchmark - Benchmarking mode (prints benchmark statistics)

libCEED example to compute surface area using DMPlex - area

This example uses the mass matrix to compute the surface area of a cube or a discrete cubed-sphere, defined via DMPlex.

To build, run make area

To run, ./area -problem cube -ceed [ceed-resource] -petscspace_degree [degree]

or

./area -problem sphere -ceed [ceed-resource] -petscspace_degree [degree]

Command line arguments

The following arguments can be specified for the area example:

  • -ceed - CEED resource specifier
  • -problem - Problem to solve, either 'cube' or 'sphere'
  • -petscspace_degree - Polynomial degree of tensor product basis
  • -q_extra - Number of extra quadrature points
  • -test - Testing mode (do not print unless error is large)
  • -mesh - Read mesh from file