partial.hpp serif/dimensions/quantities.hpp serif/utils/misc/concepts/numeric.hpp src/compile_time_checks/serif/dimensions/quantities.cpp src/include/serif/eos/models/polytropic.hpp src/include/serif/eos/models/traits.hpp src/include/serif/dimensions/partial.hpp serif/dimensions/quantities.hpp serif/utils/misc/concepts/numeric.hpp serif/utils/misc/std_helper/cleaning.hpp compare concepts tuple type_traits src/compile_time_checks/serif/backend_independent.cpp src/compile_time_checks/serif/dimensions/quantities.cpp src/compile_time_checks/serif/eos/consumers_meanfield.cpp src/compile_time_checks/serif/eos/evaluation.cpp src/compile_time_checks/serif/eos/models/traits.cpp src/compile_time_checks/serif/geometry/surface.cpp src/compile_time_checks/serif/geometry/surface_boundary.cpp src/compile_time_checks/serif/model/barotropic_density.cpp src/compile_time_checks/serif/model/density_primal.cpp src/compile_time_checks/serif/model/execution_validation.cpp src/compile_time_checks/serif/model/finite_elements.cpp src/compile_time_checks/serif/model/normalization.cpp src/compile_time_checks/serif/model/state.cpp src/compile_time_checks/serif/model/validation.cpp src/compile_time_checks/serif/operators/coordinate_identity_meanfield.cpp src/compile_time_checks/serif/operators/plan.cpp src/compile_time_checks/serif/solver/continuation.cpp src/compile_time_checks/serif/solver/normalization/plan.cpp src/compile_time_checks/serif/solver/public_contracts.cpp src/include/serif/backend/mfem/geometry/surface/boundary_computations.hpp src/include/serif/backend/mfem/geometry/surface/boundary_data.hpp src/include/serif/backend/mfem/geometry/surface/boundary_inputs.hpp src/include/serif/backend/mfem/geometry/surface/boundary_rows.hpp src/include/serif/backend/mfem/geometry/surface/pressure_boundary.hpp src/include/serif/backend/mfem/geometry/surface/pressure_field_block.hpp src/include/serif/backend/mfem/model/binding.hpp src/include/serif/backend/mfem/model/computations/pressure_boundary.hpp src/include/serif/backend/mfem/model/equilibrium.hpp src/include/serif/backend/mfem/model/equilibrium_evaluator.hpp src/include/serif/backend/mfem/model/failure.hpp src/include/serif/backend/mfem/model/finite_elements.hpp src/include/serif/backend/mfem/model/graph.hpp src/include/serif/backend/mfem/model/normalization/computations.hpp src/include/serif/backend/mfem/serialize/configuration.hpp src/include/serif/backend/mfem/serialize/discretization.hpp src/include/serif/backend/mfem/serialize/options.hpp src/include/serif/backend/mfem/serialize/top.hpp src/include/serif/backend/mfem/solver/admissibility/geometry.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/seed/fitted.hpp src/include/serif/backend/mfem/solver/seed/state.hpp src/include/serif/dimensions/partial.hpp src/include/serif/eos/evaluation.hpp src/include/serif/eos/models/concepts.hpp src/include/serif/eos/models/polytropic.hpp src/include/serif/eos/models/traits.hpp src/include/serif/equations/barotropic_density.hpp src/include/serif/equations/equations.hpp src/include/serif/geometry/surface/automatic.hpp src/include/serif/geometry/surface/compilation.hpp src/include/serif/geometry/surface/computations.hpp src/include/serif/geometry/surface/detail/automatic_relations.hpp src/include/serif/geometry/surface/detail/batches.hpp src/include/serif/geometry/surface/detail/evaluation.hpp src/include/serif/geometry/surface/detail/relations.hpp src/include/serif/geometry/surface/evaluation.hpp src/include/serif/geometry/surface/linearization.hpp src/include/serif/geometry/surface/samples.hpp src/include/serif/geometry/surface/surface.hpp src/include/serif/model/detail/surface_plan.hpp src/include/serif/model/equilibrium_plan.hpp src/include/serif/model/execution_validation.hpp src/include/serif/model/finite_elements.hpp src/include/serif/model/model.hpp src/include/serif/model/normalization/plan.hpp src/include/serif/solver/continuation/continuation.hpp src/include/serif/solver/continuation/equilibrium.hpp src/include/serif/solver/continuation/polytropic.hpp src/include/serif/solver/continuation/stellar_parameters.hpp src/include/serif/solver/equilibrium.hpp src/include/serif/solver/seed/fitted.hpp src/include/serif/solver/seed/polytrope.hpp src/include/serif/solver/seed/state.hpp src/include/serif/solver/solver.hpp src/include/serif/utils/serialize/equations.hpp src/include/serif/utils/serialize/model.hpp src/include/serif/utils/serialize/numerics.hpp src/include/serif/utils/serialize/physics.hpp src/include/serif/utils/serialize/serialize.hpp src/include/serif/utils/serialize/solver.hpp src/lib/serif/backend/mfem/model/finite_elements.cpp src/lib/serif/backend/mfem/serialize/top.cpp src/lib/serif/eos/models/polytropic.cpp src/lib/serif/solver/seed/fitted.cpp src/lib/serif/solver/seed/polytrope.cpp serif::dimensions::WithRespectTo serif::dimensions::SecondPartialDerivative serif::dimensions::PartialDerivative serif serif::dimensions Compile time tools to represent a partial derivative of one quantity wrt. another quantity. Note that this system is not responsible for computing partials or in fact derivatives of any kind. Rather, it is responsible for forcing the developer or compiler at least to ask for a specific, explicit, partial and have to spell out what that is. That in turn gives the compiler enough information to do a lot of downstream things. A motivating example would be automatic verification of the provided EOS against all partials that might be needed to ensure they are all exported by the EOS. Emily M. Boudreaux #pragmaonce #include"serif/dimensions/quantities.hpp" #include"serif/utils/misc/concepts/numeric.hpp" namespaceserif::dimensions{ template<ThermodynamicQuantityTypeQuantity> structWithRespectTofinal{}; template<ThermodynamicQuantityTypeOutputQuantity,ThermodynamicQuantityTypeFirstInput,ThermodynamicQuantityTypeSecondInput> classSecondPartialDerivativefinal{ public: explicitconstexprSecondPartialDerivative(doublevalue)noexcept:m_value(value){} [[nodiscard]]constexprdoublevalue()constnoexcept{returnm_value;} private: doublem_value; }; template<ThermodynamicQuantityTypeOutputQuantity,ThermodynamicQuantityTypeInputQuantity> classPartialDerivativefinal{ public: explicitconstexprPartialDerivative(constdoublevalue)noexcept:m_value(value){} [[nodiscard]]constexprdoublevalue()constnoexcept{ returnm_value; } friendconstexprPartialDerivative<OutputQuantity, InputQuantity>operator+( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, constPartialDerivative<OutputQuantity, InputQuantity>&rhs )noexcept{ returnPartialDerivative<OutputQuantity, InputQuantity>{lhs.m_value+rhs.m_value}; } friendconstexprPartialDerivative<OutputQuantity, InputQuantity>operator-( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, constPartialDerivative<OutputQuantity, InputQuantity>&rhs )noexcept{ returnPartialDerivative<OutputQuantity, InputQuantity>{lhs.m_value-rhs.m_value}; } template<utils::misc::concepts::IsNumericScalar> friendconstexprPartialDerivative<OutputQuantity, InputQuantity>operator*( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, constScalar&rhs )noexcept{ returnPartialDerivative<OutputQuantity, InputQuantity>{lhs.m_value*rhs}; } template<utils::misc::concepts::IsNumericScalar> friendconstexprPartialDerivative<OutputQuantity, InputQuantity>operator*( constScalar&lhs, constPartialDerivative<OutputQuantity, InputQuantity>&rhs )noexcept{ returnPartialDerivative<OutputQuantity, InputQuantity>{lhs*rhs.m_value}; } template<utils::misc::concepts::IsNumericScalar> friendconstexprPartialDerivative<OutputQuantity, InputQuantity>operator/( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, constScalar&rhs )noexcept{ returnPartialDerivative<OutputQuantity, InputQuantity>{lhs.m_value/rhs}; } template<utils::misc::concepts::IsNumericScalar> friendconstexprstd::partial_orderingoperator<=>( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, Scalarrhs )noexcept{ returnlhs.m_value<=>rhs; } template<utils::misc::concepts::IsNumericScalar> friendconstexprstd::partial_orderingoperator<=>( Scalarlhs, constPartialDerivative<OutputQuantity, InputQuantity>&rhs )noexcept{ returnlhs<=>rhs.m_value; } friendconstexprstd::partial_orderingoperator<=>( constPartialDerivative<OutputQuantity, InputQuantity>&lhs, constPartialDerivative<OutputQuantity, InputQuantity>&rhs )noexcept{ returnlhs.m_value<=>rhs.m_value; } private: doublem_value; }; }