Collaboration diagram for MixedPoisson:

[legend]

Classes
struct	CommonData

struct	OpError

Public Member Functions
	MixedPoisson (MoFEM::Interface &m_field)

MoFEMErrorCode	runProblem ()
	[Run programme] More...

Private Member Functions
MoFEMErrorCode	readMesh ()
	[Run programme] More...

MoFEMErrorCode	setupProblem ()
	[Read mesh] More...

MoFEMErrorCode	setIntegrationRules ()
	[Set up problem] More...

MoFEMErrorCode	createCommonData ()
	[Set integration rule] More...

MoFEMErrorCode	assembleSystem ()
	[Create common data] More...

MoFEMErrorCode	solveSystem ()
	[Assemble system] More...

MoFEMErrorCode	checkError (int iter_num=0)
	[Solve and refine loop] More...

MoFEMErrorCode	refineOrder ()
	[Solve] More...

MoFEMErrorCode	solveRefineLoop ()
	[Refine] More...

MoFEMErrorCode	outputResults (int iter_num=0)
	[Check error] More...

Static Private Member Functions
static double	analyticalFunction (const double x, const double y, const double z)
	[Analytical function] More...

static VectorDouble	analyticalFunctionGrad (const double x, const double y, const double z)
	[Analytical function] More...

static double	sourceFunction (const double x, const double y, const double z)
	[Analytical function gradient] More...

static MoFEMErrorCode	getTagHandle (MoFEM::Interface &m_field, const char *name, DataType type, Tag &tag_handle)
	[Source function] More...

Private Attributes
MoFEM::Interface &	mField

Simple *	simpleInterface

Range	domainEntities

double	errorIndicatorIntegral

int	totalElementNumber

int	baseOrder

int	refIterNum

boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Examples: mixed_poisson.cpp.

Definition at line 31 of file mixed_poisson.cpp.

Constructor & Destructor Documentation

◆ MixedPoisson()

MixedPoisson::MixedPoisson ( MoFEM::Interface & m_field )

inline

Definition at line 33 of file mixed_poisson.cpp.

33: mField(m_field) {}

MixedPoisson::mField

MoFEM::Interface & mField

Definition: mixed_poisson.cpp:37

Member Function Documentation

◆ analyticalFunction()

static double MixedPoisson::analyticalFunction	(	const double	x,
		const double	y,
		const double	z
	)

inlinestaticprivate

[Analytical function]

Examples: mixed_poisson.cpp.

Definition at line 47 of file mixed_poisson.cpp.

                                                   {
    return exp(-100. * (sqr(x) + sqr(y))) * cos(M_PI * x) * cos(M_PI * y);
  }

◆ analyticalFunctionGrad()

static VectorDouble MixedPoisson::analyticalFunctionGrad	(	const double	x,
		const double	y,
		const double	z
	)

inlinestaticprivate

[Analytical function]

[Analytical function gradient]

Examples: mixed_poisson.cpp.

Definition at line 54 of file mixed_poisson.cpp.

                                                             {
    VectorDouble res;
    res.resize(2);
    res[0] = -exp(-100. * (sqr(x) + sqr(y))) *
             (200. * x * cos(M_PI * x) + M_PI * sin(M_PI * x)) * cos(M_PI * y);
    res[1] = -exp(-100. * (sqr(x) + sqr(y))) *
             (200. * y * cos(M_PI * y) + M_PI * sin(M_PI * y)) * cos(M_PI * x);
    return res;
  }

◆ assembleSystem()

MoFEMErrorCode MixedPoisson::assembleSystem ( )

private

[Create common data]

[Assemble system]

Examples: mixed_poisson.cpp.

