SetUpSchurImpl Struct Reference

Inheritance diagram for SetUpSchurImpl:

Collaboration diagram for SetUpSchurImpl:

Public Member Functions
	SetUpSchurImpl (MoFEM::Interface &m_field)
virtual	~SetUpSchurImpl ()=default
MoFEMErrorCode	setUp (SmartPetscObj< KSP >)
	SetUpSchurImpl (MoFEM::Interface &m_field, SmartPetscObj< Mat > m, SmartPetscObj< Mat > p, EshelbianCore *ep_core_ptr)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (KSP solver)
MoFEMErrorCode	preProc ()
MoFEMErrorCode	postProc ()
	SetUpSchurImpl (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_is, SmartPetscObj< AO > ao_up)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (TS solver)
MoFEMErrorCode	preProc ()
MoFEMErrorCode	postProc ()
	SetUpSchurImpl (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_is, SmartPetscObj< AO > ao_up)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (TS solver)
MoFEMErrorCode	preProc ()
MoFEMErrorCode	postProc ()
	SetUpSchurImpl (MoFEM::Interface &m_field)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (SmartPetscObj< TS > solver)
	SetUpSchurImpl (MoFEM::Interface &m_field)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (SmartPetscObj< TS > solver)
	SetUpSchurImpl (MoFEM::Interface &m_field)
virtual	~SetUpSchurImpl ()=default
MoFEMErrorCode	setUp (SmartPetscObj< KSP > solver)

Private Member Functions
MoFEMErrorCode	createSubDM ()
MoFEMErrorCode	setOperator ()
MoFEMErrorCode	setPC (PC pc)
MoFEMErrorCode	setDiagonalPC (PC pc)
MoFEMErrorCode	setOperator ()
MoFEMErrorCode	setPC (PC pc)
SmartPetscObj< DM >	createSubDM ()
MoFEMErrorCode	setEntities ()
MoFEMErrorCode	createSubDM ()
MoFEMErrorCode	setOperator ()
MoFEMErrorCode	setPC (PC pc)
MoFEMErrorCode	setDiagonalPC (PC pc)

Private Attributes
MoFEM::Interface &	mField
SmartPetscObj< Mat >	S
SmartPetscObj< Mat >	M
SmartPetscObj< Mat >	P
SmartPetscObj< AO >	aoUp
SmartPetscObj< IS >	isA00
SmartPetscObj< IS >	isShift
EshelbianCore *	epCorePtr
SmartPetscObj< DM >	subDM
	field split sub dm More...
SmartPetscObj< IS >	fieldSplitIS
	IS for split Schur block. More...
SmartPetscObj< AO >	aoSchur
SmartPetscObj< DM >	schurDM
SmartPetscObj< DM >	blockDM
Range	volEnts
Range	subEnts

Additional Inherited Members
Static Public Member Functions inherited from SetUpSchur
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field)
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_it, SmartPetscObj< AO > ao_map)
	Create data structure for handling Schur complement. More...
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_it, SmartPetscObj< AO > ao_map)
	Create data structure for handling Schur complement. More...
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field)
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field)
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field)
Protected Member Functions inherited from SetUpSchur
	SetUpSchur ()=default
virtual	~SetUpSchur ()=default
	SetUpSchur ()=default
	SetUpSchur ()=default
	SetUpSchur ()=default
	SetUpSchur ()=default
	SetUpSchur ()=default

Detailed Description

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 512 of file test_broken_space.cpp.

Constructor & Destructor Documentation

SetUpSchurImpl() [1/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field )

inline

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 514 of file test_broken_space.cpp.

: SetUpSchur(), mField(m_field) {}

