normalized_system.hpp algorithm cmath expected memory string vector serif/backend/mfem/model/system.hpp serif/backend/mfem/solver/normalization/operators.hpp serif/solver/normalization/scales.hpp serif/utils/error/errors.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/model/preconditioned_system.hpp src/include/serif/backend/mfem/model/normalized_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/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/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/normalization/computations.hpp src/include/serif/backend/mfem/model/normalized_system.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::NormalizedSystem serif serif::model serif::model::backend serif::model::backend::mfem #pragmaonce #include<algorithm> #include<cmath> #include<expected> #include<memory> #include<string> #include<vector> #include"serif/backend/mfem/model/system.hpp" #include"serif/backend/mfem/solver/normalization/operators.hpp" #include"serif/solver/normalization/scales.hpp" #include"serif/utils/error/errors.hpp" namespaceserif::model::backend::mfem{ classNormalizedSystemfinal{ public: NormalizedSystem(std::shared_ptr<constPhysicalSystem>physical, solver::normalization::DiagonalNormalizationcoordinates,solver::normalization::StellarReferenceScalesscales): m_physical(checked_physical(std::move(physical))),m_scales(std::move(scales)), m_jacobian(*m_physical->jacobian().action(),std::move(coordinates)),m_residual(m_physical->residual().size()){ m_jacobian.normalization().normalize_residual(m_physical->residual(),m_residual); } [[nodiscard]]conststd::shared_ptr<constPhysicalSystem>&physical()constnoexcept{ returnm_physical; } [[nodiscard]]constsolver::normalization::DiagonalNormalization&coordinates()constnoexcept{ returnm_jacobian.normalization(); } [[nodiscard]]constsolver::normalization::StellarReferenceScales&reference_scales()constnoexcept{ returnm_scales; } [[nodiscard]]conststd::vector<double>&residual()constnoexcept{ returnm_residual; } [[nodiscard]]constsolver::normalization::backend::mfem::ScaledJacobianOperator&jacobian()constnoexcept{ returnm_jacobian; } private: [[nodiscard]]staticstd::shared_ptr<constPhysicalSystem>checked_physical(std::shared_ptr<constPhysicalSystem>physical){ usingnamespaceutils::errors; serif_ensure<SERiFErrorCode::INVALID_ARGUMENT>(physical!=nullptr,"Normalizedcoordinatesrequireanfrozenphysicalsystem.Preparethestellarmodelfirst."); returnphysical; } std::shared_ptr<constPhysicalSystem>m_physical; solver::normalization::StellarReferenceScalesm_scales; solver::normalization::backend::mfem::ScaledJacobianOperatorm_jacobian; std::vector<double>m_residual; }; [[nodiscard]]inlinestd::expected<std::shared_ptr<constNormalizedSystem>,std::string>make_normalized_system( std::shared_ptr<constPhysicalSystem>physical,solver::normalization::DiagonalNormalizationcoordinates, solver::normalization::StellarReferenceScalesreference){ autoresult=std::make_shared<constNormalizedSystem>(std::move(physical),std::move(coordinates),std::move(reference)); if(!std::ranges::all_of(result->residual(),[](doublevalue){returnstd::isfinite(value);})){ returnstd::unexpected("Thenormalizedstellarresidualisnotfinite.Checkthereferencescalesandthephysicalcandidate."); } returnresult; } }