nodal_computations.hpp expected utility tuple serif/operators/backend/mfem/nodal_evaluation.hpp serif/operators/backend/mfem/nodal_inputs.hpp serif/operators/backend/mfem/preparation.hpp serif/utils/dependencies/dependencies.hpp src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp src/include/serif/model/computations/point_constraints.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/backend/mfem/coefficient_tracking.hpp serif/discretization/backend/mfem/dof_map.hpp serif/discretization/backend/mfem/dof_selection.hpp serif/discretization/backend/mfem/point_dofs.hpp serif/discretization/backend/mfem/term_binding.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/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/operands.hpp serif/discretization/forms/value.hpp serif/discretization/functional_space/space.hpp serif/equations/concepts.hpp serif/equations/contributions.hpp serif/equations/point.hpp serif/mapping/status.hpp serif/meta/enums.hpp serif/meta/names.hpp src/include/serif/operators/backend/mfem/nodal_computations.hpp serif/operators/backend/mfem/nodal_data.hpp serif/operators/backend/mfem/nodal_evaluation.hpp serif/operators/backend/mfem/nodal_inputs.hpp serif/operators/backend/mfem/preparation.hpp serif/operators/backend/mfem/rejection.hpp serif/operators/backend/mfem/term_coefficients.hpp serif/operators/concepts.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/computation_store.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/incremental_computation.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/keys.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/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 algorithm array cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional iostream magic_enum/magic_enum_all.hpp memory mfem/fem/fespace.hpp mfem/fem/pfespace.hpp mfem/linalg/vector.hpp nameof.hpp numeric optional span stdexcept string string_view tuple type_traits unistd.h unordered_map utility variant version src/compile_time_checks/serif/operators/computation_owners.cpp src/compile_time_checks/serif/operators/serial_backend.cpp src/include/serif/model/backend/mfem/binding.hpp src/include/serif/model/backend/mfem/equilibrium.hpp src/include/serif/model/computations/point_constraints.hpp src/include/serif/model/graph.hpp src/include/serif/operators/backend/mfem/nodal_computations.hpp serif serif::operators serif::operators::backend serif::operators::backend::mfem #pragmaonce #include<expected> #include<utility> #include<tuple> #include"serif/operators/backend/mfem/nodal_evaluation.hpp" #include"serif/operators/backend/mfem/nodal_inputs.hpp" #include"serif/operators/backend/mfem/preparation.hpp" #include"serif/utils/dependencies/dependencies.hpp" namespaceserif::operators::backend::mfem { template<typenameScope,typenameInputKey,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_nodal_evaluation_computation(EquationTequation,NodalDiscretization<EquationT, SchemaT>discretization) { usingnamespaceutils::dependencies; validate_nodal_discretization(discretization); returnmake_computation<typenamenodal_keys::Scoped<Scope>::Linearization,InputList<InputKey>>( [equation=std::move(equation),discretization](constSelectedNodalValues&input) { usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::invalid_argument>(input.selected_rows==discretization.selected_rows, "Preparethisnodalgraphwiththeselectionusedtoconstructit."); returnevaluate_nodal_values(equation,discretization,input.values); }); } template<typenameScope=void,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_nodal_equation_computation(EquationTequation,NodalDiscretization<EquationT, SchemaT>discretization) { returnmake_nodal_evaluation_computation<Scope, nodal_keys::SelectedInput<Scope, typename EquationT::Test>>(std::move(equation),discretization); } template<typenameScope=void,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_nodal_sampling_computation(NodalDiscretization<EquationT, SchemaT>discretization) { usingnamespaceutils::dependencies; usingBinding=discretization::backend::mfem::TermBinding<typenameEquationT::Test,SchemaT>; validate_nodal_discretization(discretization); returnmake_computation<typenamenodal_keys::Scoped<Scope>::Samples,InputList<typenamenodal_keys::Scoped<Scope>::Bindings,typenameEquationT::Test>>( [discretization](constNodalBindings&bindings,constCoefficientView<Binding>&input) { usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::invalid_argument>(bindings.selection==discretization.selected_rows, "Preparenodalcomputationswiththeiroriginalequationselection."); returnselect_nodal_values(discretization,make_term_coefficients(input)); }); } template<typenameScope=void,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_nodal_equation_computations(EquationTequation,NodalDiscretization<EquationT, SchemaT>discretization) { returnstd::tuple{make_nodal_sampling_computation<Scope>(discretization), make_nodal_evaluation_computation<Scope, typename nodal_keys::Scoped<Scope>::Samples>(std::move(equation),discretization)}; } template<typenameScope=void,AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_point_border_computation(constPointBorderDiscretization<EquationT, SchemaT>&discretization) { usingnamespaceutils::dependencies; returnmake_computation<typenamenodal_keys::Scoped<Scope>::BorderLinearization,InputList<typenamenodal_keys::Scoped<Scope>::Linearization,typenameEquationT::Border>>( [discretization](constNodalData<EquationT, SchemaT>&nodal,constauto&border) { usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::CACHE_FE_BINDING_ERROR,std::invalid_argument>(&border.binding==&discretization.border_binding, "Pointbordercoefficientsmustusetheequation'sscalarbinding."); discretization.border_binding.validate_values(border.values); returnPointBorderData<EquationT, SchemaT>{.nodal=nodal,.border_binding=discretization.border_binding,.border_coefficient=border.values[0]}; }); } template<typenameScope=void,AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT> [[nodiscard]]automake_point_border_computations(EquationTequation,constPointBorderDiscretization<EquationT, SchemaT>&discretization) { returnstd::tuple_cat(make_nodal_equation_computations<Scope>(std::move(equation),discretization.nodal),std::tuple{make_point_border_computation<Scope>(discretization)}); } template<typenameScope=void,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT,typenameCoefficientsT> [[nodiscard]]automake_nodal_equation_cache(EquationTequation,constNodalDiscretization<EquationT, SchemaT>&discretization,constCoefficientsT&term_coefficients) { usingnamespaceutils::dependencies; constautocontext=make_nodal_dependency_context<Scope>(discretization,term_coefficients); autonodal_equation=make_nodal_equation_computation<Scope>(std::move(equation),discretization); returnmake_computation_cache<PreparationRejection>(context,std::move(nodal_equation)); } template<typenameScope=void,AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT,typenameCoefficientsT> [[nodiscard]]automake_point_border_cache(EquationTequation,constPointBorderDiscretization<EquationT, SchemaT>&discretization,constCoefficientsT&term_coefficients) { usingnamespaceutils::dependencies; constautocontext=make_point_border_dependency_context<Scope>(discretization,term_coefficients); autonodal_equation=make_nodal_equation_computation<Scope>(std::move(equation),discretization.nodal); constautopoint_border=make_point_border_computation<Scope>(discretization); returnmake_computation_cache<PreparationRejection>(context,std::move(nodal_equation),point_border); } template<typenameScope=void,typenameCacheT,NodalValueEquationEquationT,discretization::domain::schema::IsSchemaSchemaT,typenameCoefficientsT> [[nodiscard]]std::expected<void,PreparationRejection>try_prepare_nodal(CacheT&cache,constNodalDiscretization<EquationT, SchemaT>&discretization, constCoefficientsT&term_coefficients,constPreparationModemode=PreparationMode::reuse_unchanged) { returntry_prepare_validated(cache, [&]() { constautocontext=make_nodal_dependency_context<Scope>(discretization,term_coefficients); returncache.try_prepare(context,{},mode); }); } template<typenameScope=void,typenameCacheT,AffinePointEquationEquationT,discretization::domain::schema::IsSchemaSchemaT,typenameCoefficientsT> [[nodiscard]]std::expected<void,PreparationRejection>try_prepare_point_border(CacheT&cache,constPointBorderDiscretization<EquationT, SchemaT>&discretization, constCoefficientsT&term_coefficients, constPreparationModemode=PreparationMode::reuse_unchanged) { returntry_prepare_validated(cache, [&]() { constautocontext=make_point_border_dependency_context<Scope>(discretization,term_coefficients); returncache.try_prepare(context,{},mode); }); } template<typenameScope=void,typenameCacheT,typenameDiscretizationT,typenameCoefficientsT> voidprepare_nodal(CacheT&cache,constDiscretizationT&discretization,constCoefficientsT&term_coefficients,constPreparationModemode=PreparationMode::reuse_unchanged) { usingnamespaceutils::errors; constautoaccepted=try_prepare_nodal<Scope>(cache,discretization,term_coefficients,mode); if(!accepted) { serif_error<std::domain_error,SERiFErrorCode::INVALID_ARGUMENT>("{}preparationrejected:{}",DiscretizationT::Equation::name,accepted.error().message); } } template<typenameScope=void,typenameCacheT,typenameDiscretizationT,typenameCoefficientsT> voidprepare_point_border(CacheT&cache,constDiscretizationT&discretization,constCoefficientsT&term_coefficients, constPreparationModemode=PreparationMode::reuse_unchanged) { usingnamespaceutils::errors; constautoaccepted=try_prepare_point_border<Scope>(cache,discretization,term_coefficients,mode); if(!accepted) { serif_error<std::domain_error,SERiFErrorCode::INVALID_ARGUMENT>("{}preparationrejected:{}",DiscretizationT::Equation::name,accepted.error().message); } } template<typenameScope=void,typenameCacheT> [[nodiscard]]NumericalWorknodal_work_report(constCacheT&cache) { NumericalWorkwork; if(!cache.ready()) { returnwork; } if(cache.report().templatestatus<typenamenodal_keys::Scoped<Scope>::Linearization>()==utils::dependencies::PreparationStatus::computed) { constauto&nodal=cache.templateget<typenamenodal_keys::Scoped<Scope>::Linearization>(); work.physics_points=nodal.residual.Size()/nodal.derivative.Width(); } ifconstexpr(requires{cache.templateget<typenamenodal_keys::Scoped<Scope>::BorderLinearization>();}) { usingBorderData=std::remove_cvref_t<decltype(cache.templateget<typenamenodal_keys::Scoped<Scope>::BorderLinearization>())>; work.refreshed_scalars=cache.report().templatesource_changed<typenameBorderData::Border>(); } returnwork; } }//namespaceserif::operators::backend::mfem