22#ifndef __RUNADAPTIVITY_HPP__
23#define __RUNADAPTIVITY_HPP__
106 double error,
int max_order)
131 PetscErrorCode
setOrder(
int set,
int up,
double error,
int max_order,
137 PetscErrorCode
runKsp(DM dm, Mat Aij, Vec
F, Vec
D);
139 PetscErrorCode
runKspAtLevel(DM dm, Mat Aij, Vec
F, Vec
D,
bool add_level);
156 PetscErrorCode
solveProblem(DM dm, Mat Aij, Vec
F, Vec
D,
int nb_ref_cycles,
157 int max_order,
double procent = 0.33,
Deprecated interface functions.
Set data structures of MG pre-conditioner via approximation orders.
Set order of finite element.
PetscErrorCode operator()()
SetOrderToElement(MoFEM::Interface &m_field, Tag th_order, int set, int up, double error, int max_order)
PetscErrorCode postProcess()
MoFEM::Interface & mField
PetscErrorCode preProcess()
PetscErrorCode postProcess()
PetscErrorCode operator()()
SetOrderToEntities(MoFEM::Interface &m_field, Tag th_order)
MoFEM::Interface & mField
PetscErrorCode preProcess()
PetscErrorCode getOptions()
get options from line command
PetscErrorCode createIsAtLevel(int kk, IS *is)
Set IS for levels.
PetscErrorCode buildProlongationOperator(PC pc, int verb=0)
PetscErrorCode destroyIsAtLevel(int kk, IS *is)
Destroy IS if internally created.
ShellPCMGSetUpViaApproxOrdersCtx(DM dm, Mat a, vector< IS > &mg_levels)
PetscErrorCode sortErrorVector(const double procent=0.33)
vector< int > vecOrderDofs
PetscErrorCode postProcFatPrims(const string &name)
Save post-processed data.
PetscErrorCode setOrder(int set, int up, double error, int max_order, int verb=0)
Set approximation order to entities.
PetscErrorCode printDisplacements()
Print displacements at selected points.
PetscErrorCode runKspAtLevel(DM dm, Mat Aij, Vec F, Vec D, bool add_level)
PetscErrorCode createIS(IS *is)
vector< Vec > solutionAtLevel
SolidShellPrismElement::CommonData & commonData
MoFEM::Interface & mField
PostProcFatPrismOnRefinedMesh postProcFatPrism
double procentErrorMax
max value at given precent
PetscErrorCode setUpOperators()
Set up operators for post-processing shell elements.
RunAdaptivity(MoFEM::Interface &m_field, DM dm, SolidShellPrismElement::CommonData &common_data, FEMethod &shell_lhs, FEMethod &shell_error, FEMethod &dirichlet_bc, FEMethod &fix_ents)
PetscErrorCode runKsp(DM dm, Mat Aij, Vec F, Vec D)
PetscErrorCode solveProblem(DM dm, Mat Aij, Vec F, Vec D, int nb_ref_cycles, int max_order, double procent=0.33, int verb=0)
Solve linear problem with p-adaptivity.
Common data for solid shell element.