~SetUpSchurImpl() [1/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

virtualdefault

Examples

plastic.cpp, and test_broken_space.cpp.

SetUpSchurImpl() [2/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field, SmartPetscObj< Mat >  m, SmartPetscObj< Mat >  p, EshelbianCore *  ep_core_ptr  )

inline

Definition at line 11 of file SetUpSchurImpl.cpp.

      : SetUpSchur(), mField(m_field), M(m), P(p), epCorePtr(ep_core_ptr) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [2/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 21 of file SetUpSchurImpl.cpp.

{ S.reset(); }

SetUpSchurImpl() [3/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field, SmartPetscObj< DM >  sub_dm, SmartPetscObj< IS >  field_split_is, SmartPetscObj< AO >  ao_up  )

inline

Definition at line 1285 of file plastic.cpp.

      : SetUpSchur(), mField(m_field), subDM(sub_dm),
        fieldSplitIS(field_split_is), aoSchur(ao_up) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [3/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 1296 of file plastic.cpp.

{ S.reset(); }

SetUpSchurImpl() [4/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field, SmartPetscObj< DM >  sub_dm, SmartPetscObj< IS >  field_split_is, SmartPetscObj< AO >  ao_up  )

inline

Definition at line 1273 of file plasticity_old_schur.cpp.

      : SetUpSchur(), mField(m_field), subDM(sub_dm),
        fieldSplitIS(field_split_is), aoUp(ao_up) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [4/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 1284 of file plasticity_old_schur.cpp.

                            {
#ifdef ADD_CONTACT
    A.reset();
    P.reset();
#endif // ADD_CONTACT
    S.reset();
  }

SetUpSchurImpl() [5/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field )

inline

Definition at line 960 of file contact.cpp.

: SetUpSchur(), mField(m_field) {}

~SetUpSchurImpl() [5/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 962 of file contact.cpp.

{}

SetUpSchurImpl() [6/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field )

inline

Definition at line 935 of file incompressible_elasticity.cpp.

: SetUpSchur(), mField(m_field) {}

~SetUpSchurImpl() [6/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 936 of file incompressible_elasticity.cpp.

{ S.reset(); }

SetUpSchurImpl() [7/7]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field )

inline

Definition at line 884 of file elastic.cpp.

                                          : SetUpSchur(), mField(m_field) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [7/7]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

virtualdefault

Member Function Documentation

createSubDM() [1/3]

MoFEMErrorCode SetUpSchurImpl::createSubDM ( )

private

createSubDM() [2/3]

SmartPetscObj<DM> SetUpSchurImpl::createSubDM ( )

private

createSubDM() [3/3]

MoFEMErrorCode SetUpSchurImpl::createSubDM ( )

private

Definition at line 1056 of file contact.cpp.

                                           {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
 
  auto create_dm = [&](const char *name, const char *field_name) {
    auto dm = createDM(mField.get_comm(), "DMMOFEM");
    auto create_dm_imp = [&]() {
      MoFEMFunctionBegin;
      CHKERR DMMoFEMCreateSubDM(dm, simple->getDM(), name);
      CHKERR DMMoFEMSetSquareProblem(dm, PETSC_TRUE);
      CHKERR DMMoFEMAddElement(dm, simple->getDomainFEName());
      CHKERR DMMoFEMAddSubFieldRow(dm, field_name);
      CHKERR DMMoFEMAddSubFieldCol(dm, field_name);
      CHKERR DMSetUp(dm);
      MoFEMFunctionReturn(0);
    };
    CHK_THROW_MESSAGE(
        create_dm_imp(),
        "Error in creating schurDM. It is possible that schurDM is "
        "already created");
    return dm;
  };
 
  // Note: here we can make block with bubbles of "U" and "SIGMA" fields. See
  // vec-0 where bubbles are added.
 
  schurDM = create_dm("SCHUR", "U");
  blockDM = create_dm("BLOCK", "SIGMA");
 
  if constexpr (AT == AssemblyType::BLOCK_SCHUR) {
 
    auto get_nested_mat_data = [&](auto schur_dm, auto block_dm) {
      auto block_mat_data = createBlockMatStructure(
          simple->getDM(),
 
          {{
 
              simple->getDomainFEName(),
 
              {
 
                  {"U", "U"}, {"SIGMA", "U"}, {"U", "SIGMA"}, {"SIGMA", "SIGMA"}
 
              }}}
 
      );
 
      return getNestSchurData(
 
          {schurDM, blockDM}, block_mat_data,
 
          {"SIGMA"}, {nullptr}, true
 
      );
    };
 
    auto nested_mat_data = get_nested_mat_data(schurDM, blockDM);
    CHKERR DMMoFEMSetNestSchurData(simple->getDM(), nested_mat_data);
 
  } else {
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
            "Only BLOCK_SCHUR is implemented");
  }
 
  MoFEMFunctionReturn(0);
}

postProc() [1/3]

MoFEMErrorCode SetUpSchurImpl::postProc

(

)

Examples

plastic.cpp.

Definition at line 213 of file SetUpSchurImpl.cpp.

                                        {
  MoFEMFunctionBegin;
  if (S) {
    MOFEM_LOG("EP", Sev::verbose) << "Assemble Schur";
    CHKERR MatAssemblyBegin(S, MAT_FINAL_ASSEMBLY);
    CHKERR MatAssemblyEnd(S, MAT_FINAL_ASSEMBLY);
  }
  CHKERR MatAssemblyBegin(M, MAT_FINAL_ASSEMBLY);
  CHKERR MatAssemblyEnd(M, MAT_FINAL_ASSEMBLY);
  CHKERR MatAssemblyBegin(P, MAT_FINAL_ASSEMBLY);
  CHKERR MatAssemblyEnd(P, MAT_FINAL_ASSEMBLY);
  CHKERR MatAXPY(P, 1, M, SAME_NONZERO_PATTERN);
 
  MoFEMFunctionReturn(0);
}

postProc() [2/3]

MoFEMErrorCode SetUpSchurImpl::postProc

(

)

postProc() [3/3]

MoFEMErrorCode SetUpSchurImpl::postProc

(

)

preProc() [1/3]

MoFEMErrorCode SetUpSchurImpl::preProc

(

)

Examples

plastic.cpp.

Definition at line 202 of file SetUpSchurImpl.cpp.

                                       {
  MoFEMFunctionBegin;
  MOFEM_LOG("EP", Sev::verbose) << "Zero Schur";
  if (SetUpSchurImpl::S) {
    CHKERR MatZeroEntries(S);
  }
  CHKERR MatZeroEntries(M);
  CHKERR MatZeroEntries(P);
  MoFEMFunctionReturn(0);
}

preProc() [2/3]

MoFEMErrorCode SetUpSchurImpl::preProc

(

)

preProc() [3/3]

MoFEMErrorCode SetUpSchurImpl::preProc

(

)

setDiagonalPC() [1/2]

MoFEMErrorCode SetUpSchurImpl::setDiagonalPC ( PC  pc )

private

setDiagonalPC() [2/2]

MoFEMErrorCode SetUpSchurImpl::setDiagonalPC ( PC  pc )

private

Definition at line 1213 of file contact.cpp.

                                                  {
  MoFEMFunctionBegin;
  KSP *subksp;
  CHKERR PCFieldSplitSchurGetSubKSP(pc, PETSC_NULL, &subksp);
  auto get_pc = [](auto ksp) {
    PC pc_raw;
    CHKERR KSPGetPC(ksp, &pc_raw);
    return SmartPetscObj<PC>(pc_raw, true); // bump reference
  };
  CHKERR setSchurMatSolvePC(get_pc(subksp[0]));
  CHKERR PetscFree(subksp);
  MoFEMFunctionReturn(0);
}

setEntities()

MoFEMErrorCode SetUpSchurImpl::setEntities ( )

private

Definition at line 960 of file elastic.cpp.

                                           {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  CHKERR mField.get_moab().get_entities_by_dimension(simple->getMeshset(),
                                                     SPACE_DIM, volEnts);
  CHKERR mField.get_moab().get_entities_by_handle(simple->getMeshset(),
                                                  subEnts);
  subEnts = subtract(subEnts, volEnts);
  MoFEMFunctionReturn(0);
};

setOperator() [1/3]

MoFEMErrorCode SetUpSchurImpl::setOperator ( )

private

setOperator() [2/3]

MoFEMErrorCode SetUpSchurImpl::setOperator ( )

private

setOperator() [3/3]

MoFEMErrorCode SetUpSchurImpl::setOperator ( )

private

Definition at line 1123 of file contact.cpp.

                                           {
  MoFEMFunctionBegin;
 
  double eps_stab = 1e-4;
  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-eps_stab", &eps_stab,
                               PETSC_NULL);
 
  using B = FormsIntegrators<BoundaryEleOp>::Assembly<
      BLOCK_PRECONDITIONER_SCHUR>::BiLinearForm<IT>;
  using OpMassStab = B::OpMass<3, SPACE_DIM * SPACE_DIM>;
 
  auto simple = mField.getInterface<Simple>();
  auto pip = mField.getInterface<PipelineManager>();
 
  // block data structure
  boost::shared_ptr<BlockStructure> block_data;
  CHKERR DMMoFEMGetBlocMatData(simple->getDM(), block_data);
 
  // Boundary
  auto dm_is = getDMSubData(schurDM)->getSmartRowIs();
  auto ao_up = createAOMappingIS(dm_is, PETSC_NULL);
 
  pip->getOpBoundaryLhsPipeline().push_front(createOpSchurAssembleBegin());
  pip->getOpBoundaryLhsPipeline().push_back(
      new OpMassStab("SIGMA", "SIGMA",
                     [eps_stab](double, double, double) { return eps_stab; }));
  pip->getOpBoundaryLhsPipeline().push_back(
 
      createOpSchurAssembleEnd({"SIGMA"}, {nullptr}, {ao_up}, {S}, {false},
                               false, block_data)
 
  );
 
  // Domain
  pip->getOpDomainLhsPipeline().push_front(createOpSchurAssembleBegin());
  pip->getOpDomainLhsPipeline().push_back(
 
      createOpSchurAssembleEnd({"SIGMA"}, {nullptr}, {ao_up}, {S}, {false},
                               false, block_data)
 
  );
 
  auto pre_proc_schur_lhs_ptr = boost::make_shared<FEMethod>();
  auto post_proc_schur_lhs_ptr = boost::make_shared<FEMethod>();
 
  pre_proc_schur_lhs_ptr->preProcessHook = [this]() {
    MoFEMFunctionBegin;
    CHKERR MatZeroEntries(S);
    MOFEM_LOG("CONTACT", Sev::verbose) << "Lhs Assemble Begin";
    MoFEMFunctionReturn(0);
  };
 
  post_proc_schur_lhs_ptr->postProcessHook = [this, ao_up,
                                              post_proc_schur_lhs_ptr]() {
    MoFEMFunctionBegin;
    MOFEM_LOG("CONTACT", Sev::verbose) << "Lhs Assemble End";
    auto print_mat_norm = [this](auto a, std::string prefix) {
      MoFEMFunctionBegin;
      double nrm;
      CHKERR MatNorm(a, NORM_FROBENIUS, &nrm);
      MOFEM_LOG("CONTACT", Sev::noisy) << prefix << " norm = " << nrm;
      MoFEMFunctionReturn(0);
    };
    CHKERR MatAssemblyBegin(S, MAT_FINAL_ASSEMBLY);
    CHKERR MatAssemblyEnd(S, MAT_FINAL_ASSEMBLY);
    CHKERR EssentialPostProcLhs<DisplacementCubitBcData>(
        mField, post_proc_schur_lhs_ptr, 1, S, ao_up)();
#ifndef NDEBUG
    CHKERR print_mat_norm(S, "S");
#endif // NDEBUG
    MOFEM_LOG("CONTACT", Sev::verbose) << "Lhs Assemble Finish";
    MoFEMFunctionReturn(0);
  };
 
  auto ts_ctx_ptr = getDMTsCtx(simple->getDM());
  ts_ctx_ptr->getPreProcessIJacobian().push_front(pre_proc_schur_lhs_ptr);
  ts_ctx_ptr->getPostProcessIJacobian().push_back(post_proc_schur_lhs_ptr);
 
  MoFEMFunctionReturn(0);
}

setPC() [1/3]

MoFEMErrorCode SetUpSchurImpl::setPC ( PC  pc )

private

setPC() [2/3]

MoFEMErrorCode SetUpSchurImpl::setPC ( PC  pc )

private

setPC() [3/3]

MoFEMErrorCode SetUpSchurImpl::setPC ( PC  pc )

private

Definition at line 1204 of file contact.cpp.

                                          {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  auto block_is = getDMSubData(blockDM)->getSmartRowIs();
  CHKERR PCFieldSplitSetIS(pc, NULL, block_is);
  CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
  MoFEMFunctionReturn(0);
}

setUp() [1/7]

MoFEMErrorCode SetUpSchurImpl::setUp

(

KSP

solver

)

Definition at line 39 of file SetUpSchurImpl.cpp.

                                               {
  MoFEMFunctionBegin;
 
  auto get_ents_by_dim = [&](int dim) {
    Range ents;
    CHKERR mField.get_moab().get_entities_by_dimension(0, dim, ents);
    return ents;
  };
 
  std::vector<std::string> schur_field_list{
 
      epCorePtr->stretchTensor, 
      
      epCorePtr->bubbleField,
 
      epCorePtr->rotAxis,       
 
      epCorePtr->spatialL2Disp,
 
      epCorePtr->piolaStress
 
  };
  std::vector<std::string> field_list{epCorePtr->piolaStress,
                                      epCorePtr->contactDisp};
 
  auto create_schur_dm = [&](SmartPetscObj<DM> &dm_sub) {
    MoFEMFunctionBegin;
 
    dm_sub = createDM(mField.get_comm(), "DMMOFEM");
    CHKERR DMMoFEMCreateSubDM(dm_sub, epCorePtr->dmElastic, "SUB_SCHUR");
    CHKERR DMMoFEMSetSquareProblem(dm_sub, PETSC_TRUE);
    CHKERR DMMoFEMSetIsPartitioned(dm_sub, PETSC_TRUE);
    CHKERR DMMoFEMAddElement(dm_sub, epCorePtr->elementVolumeName);
    CHKERR DMMoFEMAddElement(dm_sub, epCorePtr->naturalBcElement);
    CHKERR DMMoFEMAddElement(dm_sub, epCorePtr->contactElement);
 
    auto faces_ptr = boost::make_shared<Range>(get_ents_by_dim(SPACE_DIM - 1));
    std::vector<boost::shared_ptr<Range>> dm_range_list{faces_ptr, nullptr};
 
    if (field_list.size() != dm_range_list.size()) {
      SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
              "Both ranges should have the same size");
    }
 
    int r_idx = 0;
    for (auto f : field_list) {
      MOFEM_LOG("EP", Sev::inform) << "Add schur field: " << f;
      CHKERR DMMoFEMAddSubFieldRow(dm_sub, f, dm_range_list[r_idx]);
      CHKERR DMMoFEMAddSubFieldCol(dm_sub, f, dm_range_list[r_idx]);
      ++r_idx;
    }
    CHKERR DMSetUp(dm_sub);
    MoFEMFunctionReturn(0);
  };
 
  auto create_a00_is = [&](SmartPetscObj<IS> &is_a00) {
    MoFEMFunctionBegin;
    auto vols = get_ents_by_dim(SPACE_DIM);
    std::vector<SmartPetscObj<IS>> is_vec;
    std::vector<Range *> range_list_ptr(schur_field_list.size(), nullptr);
    range_list_ptr.back() = &vols;
 
    int r_idx = 0;
    for (auto f : schur_field_list) {
      SmartPetscObj<IS> is;
      CHKERR mField.getInterface<ISManager>()->isCreateProblemFieldAndRank(
          "ELASTIC_PROBLEM", ROW, f, 0, MAX_DOFS_ON_ENTITY, is,
          range_list_ptr[r_idx]);
      is_vec.push_back(is);
      ++r_idx;
    }
    if (is_vec.size()) {
      is_a00 = is_vec[0];
      for (auto i = 1; i < is_vec.size(); ++i) {
        IS is_union_raw;
        CHKERR ISExpand(is_a00, is_vec[i], &is_union_raw);
        is_a00 = SmartPetscObj<IS>(is_union_raw);
      }
      CHKERR ISSort(is_a00);
    }
    MoFEMFunctionReturn(0);
  };
 
  PC pc;
  CHKERR KSPSetFromOptions(solver);
  CHKERR KSPGetPC(solver, &pc);
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
 
    SmartPetscObj<DM> schur_dm, a00_dm;
    CHKERR create_schur_dm(schur_dm);
    CHKERR create_a00_is(isA00);
 
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
 
    S = createDMMatrix(schur_dm);
 
    auto set_ops = [&]() {
      MoFEMFunctionBegin;
      auto range_ptr = boost::make_shared<Range>(get_ents_by_dim(SPACE_DIM));
      std::vector<boost::shared_ptr<Range>> ranges_list(schur_field_list.size(),
                                                        nullptr);
      ranges_list.back() = range_ptr;
      std::vector<SmartPetscObj<AO>> ao_list(schur_field_list.size(),
                                             SmartPetscObj<AO>());
      std::vector<SmartPetscObj<Mat>> mat_list(schur_field_list.size(),
                                               SmartPetscObj<Mat>());
      std::vector<bool> symm_list(schur_field_list.size(), false);
 
      auto dm_is = getDMSubData(schur_dm)->getSmartRowIs();
      aoUp = createAOMappingIS(dm_is, PETSC_NULL);
      ao_list.back() = aoUp;
      mat_list.back() = S;
 
      const bool symmetric_system =
          EshelbianCore::gradApperoximator <= MODERATE_ROT &&
          EshelbianCore::rotSelector == SMALL_ROT;
      epCorePtr->elasticFeLhs->getOpPtrVector().push_front(
          new OpSchurAssembleBegin());
      epCorePtr->elasticFeLhs->getOpPtrVector().push_back(
          new OpSchurAssembleEnd<SCHUR_DGESV>(schur_field_list, ranges_list,
                                              ao_list, mat_list, symm_list,
                                              symmetric_system));
      epCorePtr->elasticBcLhs->getOpPtrVector().push_front(
          new OpSchurAssembleBegin());
      epCorePtr->elasticBcLhs->getOpPtrVector().push_back(
          new OpSchurAssembleEnd<SCHUR_DGESV>(schur_field_list, ranges_list,
                                              ao_list, mat_list, symm_list,
                                              symmetric_system));
 
      MoFEMFunctionReturn(0);
    };
 
    auto set_pc = [&]() {
      MoFEMFunctionBegin;
      CHKERR PCFieldSplitSetIS(pc, NULL, isA00);
      CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
      MoFEMFunctionReturn(0);
    };
 
    CHKERR set_ops();
    CHKERR set_pc();
 
  } else {
    epCorePtr->elasticFeLhs->getOpPtrVector().push_front(
        new OpSchurAssembleBegin());
    epCorePtr->elasticFeLhs->getOpPtrVector().push_back(
        new OpSchurAssembleEnd<SCHUR_DGESV>({}, {}, {}, {}, {}));
    epCorePtr->elasticBcLhs->getOpPtrVector().push_front(
        new OpSchurAssembleBegin());
    epCorePtr->elasticBcLhs->getOpPtrVector().push_back(
        new OpSchurAssembleEnd<SCHUR_DGESV>({}, {}, {}, {}, {}));
  }
 
  MoFEMFunctionReturn(0);
}

setUp() [2/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( SmartPetscObj< KSP >  solver )

virtual

Implements SetUpSchur.

setUp() [3/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( SmartPetscObj< KSP >  ksp )

virtual

Implements SetUpSchur.

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 525 of file test_broken_space.cpp.

                                                           {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  auto pip_mng = mField.getInterface<PipelineManager>();
 
  CHKERR KSPSetFromOptions(ksp);
  PC pc;
  CHKERR KSPGetPC(ksp, &pc);
 
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
 
    MOFEM_LOG("AT", Sev::inform) << "Setup Schur pc";
 
    auto create_sub_dm = [&]() {
      auto simple = mField.getInterface<Simple>();
 
      auto create_dm = [&](
 
                           std::string problem_name,
                           std::vector<std::string> fe_names,
                           std::vector<std::string> fields,
 
                           auto dm_type
 
                       ) {
        auto dm = createDM(mField.get_comm(), dm_type);
        auto create_dm_imp = [&]() {
          MoFEMFunctionBegin;
          CHKERR DMMoFEMCreateSubDM(dm, simple->getDM(), problem_name.c_str());
          CHKERR DMMoFEMSetSquareProblem(dm, PETSC_TRUE);
          for (auto fe : fe_names) {
            CHKERR DMMoFEMAddElement(dm, fe);
          }
          CHKERR DMMoFEMAddElement(dm, simple->getSkeletonFEName());
          for (auto field : fields) {
            CHKERR DMMoFEMAddSubFieldRow(dm, field);
            CHKERR DMMoFEMAddSubFieldCol(dm, field);
          }
          CHKERR DMSetUp(dm);
          MoFEMFunctionReturn(0);
        };
        CHK_THROW_MESSAGE(
            create_dm_imp(),
            "Error in creating schurDM. It is possible that schurDM is "
            "already created");
        return dm;
      };
 
      auto schur_dm = create_dm(
 
          "SCHUR",
 
          {simple->getDomainFEName(), simple->getSkeletonFEName()},
 
          {"HYBRID"},
 
          "DMMOFEM_MG");
 
      auto block_dm = create_dm(
 
          "BLOCK",
 
          {simple->getDomainFEName(), simple->getSkeletonFEName()},
 
          {"BROKEN", "U"},
 
          "DMMOFEM");
 
      return std::make_tuple(schur_dm, block_dm);
    };
 
    auto get_nested_mat_data = [&](auto schur_dm, auto block_dm) {
      auto block_mat_data = createBlockMatStructure(
          simple->getDM(),
 
          {
 
              {
 
                  simple->getDomainFEName(),
 
                  {
 
                      {"BROKEN", "BROKEN"},
                      {"U", "U"},
                      {"BROKEN", "U"},
                      {"U", "BROKEN"}
 
                  }
 
              },
 
              {
 
                  simple->getSkeletonFEName(),
 
                  {
 
                      {"BROKEN", "HYBRID"}, {"HYBRID", "BROKEN"}
 
                  }
 
              }
 
          }
 
      );
 
      return getNestSchurData(
 
          {schur_dm, block_dm}, block_mat_data,
 
          {"BROKEN", "U"}, {nullptr, nullptr}, true
 
      );
    };
 
    auto set_ops = [&](auto schur_dm) {
      MoFEMFunctionBegin;
      auto dm_is = getDMSubData(schur_dm)->getSmartRowIs();
      auto ao_up = createAOMappingIS(dm_is, PETSC_NULL);
 
      boost::shared_ptr<BlockStructure> block_data;
      CHKERR DMMoFEMGetBlocMatData(simple->getDM(), block_data);
 
      pip_mng->getOpDomainLhsPipeline().push_front(
          createOpSchurAssembleBegin());
      pip_mng->getOpDomainLhsPipeline().push_back(
 
          createOpSchurAssembleEnd(
              {"BROKEN", "U"}, {nullptr, nullptr}, {SmartPetscObj<AO>(), ao_up},
              {SmartPetscObj<Mat>(), S}, {true, true}, true, block_data)
 
      );
 
      auto pre_proc_schur_lhs_ptr = boost::make_shared<FEMethod>();
      auto post_proc_schur_lhs_ptr = boost::make_shared<FEMethod>();
 
      pre_proc_schur_lhs_ptr->preProcessHook = [this]() {
        MoFEMFunctionBegin;
        CHKERR MatZeroEntries(S);
        MOFEM_LOG("AT", Sev::verbose) << "Lhs Assemble Begin";
        MoFEMFunctionReturn(0);
      };
 
      post_proc_schur_lhs_ptr->postProcessHook = [this, ao_up,
                                                  post_proc_schur_lhs_ptr]() {
        MoFEMFunctionBegin;
        MOFEM_LOG("AT", Sev::verbose) << "Lhs Assemble End";
        CHKERR MatAssemblyBegin(S, MAT_FINAL_ASSEMBLY);
        CHKERR MatAssemblyEnd(S, MAT_FINAL_ASSEMBLY);
        MOFEM_LOG("AT", Sev::verbose) << "Lhs Assemble Finish";
        MoFEMFunctionReturn(0);
      };
 
      auto ksp_ctx_ptr = getDMKspCtx(simple->getDM());
      ksp_ctx_ptr->getPreProcSetOperators().push_front(pre_proc_schur_lhs_ptr);
      ksp_ctx_ptr->getPostProcSetOperators().push_back(post_proc_schur_lhs_ptr);
 
      MoFEMFunctionReturn(0);
    };
 
    auto set_pc = [&](auto pc, auto block_dm) {
      MoFEMFunctionBegin;
      auto block_is = getDMSubData(block_dm)->getSmartRowIs();
      CHKERR PCFieldSplitSetIS(pc, NULL, block_is);
      CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
      MoFEMFunctionReturn(0);
    };
 
    auto set_diagonal_pc = [&](auto pc, auto schur_dm) {
      MoFEMFunctionBegin;
 
      auto A = createDMBlockMat(simple->getDM());
      auto P = createDMNestSchurMat(simple->getDM());
      CHKERR PCSetOperators(pc, A, P);
 
      KSP *subksp;
      CHKERR PCFieldSplitSchurGetSubKSP(pc, PETSC_NULL, &subksp);
      auto get_pc = [](auto ksp) {
        PC pc_raw;
        CHKERR KSPGetPC(ksp, &pc_raw);
        return pc_raw;
      };
      CHKERR setSchurA00MatSolvePC(SmartPetscObj<PC>(get_pc(subksp[0]), true));
 
      auto set_pc_p_mg = [&](auto dm, auto pc) {
        MoFEMFunctionBegin;
 
        CHKERR PCSetDM(pc, dm);
        PetscBool same = PETSC_FALSE;
        PetscObjectTypeCompare((PetscObject)pc, PCMG, &same);
        if (same) {
          // By default do not use shell mg mat. Implementation of SOR is slow.
          CHKERR PCMGSetUpViaApproxOrders(
              pc, createPCMGSetUpViaApproxOrdersCtx(dm, S, false));
          CHKERR PCSetFromOptions(pc);
        }
        MoFEMFunctionReturn(0);
      };
 
      CHKERR set_pc_p_mg(schur_dm, get_pc(subksp[1]));
 
      CHKERR PetscFree(subksp);
      MoFEMFunctionReturn(0);
    };
 
    auto [schur_dm, block_dm] = create_sub_dm();
    auto nested_mat_data = get_nested_mat_data(schur_dm, block_dm);
    CHKERR DMMoFEMSetNestSchurData(simple->getDM(), nested_mat_data);
    S = createDMHybridisedL2Matrix(schur_dm);
    CHKERR MatSetDM(S, PETSC_NULL);
    int bs = (SPACE_DIM == 2) ? NBEDGE_L2(approx_order - 1)
                              : NBFACETRI_L2(approx_order - 1);
    CHKERR MatSetBlockSize(S, bs);
 
    CHKERR set_ops(schur_dm);
    CHKERR set_pc(pc, block_dm);
    DM solver_dm;
    CHKERR KSPGetDM(ksp, &solver_dm);
    CHKERR DMSetMatType(solver_dm, MATSHELL);
 
    CHKERR KSPSetUp(ksp);
    CHKERR set_diagonal_pc(pc, schur_dm);
 
  } else {
    SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
            "PC is not set to PCFIELDSPLIT");
  }
  MoFEMFunctionReturn(0);
}

setUp() [4/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( SmartPetscObj< TS >  solver )

virtual

Implements SetUpSchur.

setUp() [5/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( SmartPetscObj< TS >  solver )

virtual

Implements SetUpSchur.

Definition at line 979 of file contact.cpp.

                                                             {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  auto pip = mField.getInterface<PipelineManager>();
 
  SNES snes;
  CHKERR TSGetSNES(solver, &snes);
  KSP ksp;
  CHKERR SNESGetKSP(snes, &ksp);
  CHKERR KSPSetFromOptions(ksp);
 
  PC pc;
  CHKERR KSPGetPC(ksp, &pc);
 
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
 
    MOFEM_LOG("CONTACT", Sev::inform) << "Setup Schur pc";
 
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "It is expected that Schur matrix is not allocated. This is "
          "possible only if PC is set up twice");
    }
 
    CHKERR createSubDM();
 
    // Add data to DM storage
    S = createDMMatrix(schurDM);
    CHKERR MatSetBlockSize(S, SPACE_DIM);
    // CHKERR MatSetOption(S, MAT_SYMMETRIC, PETSC_TRUE);
 
    // Set DM to use shell block matrix
    DM solver_dm;
    CHKERR TSGetDM(solver, &solver_dm);
    CHKERR DMSetMatType(solver_dm, MATSHELL);
 
    auto ts_ctx_ptr = getDMTsCtx(solver_dm);
    auto A = createDMBlockMat(simple->getDM());
    auto P = createDMNestSchurMat(simple->getDM());
 
    if (is_quasi_static == PETSC_TRUE) {
      auto swap_assemble = [](TS ts, PetscReal t, Vec u, Vec u_t, PetscReal a,
                              Mat A, Mat B, void *ctx) {
        return TsSetIJacobian(ts, t, u, u_t, a, B, A, ctx);
      };
      CHKERR TSSetIJacobian(solver, A, P, swap_assemble, ts_ctx_ptr.get());
    } else {
      auto swap_assemble = [](TS ts, PetscReal t, Vec u, Vec u_t, Vec utt,
                              PetscReal a, PetscReal aa, Mat A, Mat B,
                              void *ctx) {
        return TsSetI2Jacobian(ts, t, u, u_t, utt, a, aa, B, A, ctx);
      };
      CHKERR TSSetI2Jacobian(solver, A, P, swap_assemble, ts_ctx_ptr.get());
    }
    CHKERR KSPSetOperators(ksp, A, P);
 
    CHKERR setOperator();
    CHKERR setPC(pc);
    CHKERR TSSetUp(solver);
    CHKERR KSPSetUp(ksp);
    CHKERR setDiagonalPC(pc);
 
  } else {
    MOFEM_LOG("CONTACT", Sev::inform) << "No Schur PC";
    pip->getOpBoundaryLhsPipeline().push_front(createOpSchurAssembleBegin());
    pip->getOpBoundaryLhsPipeline().push_back(
        createOpSchurAssembleEnd({}, {}, {}, {}, {}, false));
    pip->getOpDomainLhsPipeline().push_front(createOpSchurAssembleBegin());
    pip->getOpDomainLhsPipeline().push_back(
        createOpSchurAssembleEnd({}, {}, {}, {}, {}, false));
  }
  MoFEMFunctionReturn(0);
}

