structure for projection matrices More...

#include "src/finite_elements/ConstrainMatrixCtx.hpp"

Collaboration diagram for MoFEM::ConstrainMatrixCtx:

Public Member Functions
	ConstrainMatrixCtx (MoFEM::Interface &m_field, string x_problem, string y_problem, bool create_ksp=true, bool own_contrain_matrix=false)

	ConstrainMatrixCtx (MoFEM::Interface &m_field, VecScatter scatter, bool create_ksp=true, bool own_contrain_matrix=false)

virtual	~ConstrainMatrixCtx ()

MoFEMErrorCode	initializeQorP (Vec x)
	initialize vectors and matrices for Q and P shell matrices, scattering is set based on x_problem and y_problem

MoFEMErrorCode	initializeQTKQ ()
	initialize vectors and matrices for CTC+QTKQ shell matrices, scattering is set based on x_problem and y_problem

MoFEMErrorCode	recalculateCTandCCT ()
	re-calculate CT and CCT if C matrix has been changed since initialization

MoFEMErrorCode	recalculateCTC ()
	re-calculate CTC matrix has been changed since initialization

MoFEMErrorCode	destroyQorP ()
	destroy sub-matrices used for shell matrices P, Q, R, RT

MoFEMErrorCode	destroyQTKQ ()
	destroy sub-matrices used for shell matrix QTKQ

Public Attributes
MoFEM::Interface &	mField

KSP	kSP

Mat	C

Mat	CT

Mat	CCT

Mat	CTC

Mat	K

Vec	Cx

Vec	CCTm1_Cx

Vec	CT_CCTm1_Cx

Vec	CTCx

Vec	X

Vec	Qx

Vec	KQx

bool	initQorP

bool	initQTKQ

bool	createKSP

bool	createScatter

bool	cancelKSPMonitor

bool	ownConstrainMatrix

VecScatter	sCatter

string	xProblem

string	yProblem

PetscLogEvent	MOFEM_EVENT_projInit

PetscLogEvent	MOFEM_EVENT_projQ

PetscLogEvent	MOFEM_EVENT_projP

PetscLogEvent	MOFEM_EVENT_projR

PetscLogEvent	MOFEM_EVENT_projRT

PetscLogEvent	MOFEM_EVENT_projCTC_QTKQ

PetscReal	rTol

PetscReal	absTol

PetscReal	dTol

PetscInt	maxIts

Friends
MoFEMErrorCode	ProjectionMatrixMultOpQ (Mat Q, Vec x, Vec f)
	Multiplication operator for Q = I-CTC(CCT)^-1C.

MoFEMErrorCode	ConstrainMatrixMultOpP (Mat P, Vec x, Vec f)
	Multiplication operator for P = CT(CCT)^-1C.

MoFEMErrorCode	ConstrainMatrixMultOpR (Mat R, Vec x, Vec f)
	Multiplication operator for R = CT(CCT)^-1.

MoFEMErrorCode	ConstrainMatrixMultOpRT (Mat RT, Vec x, Vec f)
	Multiplication operator for RT = (CCT)^-TC.

MoFEMErrorCode	ConstrainMatrixMultOpCTC_QTKQ (Mat CTC_QTKQ, Vec x, Vec f)
	Multiplication operator for RT = (CCT)^-TC.

MoFEMErrorCode	ConstrainMatrixDestroyOpPorQ ()

MoFEMErrorCode	ConstrainMatrixDestroyOpQTKQ ()

Detailed Description

structure for projection matrices

Definition at line 18 of file ConstrainMatrixCtx.hpp.

Constructor & Destructor Documentation

◆ ConstrainMatrixCtx() [1/2]

MoFEM::ConstrainMatrixCtx::ConstrainMatrixCtx	(	MoFEM::Interface &	m_field,
		string	x_problem,
		string	y_problem,
		bool	create_ksp = true,
		bool	own_contrain_matrix = false )

Construct data structure to build operators for projection matrices

