| 435 | } |
| 436 | |
| 437 | std::pair<double, double> elliptic2cart(double theta, double rho_w, |
| 438 | double rho_h) { |
| 439 | return std::make_pair(rho_w * cos(theta), rho_h * sin(theta)); |
| 440 | } |
| 441 | |
| 442 | std::pair<vector_1d, vector_1d> |
| 443 | elliptic2cart(vector_1d theta, vector_1d rho_w, vector_1d rho_h) { |