Interface for linear (KSP) solver. More...

#include <src/petsc/KspCtx.hpp>

Collaboration diagram for MoFEM::KspCtx:

Public Types
typedef MoFEM::PairNameFEMethodPtr	PairNameFEMethodPtr

typedef MoFEM::FEMethodsSequence	FEMethodsSequence

typedef MoFEM::BasicMethodsSequence	BasicMethodsSequence

Public Member Functions
	KspCtx (MoFEM::Interface &m_field, const std::string &_problem_name)

virtual	~KspCtx ()

FEMethodsSequence &	get_loops_to_do_Mat ()

FEMethodsSequence &	get_loops_to_do_Rhs ()

BasicMethodsSequence &	get_preProcess_to_do_Rhs ()

BasicMethodsSequence &	get_postProcess_to_do_Rhs ()

BasicMethodsSequence &	get_preProcess_to_do_Mat ()

BasicMethodsSequence &	get_postProcess_to_do_Mat ()

Public Attributes
DEPRECATED typedef MoFEM::PairNameFEMethodPtr	loop_pair_type

DEPRECATED typedef MoFEM::FEMethodsSequence	loops_to_do_type

DEPRECATED typedef MoFEM::BasicMethodsSequence	basic_method_to_do

MoFEM::Interface &	mField

moab::Interface &	moab

std::string	problemName
	Problem name. More...

MoFEMTypes	bH
	If set to MF_EXIST check if element exist. More...

FEMethodsSequence	loops_to_do_Mat

FEMethodsSequence	loops_to_do_Rhs

BasicMethodsSequence	preProcess_Mat

BasicMethodsSequence	postProcess_Mat

BasicMethodsSequence	preProcess_Rhs
	Sequence of methods run before residual is assembled. More...

BasicMethodsSequence	postProcess_Rhs
	Sequence of methods run after residual is assembled. More...

PetscLogEvent	MOFEM_EVENT_KspRhs

PetscLogEvent	MOFEM_EVENT_KspMat

Friends
PetscErrorCode	KspRhs (KSP ksp, Vec f, void *ctx)
	Run over elements in the lists. More...

PetscErrorCode	KspMat (KSP ksp, Mat A, Mat B, void *ctx)
	Run over elenents in the list. More...

Detailed Description

Interface for linear (KSP) solver.

Definition at line 27 of file KspCtx.hpp.

Member Typedef Documentation

◆ BasicMethodsSequence

typedef MoFEM::BasicMethodsSequence MoFEM::KspCtx::BasicMethodsSequence

Definition at line 46 of file KspCtx.hpp.

◆ FEMethodsSequence

typedef MoFEM::FEMethodsSequence MoFEM::KspCtx::FEMethodsSequence

Definition at line 45 of file KspCtx.hpp.

◆ PairNameFEMethodPtr

typedef MoFEM::PairNameFEMethodPtr MoFEM::KspCtx::PairNameFEMethodPtr

Definition at line 44 of file KspCtx.hpp.

Constructor & Destructor Documentation

◆ KspCtx()

MoFEM::KspCtx::KspCtx	(	MoFEM::Interface &	m_field,
		const std::string &	_problem_name
	)

Definition at line 64 of file KspCtx.hpp.

       : mField(m_field), moab(m_field.get_moab()), problemName(_problem_name),
         bH(MF_EXIST) {
     PetscLogEventRegister("LoopKSPRhs", 0, &MOFEM_EVENT_KspRhs);
     PetscLogEventRegister("LoopKSPMat", 0, &MOFEM_EVENT_KspMat);
   }

◆ ~KspCtx()

virtual MoFEM::KspCtx::~KspCtx ( )

virtual

Definition at line 70 of file KspCtx.hpp.

70 {}

Member Function Documentation

◆ get_loops_to_do_Mat()

FEMethodsSequence& MoFEM::KspCtx::get_loops_to_do_Mat ( )

Returns: return reference to vector with FEMethod to calculate matrix

Definition at line 75 of file KspCtx.hpp.

75 { return loops_to_do_Mat; }

MoFEM::KspCtx::loops_to_do_Mat

FEMethodsSequence loops_to_do_Mat

