A_mul_B! to mul!

Author: ChrisRackauckas

src/ode/brusselator_prob.jl

@@ -109,11 +109,11 @@ function brusselator_1d(du, u_, p, t)
     A, B, α, buffer = p
     u = @view(u_[:, 1])
     v = @view(u_[:, 2])
-    A_mul_B!(buffer, D_brusselator_u, u)
+    mul!(buffer, D_brusselator_u, u)
     Du = buffer
     @. du[:, 1] = A + u^2*v - (B+1)*u + α*Du
 
-    A_mul_B!(buffer, D_brusselator_v, v)
+    mul!(buffer, D_brusselator_v, v)
     Dv = buffer
     @. du[:, 2] = B*u - u^2*v + α*Dv
     nothing

src/ode/filament_prob.jl

@@ -155,15 +155,15 @@ function (f::FilamentCache)(dr, r, p, t)
     A, P, F, S1, S2 = f.A, f.P.P, f.F.F, f.Sc.S1, f.Sc.S2
 
     # implement dr = P * (A*r + F) in an optimized way to avoid temporaries
-    A_mul_B!(S1, A, r)
+    mul!(S1, A, r)
     S1 .+= F
-    A_mul_B!(S2, P, S1)
+    mul!(S2, P, S1)
     copy!(dr, S2)
     return dr
 end
 
 function (f::FilamentCache)(::Type{Val{:jac}}, J, r, p, t)
-    A_mul_B!(J, f.P.P, f.A)
+    mul!(J, f.P.P, f.A)
     nothing
 end
 

src/ode/ode_simple_nonlinear_prob.jl

@@ -244,7 +244,7 @@ prob_ode_pleiades = ODEProblem(pleiades,[3.0,3.0,-1.0,-3.0,2.0,-2.0,2.0,3.0,-3.0
 srand(100)
 const mm_A = rand(4,4)
 mm_linear = function (du,u,p,t)
-  A_mul_B!(du,mm_A,u)
+  mul!(du,mm_A,u)
 end
 const MM_linear =full(Diagonal(0.5ones(4)))
 (::typeof(mm_linear))(::Type{Val{:analytic}},u0,p,t) = expm(inv(MM_linear)*mm_A*t)*u0

src/sde_premade_problems.jl

@@ -389,7 +389,7 @@ function generate_stiff_stoch_heat(D=1,k=1;N = 100)
     dx = 1/N
     A = D/(dx^2) * A
     function f(du,u,p,t)
-        A_mul_B!(du,A,u)
+        mul!(du,A,u)
     end
     #=
     function f(::Type{Val{:analytic}},u0,p,t,W)