geometry_derivatives.hpp array cmath expected serif/backend/mfem/operators/element_data.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/rejection.hpp serif/backend/mfem/operators/workspace.hpp serif/operators/plan.hpp serif/backend/mfem/mapping/domain_mapper.hpp serif/backend/mfem/shared/parallel.hpp src/include/serif/backend/mfem/operators/volume_stages.hpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/tensor_basis.hpp serif/backend/mfem/mapping/compactification/compactification.hpp serif/backend/mfem/mapping/compactification/input_result.hpp serif/backend/mfem/mapping/compactification/kelvin.hpp serif/backend/mfem/mapping/compactification/variation.hpp serif/backend/mfem/mapping/context.hpp serif/backend/mfem/mapping/data.hpp serif/backend/mfem/mapping/domain_mapper.hpp serif/backend/mfem/mapping/status.hpp serif/backend/mfem/mapping/transformations.hpp serif/backend/mfem/mapping/variation.hpp serif/backend/mfem/mapping/volume_cache.hpp serif/backend/mfem/operators/element_data.hpp src/include/serif/backend/mfem/operators/geometry_derivatives.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/rejection.hpp serif/backend/mfem/operators/workspace.hpp serif/backend/mfem/shared/parallel.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/concepts.hpp serif/discretization/blocks/fields.hpp serif/discretization/blocks/traits.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/types.hpp serif/discretization/forms/forms.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/reconstruction/surface_power.hpp serif/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/meta/enums.hpp serif/meta/names.hpp serif/operators/concepts.hpp serif/operators/plan.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/error/nonfinite.hpp serif/utils/misc/compensated_sum.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/std_helper/variadic_cleaning.hpp serif/utils/misc/terminal/colors.hpp serif/utils/misc/traits/type_uniqueness.hpp serif/utils/types/strings.hpp serif/utils/types/type_list.hpp serif/utils/user_options/mapping/compactification/kelvin.hpp serif/utils/user_options/mapping/domain_mapper.hpp algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format iostream limits magic_enum/magic_enum_all.hpp memory mfem/fem/eltrans.hpp mfem/fem/fe/fe_base.hpp mfem/fem/fe/fe_h1.hpp mfem/fem/fe/fe_l2.hpp mfem/fem/fespace.hpp mfem/fem/intrules.hpp mfem/fem/pfespace.hpp mfem/linalg/densemat.hpp mfem/linalg/vector.hpp mfem/mesh/mesh.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility vector version src/compile_time_checks/serif/model/scalar_row.cpp src/compile_time_checks/serif/model/state.cpp src/compile_time_checks/serif/operators/computation_owners.cpp src/include/serif/backend/mfem/model/computations/interior.hpp src/include/serif/backend/mfem/model/computations/potential_infinity.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp src/include/serif/backend/mfem/model/computations/scalar_row.hpp src/include/serif/backend/mfem/model/computations/surface_potential.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/operators/geometry_derivatives.hpp src/include/serif/backend/mfem/operators/volume_computations.hpp src/include/serif/backend/mfem/operators/volume_stages.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.hpp src/include/serif/backend/mfem/solver/equilibrium.hpp src/include/serif/backend/mfem/solver/nonlinear/stellar_problem.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/solver.hpp serif serif::operators serif::operators::backend serif::operators::backend::mfem #pragmaonce #include<array> #include<cmath> #include<expected> #include"serif/backend/mfem/operators/element_data.hpp" #include"serif/backend/mfem/operators/operand_evaluation.hpp" #include"serif/backend/mfem/operators/rejection.hpp" #include"serif/backend/mfem/operators/workspace.hpp" #include"serif/operators/plan.hpp" #include"serif/backend/mfem/mapping/domain_mapper.hpp" #include"serif/backend/mfem/shared/parallel.hpp" namespaceserif::operators::backend::mfem{ template<VolumeEquationEquation,typenameContribution> [[nodiscard]]autogeometry_integrand_variation( constElementData<EquationPlan<Equation>::operand_count>&data, intpoint, constmapping::VolumeMappingContext&context, constmapping::VolumeMappingVariation&variation ){ usingPlan=EquationPlan<Equation>; usingTest=OperandEvaluator<equations::TestOperandOf<Contribution>>; typenameTest::Valueresult{}; for(intoutput=0;output<Test::components;++output){ constintrow=point*Test::components+output; doublevalue=data.weighted_residual[row]*variation.weight_variation/context.quadrature.physical_integration_weight; for(intaxis=0;axis<3;++axis){ value+=data.weighted_position_derivative(row,axis)*variation.mapping.physical_position_variation[axis]; } discretization::forms::component(result,output)=value; } Plan::templatefor_each_input<Contribution>([&]<typenameOperand>(Operand,std::size_tindex){ usingEvaluator=OperandEvaluator<Operand>; ifconstexpr(Evaluator::geometry_dependent){ constautodelta=Evaluator::geometry_variation( Evaluator::reference_value(*data.reference_operand_values[index],point),context.mapping,variation.mapping ); for(intoutput=0;output<Test::components;++output){ for(intinput=0;input<Evaluator::components;++input){ discretization::forms::component(result,output)+=data.weighted_field_derivatives( point*Test::components+output,Plan::templatecomponent_offset<Operand>()+input )*discretization::forms::component(delta,input); } } } }); returnresult; } template<VolumeEquationEquation,typenameContribution> [[nodiscard]]std::expected<::mfem::DenseMatrix,PreparationRejection>element_geometry_derivatives( constElementData<EquationPlan<Equation>::operand_count>&data, constdiscretization::backend::mfem::DofMap&displacement_binding, const::mfem::FiniteElementSpace&compactification_space, constmapping::DomainMapper&mapper, Workspace&workspace ){ usingTest=OperandEvaluator<equations::TestOperandOf<Contribution>>; constauto&element=*displacement_binding.space().GetFE(data.element); workspace.element_values=data.geometry->displacement; constautobase=mapping::element_displacement_data_from_element_vdofs(element,workspace.element_values); constmapping::ElementCompactificationDatacompactification( *compactification_space.GetFE(data.element),data.geometry->compactification ); constmapping::ElementMappingDatamapping_data{base,compactification}; ::mfem::IsoparametricTransformationtransformation; ::serif::utils::parallel::element_transformation(*displacement_binding.space().GetMesh(),data.element,transformation); constintpoints=data.basis->rule->GetNPoints(); ::mfem::DenseMatrixresult(points*Test::reference_components,12); constmapping::DMReferenceDisplacementSamplesample{workspace.displacement_value,workspace.displacement_gradient}; mapping::CompactificationPointDatacompactification_sample; ::mfem::DenseMatrixinverse_mesh_jacobian(3); for(intpoint=0;point<points;++point){ data.geometry->mapping.load(point,workspace.mapping_context); //Alltwelvedirectionssharetheinversefixed-meshJacobian,thecompactificationsample //andthephysicaltestvaluesatthispoint. ::mfem::Mult(workspace.mapping_context.quadrature.J_element_to_physical_inverse, workspace.mapping_context.mapping.mapping_jacobian,inverse_mesh_jacobian); constmapping::CompactificationPointData*fixed_sample=nullptr; if(workspace.mapping_context.mapping.compactified){ constautostatus=mapper.evaluate_compactification_coordinate( compactification,transformation,data.basis->rule->IntPoint(point), workspace.mapping_workspace,compactification_sample,&inverse_mesh_jacobian ); if(status!=mapping::MappingStatus::VALID){ returnstd::unexpected(PreparationRejection{ .reason=PreparationRejectionReason::MAPPING,.mapping_status=status, .element=data.element,.quadrature_point=point, .message="Invalidfixedcompactificationsampleforanequationderivative." }); } fixed_sample=&compactification_sample; } std::array<typenameTest::Value,Test::reference_components>tests; for(inttest_component=0;test_component<Test::reference_components;++test_component){ typenameTest::ReferenceValuereference{}; discretization::forms::component(reference,test_component)=1.0; tests[test_component]=Test::physical_value(reference,workspace.mapping_context.mapping); } for(intcomponent=0;component<12;++component){ workspace.displacement_value=0.0; workspace.displacement_gradient=0.0; if(component%4==0){ workspace.displacement_value[component/4]=1.0; }else{ workspace.displacement_gradient(component/4,component%4-1)=1.0; } constautostatus=mapper.evaluate_volume_direction_variation( mapping_data,base,transformation,data.basis->rule->IntPoint(point), workspace.mapping_context,workspace.mapping_workspace,workspace.mapping_variation,sample,fixed_sample, inverse_mesh_jacobian ); if(status!=mapping::MappingStatus::VALID){ returnstd::unexpected(PreparationRejection{ .reason=PreparationRejectionReason::MAPPING,.mapping_status=status, .element=data.element,.quadrature_point=point, .message="Invalidlocaldisplacementderivativeofanequation." }); } constautophysics=geometry_integrand_variation<Equation, Contribution>( data,point,workspace.mapping_context,workspace.mapping_variation ); for(inttest_component=0;test_component<Test::reference_components;++test_component){ typenameTest::ReferenceValuereference{}; discretization::forms::component(reference,test_component)=1.0; constauto&test=tests[test_component]; constautotest_variation=Test::geometry_variation(reference,workspace.mapping_context.mapping,workspace.mapping_variation.mapping); doublevalue=0.0; for(intoutput=0;output<Test::components;++output){ value+=discretization::forms::component(test,output)*discretization::forms::component(physics,output) +data.weighted_residual[point*Test::components+output]*discretization::forms::component(test_variation,output); } if(!std::isfinite(value)){ returnstd::unexpected(PreparationRejection{ .reason=PreparationRejectionReason::NON_FINITE,.element=data.element, .quadrature_point=point,.message="Non-finitelocaldisplacementderivativeofanequation." }); } result(point*Test::reference_components+test_component,component)=value; } } } returnresult; } }//namespaceserif::operators::backend::mfem