Definition: KspCtx.hpp:48

◆ get_loops_to_do_Rhs()

FEMethodsSequence& MoFEM::KspCtx::get_loops_to_do_Rhs ( )

Returns: return vector to vector with FEMethod to vector

Definition at line 80 of file KspCtx.hpp.

80 { return loops_to_do_Rhs; }

MoFEM::KspCtx::loops_to_do_Rhs

FEMethodsSequence loops_to_do_Rhs

Definition: KspCtx.hpp:50

◆ get_postProcess_to_do_Mat()

BasicMethodsSequence& MoFEM::KspCtx::get_postProcess_to_do_Mat ( )

The sequence of BasicMethod is executed after tangent matrix is calculated. It can be used to setup data structures, e.g. aggregate data from processors or to apply essential boundary conditions.

Returns: reference to BasicMethod for postprocessing

Definition at line 112 of file KspCtx.hpp.

112 { return postProcess_Mat; }

MoFEM::KspCtx::postProcess_Mat

BasicMethodsSequence postProcess_Mat

Definition: KspCtx.hpp:54

◆ get_postProcess_to_do_Rhs()

BasicMethodsSequence& MoFEM::KspCtx::get_postProcess_to_do_Rhs ( )

The sequence of BasicMethod is executed after residual is calculated. It can be used to setup data structures, e.g. aggregate data from processors or to apply essential boundary conditions.

Returns: reference to BasicMethod for postprocessing

Definition at line 98 of file KspCtx.hpp.

98 { return postProcess_Rhs; }

MoFEM::KspCtx::postProcess_Rhs

BasicMethodsSequence postProcess_Rhs

Sequence of methods run after residual is assembled.

Definition: KspCtx.hpp:59

◆ get_preProcess_to_do_Mat()

BasicMethodsSequence& MoFEM::KspCtx::get_preProcess_to_do_Mat ( )

Returns: reference to BasicMethod for preprocessing

Definition at line 103 of file KspCtx.hpp.

103 { return preProcess_Mat; }

MoFEM::KspCtx::preProcess_Mat

BasicMethodsSequence preProcess_Mat

Definition: KspCtx.hpp:52

◆ get_preProcess_to_do_Rhs()

BasicMethodsSequence& MoFEM::KspCtx::get_preProcess_to_do_Rhs ( )

The sequence of BasicMethod is executed before residual is calculated. It can be used to setup data structures, e.g. zero global variable which is integrated in domain, e.g. for calculation of strain energy.

Returns: reference to BasicMethod for preprocessing

Definition at line 89 of file KspCtx.hpp.

89 { return preProcess_Rhs; }

MoFEM::KspCtx::preProcess_Rhs

BasicMethodsSequence preProcess_Rhs

Sequence of methods run before residual is assembled.

Definition: KspCtx.hpp:57

Friends And Related Function Documentation

◆ KspMat

PetscErrorCode KspMat	(	KSP	ksp,
		Mat	A,
		Mat	B,
		void *	ctx
	)

friend

Run over elenents in the list.

Parameters

ksp	KSP solver
A	matrix
B	Preconditioner matrix
ctx	data context, i.e. KspCtx

Returns: error code