setUp() [6/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( TS  solver )

virtual

Implements SetUpSchur.

setUp() [7/7]

MoFEMErrorCode SetUpSchurImpl::setUp ( TS  solver )

virtual

Implements SetUpSchur.

Definition at line 1311 of file plastic.cpp.

                                              {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  auto pip_mng = mField.getInterface<PipelineManager>();
 
  SNES snes;
  CHKERR TSGetSNES(solver, &snes);
  KSP ksp;
  CHKERR SNESGetKSP(snes, &ksp);
  CHKERR KSPSetFromOptions(ksp);
 
  PC pc;
  CHKERR KSPSetFromOptions(ksp);
  CHKERR KSPGetPC(ksp, &pc);
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
 
    S = createDMMatrix(subDM);
    CHKERR MatSetBlockSize(S, SPACE_DIM);
 
    // Set DM to use shell block matrix
    DM solver_dm;
    CHKERR TSGetDM(solver, &solver_dm);
    CHKERR DMSetMatType(solver_dm, MATSHELL);
 
    auto ts_ctx_ptr = getDMTsCtx(solver_dm);
    auto A = createDMBlockMat(simple->getDM());
    auto P = createDMNestSchurMat(simple->getDM());
 
    if (is_quasi_static == PETSC_TRUE) {
      auto swap_assemble = [](TS ts, PetscReal t, Vec u, Vec u_t, PetscReal a,
                              Mat A, Mat B, void *ctx) {
        return TsSetIJacobian(ts, t, u, u_t, a, B, A, ctx);
      };
      CHKERR TSSetIJacobian(solver, A, P, swap_assemble, ts_ctx_ptr.get());
    } else {
      auto swap_assemble = [](TS ts, PetscReal t, Vec u, Vec u_t, Vec utt,
                              PetscReal a, PetscReal aa, Mat A, Mat B,
                              void *ctx) {
        return TsSetI2Jacobian(ts, t, u, u_t, utt, a, aa, B, A, ctx);
      };
      CHKERR TSSetI2Jacobian(solver, A, P, swap_assemble, ts_ctx_ptr.get());
    }
    CHKERR KSPSetOperators(ksp, A, P);
 
    auto set_ops = [&]() {
      MoFEMFunctionBegin;
      auto pip_mng = mField.getInterface<PipelineManager>();
 
      boost::shared_ptr<BlockStructure> block_data;
      CHKERR DMMoFEMGetBlocMatData(simple->getDM(), block_data);
 
#ifndef ADD_CONTACT
      // Boundary
      pip_mng->getOpBoundaryLhsPipeline().push_front(
          createOpSchurAssembleBegin());
      pip_mng->getOpBoundaryLhsPipeline().push_back(createOpSchurAssembleEnd(
 
          {"EP", "TAU"}, {nullptr, nullptr}, {SmartPetscObj<AO>(), aoSchur},
          {SmartPetscObj<Mat>(), S}, {false, false}, false, block_data
 
          ));
      // Domain
      pip_mng->getOpDomainLhsPipeline().push_front(
          createOpSchurAssembleBegin());
      pip_mng->getOpDomainLhsPipeline().push_back(createOpSchurAssembleEnd(
 
          {"EP", "TAU"}, {nullptr, nullptr}, {SmartPetscObj<AO>(), aoSchur},
          {SmartPetscObj<Mat>(), S}, {false, false}, false, block_data
 
          ));
#else
 
      double eps_stab = 1e-4;
      CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-eps_stab", &eps_stab,
                                   PETSC_NULL);
 
      using B = FormsIntegrators<BoundaryEleOp>::Assembly<
          BLOCK_PRECONDITIONER_SCHUR>::BiLinearForm<IT>;
      using OpMassStab = B::OpMass<3, SPACE_DIM * SPACE_DIM>;
 
      // Boundary
      pip_mng->getOpBoundaryLhsPipeline().push_front(
          createOpSchurAssembleBegin());
      pip_mng->getOpBoundaryLhsPipeline().push_back(
          new OpMassStab("SIGMA", "SIGMA", [eps_stab](double, double, double) {
            return eps_stab;
          }));
      pip_mng->getOpBoundaryLhsPipeline().push_back(createOpSchurAssembleEnd(
 
          {"SIGMA", "EP", "TAU"}, {nullptr, nullptr, nullptr},
          {SmartPetscObj<AO>(), SmartPetscObj<AO>(), aoSchur},
          {SmartPetscObj<Mat>(), SmartPetscObj<Mat>(), S},
          {false, false, false}, false, block_data
 
          ));
      // Domain
      pip_mng->getOpDomainLhsPipeline().push_front(
          createOpSchurAssembleBegin());
      pip_mng->getOpDomainLhsPipeline().push_back(createOpSchurAssembleEnd(
 
          {"SIGMA", "EP", "TAU"}, {nullptr, nullptr, nullptr},
          {SmartPetscObj<AO>(), SmartPetscObj<AO>(), aoSchur},
          {SmartPetscObj<Mat>(), SmartPetscObj<Mat>(), S},
          {false, false, false}, false, block_data
 
          ));
#endif // ADD_CONTACT
      MoFEMFunctionReturn(0);
    };
 
    auto set_assemble_elems = [&]() {
      MoFEMFunctionBegin;
      auto schur_asmb_pre_proc = boost::make_shared<FEMethod>();
      schur_asmb_pre_proc->preProcessHook = [this]() {
        MoFEMFunctionBegin;
        CHKERR MatZeroEntries(S);
        MOFEM_LOG("TIMER", Sev::verbose) << "Lhs Assemble Begin";
        MoFEMFunctionReturn(0);
      };
      auto schur_asmb_post_proc = boost::make_shared<FEMethod>();
 
      schur_asmb_post_proc->postProcessHook = [this, schur_asmb_post_proc]() {
        MoFEMFunctionBegin;
        MOFEM_LOG("TIMER", Sev::verbose) << "Lhs Assemble End";
 
        // Apply essential constrains to Schur complement
        CHKERR MatAssemblyBegin(S, MAT_FINAL_ASSEMBLY);
        CHKERR MatAssemblyEnd(S, MAT_FINAL_ASSEMBLY);
        CHKERR EssentialPostProcLhs<DisplacementCubitBcData>(
            mField, schur_asmb_post_proc, 1, S, aoSchur)();
 
        MoFEMFunctionReturn(0);
      };
      auto ts_ctx_ptr = getDMTsCtx(simple->getDM());
      ts_ctx_ptr->getPreProcessIJacobian().push_front(schur_asmb_pre_proc);
      ts_ctx_ptr->getPostProcessIJacobian().push_front(schur_asmb_post_proc);
      MoFEMFunctionReturn(0);
    };
 
    auto set_pc = [&]() {
      MoFEMFunctionBegin;
      CHKERR PCFieldSplitSetIS(pc, NULL, fieldSplitIS);
      CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
      MoFEMFunctionReturn(0);
    };
 
    auto set_diagonal_pc = [&]() {
      MoFEMFunctionBegin;
      KSP *subksp;
      CHKERR PCFieldSplitSchurGetSubKSP(pc, PETSC_NULL, &subksp);
      auto get_pc = [](auto ksp) {
        PC pc_raw;
        CHKERR KSPGetPC(ksp, &pc_raw);
        return SmartPetscObj<PC>(pc_raw, true); // bump reference
      };
      CHKERR setSchurMatSolvePC(get_pc(subksp[0]));
      CHKERR PetscFree(subksp);
      MoFEMFunctionReturn(0);
    };
 
    CHKERR set_ops();
    CHKERR set_pc();
    CHKERR set_assemble_elems();
 
    CHKERR TSSetUp(solver);
    CHKERR KSPSetUp(ksp);
    CHKERR set_diagonal_pc();
 
  } else {
    pip_mng->getOpBoundaryLhsPipeline().push_front(
        createOpSchurAssembleBegin());
    pip_mng->getOpBoundaryLhsPipeline().push_back(
        createOpSchurAssembleEnd({}, {}, {}, {}, {}, false));
    pip_mng->getOpDomainLhsPipeline().push_front(createOpSchurAssembleBegin());
    pip_mng->getOpDomainLhsPipeline().push_back(
        createOpSchurAssembleEnd({}, {}, {}, {}, {}, false));
  }
 
  // fieldSplitIS.reset();
  // aoSchur.reset();
  MoFEMFunctionReturn(0);
}

