Numerical Solvers#

The RatelSolverType determines how the composite CeedOperator are build and used to set the appropriate DMSNES or DMTS options.

PetscLogEvent RATEL_Prolong[RATEL_MAX_MULTIGRID_LEVELS]#
PetscLogEvent RATEL_Restrict[RATEL_MAX_MULTIGRID_LEVELS]#
static PetscErrorCode RatelRegisterPMultigridLogEvents(PetscInt num_multigrid_levels)#

Register p-multigrid Ratel log events.

Not collective across MPI processes.

Parameters:

  • num_multigrid_levels – [in] Number of multigrid levels to register

Returns:

An error code: 0 - success, otherwise - failure

static PetscErrorCode RatelPMGProcessCommandLineOptions(RatelPMGContext pmg, PetscInt pmg_field, const char *pc_prefix)#

Get PCpMG options.

Collective across MPI processes.

Parameters:

  • pmg – [inout] PCpMG context

  • pmg_field – [in] Index of field on the DM

  • pc_prefix – [in] Command line option prefix for PC

Returns:

An error code: 0 - success, otherwise - failure

static PetscErrorCode RatelSetupMultigridLevel(Ratel ratel, DM dm_level, Vec M_loc, PetscInt level, CeedOperator op_jacobian_fine, CeedOperator op_jacobian, CeedOperator op_prolong, CeedOperator op_restrict)#

Setup CeedOperator for multigrid prolongation, restriction, and coarse grid Jacobian evaluation.

Collective across MPI processes.

Parameters:

  • ratel – [in] Ratel context

  • dm_level – [in] DMPlex for multigrid level to setup

  • M_loc – [in] PETSc local vector holding multiplicity data

  • level – [in] Multigrid level to set up

  • op_jacobian_fine – [inout] Composite CeedOperator for Jacobian

  • op_jacobian – [inout] Composite CeedOperator for Jacobian

  • op_prolong – [inout] Composite CeedOperator for prolongation

  • op_restrict – [inout] Composite CeedOperator for restriction

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCPMGCreate(PC pc)#

Setup PCpMG preconditioner from Ratel context.

Collective across MPI processes.

Parameters:

  • pc – [inout] PCpMG object

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCSetUp_PMG(PC pc)#

SetUp for PCpMG.

Reassemble coarse operator.

Collective across MPI processes.

Parameters:

  • pc – [inout] PC object to setup

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCView_PMG(PC pc, PetscViewer viewer)#

View PCpMG.

Collective across MPI processes.

Parameters:

  • pc – [in] PC object to view

  • viewer – [inout] Visualization context

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCApply_PMG(PC pc, Vec X_in, Vec X_out)#

Apply PCpMG.

Collective across MPI processes.

Parameters:

  • pc – [in] PC object to apply

  • X_in – [in] Input vector

  • X_out – [out] Output vector

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCApplyTranspose_PMG(PC pc, Vec X_in, Vec X_out)#

Apply PCpMG transpose.

Collective across MPI processes.

Parameters:

  • pc – [in] PC object to apply

  • X_in – [in] Input vector

  • X_out – [out] Output vector

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCMatApply_PMG(PC pc, Mat X_in, Mat X_out)#

Apply PCpMG to multiple vectors stored as MATDENSE.

Collective across MPI processes.

Parameters:

  • pc – [in] PC object to apply

  • X_in – [in] Input matrix

  • X_out – [out] Output matrix

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCDestroy_PMG(PC pc)#

Destroy PCpMG context data.

Collective across MPI processes.

Parameters:

  • pc – [inout] PC object to destroy

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCPMGContextDestroy(RatelPMGContext pmg)#

Destroy Ratel PMG preconditioner.

Parameters:

  • pmg – [inout] Ratel pMG context

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelPCRegisterPMG(Ratel ratel)#

Register PCpMG preconditioner.

Not collective across MPI processes.

Parameters:

  • ratel – [in] Ratel context

Returns:

An error code: 0 - success, otherwise - failure

PetscErrorCode RatelDMSetupSolver(Ratel ratel)#

Setup solver contexts.

Collective across MPI processes.

Parameters:

  • ratel – [inout] Ratel context

Returns:

An error code: 0 - success, otherwise - failure