use nameof and operation to get symbols

Author: isaacsas

src/systems/diffeqs/abstractodesystem.jl

@@ -163,7 +163,7 @@ function DiffEqBase.ODEFunction{iip}(sys::AbstractODESystem, dvs = states(sys),
                      tgrad = _tgrad === nothing ? nothing : _tgrad,
                      mass_matrix = _M,
                      jac_prototype = sparse ? similar(sys.jac[],Float64) : nothing,
-                     syms = Symbol.(states(sys)))
+                     syms = nameof.(operation.(states(sys))))
 end
 
 """
@@ -225,7 +225,7 @@ function ODEFunctionExpr{iip}(sys::AbstractODESystem, dvs = states(sys),
                          tgrad = tgrad,
                          mass_matrix = M,
                          jac_prototype = $jp_expr,
-                         syms = $(Symbol.(states(sys))))
+                         syms = $(nameof.(operation.(states(sys)))))
     end
     !linenumbers ? striplines(ex) : ex
 end

src/systems/diffeqs/sdesystem.jl

@@ -196,7 +196,7 @@ function DiffEqBase.SDEFunction{iip}(sys::SDESystem, dvs = sys.states, ps = sys.
                      Wfact = _Wfact === nothing ? nothing : _Wfact,
                      Wfact_t = _Wfact_t === nothing ? nothing : _Wfact_t,
                      mass_matrix = _M,
-                     syms = Symbol.(sys.states))
+                     syms = nameof.(operation.(states(sys))))
 end
 
 function DiffEqBase.SDEFunction(sys::SDESystem, args...; kwargs...)
@@ -269,7 +269,7 @@ function SDEFunctionExpr{iip}(sys::SDESystem, dvs = states(sys),
                          Wfact = Wfact,
                          Wfact_t = Wfact_t,
                          mass_matrix = M,
-                         syms = $(Symbol.(states(sys))),kwargs...)
+                         syms = $(nameof.(operation.(states(sys)))),kwargs...)
     end
     !linenumbers ? striplines(ex) : ex
 end

src/systems/jumps/jumpsystem.jl

@@ -197,7 +197,7 @@ function DiffEqBase.DiscreteProblem(sys::JumpSystem, u0map, tspan::Tuple,
         p = parammap
     end
     f  = DiffEqBase.DISCRETE_INPLACE_DEFAULT
-    df = DiscreteFunction(f, syms=Symbol.(states(sys)))
+    df = DiscreteFunction{true,true}(f, syms=nameof.(operation.(states(sys))))
     DiscreteProblem(df, u0, tspan, p; kwargs...)
 end
 
@@ -230,7 +230,7 @@ function DiscreteProblemExpr(sys::JumpSystem, u0map, tspan::Tuple,
         u0 = $u0
         p = $p
         tspan = $tspan
-        df = DiscreteFunction(f, syms=$(Symbol.(states(sys))))
+        df = DiscreteFunction{true,true}(f, syms=$(nameof.(operation.(states(sys)))))
         DiscreteProblem(df, u0, tspan, p; kwargs...)
     end
 end