| 62 | |
| 63 | template<typename T> |
| 64 | void gaussian1D(T* out, const int dim, double sigma = 0.0) { |
| 65 | if (!(sigma > 0)) sigma = 0.25 * dim; |
| 66 | |
| 67 | T sum = (T)0; |
| 68 | for (int i = 0; i < dim; i++) { |
| 69 | int x = i - (dim - 1) / 2; |
| 70 | T el = 1. / sqrt(2 * PI_VAL * sigma * sigma) * |
| 71 | exp(-((x * x) / (2 * (sigma * sigma)))); |
| 72 | out[i] = el; |
| 73 | sum += el; |
| 74 | } |
| 75 | |
| 76 | for (int k = 0; k < dim; k++) out[k] /= sum; |
| 77 | } |
| 78 | |
| 79 | template<typename T> |
| 80 | std::array<Kernel, 4> getOrbKernels() { |