Fix index unraveling

src/interp.rs

@@ -22,7 +22,7 @@ where
     // Interpolate over x
     for iy in 0..2 {
         for iz in 0..2 {
-            let iix = iy + (2 * iz);
+            let iix = (2 * iy) + iz;
             c[iix] = values[[0, iy, iz]] * m_x[[0]] + values[[1, iy, iz]] * x[[0]];
         }
     }