275 |
double cv = td[i].cv * 1e3; /* J/kgK */ |
double cv = td[i].cv * 1e3; /* J/kgK */ |
276 |
double w = td[i].w; /* m/s */ |
double w = td[i].w; /* m/s */ |
277 |
double s = td[i].s * 1e3; /* J/kgK */ |
double s = td[i].s * 1e3; /* J/kgK */ |
278 |
//fprintf(stderr,"T = %f, rho = %f, p = %f\n",T,rho,p); |
//fprintf(stderr,"T = %f, rho = %f, p = %f, w = %f, wcalc = %f\n",T,rho,p,w, helmholtz_w(T,rho,d)); |
279 |
ASSERT_TOL(helmholtz_s, T, rho, d, s, s*1e-8); |
ASSERT_TOL(helmholtz_s, T, rho, d, s, s*1e-8); |
280 |
ASSERT_TOL(helmholtz_p, T, rho, d, p, p*1e-8); |
ASSERT_TOL(helmholtz_p, T, rho, d, p, p*1e-8); |
281 |
ASSERT_TOL(helmholtz_cv, T, rho, d, cv, cv*1e-8); |
ASSERT_TOL(helmholtz_cv, T, rho, d, cv, cv*1e-8); |
282 |
|
ASSERT_TOL(helmholtz_w, T, rho, d, w, w*1e-8); |
283 |
} |
} |
284 |
|
|
285 |
fprintf(stderr,"Tests completed OK (maximum error = %0.2f%%)\n",maxerr); |
fprintf(stderr,"Tests completed OK (maximum error = %0.2f%%)\n",maxerr); |