format

Author: ChrisRackauckas

src/dense/interpolants.jl

@@ -3000,23 +3000,26 @@ end
     @dprkn6pre0
     halfsize = length(y₀) ÷ 2
     isfirsthalf = idxs .<= halfsize
-    secondhalf =  idxs .> halfsize
+    secondhalf = idxs .> halfsize
     firstidxs = idxs[isfirsthalf]
     secondidxs_shifted = idxs[secondhalf]
     secondidxs = secondidxs_shifted .- halfsize
 
     @views @.. broadcast=false out[secondhalf]=uprev[secondidxs] +
-                                           dt * Θ *
-                                           (duprev[secondidxs] +
-                                            dt * Θ *
-                                            (b1Θ * k1[secondidxs] +
-                                             b3Θ * k3[secondidxs] +
-                                             b4Θ * k4[secondidxs] + b5Θ * k5[secondidxs] +
-                                             b6Θ * k6[secondidxs]))
+                                               dt * Θ *
+                                               (duprev[secondidxs] +
+                                                dt * Θ *
+                                                (b1Θ * k1[secondidxs] +
+                                                 b3Θ * k3[secondidxs] +
+                                                 b4Θ * k4[secondidxs] +
+                                                 b5Θ * k5[secondidxs] +
+                                                 b6Θ * k6[secondidxs]))
     @views @.. broadcast=false out[isfirsthalf]=duprev[firstidxs] +
-                                           dt * Θ *
-                                           (bp1Θ * k1[firstidxs] + bp3Θ * k3[firstidxs] +
-                                            bp4Θ * k4[firstidxs] + bp5Θ * k5[firstidxs] +
-                                            bp6Θ * k6[firstidxs])
+                                                dt * Θ *
+                                                (bp1Θ * k1[firstidxs] +
+                                                 bp3Θ * k3[firstidxs] +
+                                                 bp4Θ * k4[firstidxs] +
+                                                 bp5Θ * k5[firstidxs] +
+                                                 bp6Θ * k6[firstidxs])
     out
 end

src/nlsolve/utils.jl

@@ -308,7 +308,8 @@ function build_nlsolver(
             end
             prob = NonlinearProblem(NonlinearFunction(nlf), copy(ztmp), nlp_params)
             cache = init(prob, nlalg.alg)
-            nlcache = NonlinearSolveCache(nothing, tstep, nothing, nothing, invγdt, prob, cache)
+            nlcache = NonlinearSolveCache(
+                nothing, tstep, nothing, nothing, invγdt, prob, cache)
         else
             nlcache = NLNewtonConstantCache(tstep, J, W, true, true, true, tType(dt), uf,
                 invγdt, tType(nlalg.new_W_dt_cutoff), t)

test/interface/interpolation_output_types.jl

@@ -63,7 +63,7 @@ sol = solve(prob, DPRKN6())
 @test sol(0.5) isa RecursiveArrayTools.ArrayPartition
 @test sol(0.5; idxs = 1) isa Number
 @test sol(0.5; idxs = [1]) isa Vector
-@test sol(0.5; idxs = [1,2]) isa Vector
-@test Vector(sol(0.5)) == sol(0.5; idxs = [1,2])
-@test reverse(Vector(sol(0.5))) == sol(0.5; idxs = [2,1])
+@test sol(0.5; idxs = [1, 2]) isa Vector
+@test Vector(sol(0.5)) == sol(0.5; idxs = [1, 2])
+@test reverse(Vector(sol(0.5))) == sol(0.5; idxs = [2, 1])
 @test Vector(sol(0.5)) == [sol(0.5; idxs = 1); sol(0.5; idxs = 2)]