Skip to content
Snippets Groups Projects
Commit 4308fa88 authored by Ian Bell's avatar Ian Bell
Browse files

Few more template

parent 3f9561f4
No related branches found
No related tags found
No related merge requests found
...@@ -107,12 +107,12 @@ auto extrapolate_from_critical(const Model& model, const Scalar Tc, const Scalar ...@@ -107,12 +107,12 @@ auto extrapolate_from_critical(const Model& model, const Scalar Tc, const Scalar
using tdx = TDXDerivatives<Model>; using tdx = TDXDerivatives<Model>;
auto z = (Eigen::ArrayXd(1) << 1.0).finished(); auto z = (Eigen::ArrayXd(1) << 1.0).finished();
auto R = model.R; auto R = model.R;
auto ders = tdx::get_Ar0n<4>(model, Tc, rhoc, z); auto ders = tdx::template get_Ar0n<4>(model, Tc, rhoc, z);
auto dpdrho = R*Tc*(1 + 2 * ders[1] + ders[2]); // Should be zero auto dpdrho = R*Tc*(1 + 2 * ders[1] + ders[2]); // Should be zero
auto d2pdrho2 = R*Tc/rhoc*(2 * ders[1] + 4 * ders[2] + ders[3]); // Should be zero auto d2pdrho2 = R*Tc/rhoc*(2 * ders[1] + 4 * ders[2] + ders[3]); // Should be zero
auto d3pdrho3 = R*Tc/(rhoc*rhoc)*(6 * ders[2] + 6 * ders[3] + ders[4]); auto d3pdrho3 = R*Tc/(rhoc*rhoc)*(6 * ders[2] + 6 * ders[3] + ders[4]);
auto Ar11 = tdx::get_Ar11(model, Tc, rhoc, z); auto Ar11 = tdx::template get_Ar11(model, Tc, rhoc, z);
auto Ar12 = tdx::get_Ar12(model, Tc, rhoc, z); auto Ar12 = tdx::template get_Ar12(model, Tc, rhoc, z);
auto d2pdrhodT = R * (1 + 2 * ders[1] + ders[2] - 2 * Ar11 - Ar12); auto d2pdrhodT = R * (1 + 2 * ders[1] + ders[2] - 2 * Ar11 - Ar12);
auto Brho = sqrt(6*d2pdrhodT*Tc/d3pdrho3); auto Brho = sqrt(6*d2pdrhodT*Tc/d3pdrho3);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment