ArcLengthMatShell Struct Reference

shell matrix for arc-length method More...

#include <users_modules/basic_finite_elements/src/ArcLengthTools.hpp>

Collaboration diagram for ArcLengthMatShell:

Public Member Functions
	ArcLengthMatShell (Mat aij, boost::shared_ptr< ArcLengthCtx > arc_ptr, string problem_name)
virtual	~ArcLengthMatShell ()=default
DEPRECATED	ArcLengthMatShell (Mat aij, ArcLengthCtx *arc_ptr_raw, string problem_name)
MoFEMErrorCode	setLambda (Vec ksp_x, double *lambda, ScatterMode scattermode)

Public Attributes
SmartPetscObj< Mat >	Aij
string	problemName
ArcLengthCtx *	arcPtrRaw

Private Attributes
boost::shared_ptr< ArcLengthCtx >	arcPtr

Friends
MoFEMErrorCode	ArcLengthMatMultShellOp (Mat A, Vec x, Vec f)

Detailed Description

shell matrix for arc-length method

Shell matrix which has structure:

\[ \left[ \begin{array}{cc} \mathbf{K} & -\mathbf{F}_\lambda \\ \textrm{d}\mathbf{b} & D \end{array} \right] \left\{ \begin{array}{c} \delta \Delta \mathbf{x} \\ \delta \Delta \lambda \end{array} \right\} = \left[ \begin{array}{c} -\mathbf{f}_\textrm{int} \\ -r_\lambda \end{array} \right] \]

Definition at line 205 of file ArcLengthTools.hpp.

Constructor & Destructor Documentation

ArcLengthMatShell() [1/2]

ArcLengthMatShell::ArcLengthMatShell	(	Mat	aij,
		boost::shared_ptr< ArcLengthCtx >	arc_ptr,
		string	problem_name
	)

Definition at line 137 of file ArcLengthTools.cpp.

    : Aij(aij, true), problemName(problem_name), arcPtrRaw(arc_ptr.get()),
      arcPtr(arc_ptr) {}

~ArcLengthMatShell()

virtual ArcLengthMatShell::~ArcLengthMatShell ( )

virtualdefault

ArcLengthMatShell() [2/2]

ArcLengthMatShell::ArcLengthMatShell	(	Mat	aij,
		ArcLengthCtx *	arc_ptr_raw,
		string	problem_name
	)

Deprecated:

use constructor with shared_ptr

Definition at line 133 of file ArcLengthTools.cpp.

    : Aij(aij, true), problemName(problem_name), arcPtrRaw(arc_ptr_raw) {}

Member Function Documentation

setLambda()

MoFEMErrorCode ArcLengthMatShell::setLambda	(	Vec	ksp_x,
		double *	lambda,
		ScatterMode	scattermode
	)

Definition at line 143 of file ArcLengthTools.cpp.

                                                                     {
  MoFEMFunctionBegin;
 
  int part = arcPtrRaw->getPart();
  int rank = arcPtrRaw->mField.get_comm_rank();
 
  switch (scattermode) {
  case SCATTER_FORWARD: {
    Vec lambda_ghost;
    if (rank == part) {
      CHKERR VecCreateGhostWithArray(arcPtrRaw->mField.get_comm(), 1, 1, 0,
                                     PETSC_NULL, lambda, &lambda_ghost);
    } else {
      int one[] = {0};
      CHKERR VecCreateGhostWithArray(arcPtrRaw->mField.get_comm(), 0, 1, 1, one,
                                     lambda, &lambda_ghost);
    }
    int idx = arcPtrRaw->getPetscGlobalDofIdx();
    if (part == rank) {
      CHKERR VecGetValues(ksp_x, 1, &idx, lambda);
    }
    CHKERR VecGhostUpdateBegin(lambda_ghost, INSERT_VALUES, SCATTER_FORWARD);
    CHKERR VecGhostUpdateEnd(lambda_ghost, INSERT_VALUES, SCATTER_FORWARD);
    CHKERR VecDestroy(&lambda_ghost);
  } break;
  case SCATTER_REVERSE: {
    if (arcPtrRaw->getPetscLocalDofIdx() != -1) {
      PetscScalar *array;
      CHKERR VecGetArray(ksp_x, &array);
      array[arcPtrRaw->getPetscLocalDofIdx()] = *lambda;
      CHKERR VecRestoreArray(ksp_x, &array);
    }
  } break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "not implemented");
  }
 
  MoFEMFunctionReturn(0);
}

Friends And Related Function Documentation

ArcLengthMatMultShellOp

MoFEMErrorCode ArcLengthMatMultShellOp ( Mat  A, Vec  x, Vec  f  )

friend

mult operator for Arc Length Shell Mat

Definition at line 184 of file ArcLengthTools.cpp.

                                                            {
  MoFEMFunctionBegin;
  void *void_ctx;
  CHKERR MatShellGetContext(A, &void_ctx);
  ArcLengthMatShell *ctx = static_cast<ArcLengthMatShell *>(void_ctx);
  CHKERR MatMult(ctx->Aij, x, f);
  double lambda;
  CHKERR ctx->setLambda(x, &lambda, SCATTER_FORWARD);
  double db_dot_x;
  CHKERR VecDot(ctx->arcPtrRaw->db, x, &db_dot_x);
  double f_lambda;
  f_lambda = ctx->arcPtrRaw->dIag * lambda + db_dot_x;
  CHKERR ctx->setLambda(f, &f_lambda, SCATTER_REVERSE);
  CHKERR VecAXPY(f, lambda, ctx->arcPtrRaw->F_lambda);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

Aij

SmartPetscObj<Mat> ArcLengthMatShell::Aij

Definition at line 207 of file ArcLengthTools.hpp.

arcPtr

boost::shared_ptr<ArcLengthCtx> ArcLengthMatShell::arcPtr

private

Definition at line 223 of file ArcLengthTools.hpp.

arcPtrRaw

ArcLengthCtx* ArcLengthMatShell::arcPtrRaw

Definition at line 209 of file ArcLengthTools.hpp.

problemName

string ArcLengthMatShell::problemName

Definition at line 208 of file ArcLengthTools.hpp.

---

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

• users_modules/basic_finite_elements/src/ArcLengthTools.hpp
• users_modules/basic_finite_elements/src/impl/ArcLengthTools.cpp