23 static const double ni[34] = {
24 1.4632971213167E-01, -8.4548187169114E-01, -3.7563603672040E+00,
25 3.3855169168385E+00, -9.5791963387872E-01, 1.5772038513228E-01,
26 -1.6616417199501E-02, 8.1214629983568E-04, 2.8319080123804E-04,
27 -6.0706301565874E-04, -1.8990068218419E-02, -3.2529748770505E-02,
28 -2.1841717175414E-02, -5.2838357969930E-05, -4.7184321073267E-04,
29 -3.0001780793026E-04, 4.7661393906987E-05, -4.4141845330846E-06,
30 -7.2694996297594E-16, -3.1679644845054E-05, -2.8270797985312E-06,
31 -8.5205128120103E-10, -2.2425281908000E-06, -6.5171222895601E-07,
32 -1.4341729937924E-13, -4.0516996860117E-07, -1.2734301741641E-09,
33 -1.7424871230634E-10, -6.8762131295531E-19, 1.4478307828521E-20,
34 2.6335781662795E-23, -1.1947622640071E-23, 1.8228094581404E-24,
35 -9.35370872924580E-26};
37 static const double li[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
38 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4,
39 4, 5, 8, 8, 21, 23, 29, 30, 31, 32};
41 static const double ji[34] = {
42 -2, -1, 0, 1, 2, 3, 4, 5, -9, -7, -1, 0, 1, 3, -3, 0, 1,
43 3, 17, -4, 0, 6, -5, -2, 10, -8, -11, -6, -29, -31, -38, -39, -40, -41};
49 for (
int i = 0; i < 34; i++)
51 val +=
ni[i] * std::pow(7.1 - pi,
li[i]) * std::pow(tau - 1.222,
ji[i]);
61 for (
int i = 0; i < 34; i++)
63 val +=
ni[i] *
ji[i] * std::pow(7.1 - pi,
li[i]) *
64 std::pow(tau - 1.222,
ji[i] - 1.0);
71 const double tau,
const double pi)
74 for (
int i = 0; i < 34; i++)
76 val +=
ni[i] *
ji[i] * (
ji[i] - 1.0) * std::pow(7.1 - pi,
li[i]) *
77 std::pow(tau - 1.222,
ji[i] - 2.0);
87 for (
int i = 0; i < 34; i++)
89 val -=
ni[i] *
li[i] * std::pow(7.1 - pi,
li[i] - 1.0) *
90 std::pow(tau - 1.222,
ji[i]);
100 for (
int i = 0; i < 34; i++)
102 val +=
ni[i] *
li[i] * (
li[i] - 1.0) * std::pow(7.1 - pi,
li[i] - 2.0) *
103 std::pow(tau - 1.222,
ji[i]);
110 const double tau,
const double pi)
113 for (
int i = 0; i < 34; i++)
115 val -=
ni[i] *
ji[i] *
li[i] * std::pow(7.1 - pi,
li[i] - 1.0) *
116 std::pow(tau - 1.222,
ji[i] - 1.0);
Declare a class for dimensionless Gibbs free energy, region1.
static const double ni[34]
static const double li[34]
static const double ji[34]
static double get_dgamma_dtau(const double tau, const double pi)
static double get_dgamma_dpi(const double tau, const double pi)
static double get_dgamma_dpi_dpi(const double tau, const double pi)
static double get_dgamma_dtau_dtau(const double tau, const double pi)
static double get_dgamma_dtau_dpi(const double tau, const double pi)
static double get_gamma(const double tau, const double pi)