From 99cfb4b325558e0b4c1c1ea53daf19c1328dbf58 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:31:56 -0600 Subject: [PATCH 01/17] Reference shift modifier for Davis Products EOS --- doc/sphinx/src/models.rst | 9 +++++++-- doc/sphinx/src/modifiers.rst | 2 ++ 2 files changed, 9 insertions(+), 2 deletions(-) diff --git a/doc/sphinx/src/models.rst b/doc/sphinx/src/models.rst index 19424ce45c..163b2d9721 100644 --- a/doc/sphinx/src/models.rst +++ b/doc/sphinx/src/models.rst @@ -1231,7 +1231,7 @@ Finally, the pressure, energy, and temperature along the isentrope are given by .. math:: - e_S(\rho) = e_{\mathrm{C}} G(\rho) \frac{1}{\rho V_{\mathrm{C}}} - e_0 + e_S(\rho) = e_{\mathrm{C}} G(\rho) \frac{1}{\rho V_{\mathrm{C}}} .. math:: @@ -1272,7 +1272,7 @@ The energy at the CJ state can be calculated as e_\mathrm{CJ} = \frac{P_0 + P_\mathrm{CJ}}{2(V_0 - V_\mathrm{CJ})}, relative to :math:`e = 0` at the reference state of the *reactants*. Therefore -the :math:`e_0` energy offset of the products EOS is given by +the energy offset of the products EOS is given by .. math:: @@ -1280,6 +1280,11 @@ the :math:`e_0` energy offset of the products EOS is given by Practically, this means :math:`e_0` should be positive for any energetic material. +To provide the energy offset to the Davis Products EOS, `the energy shift +modifier`_ should be used. Note that the convention there +is that the shift is positive, so :math:`-e_0` should be provided to the shift +modifier. + The constructor for the Davis Products EOS is .. code-block:: cpp diff --git a/doc/sphinx/src/modifiers.rst b/doc/sphinx/src/modifiers.rst index 552b16940b..75e11310af 100644 --- a/doc/sphinx/src/modifiers.rst +++ b/doc/sphinx/src/modifiers.rst @@ -22,6 +22,8 @@ We list below the available modifiers and their constructors. The Shifted EOS ----------------- +.. _modifiers shifted EOS: + The shifted equation of state modifies zero point energy of an underlying model by some shift. So for example, it transforms From 83700480eb8769959f92bfd5d5883785ab1f3f96 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:32:58 -0600 Subject: [PATCH 02/17] Remove E0 from input --- singularity-eos/eos/eos_davis.hpp | 11 +++++------ 1 file changed, 5 insertions(+), 6 deletions(-) diff --git a/singularity-eos/eos/eos_davis.hpp b/singularity-eos/eos/eos_davis.hpp index 86eedd224b..993d8d5b5b 100644 --- a/singularity-eos/eos/eos_davis.hpp +++ b/singularity-eos/eos/eos_davis.hpp @@ -177,8 +177,8 @@ class DavisProducts : public EosBase { DavisProducts() = default; PORTABLE_INLINE_FUNCTION DavisProducts(const Real a, const Real b, const Real k, const Real n, const Real vc, - const Real pc, const Real Cv, const Real E0) - : _a(a), _b(b), _k(k), _n(n), _vc(vc), _pc(pc), _Cv(Cv), _E0(E0) {} + const Real pc, const Real Cv) + : _a(a), _b(b), _k(k), _n(n), _vc(vc), _pc(pc), _Cv(Cv), {} DavisProducts GetOnDevice() { return *this; } template PORTABLE_INLINE_FUNCTION Real TemperatureFromDensityInternalEnergy( @@ -287,7 +287,7 @@ class DavisProducts : public EosBase { PORTABLE_INLINE_FUNCTION void PrintParams() const { static constexpr char s1[]{"DavisProducts Params: "}; printf("%sa:%e b:%e k:%e\nn:%e vc:%e pc:%e\nCv:%e E0:%e\n", s1, _a, _b, _k, _n, _vc, - _pc, _Cv, _E0); + _pc, _Cv); } inline void Finalize() {} static std::string EosType() { return std::string("DavisProducts"); } @@ -296,7 +296,7 @@ class DavisProducts : public EosBase { private: static constexpr Real onethird = 1.0 / 3.0; - Real _a, _b, _k, _n, _vc, _pc, _Cv, _E0; + Real _a, _b, _k, _n, _vc, _pc, _Cv; // static constexpr const char _eos_type[] = "DavisProducts"; static constexpr const unsigned long _preferred_input = thermalqs::density | thermalqs::specific_internal_energy; @@ -314,8 +314,7 @@ class DavisProducts : public EosBase { const Real ec = _pc * _vc / (_k - 1.0 + _a); // const Real de = ecj-(Es(rho0)-_E0); return ec * std::pow(0.5 * (std::pow(vvc, _n) + std::pow(vvc, -_n)), _a / _n) / - std::pow(vvc, _k - 1.0 + _a) - - _E0; + std::pow(vvc, _k - 1.0 + _a); } PORTABLE_INLINE_FUNCTION Real Ts(const Real rho) const { const Real vvc = 1 / (rho * _vc); From 08da5dde4d53c742269f1c9f70162f30c778964b Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:33:53 -0600 Subject: [PATCH 03/17] Clang format --- singularity-eos/eos/eos_davis.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/singularity-eos/eos/eos_davis.hpp b/singularity-eos/eos/eos_davis.hpp index 993d8d5b5b..6d3d413211 100644 --- a/singularity-eos/eos/eos_davis.hpp +++ b/singularity-eos/eos/eos_davis.hpp @@ -314,7 +314,7 @@ class DavisProducts : public EosBase { const Real ec = _pc * _vc / (_k - 1.0 + _a); // const Real de = ecj-(Es(rho0)-_E0); return ec * std::pow(0.5 * (std::pow(vvc, _n) + std::pow(vvc, -_n)), _a / _n) / - std::pow(vvc, _k - 1.0 + _a); + std::pow(vvc, _k - 1.0 + _a); } PORTABLE_INLINE_FUNCTION Real Ts(const Real rho) const { const Real vvc = 1 / (rho * _vc); From 257e33fcd85f9486d46fe15868f1fdbe710615fe Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:35:54 -0600 Subject: [PATCH 04/17] Typo --- singularity-eos/eos/eos_davis.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/singularity-eos/eos/eos_davis.hpp b/singularity-eos/eos/eos_davis.hpp index 6d3d413211..5691e10f39 100644 --- a/singularity-eos/eos/eos_davis.hpp +++ b/singularity-eos/eos/eos_davis.hpp @@ -178,7 +178,7 @@ class DavisProducts : public EosBase { PORTABLE_INLINE_FUNCTION DavisProducts(const Real a, const Real b, const Real k, const Real n, const Real vc, const Real pc, const Real Cv) - : _a(a), _b(b), _k(k), _n(n), _vc(vc), _pc(pc), _Cv(Cv), {} + : _a(a), _b(b), _k(k), _n(n), _vc(vc), _pc(pc), _Cv(Cv) {} DavisProducts GetOnDevice() { return *this; } template PORTABLE_INLINE_FUNCTION Real TemperatureFromDensityInternalEnergy( From 0032e6793d4c602f9032daa514d80e1926819be1 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:40:48 -0600 Subject: [PATCH 05/17] Remove E0 from Davis Products EOS --- singularity-eos/eos/singularity_eos.cpp | 10 +++++----- singularity-eos/eos/singularity_eos.f90 | 4 ++-- 2 files changed, 7 insertions(+), 7 deletions(-) diff --git a/singularity-eos/eos/singularity_eos.cpp b/singularity-eos/eos/singularity_eos.cpp index e000bfe268..c7c3769bb7 100644 --- a/singularity-eos/eos/singularity_eos.cpp +++ b/singularity-eos/eos/singularity_eos.cpp @@ -154,23 +154,23 @@ int init_sg_JWL(const int matindex, EOS *eos, const double A, const double B, int init_sg_DavisProducts(const int matindex, EOS *eos, const double a, const double b, const double k, const double n, const double vc, - const double pc, const double Cv, const double E0, + const double pc, const double Cv, int const *const enabled, double *const vals) { assert(matindex >= 0); - EOS eosi = SGAPPLYMODSIMPLE(DavisProducts(a, b, k, n, vc, pc, Cv, E0)); + EOS eosi = SGAPPLYMODSIMPLE(DavisProducts(a, b, k, n, vc, pc, Cv)); if (enabled[3] == 1) { singularity::pAlpha2BilinearRampParams(eosi, vals[2], vals[3], vals[4], vals[2], vals[3], vals[4], vals[5]); } - EOS eos_ = SGAPPLYMOD(DavisProducts(a, b, k, n, vc, pc, Cv, E0)); + EOS eos_ = SGAPPLYMOD(DavisProducts(a, b, k, n, vc, pc, Cv)); eos[matindex] = eos_.GetOnDevice(); return 0; } int init_sg_DavisProducts(const int matindex, EOS *eos, const double a, const double b, const double k, const double n, const double vc, - const double pc, const double Cv, const double E0) { - return init_sg_DavisProducts(matindex, eos, a, b, k, n, vc, pc, Cv, E0, def_en, def_v); + const double pc, const double Cv) { + return init_sg_DavisProducts(matindex, eos, a, b, k, n, vc, pc, Cv, def_en, def_v); } int init_sg_DavisReactants(const int matindex, EOS *eos, const double rho0, diff --git a/singularity-eos/eos/singularity_eos.f90 b/singularity-eos/eos/singularity_eos.f90 index 2c0ef227ec..ae131206b1 100644 --- a/singularity-eos/eos/singularity_eos.f90 +++ b/singularity-eos/eos/singularity_eos.f90 @@ -112,13 +112,13 @@ end function init_sg_JWL interface integer(kind=c_int) function & - init_sg_DavisProducts(matindex, eos, a, b, k, n, vc, pc, Cv, E0, & + init_sg_DavisProducts(matindex, eos, a, b, k, n, vc, pc, Cv, & sg_mods_enabled, sg_mods_values) & bind(C, name='init_sg_DavisProducts') import integer(c_int), value, intent(in) :: matindex type(c_ptr), value, intent(in) :: eos - real(kind=c_double), value, intent(in) :: a, b, k, n, vc, pc, Cv, E0 + real(kind=c_double), value, intent(in) :: a, b, k, n, vc, pc, Cv type(c_ptr), value, intent(in) :: sg_mods_enabled, sg_mods_values end function init_sg_DavisProducts end interface From 01627d90fc5b2d9b8bc4d1e3c83863b727c18362 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:41:50 -0600 Subject: [PATCH 06/17] Remove E0 from Davis Products signature --- doc/sphinx/src/models.rst | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/doc/sphinx/src/models.rst b/doc/sphinx/src/models.rst index 163b2d9721..fe37f8e332 100644 --- a/doc/sphinx/src/models.rst +++ b/doc/sphinx/src/models.rst @@ -1290,12 +1290,11 @@ The constructor for the Davis Products EOS is .. code-block:: cpp DavisProducts(const Real a, const Real b, const Real k, const Real n, const Real vc, - const Real pc, const Real Cv, const Real E0) + const Real pc, const Real Cv) where ``a`` is :math:`a`, ``b`` is :math:`b`, ``k`` is :math:`k`, ``n`` is :math:`n`, ``vc`` is :math:`V_\mathrm{C}`, ``pc`` is -:math:`P_\mathrm{C}`, ``Cv`` is :math:`C_{V,0}`, and ``E0`` is -:math:`e_\mathrm{C}`. +:math:`P_\mathrm{C}`, ``Cv`` is :math:`C_{V,0}`. Spiner EOS ```````````` From 234512fa199e2fff9f900055e2a279a3aca490f1 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:43:14 -0600 Subject: [PATCH 07/17] Clang format --- singularity-eos/eos/singularity_eos.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/singularity-eos/eos/singularity_eos.cpp b/singularity-eos/eos/singularity_eos.cpp index c7c3769bb7..c80ab710b5 100644 --- a/singularity-eos/eos/singularity_eos.cpp +++ b/singularity-eos/eos/singularity_eos.cpp @@ -154,8 +154,8 @@ int init_sg_JWL(const int matindex, EOS *eos, const double A, const double B, int init_sg_DavisProducts(const int matindex, EOS *eos, const double a, const double b, const double k, const double n, const double vc, - const double pc, const double Cv, - int const *const enabled, double *const vals) { + const double pc, const double Cv, int const *const enabled, + double *const vals) { assert(matindex >= 0); EOS eosi = SGAPPLYMODSIMPLE(DavisProducts(a, b, k, n, vc, pc, Cv)); if (enabled[3] == 1) { From 571fe88b88e7c911071c92a8c8d509fd9e82f29f Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 21:45:30 -0600 Subject: [PATCH 08/17] Remove final Davis Products EO --- singularity-eos/eos/singularity_eos.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/singularity-eos/eos/singularity_eos.f90 b/singularity-eos/eos/singularity_eos.f90 index ae131206b1..e325601758 100644 --- a/singularity-eos/eos/singularity_eos.f90 +++ b/singularity-eos/eos/singularity_eos.f90 @@ -460,12 +460,12 @@ integer function init_sg_JWL_f(matindex, eos, A, B, R1, R2, w, rho0, Cv, & end function init_sg_JWL_f integer function init_sg_DavisProducts_f(matindex, eos, a, b, k, n, vc, pc, & - Cv, E0, sg_mods_enabled, & + Cv, sg_mods_enabled, & sg_mods_values) & result(err) integer(c_int), value, intent(in) :: matindex type(sg_eos_ary_t), intent(in) :: eos - real(kind=8), value, intent(in) :: a, b, k, n, vc, pc, Cv, E0 + real(kind=8), value, intent(in) :: a, b, k, n, vc, pc, Cv integer(kind=c_int), dimension(:), target, optional, intent(inout) :: sg_mods_enabled real(kind=8), dimension(:), target, optional, intent(inout) :: sg_mods_values ! local vars @@ -478,7 +478,7 @@ integer function init_sg_DavisProducts_f(matindex, eos, a, b, k, n, vc, pc, & if(present(sg_mods_values)) sg_mods_values_use = sg_mods_values err = init_sg_DavisProducts(matindex-1, eos%ptr, a, b, k, n, vc, pc, Cv, & - E0, c_loc(sg_mods_enabled_use), & + c_loc(sg_mods_enabled_use), & c_loc(sg_mods_values_use)) end function init_sg_DavisProducts_f From c519d29e0f7d05ae9b739dd1460ca437074c770b Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 22:19:36 -0600 Subject: [PATCH 09/17] Remove E0 from print parameters function --- singularity-eos/eos/eos_davis.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/singularity-eos/eos/eos_davis.hpp b/singularity-eos/eos/eos_davis.hpp index 5691e10f39..a00a9bfd4e 100644 --- a/singularity-eos/eos/eos_davis.hpp +++ b/singularity-eos/eos/eos_davis.hpp @@ -286,8 +286,8 @@ class DavisProducts : public EosBase { static inline unsigned long max_scratch_size(unsigned int nelements) { return 0; } PORTABLE_INLINE_FUNCTION void PrintParams() const { static constexpr char s1[]{"DavisProducts Params: "}; - printf("%sa:%e b:%e k:%e\nn:%e vc:%e pc:%e\nCv:%e E0:%e\n", s1, _a, _b, _k, _n, _vc, - _pc, _Cv); + printf("%sa:%e b:%e k:%e\nn:%e vc:%e pc:%e\nCv:%e \n", s1, _a, _b, _k, _n, _vc, _pc, + _Cv); } inline void Finalize() {} static std::string EosType() { return std::string("DavisProducts"); } From 06289f4145182a07ec33b857086a0ad7c075e818 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Fri, 10 May 2024 22:32:52 -0600 Subject: [PATCH 10/17] Remove E0 From Davis Products in python module --- python/module.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/python/module.cpp b/python/module.cpp index a6b3ad1c79..711e0990b3 100644 --- a/python/module.cpp +++ b/python/module.cpp @@ -62,7 +62,7 @@ PYBIND11_MODULE(singularity_eos, m) { .def( py::init(), py::arg("a"), py::arg("b"), py::arg("k"), py::arg("n"), py::arg("vc"), - py::arg("pc"), py::arg("Cv"), py::arg("E0") + py::arg("pc"), py::arg("Cv") ); #ifdef SPINER_USE_HDF From cbe4c76f60a764cc793acbcb1431cfc10681000a Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 11:24:45 -0600 Subject: [PATCH 11/17] Remove E0 argument from Davis Products --- test/pte_test_utils.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/pte_test_utils.hpp b/test/pte_test_utils.hpp index d61b09747e..af2e68ad94 100644 --- a/test/pte_test_utils.hpp +++ b/test/pte_test_utils.hpp @@ -59,7 +59,7 @@ inline void set_eos(T *eos) { singularity::EOS dr = singularity::DavisReactants( 1.890, 4.115e10, 1.0e6, 297.0, 1.8e5, 4.6, 0.34, 0.56, 0.0, 0.4265, 0.001074e10); singularity::EOS dp = singularity::DavisProducts(0.798311, 0.58, 1.35, 2.66182, 0.75419, - 3.2e10, 0.001072e10, 0.0); + 3.2e10, 0.001072e10); eos[0] = gr.GetOnDevice(); eos[1] = dr.GetOnDevice(); eos[2] = dp.GetOnDevice(); From bbf23e1e843fd79ce46f75befe710bb79a60d369 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 11:26:31 -0600 Subject: [PATCH 12/17] Update changelong to reflect that Davis Products E0 is now gone instead of 'fixed' --- CHANGELOG.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/CHANGELOG.md b/CHANGELOG.md index 7a49bee9d5..6bbc4c9b0c 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -20,10 +20,10 @@ - [[PR335]](https://github.com/lanl/singularity-eos/pull/335) Fix missing hermite.hpp in CMake install required for Helmholtz EOS - [[PR356]](https://github.com/lanl/singularity-eos/pull/356) Guard against FPEs in the PTE solver - [[PR356]](https://github.com/lanl/singularity-eos/pull/356) Update CMake for proper Kokkos linking in Fortran interface -- [[PR372]](https://github.com/lanl/singularity-eos/pull/372) Fix missing utilization of E0 in Davis Products EOS ### Changed (changing behavior/API/variables/...) - [[PR363]](https://github.com/lanl/singularity-eos/pull/363) Template lambda values for scalar calls +- [[PR372]](https://github.com/lanl/singularity-eos/pull/372) Removed E0 from Davis Products EOS in favor of using the shifted EOS modifier. CHANGES API! ### Infrastructure (changes irrelevant to downstream codes) - [[PR329]](https://github.com/lanl/singularity-eos/pull/329) Move vinet tests into analytic test suite From 4eac409d506655219d480993ac2d52f353f9531e Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 16:17:46 -0600 Subject: [PATCH 13/17] Remove E0 argument from Davis Products --- test/test_f_iface.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/test_f_iface.f90 b/test/test_f_iface.f90 index 52a92e95cc..7b845a2119 100644 --- a/test/test_f_iface.f90 +++ b/test/test_f_iface.f90 @@ -43,7 +43,7 @@ program test_sg_fortran_interface 4.6d0, 0.34d0, 0.56d0,0.d0, 0.4265d0, 0.001074d10) mat = mat + 1 res = init_sg_DavisProducts_f(mat, eos, 0.798311d0, 0.58d0, 1.35d0, 2.66182d0,& - 0.75419d0, 3.2d10, 0.001072d10, 0.d0) + 0.75419d0, 3.2d10, 0.001072d10) ! cleanup res = finalize_sg_eos_f(nmat, eos) From c265ca4b0c46da4834af0df2693016868b914424 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 16:56:19 -0600 Subject: [PATCH 14/17] Remove E0 from Davis Products EOS --- singularity-eos/eos/singularity_eos.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/singularity-eos/eos/singularity_eos.hpp b/singularity-eos/eos/singularity_eos.hpp index 9a7ddd43f7..5cec1491f4 100644 --- a/singularity-eos/eos/singularity_eos.hpp +++ b/singularity-eos/eos/singularity_eos.hpp @@ -39,8 +39,8 @@ int init_sg_Gruneisen(const int matindex, EOS *eos, const double C0, const doubl int init_sg_DavisProducts(const int matindex, EOS *eos, const double a, const double b, const double k, const double n, const double vc, - const double pc, const double Cv, const double E0, - int const *const enabled, double *const vals); + const double pc, const double Cv, int const *const enabled, + double *const vals); int init_sg_DavisReactants(const int matindex, EOS *eos, const double rho0, const double e0, const double P0, const double T0, From a6f09ba8d93a58da633a07210fed151b977596cc Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 16:58:15 -0600 Subject: [PATCH 15/17] Remove Davis Products EOS E0 parameter --- singularity-eos/eos/singularity_eos.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/singularity-eos/eos/singularity_eos.hpp b/singularity-eos/eos/singularity_eos.hpp index 5cec1491f4..9a31d10f54 100644 --- a/singularity-eos/eos/singularity_eos.hpp +++ b/singularity-eos/eos/singularity_eos.hpp @@ -134,7 +134,7 @@ int init_sg_Gruneisen(const int matindex, EOS *eos, const double C0, const doubl int init_sg_DavisProducts(const int matindex, EOS *eos, const double a, const double b, const double k, const double n, const double vc, - const double pc, const double Cv, const double E0); + const double pc, const double Cv); int init_sg_DavisReactants(const int matindex, EOS *eos, const double rho0, const double e0, const double P0, const double T0, From a0b5a9ffd7de6839c6179d90a12047bec682ae2d Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 18:22:19 -0600 Subject: [PATCH 16/17] Fix number of arguments to Davis Products EOS --- test/python_bindings.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/python_bindings.py b/test/python_bindings.py index a0ff479ad7..951e90523e 100644 --- a/test/python_bindings.py +++ b/test/python_bindings.py @@ -87,7 +87,7 @@ def testDavisReactants(self): eos = singularity_eos.DavisReactants(1,1,1,1,1,1,1,1,1,1,1) def testDavisProducts(self): - eos = singularity_eos.DavisProducts(1,1,1,1,1,1,1,1) + eos = singularity_eos.DavisProducts(1,1,1,1,1,1,1) class Modifiers(unittest.TestCase, EOSTestBase): def setUp(self): From 768f629c95178b08189fc999183b8b9439471574 Mon Sep 17 00:00:00 2001 From: Jeffrey H Peterson Date: Mon, 13 May 2024 19:18:18 -0600 Subject: [PATCH 17/17] Finally remove E0 from DavisProducts EOS pybinding init --- python/module.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/python/module.cpp b/python/module.cpp index 711e0990b3..2b89795552 100644 --- a/python/module.cpp +++ b/python/module.cpp @@ -60,7 +60,7 @@ PYBIND11_MODULE(singularity_eos, m) { eos_class(m, "DavisProducts") .def(py::init()) .def( - py::init(), + py::init(), py::arg("a"), py::arg("b"), py::arg("k"), py::arg("n"), py::arg("vc"), py::arg("pc"), py::arg("Cv") );