User need to set matrix C to make it work

Parameters

x_problem	problem on which vector is projected
y_problem	problem used to construct projection matrices
create_ksp	create ksp solver otherwise user need to set it up

Definition at line 19 of file ConstrainMatrixCtx.cpp.

    : mField(m_field), INIT_DATA_CONSTRAINMATRIXCTX, xProblem(x_problem),
      yProblem(y_problem) {
  PetscLogEventRegister("ProjectionInit", 0, &MOFEM_EVENT_projInit);
  PetscLogEventRegister("ProjectionQ", 0, &MOFEM_EVENT_projQ);
  PetscLogEventRegister("ProjectionP", 0, &MOFEM_EVENT_projP);
  PetscLogEventRegister("ProjectionR", 0, &MOFEM_EVENT_projR);
  PetscLogEventRegister("ProjectionRT", 0, &MOFEM_EVENT_projRT);
  PetscLogEventRegister("ProjectionCTC_QTKQ", 0, &MOFEM_EVENT_projCTC_QTKQ);
}

◆ ConstrainMatrixCtx() [2/2]

MoFEM::ConstrainMatrixCtx::ConstrainMatrixCtx	(	MoFEM::Interface &	m_field,
		VecScatter	scatter,
		bool	create_ksp = true,
		bool	own_contrain_matrix = false )

Definition at line 33 of file ConstrainMatrixCtx.cpp.

    : mField(m_field), INIT_DATA_CONSTRAINMATRIXCTX, sCatter(scatter) {
  PetscLogEventRegister("ProjectionInit", 0, &MOFEM_EVENT_projInit);
  PetscLogEventRegister("ProjectionQ", 0, &MOFEM_EVENT_projQ);
  PetscLogEventRegister("ProjectionP", 0, &MOFEM_EVENT_projP);
  PetscLogEventRegister("ProjectionR", 0, &MOFEM_EVENT_projR);
  PetscLogEventRegister("ProjectionRT", 0, &MOFEM_EVENT_projRT);
  PetscLogEventRegister("ProjectionCTC_QTKQ", 0, &MOFEM_EVENT_projCTC_QTKQ);
}

◆ ~ConstrainMatrixCtx()

virtual MoFEM::ConstrainMatrixCtx::~ConstrainMatrixCtx ( )

inlinevirtual

Definition at line 61 of file ConstrainMatrixCtx.hpp.

                                {
    ierr = destroyQorP();
    CHKERRABORT(mField.get_comm(), ierr);
    ierr = destroyQTKQ();
    CHKERRABORT(mField.get_comm(), ierr);
    if (ownConstrainMatrix) {
      ierr = MatDestroy(&C);
      CHKERRABORT(mField.get_comm(), ierr);
    }
  };

Member Function Documentation

◆ destroyQorP()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::destroyQorP ( )

destroy sub-matrices used for shell matrices P, Q, R, RT

Definition at line 94 of file ConstrainMatrixCtx.cpp.

                                               {
  MoFEMFunctionBegin;
  if (initQorP)
    MoFEMFunctionReturnHot(0);
  CHKERR MatDestroy(&CT);
  CHKERR MatDestroy(&CCT);
  if (createKSP) {
    CHKERR KSPDestroy(&kSP);
  }
  CHKERR VecDestroy(&X);
  CHKERR VecDestroy(&Cx);
  CHKERR VecDestroy(&CCTm1_Cx);
  CHKERR VecDestroy(&CT_CCTm1_Cx);
  if (createScatter) {
    CHKERR VecScatterDestroy(&sCatter);
  }
  initQorP = true;
  MoFEMFunctionReturn(0);
}

◆ destroyQTKQ()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::destroyQTKQ ( )

destroy sub-matrices used for shell matrix QTKQ

