L2EnrichedQuadFESpace::L2EnrichedQuadFESpace (shared_ptr<MeshAccess> ama,
const Flags & flags)
: FESpace (ama, flags) {
_k = int(flags.GetNumFlag ("order", 2));
evaluator =
make_shared<T_DifferentialOperator<DiffOpId<2>>>();
flux_evaluator =
make_shared<T_DifferentialOperator<DiffOpGradient<2>>>();
boundary_evaluator =
make_shared<T_DifferentialOperator<DiffOpIdBoundary<2>>>();
integrator = GetIntegrators() .
CreateBFI("mass", ma->GetDimension(),
make_shared<ConstantCoefficientFunction>(1));
}
MyHighOrderFESpace :: MyHighOrderFESpace (shared_ptr<MeshAccess> ama, const Flags & flags)
: FESpace (ama, flags)
{
cout << "Constructor of MyHighOrderFESpace" << endl;
order = int(flags.GetNumFlag ("order", 2));
// needed to draw solution function
evaluator = make_shared<T_DifferentialOperator<DiffOpId<2>>>();
flux_evaluator = make_shared<T_DifferentialOperator<DiffOpGradient<2>>>();
boundary_evaluator = make_shared<T_DifferentialOperator<DiffOpIdBoundary<2>>>();
integrator = GetIntegrators().CreateBFI("mass", ma->GetDimension(),
make_shared<ConstantCoefficientFunction>(1));
}
