physical_blocks.hpp array limits tuple serif/backend/mfem/model/system.hpp serif/model/keys.hpp serif/operators/linearized_system_view.hpp serif/utils/dependencies/computation.hpp src/include/serif/backend/mfem/model/preconditioning/computations.hpp src/include/serif/backend/mfem/model/computations/physical_blocks.hpp serif/backend/mfem/model/system.hpp serif/backend/mfem/operators/actions.hpp serif/backend/mfem/operators/linearized_block.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/system_layout.hpp serif/model/keys.hpp serif/operators/linearized_system_view.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/type_lists.hpp serif/utils/dependencies/key_list.hpp serif/utils/dependencies/key_list_operations.hpp serif/utils/dependencies/resource.hpp serif/utils/dependencies/snapshot.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/std_helper/cleaning.hpp serif/utils/misc/traits/type_uniqueness.hpp serif/utils/types/type_list.hpp array concepts cstddef functional limits memory mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp optional span tuple type_traits utility vector src/compile_time_checks/serif/model/state.cpp src/include/serif/backend/mfem/model/computations/physical_blocks.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/preconditioning/computations.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::model serif::model::computations serif::model::computations::detail #pragmaonce #include<array> #include<limits> #include<tuple> #include"serif/backend/mfem/model/system.hpp" #include"serif/model/keys.hpp" #include"serif/operators/linearized_system_view.hpp" #include"serif/utils/dependencies/computation.hpp" namespaceserif::model::computations{ namespacedetail{ template<typenameLayout,typename...Rows,typename...Fields> [[nodiscard]]autophysical_block(constbackend::mfem::PhysicalSystem&system,constLayout&layout, utils::dependencies::KeyList<Rows...>,utils::dependencies::KeyList<Fields...>){ usingKind=discretization::blocks::RepresentedValueKind; static_assert(sizeof...(Rows)>0&&sizeof...(Fields)>0,"Selectatleastoneequationandonephysicalcoordinateforasolverblock."); conststd::arrayrows{backend::mfem::CoordinateRange{static_cast<int>(layout.templateoffset<Kind::residual,Rows>()),static_cast<int>(layout.templatesize<Kind::residual,Rows>())}...}; conststd::arraycolumns{backend::mfem::CoordinateRange{static_cast<int>(layout.templateoffset<Kind::value,Fields>()),static_cast<int>(layout.templatesize<Kind::value,Fields>())}...}; returnsystem.block(rows,columns); } } template<typenameBlock,discretization::blocks::SystemLayoutTypeLayout> [[nodiscard]]autophysical_block(constLayout&layout){ usingnamespaceutils::dependencies; usingnamespaceutils::errors; usingKind=discretization::blocks::RepresentedValueKind; serif_ensure<SERiFErrorCode::INVALID_DOF_COUNT>(layout.templatetotal_size<Kind::residual>()<=std::numeric_limits<int>::max()&& layout.templatetotal_size<Kind::value>()<=std::numeric_limits<int>::max(),"PhysicalsolverblockcoordinatesexceedMFEM'sintegerindexrange.Supplyarepresentablesystemlayout."); returnmake_computation<Block,InputList<keys::PhysicalSystem>>([layout](conststd::shared_ptr<constbackend::mfem::PhysicalSystem>&system){ serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(system!=nullptr,"Physicalsolverblocksrequireanevaluatedstellarsystem."); serif_ensure<SERiFErrorCode::MATRIX_SHAPE_ERROR>(layout.templatetotal_size<Kind::residual>()==static_cast<std::size_t>(system->jacobian().action()->Height())&& layout.templatetotal_size<Kind::value>()==static_cast<std::size_t>(system->jacobian().action()->Width()),"Therequestedsolverblocklayoutdoesnotmatchitsfrozenphysicalsystem.Usethatsystem'scoordinatedescription."); returndetail::physical_block(*system,layout,typenameBlock::Rows{},typenameBlock::Values{}); }); } template<discretization::blocks::SystemLayoutTypeLayout,typename...Blocks> [[nodiscard]]automake_physical_block_computations(constLayout&layout,utils::dependencies::KeyList<Blocks...>){ returnstd::tuple{physical_block<Blocks>(layout)...}; } }