fixes

Co-authored-by: "Yingbo Ma" <[email protected]>

Author: shashi

src/systems/diffeqs/sdesystem.jl

@@ -166,6 +166,9 @@ function DiffEqBase.SDEFunction{iip}(sys::SDESystem, dvs = states(sys), ps = par
                                      u0 = nothing;
                                      version = nothing, tgrad=false, sparse = false,
                                      jac = false, Wfact = false, eval_expression = true, kwargs...) where {iip}
+    dvs = scalarize.(dvs)
+    ps = scalarize.(ps)
+
     f_gen = generate_function(sys, dvs, ps; expression=Val{eval_expression}, kwargs...)
     f_oop,f_iip = eval_expression ? (@RuntimeGeneratedFunction(ex) for ex in f_gen) : f_gen
     g_gen = generate_diffusion_function(sys, dvs, ps; expression=Val{eval_expression}, kwargs...)

src/systems/reaction/reactionsystem.jl

@@ -461,11 +461,13 @@ This input may contain repeat parameters.
 function Base.convert(::Type{<:SDESystem}, rs::ReactionSystem; 
                       noise_scaling=nothing, name=nameof(rs), combinatoric_ratelaws=true, kwargs...)
 
-    if noise_scaling isa Vector
+    if noise_scaling isa AbstractArray
         (length(noise_scaling)!=length(equations(rs))) &&
         error("The number of elements in 'noise_scaling' must be equal " *
               "to the number of reactions in the reaction system.")
-        noise_scaling = value.(noise_scaling)
+        if !(noise_scaling isa Symbolics.Arr)
+            noise_scaling = value.(noise_scaling)
+        end
     elseif !isnothing(noise_scaling)
         noise_scaling = fill(value(noise_scaling),length(equations(rs)))
     end