Definition at line 229 of file mixed_poisson.cpp.

                                            {
  MoFEMFunctionBegin;
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
  pipeline_mng->getDomainLhsFE().reset();
  pipeline_mng->getDomainRhsFE().reset();
  pipeline_mng->getOpDomainRhsPipeline().clear();
  pipeline_mng->getOpDomainLhsPipeline().clear();
 
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpCalculateHOJacForFace(jac_ptr));
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
  pipeline_mng->getOpDomainLhsPipeline().push_back(new OpMakeHdivFromHcurl());
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpSetContravariantPiolaTransformOnFace2D(jac_ptr));
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpSetInvJacHcurlFace(inv_jac_ptr));
  pipeline_mng->getOpDomainLhsPipeline().push_back(new OpSetHOWeightsOnFace());
 
  auto beta = [](const double, const double, const double) { return 1; };
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpHdivHdiv("FLUX", "FLUX", beta));
  auto unity = []() { return 1; };
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpHdivU("FLUX", "U", unity, true));
  auto source = [&](const double x, const double y, const double z) {
    return -sourceFunction(x, y, z);
  };
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpDomainSource("U", source));
  MoFEMFunctionReturn(0);
}

◆ checkError()

MoFEMErrorCode MixedPoisson::checkError ( int iter_num = 0 )

private

[Solve and refine loop]

[Check error]

Examples: mixed_poisson.cpp.

Definition at line 355 of file mixed_poisson.cpp.

                                                    {
  MoFEMFunctionBegin;
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
  pipeline_mng->getDomainLhsFE().reset();
  pipeline_mng->getDomainRhsFE().reset();
  pipeline_mng->getOpDomainRhsPipeline().clear();
 
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
 
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpCalculateHOJacForFace(jac_ptr));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
  pipeline_mng->getOpDomainRhsPipeline().push_back(new OpMakeHdivFromHcurl());
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpSetContravariantPiolaTransformOnFace2D(jac_ptr));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpSetInvJacHcurlFace(inv_jac_ptr));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpSetInvJacL2ForFace(inv_jac_ptr));
  pipeline_mng->getOpDomainRhsPipeline().push_back(new OpSetHOWeightsOnFace());
 
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpCalculateScalarFieldValues("U", commonDataPtr->approxVals));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpCalculateScalarFieldGradient<2>("U",
                                            commonDataPtr->approxValsGrad));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpCalculateHVecVectorField<3>("FLUX", commonDataPtr->approxFlux));
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpError(commonDataPtr, mField));
 
  CHKERR VecZeroEntries(commonDataPtr->petscVec);
  CHKERR VecGhostUpdateBegin(commonDataPtr->petscVec, INSERT_VALUES,
                             SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(commonDataPtr->petscVec, INSERT_VALUES,
                           SCATTER_FORWARD);
  CHKERR pipeline_mng->loopFiniteElements();
  CHKERR VecAssemblyBegin(commonDataPtr->petscVec);
  CHKERR VecAssemblyEnd(commonDataPtr->petscVec);
  CHKERR VecGhostUpdateBegin(commonDataPtr->petscVec, ADD_VALUES,
                             SCATTER_REVERSE);
  CHKERR VecGhostUpdateEnd(commonDataPtr->petscVec, ADD_VALUES,
                           SCATTER_REVERSE);
  CHKERR VecGhostUpdateBegin(commonDataPtr->petscVec, INSERT_VALUES,
                             SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(commonDataPtr->petscVec, INSERT_VALUES,
                           SCATTER_FORWARD);
  const double *array;
  CHKERR VecGetArrayRead(commonDataPtr->petscVec, &array);
  MOFEM_LOG("EXAMPLE", Sev::inform)
      << "Global error L2 norm: " << std::sqrt(array[0]);
  MOFEM_LOG("EXAMPLE", Sev::inform)
      << "Global error H1 seminorm: " << std::sqrt(array[1]);
  MOFEM_LOG("EXAMPLE", Sev::inform)
      << "Global error indicator: " << std::sqrt(array[2]);
  MOFEM_LOG("EXAMPLE", Sev::inform)
      << "Total number of elements: " << (int)array[3];
 
  errorIndicatorIntegral = array[2];
  totalElementNumber = (int)array[3];
  CHKERR VecRestoreArrayRead(commonDataPtr->petscVec, &array);
 
  std::vector<Tag> tag_handles;
  tag_handles.resize(4);
  CHKERR getTagHandle(mField, "ERROR_L2_NORM", MB_TYPE_DOUBLE, tag_handles[0]);
  CHKERR getTagHandle(mField, "ERROR_H1_SEMINORM", MB_TYPE_DOUBLE,
                      tag_handles[1]);
  CHKERR getTagHandle(mField, "ERROR_INDICATOR", MB_TYPE_DOUBLE,
                      tag_handles[2]);
  CHKERR getTagHandle(mField, "ORDER", MB_TYPE_INTEGER, tag_handles[3]);
 
  ParallelComm *pcomm =
      ParallelComm::get_pcomm(&mField.get_moab(), MYPCOMM_INDEX);
  if (pcomm == NULL)
    SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY, "Communicator not set");
 
  tag_handles.push_back(pcomm->part_tag());
  std::ostringstream strm;
  strm << "error_" << iter_num << ".h5m";
  CHKERR mField.get_moab().write_file(strm.str().c_str(), "MOAB",
                                      "PARALLEL=WRITE_PART", 0, 0,
                                      tag_handles.data(), tag_handles.size());
  MoFEMFunctionReturn(0);
}

