| 468 | } |
| 469 | } |
| 470 | double AbstractCubic::PI_12(double delta, const std::vector<double>& x, std::size_t idelta) { |
| 471 | double bm = bm_term(x); |
| 472 | double cm = cm_term(); |
| 473 | switch (idelta) { |
| 474 | case 0: |
| 475 | return (1 + (Delta_1 * bm + cm) * rho_r * delta) * (1 + (Delta_2 * bm + cm) * rho_r * delta); |
| 476 | case 1: |
| 477 | return rho_r * (2 * (bm * Delta_1 + cm) * (bm * Delta_2 + cm) * delta * rho_r + (Delta_1 + Delta_2) * bm + 2 * cm); |
| 478 | case 2: |
| 479 | return 2 * (Delta_1 * bm + cm) * (Delta_2 * bm + cm) * rho_r * rho_r; |
| 480 | case 3: |
| 481 | return 0; |
| 482 | case 4: |
| 483 | return 0; |
| 484 | default: |
| 485 | throw -1; |
| 486 | } |
| 487 | } |
| 488 | double AbstractCubic::d_PI_12_dxi(double delta, const std::vector<double>& x, std::size_t idelta, std::size_t i, bool xN_independent) { |
| 489 | double bm = bm_term(x); |
| 490 | double cm = cm_term(); |
nothing calls this directly
no outgoing calls
no test coverage detected