Merge pull request #387 from jfeist/fix_verner_interpolation

fix interpolation for GPUVern7 and GPUVern9

Author: ChrisRackauckas

src/ensemblegpukernel/integrators/nonstiff/interpolants.jl

@@ -54,48 +54,48 @@ end
     a1512, a1513, c16, a1601, a1604, a1605, a1606, a1607, a1608, a1609,
     a1611, a1612, a1613 = integ.tab.extra
 
-    @unpack k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, uprev, f, t, p = integ
+    @unpack k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, uprev, f, tprev, p = integ
 
     k11 = f(
         uprev +
         dt * (a1101 * k1 + a1104 * k4 + a1105 * k5 + a1106 * k6 +
          a1107 * k7 + a1108 * k8 + a1109 * k9),
         p,
-        t + c11 * dt)
+        tprev + c11 * dt)
     k12 = f(
         uprev +
         dt * (a1201 * k1 + a1204 * k4 + a1205 * k5 + a1206 * k6 +
          a1207 * k7 + a1208 * k8 + a1209 * k9 + a1211 * k11),
         p,
-        t + c12 * dt)
+        tprev + c12 * dt)
     k13 = f(
         uprev +
         dt * (a1301 * k1 + a1304 * k4 + a1305 * k5 + a1306 * k6 +
          a1307 * k7 + a1308 * k8 + a1309 * k9 + a1311 * k11 +
          a1312 * k12),
         p,
-        t + c13 * dt)
+        tprev + c13 * dt)
     k14 = f(
         uprev +
         dt * (a1401 * k1 + a1404 * k4 + a1405 * k5 + a1406 * k6 +
          a1407 * k7 + a1408 * k8 + a1409 * k9 + a1411 * k11 +
          a1412 * k12 + a1413 * k13),
         p,
-        t + c14 * dt)
+        tprev + c14 * dt)
     k15 = f(
         uprev +
         dt * (a1501 * k1 + a1504 * k4 + a1505 * k5 + a1506 * k6 +
          a1507 * k7 + a1508 * k8 + a1509 * k9 + a1511 * k11 +
          a1512 * k12 + a1513 * k13),
         p,
-        t + c15 * dt)
+        tprev + c15 * dt)
     k16 = f(
         uprev +
         dt * (a1601 * k1 + a1604 * k4 + a1605 * k5 + a1606 * k6 +
          a1607 * k7 + a1608 * k8 + a1609 * k9 + a1611 * k11 +
          a1612 * k12 + a1613 * k13),
         p,
-        t + c16 * dt)
+        tprev + c16 * dt)
 
     return y₀ +
            dt * (integ.k1 * b1Θ
@@ -161,82 +161,82 @@ end
     a2514, a2515, a2517, a2518, a2519, a2520, a2521, c26, a2601, a2608, a2609, a2610,
     a2611, a2612, a2613, a2614, a2615, a2617, a2618, a2619, a2620, a2621 = integ.tab.extra
 
-    @unpack k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, uprev, f, t, p = integ
+    @unpack k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, uprev, f, tprev, p = integ
 
     k11 = f(
         uprev +
         dt * (a1701 * k1 + a1708 * k2 + a1709 * k3 + a1710 * k4 +
          a1711 * k5 + a1712 * k6 + a1713 * k7 + a1714 * k8 + a1715 * k9),
-        p, t + c17 * dt)
+        p, tprev + c17 * dt)
     k12 = f(
         uprev +
         dt * (a1801 * k1 + a1808 * k2 + a1809 * k3 + a1810 * k4 +
          a1811 * k5 + a1812 * k6 + a1813 * k7 + a1814 * k8 +
          a1815 * k9 + a1817 * k11),
         p,
-        t + c18 * dt)
+        tprev + c18 * dt)
     k13 = f(
         uprev +
         dt * (a1901 * k1 + a1908 * k2 + a1909 * k3 + a1910 * k4 +
          a1911 * k5 + a1912 * k6 + a1913 * k7 + a1914 * k8 +
          a1915 * k9 + a1917 * k11 + a1918 * k12),
         p,
-        t + c19 * dt)
+        tprev + c19 * dt)
     k14 = f(
         uprev +
         dt * (a2001 * k1 + a2008 * k2 + a2009 * k3 + a2010 * k4 +
          a2011 * k5 + a2012 * k6 + a2013 * k7 + a2014 * k8 +
          a2015 * k9 + a2017 * k11 + a2018 * k12 + a2019 * k13),
         p,
-        t + c20 * dt)
+        tprev + c20 * dt)
     k15 = f(
         uprev +
         dt * (a2101 * k1 + a2108 * k2 + a2109 * k3 + a2110 * k4 +
          a2111 * k5 + a2112 * k6 + a2113 * k7 + a2114 * k8 +
          a2115 * k9 + a2117 * k11 + a2118 * k12 + a2119 * k13 +
          a2120 * k14),
         p,
-        t + c21 * dt)
+        tprev + c21 * dt)
     k16 = f(
         uprev +
         dt * (a2201 * k1 + a2208 * k2 + a2209 * k3 + a2210 * k4 +
          a2211 * k5 + a2212 * k6 + a2213 * k7 + a2214 * k8 +
          a2215 * k9 + a2217 * k11 + a2218 * k12 + a2219 * k13 +
          a2220 * k14 + a2221 * k15),
         p,
-        t + c22 * dt)
+        tprev + c22 * dt)
     k17 = f(
         uprev +
         dt * (a2301 * k1 + a2308 * k2 + a2309 * k3 + a2310 * k4 +
          a2311 * k5 + a2312 * k6 + a2313 * k7 + a2314 * k8 +
          a2315 * k9 + a2317 * k11 + a2318 * k12 + a2319 * k13 +
          a2320 * k14 + a2321 * k15),
         p,
-        t + c23 * dt)
+        tprev + c23 * dt)
     k18 = f(
         uprev +
         dt * (a2401 * k1 + a2408 * k2 + a2409 * k3 + a2410 * k4 +
          a2411 * k5 + a2412 * k6 + a2413 * k7 + a2414 * k8 +
          a2415 * k9 + a2417 * k11 + a2418 * k12 + a2419 * k13 +
          a2420 * k14 + a2421 * k15),
         p,
-        t + c24 * dt)
+        tprev + c24 * dt)
     k19 = f(
         uprev +
         dt * (a2501 * k1 + a2508 * k2 + a2509 * k3 + a2510 * k4 +
          a2511 * k5 + a2512 * k6 + a2513 * k7 + a2514 * k8 +
          a2515 * k9 + a2517 * k11 + a2518 * k12 + a2519 * k13 +
          a2520 * k14 + a2521 * k15),
         p,
-        t + c25 * dt)
+        tprev + c25 * dt)
     k20 = f(
         uprev +
         dt * (a2601 * k1 + a2608 * k2 + a2609 * k3 + a2610 * k4 +
          a2611 * k5 + a2612 * k6 + a2613 * k7 + a2614 * k8 +
          a2615 * k9 + a2617 * k11 + a2618 * k12 + a2619 * k13 +
          a2620 * k14 + a2621 * k15),
         p,
-        t + c26 * dt)
+        tprev + c26 * dt)
 
     return y₀ +
            dt *