fix interplate module and opt attributes

Author: 9il

source/mir/interpolate/constant.d

@@ -183,6 +183,8 @@ struct Constant(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators = Rep
         this._ownsData = false;
     }
 
+@trusted:
+
     ///
     GridVectors[dimension] grid(size_t dimension = 0)() const @property
         if (dimension < N)
@@ -232,7 +234,7 @@ struct Constant(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators = Rep
                 else
                     indexes[i] = this.findInterval!i(xs[i]);
             }
-            return _data[indexes];
+            return _data[][indexes];
         }
     }
 }

source/mir/interpolate/linear.d

@@ -269,6 +269,8 @@ struct Linear(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators...)
         this._ownsData = false;
     }
 
+@trusted:
+
     ///
     GridVectors[dimension] grid(size_t dimension = 0)() const @property
         if (dimension < N)
@@ -362,7 +364,7 @@ struct Linear(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators...)
                 static if (rp2d == 1)
                     shuffle3!1(local[1][0 .. L], local[1][L .. 2 * L], local[0][0 .. L], local[0][L .. 2 * L]);
                 static if (i + 1 == N)
-                    return *cast(SplineReturnType!(F, N, rp2d)*) local[0].ptr;
+                    return *cast(SplineReturnType!(F, N, 2 ^^ rp2d)*) local[0].ptr;
             }
         }
     }

source/mir/interpolate/package.d

@@ -67,10 +67,7 @@ template findInterval(size_t dimension = 0)
             auto grid = interpolant.grid[1 .. $][0 .. len];
         }
         assert(len >= 0);
-        return len - grid
-            .assumeSorted
-            .upperBound(x)
-            .length;
+        return len - grid.assumeSorted.upperBound(x).length;
     }
 }
 
@@ -262,10 +259,10 @@ auto copyvec(F, size_t N)(ref const F[N] from, ref F[N] to)
 
 package template SplineReturnType(F, size_t N, size_t P)
 {
-    static if (P == 0 || N == 0)
+    static if (P <= 1 || N == 0)
         alias SplineReturnType = F;
     else
-        alias SplineReturnType = .SplineReturnType!(F, N - 1, P)[2 ^^ P];
+        alias SplineReturnType = .SplineReturnType!(F, N - 1, P)[P];
 }
 
 template generateShuffles3(size_t N, size_t P)

source/mir/interpolate/spline.d

@@ -566,12 +566,12 @@ struct Spline(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators...)
             // X.length == N || derivative == 0 && X.length && X.length <= N
         {
             auto d4 = this.opCall!3(xs);
-            SplineReturnType!(F, N, 2) d3;
+            SplineReturnType!(F, N, 3) d3;
             void fun(size_t d, A, B)(ref A a, ref B b)
             {
                 static if (d)
                     foreach(i; Iota!3)
-                        fun!(d - 1)(a, b);
+                        fun!(d - 1)(a[i], b[i]);
                 else
                     b = a;
             }
@@ -653,7 +653,7 @@ struct Spline(F, size_t N = 1, FirstGridIterator = F*, NextGridIterators...)
                     );
                 static if (i + 1 == N)
                 {
-                    return *cast(SplineReturnType!(F, N, rp2d)*) local[0].ptr;
+                    return *cast(SplineReturnType!(F, N, 2 ^^ rp2d)*) local[0].ptr;
                 }
                 else
                 {

source/mir/math/common.d

@@ -69,7 +69,7 @@ version(LDC)
 }
 else
 {
-    enum { fastmath, fusedmath, optmath };
+    enum { fmamath, optmath, fastmath };
 }
 
 version(LDC)