volume_computations.hpp expected memory tuple utility serif/backend/mfem/operators/volume_stages.hpp serif/backend/mfem/operators/volume_inputs.hpp serif/utils/dependencies/dependencies.hpp 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/scalar_row.hpp serif/backend/mfem/discretization/basis_evaluation.hpp serif/backend/mfem/discretization/coefficient_tracking.hpp serif/backend/mfem/discretization/dof_map.hpp serif/backend/mfem/discretization/quadrature/concepts.hpp serif/backend/mfem/discretization/quadrature/configure.hpp serif/backend/mfem/discretization/quadrature/mfem_resolver.hpp serif/backend/mfem/discretization/quadrature/mfem_static_order.hpp serif/backend/mfem/discretization/quadrature/order.hpp serif/backend/mfem/discretization/reconstruction.hpp serif/backend/mfem/discretization/tensor_basis.hpp serif/backend/mfem/discretization/term_binding.hpp serif/backend/mfem/discretization/term_coefficients.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 serif/backend/mfem/operators/element_evaluation.hpp serif/backend/mfem/operators/element_preparation.hpp serif/backend/mfem/operators/geometry_derivatives.hpp serif/backend/mfem/operators/operand_evaluation.hpp serif/backend/mfem/operators/preparation.hpp serif/backend/mfem/operators/quadrature.hpp serif/backend/mfem/operators/rejection.hpp src/include/serif/backend/mfem/operators/volume_computations.hpp serif/backend/mfem/operators/volume_data.hpp serif/backend/mfem/operators/volume_inputs.hpp serif/backend/mfem/operators/volume_stages.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/runtime.hpp serif/discretization/blocks/traits.hpp serif/discretization/domain/concepts.hpp serif/discretization/domain/ids/boundary.hpp serif/discretization/domain/ids/concepts.hpp serif/discretization/domain/ids/domain.hpp serif/discretization/domain/ids/lists/concepts.hpp serif/discretization/domain/ids/lists/lists.hpp serif/discretization/domain/ids/lists/traits.hpp serif/discretization/domain/ids/validation/concepts.hpp serif/discretization/domain/ids/validation/unique_ids.hpp serif/discretization/domain/physical_domains.hpp serif/discretization/domain/relation/concepts.hpp serif/discretization/domain/relation/lists/concepts.hpp serif/discretization/domain/relation/lists/relation_list.hpp serif/discretization/domain/relation/lists/traits.hpp serif/discretization/domain/relation/relations.hpp serif/discretization/domain/relation/validation/concepts.hpp serif/discretization/domain/relation/validation/validation.hpp serif/discretization/domain/resolver.hpp serif/discretization/domain/schema/concepts.hpp serif/discretization/domain/schema/domain_schema.hpp serif/discretization/domain/schema/traits.hpp serif/discretization/domain/types.hpp serif/discretization/forms/concept.hpp serif/discretization/forms/form_traits.hpp serif/discretization/forms/forms.hpp serif/discretization/forms/operands.hpp serif/discretization/forms/runtime.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/discretization/quadrature/policy.hpp serif/discretization/quadrature/query.hpp serif/discretization/quadrature/rules.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/contributions.hpp serif/operators/plan.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_cache.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependencies.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/cache_storage.hpp serif/utils/dependencies/detail/cache_traits.hpp serif/utils/dependencies/detail/graph.hpp serif/utils/dependencies/detail/result.hpp serif/utils/dependencies/detail/schema.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.hpp serif/utils/dependencies/observers.hpp serif/utils/dependencies/preparation.hpp serif/utils/dependencies/resource.hpp serif/utils/dependencies/snapshot.hpp serif/utils/error/codes.hpp serif/utils/error/errors.hpp serif/utils/error/nonfinite.hpp serif/utils/format/strings.hpp serif/utils/misc/compensated_sum.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/filter.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/concepts/type_uniqueness.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 atomic cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional 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/geom.hpp mfem/fem/gridfunc.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 variant 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/volume_computations.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<expected> #include<memory> #include<tuple> #include<utility> #include"serif/backend/mfem/operators/volume_stages.hpp" #include"serif/backend/mfem/operators/volume_inputs.hpp" #include"serif/utils/dependencies/dependencies.hpp" namespaceserif::operators::backend::mfem { template<std::size_tContribution,typenameOperandT,typenameScope=void> usingVolumeOperandResult= std::conditional_t<OperandEvaluator<OperandT>::geometry_dependent,typenamevolume_keys::Scoped<Scope>::templatePhysicalValues<Contribution,OperandT>, typenamevolume_keys::Scoped<Scope>::templateReferenceValues<Contribution,OperandT>>; template<std::size_tContribution,typenameOperandT,VolumeEquationEquationT,typenameSchemaT,typenameScope=void> [[nodiscard]]automake_operand_computations(conststd::shared_ptr<VolumeScratch>&scratch) { usingnamespaceutils::dependencies; usingBasis=ElementBatch<ElementBasis<EquationPlan<EquationT>::operand_count>>; usingGeometry=ElementBatch<ElementGeometry<EquationPlan<EquationT>::operand_count>>; usingBinding=discretization::backend::mfem::TermBinding<typenameOperandT::term_type,SchemaT>; usingReferenceValues=typenamevolume_keys::Scoped<Scope>::templateReferenceValues<Contribution,OperandT>; constautosample_values=make_computation<ReferenceValues,InputList<typenamevolume_keys::Scoped<Scope>::templateBasis<Contribution>,typenameOperandT::term_type>>( [scratch](constBasis&basis,constCoefficientView<Binding>&input) { returnsample_contribution_operand<EquationT, OperandT>(basis,input,scratch->workspace); }); ifconstexpr(OperandEvaluator<OperandT>::geometry_dependent) { constautomap_values=make_computation<typenamevolume_keys::Scoped<Scope>::templatePhysicalValues<Contribution,OperandT>, InputList<ReferenceValues,typenamevolume_keys::Scoped<Scope>::templateGeometry<Contribution>>>( [scratch](constOperandSamples&samples,constGeometry&geometry) { returnmap_contribution_operand<OperandT>(samples,geometry,scratch->workspace); }); returnstd::tuple{sample_values,map_values}; } else { returnstd::tuple{sample_values}; } } template<std::size_tIndex,typenameScope,VolumeEquationEquationT,typenameSchemaT,typenameContributionT,typename...Operands> [[nodiscard]]automake_contribution_computations(conststd::shared_ptr<constVolumeConfiguration<EquationT, SchemaT>>&configuration, conststd::shared_ptr<VolumeScratch>&scratch,ContributionBinding<ContributionT>contribution, utils::types::TypeList<Operands...>) { usingnamespaceutils::dependencies; usingPlan=EquationPlan<EquationT>; usingBasis=ElementBatch<ElementBasis<Plan::operand_count>>; usingGeometry=ElementBatch<ElementGeometry<Plan::operand_count>>; usingDisplacementBinding=typenameVolumeConfiguration<EquationT, SchemaT>::DisplacementBinding; constautoreference_basis=make_computation<typenamevolume_keys::Scoped<Scope>::templateBasis<Index>,InputList<typenamevolume_keys::Scoped<Scope>::Bindings>>( [configuration,scratch,contribution](constBasisBindings<Plan::column_count>&bindings) { returnsample_contribution_basis(*configuration,contribution,bindings,scratch->basis); }); constautogeometry=make_computation<typenamevolume_keys::Scoped<Scope>::templateGeometry<Index>, InputList<typenamevolume_keys::Scoped<Scope>::templateBasis<Index>,DisplacementCoefficients,CompactificationCoefficients>>( [configuration,scratch](constBasis&basis,constCoefficientView<DisplacementBinding>&displacement,constMappingCoefficients&compactification) { returnmap_contribution_geometry<EquationT, SchemaT, ContributionT>(*configuration,basis,displacement,compactification,scratch->workspace); }); constautooperand_values=std::tuple_cat(make_operand_computations<Index, Operands, EquationT, SchemaT, Scope>(scratch)...); constautoequation_values=make_computation<typenamevolume_keys::Scoped<Scope>::templateQuadratureLinearization<Index>, InputList<typenamevolume_keys::Scoped<Scope>::templateBasis<Index>, typenamevolume_keys::Scoped<Scope>::templateGeometry<Index>,VolumeOperandResult<Index, Operands, Scope>...>>( [scratch,contribution](constBasis&basis,constGeometry&geometry,constSamplesFor<Operands>&...samples) { returnevaluate_contribution_equation<EquationT,ContributionT,Operands...>(contribution,basis,geometry,scratch->workspace,samples...); }); constautogeometry_derivatives=make_computation<typenamevolume_keys::Scoped<Scope>::templateGeometryDerivatives<Index>, InputList<typenamevolume_keys::Scoped<Scope>::templateQuadratureLinearization<Index>>>( [configuration,scratch](constContributionData<Plan::operand_count>&data){ returndifferentiate_contribution_geometry<EquationT, SchemaT, ContributionT>(*configuration,data,scratch->workspace); }); returnstd::tuple_cat(std::tuple{reference_basis,geometry},operand_values,std::tuple{equation_values,geometry_derivatives}); } template<typenameScope=void,VolumeEquationEquationT,typenameSchemaT,std::size_t...Indices> [[nodiscard]]autoconnect_volume_contributions(conststd::shared_ptr<constVolumeConfiguration<EquationT, SchemaT>>&configuration,std::index_sequence<Indices...>) { usingnamespaceutils::dependencies; usingPlan=EquationPlan<EquationT>; constautoscratch=std::make_shared<VolumeScratch>(); constautocontributions=bind_contributions(configuration->equation); constautocalculations=std::tuple_cat(make_contribution_computations<Indices, Scope>( configuration,scratch,std::get<Indices>(contributions),typenamestd::tuple_element_t<Indices,equations::ContributionsOf<EquationT>>::Operands{})...); constautovolume= make_computation<typenamevolume_keys::Scoped<Scope>::Linearization, InputList<typenamevolume_keys::Scoped<Scope>::Bindings,typenamevolume_keys::Scoped<Scope>::templateQuadratureLinearization<Indices>..., typenamevolume_keys::Scoped<Scope>::templateGeometryDerivatives<Indices>...>>( [configuration,scratch](constBasisBindings<Plan::column_count>&bindings,constauto&...results) { //Thefirstgroupcontainsthephysicalequations;thesecondtheirlocalgeometryderivatives. constautovalues=std::tie(results...); autovolume=collect_volume_results(configuration,bindings,scratch->workspace,std::get<Indices>(values)...); if(volume){ volume->geometry_derivatives={std::get<sizeof...(Indices)+Indices>(values)...}; } returnvolume; }); returnstd::tuple_cat(calculations,std::tuple{volume}); } template<VolumeEquationEquationT,std::size_tIndex,typenameScope=void,typenameCacheT> voidreport_contribution_preparation(constCacheT&cache,NumericalWork&report) { usingPlan=EquationPlan<EquationT>; usingContribution=std::tuple_element_t<Index,equations::ContributionsOf<EquationT>>; constauto&basis=cache.templateget<typenamevolume_keys::Scoped<Scope>::templateBasis<Index>>(); constauto&status=cache.report(); std::size_tpoints=0; for(constauto&element:basis) { points+=element->rule->GetNPoints(); } if(status.templatestatus<typenamevolume_keys::Scoped<Scope>::templateGeometry<Index>>()==utils::dependencies::PreparationStatus::computed) { report.geometry_points+=points; } if(status.templatestatus<typenamevolume_keys::Scoped<Scope>::templateQuadratureLinearization<Index>>()==utils::dependencies::PreparationStatus::computed) { report.physics_points+=points; } else { report.reused_elements+=basis.size(); } Plan::templatefor_each_input<Contribution>( [&]<typenameOperandT>(OperandT,std::size_t) { if(status.templatestatus<typenamevolume_keys::Scoped<Scope>::templateReferenceValues<Index,OperandT>>()==utils::dependencies::PreparationStatus::computed) { report.sampled_operands+=basis.size(); } }); } template<typenameScope=void,VolumeEquationEquationT,typenameSchemaT> [[nodiscard]]automake_volume_equation_computations(std::shared_ptr<constEquationT>equation,constVolumeDiscretization<typename EquationT::Test, SchemaT>&discretization, VolumeOptionsoptions={}) { utils::errors::serif_ensure<utils::errors::SERiFErrorCode::INVALID_ARGUMENT>(equation!=nullptr,"Volumecomputationsrequireanowningequation;supplyanon-nullequationorretainedmodelmember."); validate_volume_discretization(discretization); constautoconfiguration=std::make_shared<constVolumeConfiguration<EquationT,SchemaT>>( VolumeConfiguration<EquationT, SchemaT>{.equation=std::move(equation),.discretization=discretization,.options=std::move(options)}); returnconnect_volume_contributions<Scope>(configuration,std::make_index_sequence<std::tuple_size_v<equations::ContributionsOf<EquationT>>>{}); } template<typenameScope=void,VolumeEquationEquationT,typenameSchemaT> [[nodiscard]]automake_volume_equation_computations(EquationTequation,constVolumeDiscretization<typename EquationT::Test, SchemaT>&discretization,VolumeOptionsoptions={}) { returnmake_volume_equation_computations<Scope>(std::make_shared<constEquationT>(std::move(equation)),discretization,std::move(options)); } template<typenameScope=void,VolumeEquationEquationT,typenameSchemaT,typenameInputsT> [[nodiscard]]automake_volume_equation_cache(EquationTequation,constVolumeDiscretization<typename EquationT::Test, SchemaT>&discretization,constInputsT&inputs, VolumeOptionsoptions={}) { returnstd::apply( [&](auto...calculations) { returnutils::dependencies::make_computation_cache<PreparationRejection>(inputs,std::move(calculations)...); }, make_volume_equation_computations<Scope>(std::move(equation),discretization,std::move(options))); } template<typenameCacheT,typenameScope=void> usingVolumeEquationOf=typenamestd::remove_cvref_t<decltype(std::declval<constCacheT&>().templateget<typenamevolume_keys::Scoped<Scope>::Linearization>())>::Equation; template<typenameScope=void,typenameCacheT,typenameDiscretizationT,typenameCoefficientsT> [[nodiscard]]std::expected<void,PreparationRejection>try_prepare_volume(CacheT&cache,constDiscretizationT&discretization,constCoefficientsT&term_coefficients, const::mfem::Vector&displacement_coefficients, PreparationModemode=PreparationMode::reuse_unchanged) { returntry_prepare_validated(cache, [&]() { constautoinputs= make_volume_inputs<VolumeEquationOf<CacheT, Scope>,Scope>(discretization,term_coefficients,displacement_coefficients); returncache.try_prepare(inputs,{},mode); }); } template<typenameScope=void,typenameCacheT,typenameDiscretizationT,typenameCoefficientsT> voidprepare_volume(CacheT&cache,constDiscretizationT&discretization,constCoefficientsT&term_coefficients,const::mfem::Vector&displacement_coefficients, PreparationModemode=PreparationMode::reuse_unchanged) { usingnamespaceutils::errors; constautoresult=try_prepare_volume<Scope>(cache,discretization,term_coefficients,displacement_coefficients,mode); if(!result) { serif_error<std::domain_error,SERiFErrorCode::INVALID_ARGUMENT>("{}preparationrejectedatelement{},point{}:{}",VolumeEquationOf<CacheT, Scope>::name, result.error().element,result.error().quadrature_point,result.error().message); } } template<typenameScope=void,typenameCacheT> [[nodiscard]]NumericalWorkvolume_work_report(constCacheT&cache) { usingEquation=VolumeEquationOf<CacheT, Scope>; NumericalWorkreport; if(cache.templatecontains<typenamevolume_keys::Scoped<Scope>::Linearization>()) { [&]<std::size_t...Indices>(std::index_sequence<Indices...>) { (report_contribution_preparation<Equation, Indices, Scope>(cache,report),...); }(std::make_index_sequence<std::tuple_size_v<equations::ContributionsOf<Equation>>>{}); } returnreport; } }//namespaceserif::operators::backend::mfem