Merge pull request #96 from isaacsas/discretefunction-fix

DiscreteFunction creation fix

Author: isaacsas

src/problem.jl

@@ -12,8 +12,8 @@ end
 
 ### DiscreteProblem, passing through syms
 function DiffEqBase.DiscreteProblem(rn::DiffEqBase.AbstractReactionNetwork, u0, tspan::Tuple, p=nothing; kwargs...)
-    f = DiffEqBase.DISCRETE_INPLACE_DEFAULT
-    df = DiscreteFunction{true,typeof(f),Nothing,typeof(rn.syms)}(f,nothing,rn.syms)
+    f = (du,u,p,t) -> du.=u    # identity function to make syms works
+    df = DiscreteFunction(f,syms=rn.syms)
     DiscreteProblem(df, u0, tspan, p; kwargs...)
 end