@@ -61,9 +61,9 @@ class SimpleMACAW : public EosBase<SimpleMACAW> {
61
61
PORTABLE_INLINE_FUNCTION void CheckParams () const {
62
62
PORTABLE_ALWAYS_REQUIRE (_A >= 0 , " Parameter, 'A', must be non-negative" );
63
63
PORTABLE_ALWAYS_REQUIRE (_B >= 0 , " Parameter, 'B', must be non-negative" );
64
+ PORTABLE_ALWAYS_REQUIRE (_Cvinf > 0 , " Specific heat capacity (at constant volume), `Cvinf`, must be positive" );
64
65
PORTABLE_ALWAYS_REQUIRE (_v0 > 0 , " Reference specific volume, 'v0', must be positive" );
65
66
PORTABLE_ALWAYS_REQUIRE (_T0 >= 0 , " Reference temperature, 'T0', must be non-negative" );
66
- PORTABLE_ALWAYS_REQUIRE (_Cvinf > 0 , " Specific heat capacity (at constant volume), `Cvinf`, must be positive" );
67
67
PORTABLE_ALWAYS_REQUIRE (_Gc > 0 && _Gc < 1 , " Gruneisen parameter, 'Gc', must be in the interval (0,1)" );
68
68
_AZbar.CheckParams ();
69
69
}
@@ -249,9 +249,9 @@ class SimpleMACAW : public EosBase<SimpleMACAW> {
249
249
250
250
static constexpr unsigned long PreferredInput () { return _preferred_input; }
251
251
PORTABLE_INLINE_FUNCTION void PrintParams () const {
252
- printf (" Simple MACAW Parameters:\n A = %g\n B = %g\n Cvinf = %g\n Gc = "
252
+ printf (" Simple MACAW Parameters:\n A = %g\n B = %g\n Cvinf = %g\n v0 = %g \n T0 = %g \ n Gc = "
253
253
" %g\n " ,
254
- _A, _B, _Cvinf, _Gc);
254
+ _A, _B, _Cvinf, _v0, _T0, _Gc);
255
255
_AZbar.PrintParams ();
256
256
}
257
257
template <typename Indexer_t>
@@ -287,9 +287,7 @@ class SimpleMACAW : public EosBase<SimpleMACAW> {
287
287
static std::string EosPyType () { return EosType (); }
288
288
289
289
private:
290
- Real _A, _B, _Cvinf, _Gc;
291
- // reference values
292
- Real _v0, _T0;
290
+ Real _A, _B, _Cvinf, _v0, _T0, _Gc;
293
291
MeanAtomicProperties _AZbar;
294
292
static constexpr const unsigned long _preferred_input =
295
293
thermalqs::density | thermalqs::specific_internal_energy;
0 commit comments