| 27 | } |
| 28 | |
| 29 | inline double norm(const vector<double>& a) { |
| 30 | double accum = 0.0; |
| 31 | for (size_t i = 0; i < a.size(); ++i) { |
| 32 | accum += a[i] * a[i]; |
| 33 | } |
| 34 | return sqrt(accum); |
| 35 | }; |
| 36 | |
| 37 | ExtendedSystem extend_system( |
| 38 | py::array_t<double> positions, |