fixed index_set helper function

Author: passscoed

src/pyGroupedTransforms/NFCTtools.py

@@ -77,8 +77,7 @@ def nfct_index_set(
     if d == 1:
         return np.array([i for i in range(0, bandwidths[0])], "int")
 
-    bandwidths = bandwidths[::-1]
-    tmp = tuple([list(range(0, bw)) for bw in bandwidths])
+    tmp = tuple([list(range(0, bw)) for bw in bandwidths[::-1]])
     tmp = itertools.product(*(tmp[::-1]))
 
     freq = np.empty((d, np.prod(bandwidths)), dtype=int)

src/pyGroupedTransforms/NFFTtools.py

@@ -87,8 +87,9 @@ def nfft_index_set(
         return np.array(
             list(range(int(-bandwidths[0] / 2), int(bandwidths[0] / 2))), dtype="int"
         )
-
-    tmp = [list(range(int(-bw / 2), int(bw / 2))) for bw in bandwidths]
+    
+    
+    tmp = [list(range(int(-bw / 2), int(bw / 2))) for bw in bandwidths[::-1]]
     tmp = itertools.product(*(tmp[::-1]))
 
     freq = np.empty((d, np.prod(bandwidths)), dtype=int)
@@ -134,7 +135,6 @@ def get_transform(
     # Output:
      * `F::LinearMap{ComplexF64}` ... Linear map of the Fourier-transform implemented by the NFFT
     """
-
     if bandwidths.ndim > 1 or bandwidths.dtype != "int32":
         return "Please use an zero or one-dimensional numpy.array with dtype 'int32' as input"