Definition at line 59 of file KspCtx.cpp.

                                                      {
   // PetscValidHeaderSpecific(ksp,KSP_CLASSID,1);
   MoFEMFunctionBeginHot;
   KspCtx* ksp_ctx = (KspCtx*)ctx;
   PetscLogEventBegin(ksp_ctx->MOFEM_EVENT_KspMat,0,0,0,0);
   KspCtx::BasicMethodsSequence::iterator bit = ksp_ctx->preProcess_Mat.begin();
   for(;bit!=ksp_ctx->preProcess_Mat.end();bit++) {
     ierr = (*bit)->setKsp(ksp); CHKERRG(ierr);
     (*bit)->ksp_A = A;
     (*bit)->ksp_B = B;
     ierr = (*bit)->setKspCtx(KspMethod::CTX_OPERATORS); CHKERRG(ierr);
     ierr = ksp_ctx->mField.problem_basic_method_preProcess(ksp_ctx->problemName,*(*(bit)));  CHKERRG(ierr);
     ierr = (*bit)->setKspCtx(KspMethod::CTX_KSPNONE);  CHKERRG(ierr);
   }
   KspCtx::FEMethodsSequence::iterator lit = ksp_ctx->loops_to_do_Mat.begin();
   for(;lit!=ksp_ctx->loops_to_do_Mat.end();lit++) {
     lit->second->ksp_A = A;
     lit->second->ksp_B = B;
     ierr = lit->second->setKsp(ksp); CHKERRG(ierr);
     ierr = lit->second->setKspCtx(KspMethod::CTX_OPERATORS); CHKERRG(ierr);
     ierr = ksp_ctx->mField.loop_finite_elements(ksp_ctx->problemName,lit->first,*(lit->second),ksp_ctx->bH);  CHKERRG(ierr);
     ierr = lit->second->setKspCtx(KspMethod::CTX_KSPNONE);
   }
   bit = ksp_ctx->postProcess_Mat.begin();
   for(;bit!=ksp_ctx->postProcess_Mat.end();bit++) {
     ierr = (*bit)->setKsp(ksp); CHKERRG(ierr);
     (*bit)->ksp_A = A;
     (*bit)->ksp_B = B;
     ierr = (*bit)->setKspCtx(KspMethod::CTX_OPERATORS); CHKERRG(ierr);
     ierr = ksp_ctx->mField.problem_basic_method_postProcess(ksp_ctx->problemName,*(*(bit)));  CHKERRG(ierr);
     ierr = (*bit)->setKspCtx(KspMethod::CTX_KSPNONE);  CHKERRG(ierr);
   }
   ierr = MatAssemblyBegin(B,MAT_FINAL_ASSEMBLY); CHKERRG(ierr);
   ierr = MatAssemblyEnd(B,MAT_FINAL_ASSEMBLY); CHKERRG(ierr);
   // MatView(A,PETSC_VIEWER_DRAW_WORLD);
   // std::string wait;
   // std::cin >> wait;
   PetscLogEventEnd(ksp_ctx->MOFEM_EVENT_KspMat,0,0,0,0);
   MoFEMFunctionReturnHot(0);
 }

◆ KspRhs

PetscErrorCode KspRhs	(	KSP	ksp,
		Vec	f,
		void *	ctx
	)

friend

Run over elements in the lists.

Parameters

ksp	KSP solver
f	the right hand side vector
ctx	data context, i.e. KspCtx

Returns: error code

Definition at line 23 of file KspCtx.cpp.

                                                {
   // PetscValidHeaderSpecific(ksp,KSP_CLASSID,1);
   MoFEMFunctionBeginHot;
   KspCtx* ksp_ctx = (KspCtx*)ctx;
   PetscLogEventBegin(ksp_ctx->MOFEM_EVENT_KspRhs,0,0,0,0);
   KspCtx::BasicMethodsSequence::iterator bit = ksp_ctx->preProcess_Rhs.begin();
   for(;bit!=ksp_ctx->preProcess_Rhs.end();bit++) {
     ierr = (*bit)->setKsp(ksp); CHKERRG(ierr);
     (*bit)->ksp_f = f;
     ierr = (*bit)->setKspCtx(KspMethod::CTX_SETFUNCTION);  CHKERRG(ierr);
     ierr = ksp_ctx->mField.problem_basic_method_preProcess(ksp_ctx->problemName,*(*(bit)));  CHKERRG(ierr);
     ierr = (*bit)->setKspCtx(KspMethod::CTX_KSPNONE);  CHKERRG(ierr);
   }
   KspCtx::FEMethodsSequence::iterator lit = ksp_ctx->loops_to_do_Rhs.begin();
   for(;lit!=ksp_ctx->loops_to_do_Rhs.end();lit++) {
     ierr = lit->second->setKspCtx(KspMethod::CTX_SETFUNCTION);  CHKERRG(ierr);
     ierr = lit->second->setKsp(ksp); CHKERRG(ierr);
     lit->second->ksp_f = f;
     ierr = ksp_ctx->mField.loop_finite_elements(ksp_ctx->problemName,lit->first,*(lit->second),ksp_ctx->bH);  CHKERRG(ierr);
     ierr = lit->second->setKspCtx(KspMethod::CTX_KSPNONE);  CHKERRG(ierr);
   }
   bit = ksp_ctx->postProcess_Rhs.begin();
   for(;bit!=ksp_ctx->postProcess_Rhs.end();bit++) {
     ierr = (*bit)->setKsp(ksp); CHKERRG(ierr);
     (*bit)->ksp_f = f;
     ierr = (*bit)->setKspCtx(KspMethod::CTX_SETFUNCTION);  CHKERRG(ierr);
     ierr = ksp_ctx->mField.problem_basic_method_postProcess(ksp_ctx->problemName,*(*(bit)));  CHKERRG(ierr);
     ierr = (*bit)->setKspCtx(KspMethod::CTX_KSPNONE);  CHKERRG(ierr);
   }
   ierr = VecGhostUpdateBegin(f,ADD_VALUES,SCATTER_REVERSE); CHKERRG(ierr);
   ierr = VecGhostUpdateEnd(f,ADD_VALUES,SCATTER_REVERSE); CHKERRG(ierr);
   ierr = VecAssemblyBegin(f); CHKERRG(ierr);
   ierr = VecAssemblyEnd(f); CHKERRG(ierr);
   PetscLogEventEnd(ksp_ctx->MOFEM_EVENT_KspRhs,0,0,0,0);
   MoFEMFunctionReturnHot(0);
 }

