coupling.hpp serif/backend/mfem/operators/linearized_block.hpp serif/solver/preconditioner/recipes/coupled.hpp serif/solver/preconditioner/recipes/bordered.hpp serif/backend/mfem/solver/preconditioner/coupling.hpp serif/backend/mfem/solver/preconditioner/border.hpp serif/utils/dependencies/computation.hpp src/include/serif/backend/mfem/solver/preconditioner/graph.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/solver/preconditioner/border.hpp src/include/serif/backend/mfem/solver/preconditioner/computations/coupling.hpp serif/backend/mfem/solver/preconditioner/coupling.hpp serif/backend/mfem/solver/preconditioner/inverse.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/system_layout.hpp serif/meta/enums.hpp serif/operators/linearized_system_view.hpp serif/solver/preconditioner/contract.hpp serif/solver/preconditioner/failure.hpp serif/solver/preconditioner/keys.hpp serif/solver/preconditioner/recipes/bordered.hpp serif/solver/preconditioner/recipes/coupled.hpp serif/solver/preconditioner/strategies/coupling/order.hpp serif/solver/preconditioner/strategies/coupling/policies.hpp serif/solver/preconditioner/strategies/inverse/dense.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/error/codes.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 cstdint expected functional magic_enum/magic_enum_all.hpp memory mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp span string tuple type_traits unordered_map utility vector src/compile_time_checks/serif/model/state.cpp 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/backend/mfem/solver/preconditioner/computations/coupling.hpp src/include/serif/backend/mfem/solver/preconditioner/graph.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/preconditioner/preconditioner.hpp src/include/serif/solver/solver.hpp serif serif::solver serif::solver::preconditioner serif::solver::preconditioner::computations #pragmaonce #include"serif/backend/mfem/operators/linearized_block.hpp" #include"serif/solver/preconditioner/recipes/coupled.hpp" #include"serif/solver/preconditioner/recipes/bordered.hpp" #include"serif/backend/mfem/solver/preconditioner/coupling.hpp" #include"serif/backend/mfem/solver/preconditioner/border.hpp" #include"serif/utils/dependencies/computation.hpp" namespaceserif::solver::preconditioner::computations { template<typenameFirst,typenameSecond,typenamePolicy> [[nodiscard]]autocoupling(constCoupledPreconditioner<First, Second, Policy>&) { usingnamespaceutils::dependencies; usingRecipe=CoupledPreconditioner<First, Second, Policy>; usingInverse=backend::mfem::Inverse; usingBlock=operators::backend::mfem::LinearizedBlock; static_assert(Policy::order==strategies::coupling::CouplingOrder::second_then_first,"Recipescouplesubsystemssecond-then-first."); returnmake_computation<typenameRecipe::Result,InputList<typenameFirst::Result,typenameSecond::Result,keys::Block<typenameFirst::Rows,typenameSecond::Values>>>( [](constInverse&first,constInverse&second,constBlock&first_from_second) { returnbackend::mfem::try_coupled_inverse(first,second,{first_from_second.action(),nullptr},Policy::order); }); } template<typenameFields,typenameConstraints> [[nodiscard]]autoborder(constBorderedPreconditioner<Fields, Constraints>&choice) { usingnamespaceutils::dependencies; usingRecipe=BorderedPreconditioner<Fields, Constraints>; usingBlock=operators::backend::mfem::LinearizedBlock; returnmake_computation<typenameRecipe::Result,InputList<typenameFields::Result,keys::Block<typenameFields::Rows,typenameConstraints::Values>,keys::Block<typename Constraints::Rows, typename Fields::Values>,keys::Block<typename Constraints::Rows, typename Constraints::Values>>>( [options=choice.constraint_inverse](constbackend::mfem::Inverse&fields,constBlock&columns,constBlock&rows,constBlock&diagonal) { returnbackend::mfem::try_bordered_inverse(fields,{columns.action(),rows.action()},diagonal.action(),options); }); } }//namespaceserif::solver::preconditioner::computations