◆ createCommonData()

MoFEMErrorCode MixedPoisson::createCommonData ( )

private

[Set integration rule]

[Create common data]

Examples: mixed_poisson.cpp.

Definition at line 211 of file mixed_poisson.cpp.

                                              {
  MoFEMFunctionBegin;
  commonDataPtr = boost::make_shared<CommonData>();
  PetscInt ghosts[4] = {0, 1, 2, 3};
  if (!mField.get_comm_rank())
    commonDataPtr->petscVec =
        createGhostVector(mField.get_comm(), 4, 4, 0, ghosts);
  else
    commonDataPtr->petscVec =
        createGhostVector(mField.get_comm(), 0, 4, 4, ghosts);
  commonDataPtr->approxVals = boost::make_shared<VectorDouble>();
  commonDataPtr->approxValsGrad = boost::make_shared<MatrixDouble>();
  commonDataPtr->approxFlux = boost::make_shared<MatrixDouble>();
  MoFEMFunctionReturn(0);
}

◆ getTagHandle()

static MoFEMErrorCode MixedPoisson::getTagHandle	(	MoFEM::Interface &	m_field,
		const char *	name,
		DataType	type,
		Tag &	tag_handle
	)

inlinestaticprivate

[Source function]

Examples: mixed_poisson.cpp.

Definition at line 77 of file mixed_poisson.cpp.

                                                      {
    MoFEMFunctionBegin;
    int int_val = 0;
    double double_val = 0;
    switch (type) {
    case MB_TYPE_INTEGER:
      CHKERR m_field.get_moab().tag_get_handle(
          name, 1, type, tag_handle, MB_TAG_CREAT | MB_TAG_SPARSE, &int_val);
      break;
    case MB_TYPE_DOUBLE:
      CHKERR m_field.get_moab().tag_get_handle(
          name, 1, type, tag_handle, MB_TAG_CREAT | MB_TAG_SPARSE, &double_val);
      break;
    default:
      SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "Wrong data type %d for tag", type);
    }
    MoFEMFunctionReturn(0);
  }

◆ outputResults()

MoFEMErrorCode MixedPoisson::outputResults ( int iter_num = 0 )

private

[Check error]

[Output results]

Examples: mixed_poisson.cpp.