Definition at line 151 of file ConstrainMatrixCtx.cpp.

                                               {
  MoFEMFunctionBegin;
  if (initQTKQ)
    MoFEMFunctionReturnHot(0);
  CHKERR MatDestroy(&CTC);
  CHKERR VecDestroy(&Qx);
  CHKERR VecDestroy(&KQx);
  CHKERR VecDestroy(&CTCx);
  initQTKQ = true;
  MoFEMFunctionReturn(0);
}

◆ initializeQorP()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::initializeQorP ( Vec x )

initialize vectors and matrices for Q and P shell matrices, scattering is set based on x_problem and y_problem

Parameters

x	is a vector from problem x

Definition at line 45 of file ConstrainMatrixCtx.cpp.

                                                       {
  MoFEMFunctionBegin;
  if (initQorP) {
    initQorP = false;
 
    PetscLogEventBegin(MOFEM_EVENT_projInit, 0, 0, 0, 0);
    CHKERR MatTranspose(C, MAT_INITIAL_MATRIX, &CT);
    // need to be calculated when C is changed
    CHKERR MatTransposeMatMult(CT, CT, MAT_INITIAL_MATRIX, PETSC_DEFAULT, &CCT);
    if (createKSP) {
      CHKERR KSPCreate(mField.get_comm(), &kSP);
      // need to be recalculated when C is changed
      CHKERR KSPSetOperators(kSP, CCT, CCT);
      CHKERR KSPSetFromOptions(kSP);
      CHKERR KSPSetInitialGuessKnoll(kSP, PETSC_TRUE);
      CHKERR KSPGetTolerances(kSP, &rTol, &absTol, &dTol, &maxIts);
      CHKERR KSPSetUp(kSP);
      if (cancelKSPMonitor) {
        CHKERR KSPMonitorCancel(kSP);
      }
    }
#if PETSC_VERSION_GE(3, 5, 3)
    CHKERR MatCreateVecs(C, &X, PETSC_NULLPTR);
    CHKERR MatCreateVecs(C, PETSC_NULLPTR, &Cx);
    CHKERR MatCreateVecs(CCT, PETSC_NULLPTR, &CCTm1_Cx);
#else
    CHKERR MatGetVecs(C, &X, PETSC_NULLPTR);
    CHKERR MatGetVecs(C, PETSC_NULLPTR, &Cx);
    CHKERR MatGetVecs(CCT, PETSC_NULLPTR, &CCTm1_Cx);
#endif
    CHKERR VecDuplicate(X, &CT_CCTm1_Cx);
    if (createScatter) {
      CHKERR mField.getInterface<VecManager>()->vecScatterCreate(
          x, xProblem, ROW, X, yProblem, COL, &sCatter);
    }
    PetscLogEventEnd(MOFEM_EVENT_projInit, 0, 0, 0, 0);
  }
  MoFEMFunctionReturn(0);
}

◆ initializeQTKQ()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::initializeQTKQ ( )

initialize vectors and matrices for CTC+QTKQ shell matrices, scattering is set based on x_problem and y_problem

Definition at line 114 of file ConstrainMatrixCtx.cpp.

                                                  {
  MoFEMFunctionBegin;
  if (initQTKQ) {
    initQTKQ = false;
    PetscLogEventBegin(MOFEM_EVENT_projInit, 0, 0, 0, 0);
    // need to be recalculated when C is changed
    CHKERR MatTransposeMatMult(C, C, MAT_INITIAL_MATRIX, PETSC_DEFAULT, &CTC);
    if (debug) {
      // MatView(CCT,PETSC_VIEWER_DRAW_WORLD);
      int m, n;
      MatGetSize(CCT, &m, &n);
      PetscPrintf(mField.get_comm(), "CCT size (%d,%d)\n", m, n);
      // std::string wait;
      // std::cin >> wait;
    }
#if PETSC_VERSION_GE(3, 5, 3)
    CHKERR MatCreateVecs(K, &Qx, PETSC_NULLPTR);
    CHKERR MatCreateVecs(K, PETSC_NULLPTR, &KQx);
    CHKERR MatCreateVecs(CTC, PETSC_NULLPTR, &CTCx);
#else
    CHKERR MatGetVecs(K, &Qx, PETSC_NULLPTR);
    CHKERR MatGetVecs(K, PETSC_NULLPTR, &KQx);
    CHKERR MatGetVecs(CTC, PETSC_NULLPTR, &CTCx);
#endif
    PetscLogEventEnd(MOFEM_EVENT_projInit, 0, 0, 0, 0);
  }
  MoFEMFunctionReturn(0);
}

