#include <PCBlockDiagonal.hpp>
Public Member Functions | |
PCBlockDiagonal (KSP &rKspObject) | |
Public Attributes | |
PCBlockDiagonalContext | mPCContext |
PC | mPetscPCObject |
Private Member Functions | |
void | PCBlockDiagonalCreate (KSP &rKspObject) |
void | PCBlockDiagonalSetUp () |
Classes | |
struct | PCBlockDiagonalContext |
Let A be a matrix arising in the FEM discretisation of the bidomain equations with the following block structure:
A = (A11 B') (B A22)
By creating an instance of this class, one will define the following preconditioner:
inv(M) = inv( (A11 0) = (inv(A11) 0) (0 A22) ) (0 inv(A22))
The inverses are approximate with one cycle of AMG.
Note: This class requires PETSc to be build including HYPRE library. If it's not available, it will throw the following error:
[0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Unknown type. Check for miss-spelling or missing external package needed for type! [0]PETSC ERROR: Unable to find requested PC type hypre!
and will approximate the inverse of the subblocks with PETSc's default preconditioner (bjacobi at the time of writing this).
Definition at line 79 of file PCBlockDiagonal.hpp.
PCBlockDiagonal::PCBlockDiagonal | ( | KSP & | rKspObject | ) |
Constructor
rKspObject | KSP object where we want to install the block diagonal preconditioner. |
Definition at line 35 of file PCBlockDiagonal.cpp.
References mPCContext, PCBlockDiagonalCreate(), and PCBlockDiagonalSetUp().
void PCBlockDiagonal::PCBlockDiagonalCreate | ( | KSP & | rKspObject | ) | [private] |
Creates all the state data required by the preconditioner
rKspObject | KSP object where we want to install the block diagonal preconditioner. |
Definition at line 77 of file PCBlockDiagonal.cpp.
References PCBlockDiagonal::PCBlockDiagonalContext::A11_matrix_subblock, PCBlockDiagonal::PCBlockDiagonalContext::A11_scatter_ctx, PCBlockDiagonal::PCBlockDiagonalContext::A22_matrix_subblock, PCBlockDiagonal::PCBlockDiagonalContext::A22_scatter_ctx, PetscTools::CreateVec(), mPCContext, mPetscPCObject, PetscVecTools::SetupInterleavedVectorScatterGather(), TERMINATE, PCBlockDiagonal::PCBlockDiagonalContext::x1_subvector, PCBlockDiagonal::PCBlockDiagonalContext::x2_subvector, PCBlockDiagonal::PCBlockDiagonalContext::y1_subvector, and PCBlockDiagonal::PCBlockDiagonalContext::y2_subvector.
Referenced by PCBlockDiagonal().
void PCBlockDiagonal::PCBlockDiagonalSetUp | ( | ) | [private] |
Setups preconditioner
Definition at line 181 of file PCBlockDiagonal.cpp.
References PCBlockDiagonal::PCBlockDiagonalContext::A11_matrix_subblock, PCBlockDiagonal::PCBlockDiagonalContext::A22_matrix_subblock, mPCContext, PCBlockDiagonal::PCBlockDiagonalContext::PC_amg_A11, and PCBlockDiagonal::PCBlockDiagonalContext::PC_amg_A22.
Referenced by PCBlockDiagonal().
PC context, this will be passed to PCBlockDiagonalApply when PETSc returns control to our preconditioner subroutine. See PCShellSetContext().
Definition at line 105 of file PCBlockDiagonal.hpp.
Referenced by PCBlockDiagonal(), PCBlockDiagonalCreate(), and PCBlockDiagonalSetUp().
Generic PETSc preconditioner object
Definition at line 106 of file PCBlockDiagonal.hpp.
Referenced by PCBlockDiagonalCreate().