Definition at line 445 of file mixed_poisson.cpp.

                                                       {
  MoFEMFunctionBegin;
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
  pipeline_mng->getDomainLhsFE().reset();
 
  using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>;
 
  auto post_proc_fe =
      boost::make_shared<PostProcEle>(mField);
 
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
  
  post_proc_fe->getOpPtrVector().push_back(
      new OpCalculateHOJacForFace(jac_ptr));
  post_proc_fe->getOpPtrVector().push_back(
      new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
  post_proc_fe->getOpPtrVector().push_back(new OpMakeHdivFromHcurl());
  post_proc_fe->getOpPtrVector().push_back(
      new OpSetContravariantPiolaTransformOnFace2D(jac_ptr));
  post_proc_fe->getOpPtrVector().push_back(
      new OpSetInvJacHcurlFace(inv_jac_ptr));
 
  auto u_ptr = boost::make_shared<VectorDouble>();
  auto flux_ptr = boost::make_shared<MatrixDouble>();
  post_proc_fe->getOpPtrVector().push_back(
      new OpCalculateScalarFieldValues("U", u_ptr));
  post_proc_fe->getOpPtrVector().push_back(
      new OpCalculateHVecVectorField<3>("FLUX", flux_ptr));
 
  using OpPPMap = OpPostProcMapInMoab<3, 3>;
 
  post_proc_fe->getOpPtrVector().push_back(
 
      new OpPPMap(post_proc_fe->getPostProcMesh(),
                  post_proc_fe->getMapGaussPts(),
 
                  OpPPMap::DataMapVec{{"U", u_ptr}},
 
                  OpPPMap::DataMapMat{{"FLUX", flux_ptr}},
 
                  OpPPMap::DataMapMat{},
 
                  OpPPMap::DataMapMat{}
 
                  )
 
  );
 
  pipeline_mng->getDomainRhsFE() = post_proc_fe;
  CHKERR pipeline_mng->loopFiniteElements();
 
  std::ostringstream strm;
  strm << "out_" << iter_num << ".h5m";
  CHKERR post_proc_fe->writeFile(strm.str().c_str());
  MoFEMFunctionReturn(0);
}

◆ readMesh()

MoFEMErrorCode MixedPoisson::readMesh ( )

private

[Run programme]

[Read mesh]

Examples: mixed_poisson.cpp.

Definition at line 153 of file mixed_poisson.cpp.

                                      {
  MoFEMFunctionBegin;
  CHKERR mField.getInterface(simpleInterface);
  CHKERR simpleInterface->getOptions();
  CHKERR simpleInterface->loadFile();
  MoFEMFunctionReturn(0);
}

◆ refineOrder()

MoFEMErrorCode MixedPoisson::refineOrder ( )

private

[Solve]

[Refine]

Examples: mixed_poisson.cpp.

Definition at line 288 of file mixed_poisson.cpp.

                                         {
  MoFEMFunctionBegin;
  Tag th_error_ind, th_order;
  CHKERR getTagHandle(mField, "ERROR_INDICATOR", MB_TYPE_DOUBLE, th_error_ind);
  CHKERR getTagHandle(mField, "ORDER", MB_TYPE_INTEGER, th_order);
 
  std::vector<Range> refinement_levels;
  refinement_levels.resize(refIterNum + 1);
  for (auto ent : domainEntities) {
    double err_indic = 0;
    CHKERR mField.get_moab().tag_get_data(th_error_ind, &ent, 1, &err_indic);
    int order, new_order;
    CHKERR mField.get_moab().tag_get_data(th_order, &ent, 1, &order);
    new_order = order + 1;
    Range refined_ents;
    if (err_indic > errorIndicatorIntegral / totalElementNumber) {
      refined_ents.insert(ent);
      Range adj;
      CHKERR mField.get_moab().get_adjacencies(&ent, 1, 1, false, adj,
                                               moab::Interface::UNION);
      refined_ents.merge(adj);
      refinement_levels[new_order - baseOrder].merge(refined_ents);
      CHKERR mField.get_moab().tag_set_data(th_order, &ent, 1, &new_order);
    }
  }
 
  for (int ll = 1; ll < refinement_levels.size(); ll++) {
    CHKERR mField.getInterface<CommInterface>()->synchroniseEntities(
        refinement_levels[ll]);
    CHKERR mField.set_field_order(refinement_levels[ll], "FLUX",
                                  baseOrder + ll);
    CHKERR mField.set_field_order(refinement_levels[ll], "U",
                                  baseOrder + ll - 1);
  }
 
  CHKERR mField.getInterface<CommInterface>()->synchroniseFieldEntities("FLUX");
  CHKERR mField.getInterface<CommInterface>()->synchroniseFieldEntities("U");
  CHKERR mField.build_fields();
  CHKERR mField.build_finite_elements();
  CHKERR mField.build_adjacencies(simpleInterface->getBitRefLevel());
  mField.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_TRUE;
  CHKERR DMSetUp_MoFEM(simpleInterface->getDM());
  MoFEMFunctionReturn(0);
}

◆ runProblem()

MoFEMErrorCode MixedPoisson::runProblem ( )

[Run programme]

Examples: mixed_poisson.cpp.

Definition at line 141 of file mixed_poisson.cpp.

                                        {
  MoFEMFunctionBegin;
  CHKERR readMesh();
  CHKERR setupProblem();
  CHKERR setIntegrationRules();
  CHKERR createCommonData();
  CHKERR solveRefineLoop();
  MoFEMFunctionReturn(0);
}

◆ setIntegrationRules()

MoFEMErrorCode MixedPoisson::setIntegrationRules ( )

private

[Set up problem]

[Set integration rule]

Examples: mixed_poisson.cpp.

Definition at line 197 of file mixed_poisson.cpp.

                                                 {
  MoFEMFunctionBegin;
 
  auto rule = [](int, int, int p) -> int { return 2 * p + 1; };
 
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
  CHKERR pipeline_mng->setDomainLhsIntegrationRule(rule);
  CHKERR pipeline_mng->setDomainRhsIntegrationRule(rule);
 
  MoFEMFunctionReturn(0);
}

◆ setupProblem()

MoFEMErrorCode MixedPoisson::setupProblem ( )

private

[Read mesh]

[Set up problem]

Examples: mixed_poisson.cpp.

Definition at line 163 of file mixed_poisson.cpp.

                                          {
  MoFEMFunctionBegin;
  // Note that in 2D case HDIV and HCURL spaces are isomorphic, and therefore
  // only base for HCURL has been implemented in 2D. Base vectors for HDIV space
  // are be obtained after rotation of HCURL base vectors by a right angle
  CHKERR simpleInterface->addDomainField("FLUX", HCURL, DEMKOWICZ_JACOBI_BASE,
                                         1);
  // We use AINSWORTH_LEGENDRE_BASE since DEMKOWICZ_JACOBI_BASE for triangle
  // is not yet implemented for L2 space. For quads DEMKOWICZ_JACOBI_BASE and
  // AINSWORTH_LEGENDRE_BASE are construcreed in the same way
  CHKERR simpleInterface->addDomainField("U", L2, AINSWORTH_LEGENDRE_BASE, 1);
 
  refIterNum = 0;
  CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-ref_iter_num", &refIterNum,
                            PETSC_NULL);
  baseOrder = 2;
  CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-base_order", &baseOrder,
                            PETSC_NULL);
  CHKERR simpleInterface->setFieldOrder("FLUX", baseOrder);
  CHKERR simpleInterface->setFieldOrder("U", baseOrder - 1);
  CHKERR simpleInterface->setUp();
 
  CHKERR mField.get_moab().get_entities_by_dimension(0, 2, domainEntities,
                                                     false);
  Tag th_order;
  CHKERR getTagHandle(mField, "ORDER", MB_TYPE_INTEGER, th_order);
  for (auto ent : domainEntities) {
    CHKERR mField.get_moab().tag_set_data(th_order, &ent, 1, &baseOrder);
  }
  MoFEMFunctionReturn(0);
}

