235 auto add_domain_base_ops = [&](
auto &pipeline) {
236 auto jac_ptr = boost::make_shared<MatrixDouble>();
237 auto det_ptr = boost::make_shared<VectorDouble>();
243 auto add_domain_lhs_ops = [&](
auto &pipeline) {
245 "U",
"U", [](
const double,
const double,
const double) {
return 1; }));
248 auto add_domain_rhs_ops = [&](
auto &pipeline) {
253 add_domain_base_ops(pipeline_mng->getOpDomainLhsPipeline());
254 add_domain_base_ops(pipeline_mng->getOpDomainRhsPipeline());
255 add_domain_lhs_ops(pipeline_mng->getOpDomainLhsPipeline());
256 add_domain_rhs_ops(pipeline_mng->getOpDomainRhsPipeline());
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, 1 > OpDomainMass
FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< 1, 1 > OpDomainSource
Simple interface for fast problem set-up.
BcMapByBlockName & getBcMapByBlockName()
Get the bc map.
Calculate jacobian on Hex or other volume which is not simplex.
Set inverse jacobian to base functions.
PipelineManager interface.
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.
static double sourceFunction(const double x, const double y, const double z)
Pulse is infinitely short.