Commit 2b68adfd authored by Martin Reinecke's avatar Martin Reinecke
Browse files

shiny new kernels

parent e03057f4
Pipeline #78850 canceled with stages
in 3 minutes and 45 seconds
...@@ -207,7 +207,6 @@ template<typename T> class Interpolator ...@@ -207,7 +207,6 @@ template<typename T> class Interpolator
vector<size_t> getIdx(const mav<T,2> &ptg) const vector<size_t> getIdx(const mav<T,2> &ptg) const
{ {
vector<size_t> idx(ptg.shape(0));
constexpr size_t cellsize=16; constexpr size_t cellsize=16;
size_t nct = ntheta/cellsize+1, size_t nct = ntheta/cellsize+1,
ncp = nphi/cellsize+1; ncp = nphi/cellsize+1;
...@@ -219,6 +218,7 @@ template<typename T> class Interpolator ...@@ -219,6 +218,7 @@ template<typename T> class Interpolator
// MR_assert((itheta<nct)&&(iphi<ncp), "oops"); // MR_assert((itheta<nct)&&(iphi<ncp), "oops");
mapper[itheta*ncp+iphi].push_back(i); mapper[itheta*ncp+iphi].push_back(i);
} }
vector<size_t> idx(ptg.shape(0));
size_t cnt=0; size_t cnt=0;
for (auto &vec: mapper) for (auto &vec: mapper)
{ {
......
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment