element_transpose.hpp serif/backend/mfem/operators/element_evaluation.hpp src/include/serif/backend/mfem/operators/volume_application.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/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 serif/backend/mfem/operators/element_evaluation.hpp src/include/serif/backend/mfem/operators/element_transpose.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/workspace.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 algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception 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/model/volume.hpp src/include/serif/backend/mfem/model/volume_actions.hpp src/include/serif/backend/mfem/operators/element_transpose.hpp src/include/serif/backend/mfem/operators/jacobian_operator.hpp src/include/serif/backend/mfem/operators/volume_application.hpp src/include/serif/backend/mfem/operators/volume_field_block.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 src/lib/serif/backend/mfem/model/system_selection.cpp serif serif::operators serif::operators::backend serif::operators::backend::mfem #pragmaonce #include"serif/backend/mfem/operators/element_evaluation.hpp" namespaceserif::operators::backend::mfem{ template<VolumeEquationEquationT,typenameContributionT> voidelement_jacobian_transpose( constElementData<EquationPlan<EquationT>::operand_count>&data, const::mfem::DenseMatrix&geometry_derivatives, constdiscretization::backend::mfem::DofMap::Elements&row, const::mfem::Vector&residual_dual, Workspace&workspace, std::array<::mfem::Vector,EquationPlan<EquationT>::column_count>&column_duals, ::mfem::Vector&displacement_dual ){ usingnamespaceutils::errors; usingPlan=EquationPlan<EquationT>; usingTest=OperandEvaluator<equations::TestOperandOf<ContributionT>>; for(::mfem::Vector&column:column_duals){ column.SetSize(0); } row.element_values(data.element,residual_dual,workspace.element_values); Test::evaluate(data.geometry->test_basis,workspace.element_values,workspace.test_values,workspace.basis); data.basis->reference_test_basis.evaluate(workspace.element_values,workspace.reference_test_values,workspace.basis); constintpoints=data.basis->rule->GetNPoints(); Plan::templatefor_each_input<ContributionT>([&]<typenameOperandT>(OperandT,std::size_tindex){ usingEvaluator=OperandEvaluator<OperandT>; workspace.operand_dual.SetSize(points*Evaluator::components); for(intpoint=0;point<points;++point){ for(intinput=0;input<Evaluator::components;++input){ doublevalue=0.0; for(intoutput=0;output<Test::components;++output){ constintrow_index=point*Test::components+output; value+=data.weighted_field_derivatives(row_index,Plan::templatecomponent_offset<OperandT>()+input)*workspace.test_values[row_index]; } workspace.operand_dual[point*Evaluator::components+input]=value; } } Evaluator::integrate(data.geometry->operand_basis[index],workspace.operand_dual,workspace.element_result,workspace.basis); auto&column=column_duals[Plan::templatecolumn_index<OperandT>()]; if(column.Size()==0){ column=workspace.element_result; }else{ column+=workspace.element_result; } }); workspace.displacement_dual.SetSize(12*points); //Usethesamestoredcolumnsandcomponentorderastheforwardgeometryaction. constintderivative_stride=geometry_derivatives.Height(); constdouble*derivatives=geometry_derivatives.GetData(); constdouble*test_values=workspace.reference_test_values.GetData(); double*displacement=workspace.displacement_dual.GetData(); for(intpoint=0;point<points;++point){ for(intinput=0;input<12;++input){ doublevalue=0.0; for(intoutput=0;output<Test::reference_components;++output){ constintrow_index=point*Test::reference_components+output; value+=derivatives[input*derivative_stride+row_index]*test_values[row_index]; } serif_ensure<SERiFErrorCode::NON_FINITE_VALUE_ERROR,NonFiniteValue>( std::isfinite(value),"Non-finitegeometryJacobiantransposeaction." ); displacement[12*point+input]=value; } } data.basis->displacement_basis.integrate(workspace.displacement_dual,displacement_dual,workspace.basis); } }//namespaceserif::operators::backend::mfem