Member Data Documentation

aoSchur

SmartPetscObj<AO> SetUpSchurImpl::aoSchur

private

Examples

plastic.cpp.

Definition at line 1308 of file plastic.cpp.

aoUp

SmartPetscObj< AO > SetUpSchurImpl::aoUp

private

Definition at line 31 of file SetUpSchurImpl.cpp.

blockDM

SmartPetscObj< DM > SetUpSchurImpl::blockDM

private

Definition at line 976 of file contact.cpp.

epCorePtr

EshelbianCore* SetUpSchurImpl::epCorePtr

private

Definition at line 36 of file SetUpSchurImpl.cpp.

fieldSplitIS

SmartPetscObj< IS > SetUpSchurImpl::fieldSplitIS

private

IS for split Schur block.

Examples

plastic.cpp.

Definition at line 1307 of file plastic.cpp.

isA00

SmartPetscObj<IS> SetUpSchurImpl::isA00

private

Definition at line 32 of file SetUpSchurImpl.cpp.

isShift

SmartPetscObj<IS> SetUpSchurImpl::isShift

private

Definition at line 33 of file SetUpSchurImpl.cpp.

M

SmartPetscObj<Mat> SetUpSchurImpl::M

private

Definition at line 28 of file SetUpSchurImpl.cpp.

