diff --git a/src/fourier/details/euclidian_maps.hpp b/src/fourier/details/euclidian_maps.hpp index e0bd939..b42320c 100644 --- a/src/fourier/details/euclidian_maps.hpp +++ b/src/fourier/details/euclidian_maps.hpp @@ -131,6 +131,20 @@ namespace CosmoTool } } + template + void get_Kvec(long p, Array2& kvec) + { + const DimArray& dims = this->getDims(); + DimArray d(delta_k.size()); + for (int q = 0; q < d.size(); q++) + { + d[q] = p%dims[q]; + p = (p-d[q])/dims[q]; + } + return get_Kvec(d, kvec); + } + + template double get_K(const Array& ik) {