◆ recalculateCTandCCT()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::recalculateCTandCCT ( )

re-calculate CT and CCT if C matrix has been changed since initialization

Definition at line 85 of file ConstrainMatrixCtx.cpp.

                                                       {
  MoFEMFunctionBegin;
  if (initQorP)
    MoFEMFunctionReturnHot(0);
  CHKERR MatTranspose(C, MAT_REUSE_MATRIX, &CT);
  CHKERR MatTransposeMatMult(CT, CT, MAT_REUSE_MATRIX, PETSC_DEFAULT, &CCT);
  MoFEMFunctionReturn(0);
}

◆ recalculateCTC()

MoFEMErrorCode MoFEM::ConstrainMatrixCtx::recalculateCTC ( )

re-calculate CTC matrix has been changed since initialization

Definition at line 143 of file ConstrainMatrixCtx.cpp.

                                                  {
  MoFEMFunctionBegin;
  if (initQTKQ)
    MoFEMFunctionReturnHot(0);
  CHKERR MatTransposeMatMult(C, C, MAT_REUSE_MATRIX, PETSC_DEFAULT, &CTC);
  MoFEMFunctionReturn(0);
}

Friends And Related Symbol Documentation

◆ ConstrainMatrixDestroyOpPorQ

MoFEMErrorCode ConstrainMatrixDestroyOpPorQ ( )

friend

◆ ConstrainMatrixDestroyOpQTKQ

MoFEMErrorCode ConstrainMatrixDestroyOpQTKQ ( )

friend

◆ ConstrainMatrixMultOpCTC_QTKQ

MoFEMErrorCode ConstrainMatrixMultOpCTC_QTKQ	(	Mat	CTC_QTKQ,
		Vec	x,
		Vec	f )

friend

Multiplication operator for RT = (CCT)^-TC.

  Mat CTC_QTKQ; //for problem
  ConstrainMatrixCtx
projection_matrix_ctx(m_field,problem_name,contrains_problem_name);
  CHKERR
MatCreateShell(PETSC_COMM_WORLD,m,m,M,M,&projection_matrix_ctx,&CTC_QTKQ);
  CHKERR
MatShellSetOperation(CTC_QTKQ,MATOP_MULT,(void(*)(void))ConstrainMatrixMultOpCTC_QTKQ);
  CHKERR
MatShellSetOperation(CTC_QTKQ,MATOP_DESTROY,(void(*)(void))ConstrainMatrixDestroyOpQTKQ);

