graph.hpp tuple utility serif/backend/mfem/solver/preconditioner/computations/inverse.hpp serif/backend/mfem/solver/preconditioner/computations/coupling.hpp serif/backend/mfem/solver/preconditioner/computations/ordering.hpp serif/utils/dependencies/computation_graph.hpp src/include/serif/backend/mfem/model/preconditioning/computations.hpp src/include/serif/solver/preconditioner/preconditioner.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/solver/preconditioner/border.hpp serif/backend/mfem/solver/preconditioner/computations/coupling.hpp serif/backend/mfem/solver/preconditioner/computations/inverse.hpp serif/backend/mfem/solver/preconditioner/computations/ordering.hpp serif/backend/mfem/solver/preconditioner/computations/physical_blocks.hpp serif/backend/mfem/solver/preconditioner/coupling.hpp src/include/serif/backend/mfem/solver/preconditioner/graph.hpp serif/backend/mfem/solver/preconditioner/inverse.hpp serif/backend/mfem/solver/preconditioner/inverse/amg.hpp serif/backend/mfem/solver/preconditioner/inverse/dense.hpp serif/backend/mfem/solver/preconditioner/inverse/diagonal.hpp serif/backend/mfem/solver/preconditioner/inverse/sparse_direct.hpp serif/discretization/blocks/base.hpp serif/discretization/blocks/keys.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/block.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/solver/preconditioner/strategies/inverse/diagonal.hpp serif/solver/preconditioner/strategies/inverse/fixed_amg.hpp serif/solver/preconditioner/strategies/inverse/sparse_direct.hpp serif/utils/dependencies/binding.hpp serif/utils/dependencies/computation.hpp serif/utils/dependencies/computation_graph.hpp serif/utils/dependencies/dependency_context.hpp serif/utils/dependencies/detail/graph.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/error/errors.hpp serif/utils/misc/concepts/concat.hpp serif/utils/misc/concepts/exceptions.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/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 concepts cstddef cstdint cstdio cstdlib exception expected format functional iostream limits magic_enum/magic_enum_all.hpp memory mfem/linalg/densemat.hpp mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp span string string_view tuple type_traits unistd.h unordered_map utility vector version 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/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::detail #pragmaonce #include<tuple> #include<utility> #include"serif/backend/mfem/solver/preconditioner/computations/inverse.hpp" #include"serif/backend/mfem/solver/preconditioner/computations/coupling.hpp" #include"serif/backend/mfem/solver/preconditioner/computations/ordering.hpp" #include"serif/utils/dependencies/computation_graph.hpp" namespaceserif::solver::preconditioner { namespacedetail { template<typenameRow,typenameField,typenameStrategy> [[nodiscard]]autopreconditioner_computations(constBlockPreconditioner<Row, Field, Strategy>&choice); template<typenameFirst,typenameSecond,typenamePolicy> [[nodiscard]]autopreconditioner_computations(constCoupledPreconditioner<First, Second, Policy>&choice); template<typenameFields,typenameConstraints> [[nodiscard]]autopreconditioner_computations(constBorderedPreconditioner<Fields, Constraints>&choice); template<typenameRow,typenameField,typenameStrategy> [[nodiscard]]autopreconditioner_computations(constBlockPreconditioner<Row, Field, Strategy>&choice) { usingRecipe=BlockPreconditioner<Row, Field, Strategy>; ifconstexpr(computations::needs_diagonal<Strategy>) { returnstd::tuple{computations::physical_diagonal<typename Recipe::Rows, typename Recipe::Values>(),computations::inverse(choice)}; } elseifconstexpr(computations::needs_matrix<Strategy>) { returnstd::tuple{computations::physical_matrix<typename Recipe::Rows, typename Recipe::Values>(),computations::inverse(choice)}; } else { returnstd::tuple{computations::inverse(choice)}; } } template<typenameFirst,typenameSecond,typenamePolicy> [[nodiscard]]autopreconditioner_computations(constCoupledPreconditioner<First, Second, Policy>&choice) { returnstd::tuple_cat(preconditioner_computations(choice.first),preconditioner_computations(choice.second),std::tuple{computations::coupling(choice)}); } template<typenameFields,typenameConstraints> [[nodiscard]]autopreconditioner_computations(constBorderedPreconditioner<Fields, Constraints>&choice) { returnstd::tuple_cat(preconditioner_computations(choice.fields),std::tuple{computations::border(choice)}); } }//namespacedetail template<typenameRecipe> [[nodiscard]]automake_preconditioner_computations(constRecipe&choice) { returndetail::preconditioner_computations(choice); } template<typenameRecipe> [[nodiscard]]automake_preconditioner_graph(constRecipe&choice) { returnstd::apply([](auto...nodes) { returnutils::dependencies::make_computation_graph(std::move(nodes)...); },make_preconditioner_computations(choice)); } template<typenameRecipe,discretization::blocks::SystemLayoutTypeLayout,utils::dependencies::IsComputation...Additional> [[nodiscard]]automake_preconditioner_computations(constRecipe&choice,constLayout&,Additional...additional){ returnstd::apply([](auto...nodes){ usingGraph=decltype(utils::dependencies::make_computation_graph(std::move(nodes)...)); usingPhysicalKeys=typenamecomputations::detail::PhysicalSources<typenameGraph::Inputs>::Type; returnstd::tuple{std::move(nodes)...,computations::ordering<Recipe, Layout>(PhysicalKeys{})}; },std::tuple_cat(std::tuple{std::move(additional)...},make_preconditioner_computations(choice))); } template<typenameRecipe,discretization::blocks::SystemLayoutTypeLayout,utils::dependencies::IsComputation...Additional> [[nodiscard]]automake_preconditioner_graph(constRecipe&choice,constLayout&layout,Additional...additional){ returnstd::apply([](auto...nodes){ returnutils::dependencies::make_computation_graph(std::move(nodes)...); },make_preconditioner_computations(choice,layout,std::move(additional)...)); } }//namespaceserif::solver::preconditioner