| 21 | } |
| 22 | |
| 23 | PressureDataPoint::PressureDataPoint(EOSFitter* EOS, double T, double rho, double p, double variance) { |
| 24 | this->EOS = EOS; |
| 25 | this->T = T; |
| 26 | this->rho = rho; |
| 27 | this->p = p; |
| 28 | this->tau = EOS->Tr / this->T; |
| 29 | this->delta = this->rho / EOS->rhor; |
| 30 | this->log_tau = log(tau); |
| 31 | this->log_delta = log(delta); |
| 32 | this->variance = variance; |
| 33 | } |
| 34 | /// The part that does not depend on the coefficients |
| 35 | double PressureDataPoint::a_0() { |
| 36 | double rhoRT = this->rho * EOS->R * this->T; |