inverse.hpp concepts serif/utils/error/errors.hpp serif/solver/preconditioner/recipes/block.hpp serif/backend/mfem/solver/preconditioner/computations/physical_blocks.hpp serif/backend/mfem/solver/preconditioner/inverse/diagonal.hpp serif/backend/mfem/solver/preconditioner/inverse/dense.hpp serif/backend/mfem/solver/preconditioner/inverse/amg.hpp serif/backend/mfem/solver/preconditioner/inverse/sparse_direct.hpp serif/solver/preconditioner/strategies/inverse/diagonal.hpp src/include/serif/backend/mfem/solver/preconditioner/graph.hpp serif/backend/mfem/operators/linearized_block.hpp src/include/serif/backend/mfem/solver/preconditioner/computations/inverse.hpp serif/backend/mfem/solver/preconditioner/computations/physical_blocks.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/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/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/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/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 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/computations/inverse.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<concepts> #include"serif/utils/error/errors.hpp" #include"serif/solver/preconditioner/recipes/block.hpp" #include"serif/backend/mfem/solver/preconditioner/computations/physical_blocks.hpp" #include"serif/backend/mfem/solver/preconditioner/inverse/diagonal.hpp" #include"serif/backend/mfem/solver/preconditioner/inverse/dense.hpp" #include"serif/backend/mfem/solver/preconditioner/inverse/amg.hpp" #include"serif/backend/mfem/solver/preconditioner/inverse/sparse_direct.hpp" #include"serif/solver/preconditioner/strategies/inverse/diagonal.hpp" namespaceserif::solver::preconditioner::computations { template<typenameStrategy> inlineconstexprboolneeds_diagonal=std::same_as<Strategy,strategies::inverse::Diagonal>; template<typenameStrategy> inlineconstexprboolneeds_matrix=std::same_as<Strategy,strategies::inverse::DenseDirect>||std::same_as<Strategy,strategies::inverse::FixedAMG>|| std::same_as<Strategy,strategies::inverse::SparseDirect>; template<typenameStrategy> [[nodiscard]]std::expected<backend::mfem::Inverse,Failure>invert_matrix(std::shared_ptr<const::mfem::SparseMatrix>matrix,constStrategy&strategy) { usingnamespacebackend::mfem; ifconstexpr(std::same_as<Strategy,strategies::inverse::DenseDirect>) { ::mfem::DenseMatrixdense; matrix->ToDenseMatrix(dense); returntry_dense_inverse(dense,strategy); } elseifconstexpr(std::same_as<Strategy,strategies::inverse::SparseDirect>) { returntry_sparse_direct_inverse(std::move(matrix),strategy); } else { returntry_amg_inverse(std::move(matrix),strategy); } } template<typenameRow,typenameField,typenameStrategy> [[nodiscard]]autoinverse(constBlockPreconditioner<Row,Field,Strategy>&choice) { usingnamespaceutils::dependencies; usingRecipe=BlockPreconditioner<Row,Field,Strategy>; usingRows=typenameRecipe::Rows; usingValues=typenameRecipe::Values; usingBlock=operators::backend::mfem::LinearizedBlock; usingResult=backend::mfem::Inverse; ifconstexpr(std::same_as<Strategy,strategies::inverse::Diagonal>) { returnmake_computation<typenameRecipe::Result,InputList<keys::Diagonal<Rows,Values>>>( [](conststd::vector<double>&diagonal) { returnbackend::mfem::try_diagonal_inverse(diagonal); }); } elseifconstexpr(std::same_as<Strategy,strategies::inverse::SparseDirect>) { //Thecomputationlivesaslongasthepreconditioner,soNewtonstepsshareonesymbolicanalysis. returnmake_computation<typenameRecipe::Result,InputList<keys::Matrix<Rows,Values>>>( [strategy=choice.inverse,cache=std::make_shared<backend::mfem::SparseDirectCache>()](conststd::shared_ptr<const::mfem::SparseMatrix>&matrix) { returnbackend::mfem::try_sparse_direct_inverse(matrix,strategy,cache.get()); }); } elseifconstexpr(needs_matrix<Strategy>) { returnmake_computation<typenameRecipe::Result,InputList<keys::Matrix<Rows,Values>>>( [strategy=choice.inverse](conststd::shared_ptr<const::mfem::SparseMatrix>&matrix) { returninvert_matrix(matrix,strategy); }); } else { returnmake_computation<typenameRecipe::Result,InputList<keys::Block<Rows,Values>>>( [strategy=choice.inverse](constBlock&block)->std::expected<Result,Failure> { usingnamespaceutils::errors; autoresult=strategy.try_inverse(block); if(!result) { returnresult; } if(!result->action) { returnstd::unexpected(Failure{SERiFErrorCode::PRECONDITIONER_INVALID_OPERATOR, format_error_message<SERiFErrorCode::PRECONDITIONER_INVALID_OPERATOR>("Acustominversereturnednoaction.Returnanownedinversefortheselectedphysicalblockorapreparationfailure.")}); } if(result->action->Height()!=block.action()->Width()||result->action->Width()!=block.action()->Height()) { returnstd::unexpected(Failure{SERiFErrorCode::PRECONDITIONER_SIZE_MISMATCH, format_error_message<SERiFErrorCode::PRECONDITIONER_SIZE_MISMATCH>("Acustominversemustreversethedimensionsofitsselectedphysicalblock.Returntheinverseinthatblock'sdeclaredequationandunknowncoordinates.")}); } if(result->contract!=ApplicationContract::fixed_linear&&result->contract!=ApplicationContract::variable) { returnstd::unexpected(Failure{SERiFErrorCode::PRECONDITIONER_INVALID_OPTIONS, format_error_message<SERiFErrorCode::PRECONDITIONER_INVALID_OPTIONS>("Acustominversemustdeclarefixed-linearorvariableapplication.Selectthecontractrequiredbyitsnumericalalgorithm.")}); } returnresult; }); } } }//namespaceserif::solver::preconditioner::computations