◆ solveRefineLoop()

MoFEMErrorCode MixedPoisson::solveRefineLoop ( )

private

[Refine]

[Solve and refine loop]

Examples: mixed_poisson.cpp.

Definition at line 335 of file mixed_poisson.cpp.

                                             {
  MoFEMFunctionBegin;
  CHKERR assembleSystem();
  CHKERR solveSystem();
  CHKERR checkError();
  CHKERR outputResults();
 
  for (int ii = 0; ii < refIterNum; ii++) {
    CHKERR refineOrder();
 
    CHKERR assembleSystem();
    CHKERR solveSystem();
    CHKERR checkError(ii + 1);
    CHKERR outputResults(ii + 1);
  }
  MoFEMFunctionReturn(0);
}

◆ solveSystem()

MoFEMErrorCode MixedPoisson::solveSystem ( )

private

[Assemble system]

[Solve]

Examples: mixed_poisson.cpp.

Definition at line 267 of file mixed_poisson.cpp.

                                         {
  MoFEMFunctionBegin;
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
  auto solver = pipeline_mng->createKSP();
  CHKERR KSPSetFromOptions(solver);
  CHKERR KSPSetUp(solver);
 
  auto dm = simpleInterface->getDM();
  auto D = createDMVector(dm);
  auto F = vectorDuplicate(D);
  CHKERR VecZeroEntries(D);
 
  CHKERR KSPSolve(solver, F, D);
  CHKERR VecGhostUpdateBegin(D, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(D, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR DMoFEMMeshToLocalVector(dm, D, INSERT_VALUES, SCATTER_REVERSE);
  MoFEMFunctionReturn(0);
}

◆ sourceFunction()

static double MixedPoisson::sourceFunction	(	const double	x,
		const double	y,
		const double	z
	)

inlinestaticprivate

[Analytical function gradient]

[Source function]

Examples: mixed_poisson.cpp.

Definition at line 67 of file mixed_poisson.cpp.

                                                                               {
    return -exp(-100. * (sqr(x) + sqr(y))) *
           (400. * M_PI *
                (x * cos(M_PI * y) * sin(M_PI * x) +
                 y * cos(M_PI * x) * sin(M_PI * y)) +
            2. * (20000. * (sqr(x) + sqr(y)) - 200. - sqr(M_PI)) *
                cos(M_PI * x) * cos(M_PI * y));
  }

Member Data Documentation

◆ baseOrder

int MixedPoisson::baseOrder

private

Examples: mixed_poisson.cpp.

Definition at line 43 of file mixed_poisson.cpp.

◆ commonDataPtr

boost::shared_ptr<CommonData> MixedPoisson::commonDataPtr

private

Examples: mixed_poisson.cpp.

Definition at line 125 of file mixed_poisson.cpp.

◆ domainEntities

Range MixedPoisson::domainEntities

private

Examples: mixed_poisson.cpp.

Definition at line 40 of file mixed_poisson.cpp.

◆ errorIndicatorIntegral

double MixedPoisson::errorIndicatorIntegral

private

Examples: mixed_poisson.cpp.

Definition at line 41 of file mixed_poisson.cpp.

◆ mField

MoFEM::Interface& MixedPoisson::mField

private

Examples: mixed_poisson.cpp.

Definition at line 37 of file mixed_poisson.cpp.

◆ refIterNum

int MixedPoisson::refIterNum

private

Examples: mixed_poisson.cpp.

Definition at line 44 of file mixed_poisson.cpp.

◆ simpleInterface

Simple* MixedPoisson::simpleInterface

private

Examples: mixed_poisson.cpp.

Definition at line 38 of file mixed_poisson.cpp.

◆ totalElementNumber

int MixedPoisson::totalElementNumber

private

Examples: mixed_poisson.cpp.

Definition at line 42 of file mixed_poisson.cpp.

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

users_modules/tutorials/mix-0/mixed_poisson.cpp

Classes

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MixedPoisson()

Member Function Documentation

◆ analyticalFunction()

◆ analyticalFunctionGrad()

◆ assembleSystem()

◆ checkError()

◆ createCommonData()

◆ getTagHandle()

◆ outputResults()

◆ readMesh()

◆ refineOrder()

◆ runProblem()

◆ setIntegrationRules()

◆ setupProblem()

◆ solveRefineLoop()

◆ solveSystem()

◆ sourceFunction()

Member Data Documentation

◆ baseOrder

◆ commonDataPtr

◆ domainEntities

◆ errorIndicatorIntegral

◆ mField

◆ refIterNum

◆ simpleInterface

◆ totalElementNumber