Member Data Documentation

◆ basic_method_to_do

DEPRECATED typedef MoFEM::BasicMethodsSequence MoFEM::KspCtx::basic_method_to_do

Deprecated:: use BasicMethodsSequence

Definition at line 42 of file KspCtx.hpp.

◆ bH

MoFEMTypes MoFEM::KspCtx::bH

If set to MF_EXIST check if element exist.

Definition at line 33 of file KspCtx.hpp.

◆ loop_pair_type

DEPRECATED typedef MoFEM::PairNameFEMethodPtr MoFEM::KspCtx::loop_pair_type

Deprecated:: use PairNameFEMethodPtr

Definition at line 36 of file KspCtx.hpp.

◆ loops_to_do_Mat

FEMethodsSequence MoFEM::KspCtx::loops_to_do_Mat

Sequence of finite elements instances assembling tangent matrix

Definition at line 48 of file KspCtx.hpp.

◆ loops_to_do_Rhs

FEMethodsSequence MoFEM::KspCtx::loops_to_do_Rhs

Sequence of finite elements instances assembling residual vector

Definition at line 50 of file KspCtx.hpp.

◆ loops_to_do_type

DEPRECATED typedef MoFEM::FEMethodsSequence MoFEM::KspCtx::loops_to_do_type

Deprecated:: use FEMethodsSequence

Definition at line 39 of file KspCtx.hpp.

◆ mField

MoFEM::Interface& MoFEM::KspCtx::mField

Definition at line 29 of file KspCtx.hpp.

◆ moab

moab::Interface& MoFEM::KspCtx::moab

Definition at line 30 of file KspCtx.hpp.

◆ MOFEM_EVENT_KspMat

PetscLogEvent MoFEM::KspCtx::MOFEM_EVENT_KspMat

Definition at line 62 of file KspCtx.hpp.

◆ MOFEM_EVENT_KspRhs

PetscLogEvent MoFEM::KspCtx::MOFEM_EVENT_KspRhs

Definition at line 61 of file KspCtx.hpp.

◆ postProcess_Mat

BasicMethodsSequence MoFEM::KspCtx::postProcess_Mat

Sequence of methods run after tangent matrix is assembled

Definition at line 54 of file KspCtx.hpp.

◆ postProcess_Rhs

BasicMethodsSequence MoFEM::KspCtx::postProcess_Rhs

Sequence of methods run after residual is assembled.

Definition at line 59 of file KspCtx.hpp.

◆ preProcess_Mat

BasicMethodsSequence MoFEM::KspCtx::preProcess_Mat

Sequence of methods run before tangent matrix is assembled

Definition at line 52 of file KspCtx.hpp.

◆ preProcess_Rhs

BasicMethodsSequence MoFEM::KspCtx::preProcess_Rhs