| 33 | inline T square(const T & x) { return x * x; } |
| 34 | |
| 35 | inline double blep(const double t, const double dt) |
| 36 | { |
| 37 | if (t < dt) return -square(t / dt - 1); |
| 38 | else if (t > 1 - dt) return square((t - 1) / dt + 1); |
| 39 | else return 0.0; |
| 40 | } |
| 41 | |
| 42 | inline double blamp(double t, const double dt) |
| 43 | { |