Definition at line 261 of file ConstrainMatrixCtx.cpp.

                                                                         {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(CTC_QTKQ, &void_ctx);
  ConstrainMatrixCtx *ctx = (ConstrainMatrixCtx *)void_ctx;
  PetscLogEventBegin(ctx->MOFEM_EVENT_projCTC_QTKQ, 0, 0, 0, 0);
  Mat Q;
  int M, N, m, n;
  CHKERR MatGetSize(ctx->K, &M, &N);
  CHKERR MatGetLocalSize(ctx->K, &m, &n);
  CHKERR MatCreateShell(ctx->mField.get_comm(), m, n, M, N, ctx, &Q);
  CHKERR MatShellSetOperation(Q, MATOP_MULT,
                              (void (*)(void))ProjectionMatrixMultOpQ);
  CHKERR ctx->initializeQTKQ();
  CHKERR MatMult(Q, x, ctx->Qx);
  CHKERR MatMult(ctx->K, ctx->Qx, ctx->KQx);
  CHKERR MatMult(Q, ctx->KQx, f);
  CHKERR VecScatterBegin(ctx->sCatter, x, ctx->X, INSERT_VALUES,
                         SCATTER_FORWARD);
  CHKERR VecScatterEnd(ctx->sCatter, x, ctx->X, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR MatMult(ctx->CTC, ctx->X, ctx->CTCx);
  CHKERR VecScatterBegin(ctx->sCatter, ctx->CTCx, f, ADD_VALUES,
                         SCATTER_REVERSE);
  CHKERR VecScatterEnd(ctx->sCatter, ctx->CTCx, f, ADD_VALUES, SCATTER_REVERSE);
  CHKERR MatDestroy(&Q);
  PetscLogEventEnd(ctx->MOFEM_EVENT_projCTC_QTKQ, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

◆ ConstrainMatrixMultOpP

MoFEMErrorCode ConstrainMatrixMultOpP	(	Mat	P,
		Vec	x,
		Vec	f )

friend

Multiplication operator for P = CT(CCT)^-1C.

  Mat P; //for problem
  ConstrainMatrixCtx
projection_matrix_ctx(m_field,problem_name,contrains_problem_name);
  CHKERR MatCreateShell(PETSC_COMM_WORLD,m,m,M,M,&projection_matrix_ctx,&P);
  CHKERR
MatShellSetOperation(P,MATOP_MULT,(void(*)(void))ConstrainMatrixMultOpP);

Definition at line 197 of file ConstrainMatrixCtx.cpp.

                                                           {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(P, &void_ctx);
  ConstrainMatrixCtx *ctx = (ConstrainMatrixCtx *)void_ctx;
  PetscLogEventBegin(ctx->MOFEM_EVENT_projP, 0, 0, 0, 0);
  CHKERR ctx->initializeQorP(x);
  CHKERR VecScatterBegin(ctx->sCatter, x, ctx->X, INSERT_VALUES,
                         SCATTER_FORWARD);
  CHKERR VecScatterEnd(ctx->sCatter, x, ctx->X, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR MatMult(ctx->C, ctx->X, ctx->Cx);
  CHKERR KSPSolve(ctx->kSP, ctx->Cx, ctx->CCTm1_Cx);
  CHKERR MatMult(ctx->CT, ctx->CCTm1_Cx, ctx->CT_CCTm1_Cx);
  CHKERR VecZeroEntries(f);
  CHKERR VecGhostUpdateBegin(f, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(f, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecScatterBegin(ctx->sCatter, ctx->CT_CCTm1_Cx, f, INSERT_VALUES,
                         SCATTER_REVERSE);
  CHKERR VecScatterEnd(ctx->sCatter, ctx->CT_CCTm1_Cx, f, INSERT_VALUES,
                       SCATTER_REVERSE);
  PetscLogEventEnd(ctx->MOFEM_EVENT_projP, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

◆ ConstrainMatrixMultOpR

MoFEMErrorCode ConstrainMatrixMultOpR	(	Mat	R,
		Vec	x,
		Vec	f )

friend

Multiplication operator for R = CT(CCT)^-1.

  Mat R; //for problem
  ConstrainMatrixCtx
projection_matrix_ctx(m_field,problem_name,contrains_problem_name);
  CHKERR MatCreateShell(PETSC_COMM_WORLD,m,m,M,M,&projection_matrix_ctx,&R);
  CHKERR
MatShellSetOperation(R,MATOP_MULT,(void(*)(void))ConstrainMatrixMultOpR);

Definition at line 221 of file ConstrainMatrixCtx.cpp.

                                                           {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(R, &void_ctx);
  ConstrainMatrixCtx *ctx = (ConstrainMatrixCtx *)void_ctx;
  PetscLogEventBegin(ctx->MOFEM_EVENT_projR, 0, 0, 0, 0);
  if (ctx->initQorP)
    SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
            "you have to call first initQorP or use Q matrix");
  CHKERR KSPSolve(ctx->kSP, x, ctx->CCTm1_Cx);
  CHKERR MatMult(ctx->CT, ctx->CCTm1_Cx, ctx->CT_CCTm1_Cx);
  CHKERR VecZeroEntries(f);
  CHKERR VecGhostUpdateBegin(f, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(f, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecScatterBegin(ctx->sCatter, ctx->CT_CCTm1_Cx, f, INSERT_VALUES,
                         SCATTER_REVERSE);
  CHKERR VecScatterEnd(ctx->sCatter, ctx->CT_CCTm1_Cx, f, INSERT_VALUES,
                       SCATTER_REVERSE);
  PetscLogEventEnd(ctx->MOFEM_EVENT_projR, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

◆ ConstrainMatrixMultOpRT

MoFEMErrorCode ConstrainMatrixMultOpRT	(	Mat	RT,
		Vec	x,
		Vec	f )

friend

Multiplication operator for RT = (CCT)^-TC.

  Mat RT; //for problem
  ConstrainMatrixCtx
projection_matrix_ctx(m_field,problem_name,contrains_problem_name);
  CHKERR MatCreateShell(PETSC_COMM_WORLD,m,m,M,M,&projection_matrix_ctx,&RT);
  CHKERR
MatShellSetOperation(RT,MATOP_MULT,(void(*)(void))ConstrainMatrixMultOpRT);

Definition at line 243 of file ConstrainMatrixCtx.cpp.

                                                             {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(RT, &void_ctx);
  ConstrainMatrixCtx *ctx = (ConstrainMatrixCtx *)void_ctx;
  PetscLogEventBegin(ctx->MOFEM_EVENT_projRT, 0, 0, 0, 0);
  if (ctx->initQorP)
    SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
            "you have to call first initQorP or use Q matrix");
  CHKERR VecScatterBegin(ctx->sCatter, x, ctx->X, INSERT_VALUES,
                         SCATTER_FORWARD);
  CHKERR VecScatterEnd(ctx->sCatter, x, ctx->X, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR MatMult(ctx->C, ctx->X, ctx->Cx);
  CHKERR KSPSolve(ctx->kSP, ctx->Cx, f);
  PetscLogEventEnd(ctx->MOFEM_EVENT_projRT, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

◆ ProjectionMatrixMultOpQ

MoFEMErrorCode ProjectionMatrixMultOpQ	(	Mat	Q,
		Vec	x,
		Vec	f )

friend

Multiplication operator for Q = I-CTC(CCT)^-1C.

  Mat Q; //for problem
  ConstrainMatrixCtx
projection_matrix_ctx(m_field,problem_name,contrains_problem_name);
  CHKERR MatCreateShell(PETSC_COMM_WORLD,m,m,M,M,&projection_matrix_ctx,&Q);
  CHKERR
MatShellSetOperation(Q,MATOP_MULT,(void(*)(void))ProjectionMatrixMultOpQ);
  CHKERR
MatShellSetOperation(Q,MATOP_DESTROY,(void(*)(void))ConstrainMatrixDestroyOpPorQ);

Definition at line 163 of file ConstrainMatrixCtx.cpp.

                                                            {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(Q, &void_ctx);
  ConstrainMatrixCtx *ctx = (ConstrainMatrixCtx *)void_ctx;
  PetscLogEventBegin(ctx->MOFEM_EVENT_projQ, 0, 0, 0, 0);
  CHKERR ctx->initializeQorP(x);
  CHKERR VecCopy(x, f);
  CHKERR VecScatterBegin(ctx->sCatter, x, ctx->X, INSERT_VALUES,
                         SCATTER_FORWARD);
  CHKERR VecScatterEnd(ctx->sCatter, x, ctx->X, INSERT_VALUES, SCATTER_FORWARD);
  if (debug) {
    // CHKERR VecView(ctx->X,PETSC_VIEWER_STDOUT_WORLD);
    CHKERR VecScatterBegin(ctx->sCatter, ctx->X, f, INSERT_VALUES,
                           SCATTER_REVERSE);
    CHKERR VecScatterEnd(ctx->sCatter, ctx->X, f, INSERT_VALUES,
                         SCATTER_REVERSE);
    PetscBool flg;
    CHKERR VecEqual(x, f, &flg);
    if (flg == PETSC_FALSE)
      SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE, "scatter is not working");
  }
  CHKERR MatMult(ctx->C, ctx->X, ctx->Cx);
  CHKERR KSPSolve(ctx->kSP, ctx->Cx, ctx->CCTm1_Cx);
  CHKERR MatMult(ctx->CT, ctx->CCTm1_Cx, ctx->CT_CCTm1_Cx);
  CHKERR VecScale(ctx->CT_CCTm1_Cx, -1);
  CHKERR VecScatterBegin(ctx->sCatter, ctx->CT_CCTm1_Cx, f, ADD_VALUES,
                         SCATTER_REVERSE);
  CHKERR VecScatterEnd(ctx->sCatter, ctx->CT_CCTm1_Cx, f, ADD_VALUES,
                       SCATTER_REVERSE);
  PetscLogEventEnd(ctx->MOFEM_EVENT_projQ, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

◆ absTol

PetscReal MoFEM::ConstrainMatrixCtx::absTol

Definition at line 72 of file ConstrainMatrixCtx.hpp.

◆ C

Mat MoFEM::ConstrainMatrixCtx::C

Definition at line 23 of file ConstrainMatrixCtx.hpp.

◆ cancelKSPMonitor

bool MoFEM::ConstrainMatrixCtx::cancelKSPMonitor

Definition at line 29 of file ConstrainMatrixCtx.hpp.

◆ CCT

Mat MoFEM::ConstrainMatrixCtx::CCT

Definition at line 23 of file ConstrainMatrixCtx.hpp.

◆ CCTm1_Cx

Vec MoFEM::ConstrainMatrixCtx::CCTm1_Cx

Definition at line 24 of file ConstrainMatrixCtx.hpp.

◆ createKSP

bool MoFEM::ConstrainMatrixCtx::createKSP

Definition at line 27 of file ConstrainMatrixCtx.hpp.

◆ createScatter

bool MoFEM::ConstrainMatrixCtx::createScatter

Definition at line 28 of file ConstrainMatrixCtx.hpp.

◆ CT

Mat MoFEM::ConstrainMatrixCtx::CT

Definition at line 23 of file ConstrainMatrixCtx.hpp.

◆ CT_CCTm1_Cx

Vec MoFEM::ConstrainMatrixCtx::CT_CCTm1_Cx

Definition at line 24 of file ConstrainMatrixCtx.hpp.

◆ CTC

Mat MoFEM::ConstrainMatrixCtx::CTC

Definition at line 23 of file ConstrainMatrixCtx.hpp.

◆ CTCx

Vec MoFEM::ConstrainMatrixCtx::CTCx

Definition at line 24 of file ConstrainMatrixCtx.hpp.

◆ Cx

Vec MoFEM::ConstrainMatrixCtx::Cx

Definition at line 24 of file ConstrainMatrixCtx.hpp.

◆ dTol

PetscReal MoFEM::ConstrainMatrixCtx::dTol

Definition at line 72 of file ConstrainMatrixCtx.hpp.

◆ initQorP

bool MoFEM::ConstrainMatrixCtx::initQorP

Definition at line 26 of file ConstrainMatrixCtx.hpp.

◆ initQTKQ

bool MoFEM::ConstrainMatrixCtx::initQTKQ

Definition at line 26 of file ConstrainMatrixCtx.hpp.

◆ K

Mat MoFEM::ConstrainMatrixCtx::K

Definition at line 23 of file ConstrainMatrixCtx.hpp.

◆ KQx

Vec MoFEM::ConstrainMatrixCtx::KQx

Definition at line 25 of file ConstrainMatrixCtx.hpp.

◆ kSP

KSP MoFEM::ConstrainMatrixCtx::kSP

Definition at line 22 of file ConstrainMatrixCtx.hpp.

◆ maxIts

PetscInt MoFEM::ConstrainMatrixCtx::maxIts

Definition at line 73 of file ConstrainMatrixCtx.hpp.

◆ mField

MoFEM::Interface& MoFEM::ConstrainMatrixCtx::mField

Definition at line 20 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projCTC_QTKQ

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projCTC_QTKQ

Definition at line 43 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projInit

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projInit

Definition at line 38 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projP

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projP

Definition at line 40 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projQ

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projQ

Definition at line 39 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projR

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projR

Definition at line 41 of file ConstrainMatrixCtx.hpp.

◆ MOFEM_EVENT_projRT

PetscLogEvent MoFEM::ConstrainMatrixCtx::MOFEM_EVENT_projRT

Definition at line 42 of file ConstrainMatrixCtx.hpp.

◆ ownConstrainMatrix

bool MoFEM::ConstrainMatrixCtx::ownConstrainMatrix

Definition at line 30 of file ConstrainMatrixCtx.hpp.

◆ Qx

Vec MoFEM::ConstrainMatrixCtx::Qx

Definition at line 25 of file ConstrainMatrixCtx.hpp.

◆ rTol

PetscReal MoFEM::ConstrainMatrixCtx::rTol

Definition at line 72 of file ConstrainMatrixCtx.hpp.

◆ sCatter

VecScatter MoFEM::ConstrainMatrixCtx::sCatter

Definition at line 35 of file ConstrainMatrixCtx.hpp.

◆ X

Vec MoFEM::ConstrainMatrixCtx::X

Definition at line 25 of file ConstrainMatrixCtx.hpp.

◆ xProblem

string MoFEM::ConstrainMatrixCtx::xProblem

Definition at line 36 of file ConstrainMatrixCtx.hpp.

◆ yProblem

string MoFEM::ConstrainMatrixCtx::yProblem

Definition at line 36 of file ConstrainMatrixCtx.hpp.

The documentation for this struct was generated from the following files:

src/finite_elements/ConstrainMatrixCtx.hpp
src/finite_elements/impl/ConstrainMatrixCtx.cpp

Public Member Functions

Public Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ ConstrainMatrixCtx() [1/2]

◆ ConstrainMatrixCtx() [2/2]

◆ ~ConstrainMatrixCtx()

Member Function Documentation

◆ destroyQorP()

◆ destroyQTKQ()

◆ initializeQorP()

◆ initializeQTKQ()

◆ recalculateCTandCCT()

◆ recalculateCTC()

Friends And Related Symbol Documentation

◆ ConstrainMatrixDestroyOpPorQ

◆ ConstrainMatrixDestroyOpQTKQ

◆ ConstrainMatrixMultOpCTC_QTKQ

◆ ConstrainMatrixMultOpP

◆ ConstrainMatrixMultOpR

◆ ConstrainMatrixMultOpRT

◆ ProjectionMatrixMultOpQ

Member Data Documentation

◆ absTol

◆ C

◆ cancelKSPMonitor

◆ CCT

◆ CCTm1_Cx

◆ createKSP

◆ createScatter

◆ CT

◆ CT_CCTm1_Cx

◆ CTC

◆ CTCx

◆ Cx

◆ dTol

◆ initQorP

◆ initQTKQ

◆ K

◆ KQx

◆ kSP

◆ maxIts

◆ mField

◆ MOFEM_EVENT_projCTC_QTKQ

◆ MOFEM_EVENT_projInit

◆ MOFEM_EVENT_projP

◆ MOFEM_EVENT_projQ

◆ MOFEM_EVENT_projR

◆ MOFEM_EVENT_projRT

◆ ownConstrainMatrix

◆ Qx

◆ rTol

◆ sCatter

◆ X

◆ xProblem

◆ yProblem