Merge pull request #151 from YingboMa/myb/array

Make oop ODEFunction Array compatible

Author: ChrisRackauckas

src/utils.jl

@@ -54,7 +54,7 @@ function build_function(rhss, vs, ps, args = (), conv = rhs -> convert(Expr, rhs
         end
         function $fname(u,p,$(args...))
             X = $let_expr
-            T = $(constructor === nothing ? :(StaticArrays.similar_type(typeof(u), eltype(X))) : constructor)
+            T = $(constructor === nothing ? :(u isa StaticArrays.StaticArray ? StaticArrays.similar_type(typeof(u), eltype(X)) : x->(du=similar(u, eltype(X)); du .= x)) : constructor)
             T(X)
         end
     end

test/system_construction.jl

@@ -161,6 +161,7 @@ de = ODESystem(eqs)
     du = [0.0,0.0]
     f(du, [1.0,2.0], [1,2,3], 0.0)
     du ≈ [-1, -1/3]
+    du ≈ f([1.0,2.0], [1,2,3], 0.0)
 end
 
 # Now nonlinear system with only variables