mField

MoFEM::Interface & SetUpSchurImpl::mField

private

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 521 of file test_broken_space.cpp.

P

SmartPetscObj<Mat> SetUpSchurImpl::P

private

Examples

plastic.cpp.

Definition at line 29 of file SetUpSchurImpl.cpp.

S

SmartPetscObj< Mat > SetUpSchurImpl::S

private

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 522 of file test_broken_space.cpp.

schurDM

SmartPetscObj< DM > SetUpSchurImpl::schurDM

private

Definition at line 975 of file contact.cpp.

subDM

SmartPetscObj< DM > SetUpSchurImpl::subDM

private

field split sub dm

Examples

plastic.cpp.

Definition at line 1306 of file plastic.cpp.

subEnts

Range SetUpSchurImpl::subEnts

private

Definition at line 910 of file elastic.cpp.

volEnts

Range SetUpSchurImpl::volEnts

private

Definition at line 909 of file elastic.cpp.

---

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

• atom_tests/test_broken_space.cpp
• users_modules/eshelbian_plasticit/src/impl/SetUpSchurImpl.cpp
• users_modules/tutorials/adv-0/plastic.cpp
• users_modules/tutorials/adv-0/plasticity_old_schur.cpp
• users_modules/tutorials/adv-1/contact.cpp
• users_modules/tutorials/mix-2/incompressible_elasticity.cpp
• users_modules/tutorials/vec-0/elastic.cpp