| 66 | } |
| 67 | |
| 68 | std::vector<double> linspace(double d1, double d2, size_t n) { |
| 69 | std::vector<double> x(n); |
| 70 | // avoid division by zero |
| 71 | double step = n > 1? (d2 - d1) / static_cast<double>(n - 1) : 0.; |
| 72 | for (size_t i = 0; i < x.size(); ++i) { |
| 73 | x[i] = d1 + static_cast<double>(i) * step; |
| 74 | } |
| 75 | return x; |
| 76 | } |
| 77 | |
| 78 | std::vector<double> linspace(double d1, double d2) { |
| 79 | return linspace(d1, d2, 100); |
no test coverage detected