| 42 | } |
| 43 | |
| 44 | SpecificHeatCPDataPoint::SpecificHeatCPDataPoint(EOSFitter* EOS, double T, double rho, double cp, double variance) { |
| 45 | this->EOS = EOS; |
| 46 | this->T = T; |
| 47 | this->rho = rho; |
| 48 | this->cp = cp; |
| 49 | this->cp_over_R = this->cp / EOS->R; |
| 50 | this->tau = EOS->Tr / this->T; |
| 51 | this->delta = this->rho / EOS->rhor; |
| 52 | this->log_tau = log(tau); |
| 53 | this->log_delta = log(delta); |
| 54 | this->variance = variance; |
| 55 | } |
| 56 | /// The part that does not depend on the coefficients |
| 57 | /// Here it requires that the coefficients be passed in to calculate the precorrelation factor |
| 58 | double SpecificHeatCPDataPoint::a_0(const std::vector<double>& n) { |