preconditioned_system.hpp concepts memory utility serif/backend/mfem/model/normalized_system.hpp serif/backend/mfem/solver/preconditioner/solver.hpp src/include/serif/backend/mfem/model/preconditioning/computations.hpp serif/backend/mfem/model/normalized_system.hpp src/include/serif/backend/mfem/model/preconditioned_system.hpp serif/backend/mfem/model/system.hpp serif/backend/mfem/operators/actions.hpp serif/backend/mfem/operators/linearized_block.hpp serif/backend/mfem/solver/normalization/operators.hpp serif/backend/mfem/solver/preconditioner/inverse.hpp serif/backend/mfem/solver/preconditioner/solver.hpp serif/dimensions/dimensions.hpp serif/dimensions/quantities.hpp serif/discretization/blocks/base.hpp serif/meta/enums.hpp serif/solver/normalization/diagonal.hpp serif/solver/normalization/scales.hpp serif/solver/preconditioner/contract.hpp serif/solver/preconditioner/failure.hpp serif/utils/error/codes.hpp serif/utils/error/errors.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/types/strings.hpp algorithm cmath compare concepts cstddef cstdint cstdio cstdlib exception expected format functional iostream magic_enum/magic_enum_all.hpp memory mfem/linalg/operator.hpp mfem/linalg/sparsemat.hpp mfem/linalg/vector.hpp optional 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/preconditioned_system.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::model::backend::mfem::PreconditionedSystem serif serif::model serif::model::backend serif::model::backend::mfem #pragmaonce #include<concepts> #include<memory> #include<utility> #include"serif/backend/mfem/model/normalized_system.hpp" #include"serif/backend/mfem/solver/preconditioner/solver.hpp" namespaceserif::model::backend::mfem{ template<typenameSystem> requires(std::same_as<System,PhysicalSystem>||std::same_as<System,NormalizedSystem>) classPreconditionedSystemfinal{ public: PreconditionedSystem(std::shared_ptr<constSystem>system,solver::preconditioner::backend::mfem::Inverseinverse): m_system(checked_system(std::move(system))), m_physical_inverse(std::make_shared<solver::preconditioner::backend::mfem::Solver>(std::move(inverse),physical_jacobian())), m_inverse(make_solver()){ } [[nodiscard]]conststd::shared_ptr<constSystem>&system()constnoexcept{ returnm_system; } [[nodiscard]]::mfem::Solver&preconditioner()constnoexcept{ return*m_inverse; } [[nodiscard]]constsolver::preconditioner::backend::mfem::Solver&physical_preconditioner()constnoexcept{ return*m_physical_inverse; } [[nodiscard]]solver::preconditioner::ApplicationContractcontract()constnoexcept{ returnm_physical_inverse->contract(); } private: [[nodiscard]]staticstd::shared_ptr<constSystem>checked_system(std::shared_ptr<constSystem>system){ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(system!=nullptr, "Apreconditionedstellarsystemrequiresanevaluatedstar.PreparethematchingresidualandJacobianbeforebindingitsinverse."); returnsystem; } [[nodiscard]]autophysical_jacobian()const{ ifconstexpr(std::same_as<System,NormalizedSystem>){ returnm_system->physical()->jacobian().action(); }else{ returnm_system->jacobian().action(); } } [[nodiscard]]std::shared_ptr<::mfem::Solver>make_solver()const{ ifconstexpr(std::same_as<System,NormalizedSystem>){ returnstd::make_shared<solver::normalization::backend::mfem::ScaledPreconditioner>(*m_physical_inverse,m_system->jacobian()); }else{ returnm_physical_inverse; } } std::shared_ptr<constSystem>m_system; std::shared_ptr<solver::preconditioner::backend::mfem::Solver>m_physical_inverse; std::shared_ptr<::mfem::Solver>m_inverse; }; }