physical_blocks.hpp expected serif/backend/mfem/operators/linearized_block.hpp serif/solver/preconditioner/failure.hpp serif/solver/preconditioner/keys.hpp serif/utils/dependencies/computation.hpp src/include/serif/backend/mfem/solver/preconditioner/computations/inverse.hpp serif/backend/mfem/operators/linearized_block.hpp src/include/serif/backend/mfem/solver/preconditioner/computations/physical_blocks.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/failure.hpp serif/solver/preconditioner/keys.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/inverse.hpp src/include/serif/backend/mfem/solver/preconditioner/computations/physical_blocks.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<expected> #include"serif/backend/mfem/operators/linearized_block.hpp" #include"serif/solver/preconditioner/failure.hpp" #include"serif/solver/preconditioner/keys.hpp" #include"serif/utils/dependencies/computation.hpp" namespaceserif::solver::preconditioner::computations { template<typenameRows,typenameValues> [[nodiscard]]autophysical_diagonal() { usingnamespaceutils::dependencies; returnmake_computation<keys::Diagonal<Rows,Values>,InputList<keys::Block<Rows,Values>>>( [](constoperators::backend::mfem::LinearizedBlock&block)->std::expected<std::vector<double>,Failure> { if(!block.has_diagonal()||block.action()->Height()!=block.action()->Width()) { returnstd::unexpected(Failure{utils::errors::SERiFErrorCode::PRECONDITIONER_INCOMPATIBLE_STRATEGY, "Thisinverseneedsasquarephysicalblockwithanefficientdiagonalextractor.Supplythematchingfieldblock."}); } std::vector<double>diagonal=block.diagonal(); if(diagonal.size()!=static_cast<std::size_t>(block.action()->Height())) { returnstd::unexpected(Failure{utils::errors::SERiFErrorCode::PRECONDITIONER_SIZE_MISMATCH, "Thephysicaldiagonaldoesnotmatchitsacceptedoperator.Correcttheselectedblockextractor."}); } returndiagonal; }); } template<typenameRows,typenameValues> [[nodiscard]]autophysical_matrix() { usingnamespaceutils::dependencies; returnmake_computation<keys::Matrix<Rows,Values>,InputList<keys::Block<Rows,Values>>>( [](constoperators::backend::mfem::LinearizedBlock&block)->std::expected<std::shared_ptr<const::mfem::SparseMatrix>,Failure> { if(!block.has_assembly()) { returnstd::unexpected(Failure{utils::errors::SERiFErrorCode::PRECONDITIONER_INCOMPATIBLE_STRATEGY, "Thisinverseneedsaselectedsparsematrix.Supplythephysicalblock'sassemblerorchooseamatrix-freeinverse."}); } automatrix=block.assemble(); if(!matrix||matrix->Height()!=block.action()->Height()||matrix->Width()!=block.action()->Width()) { returnstd::unexpected(Failure{utils::errors::SERiFErrorCode::PRECONDITIONER_SIZE_MISMATCH, "Selectedsparseassemblydoesnotmatchitsacceptedaction.Correcttheblockassembler."}); } returnmatrix; }); } }//namespaceserif::solver::preconditioner::computations namespaceserif::solver::preconditioner{ template<typenameRow,typenameField> [[nodiscard]]autophysical_block(operators::backend::mfem::LinearizedBlockblock) { usingnamespaceutils::dependencies; returnown<keys::Block<KeyList<Row>,KeyList<Field>>>(std::move(block)); } }