Merge pull request #355 from SciML/bgc/fix_sampleddata

fixed SampledData initializaiton

Author: ChrisRackauckas

src/Blocks/sources.jl

@@ -672,10 +672,7 @@ data input component.
     eqs = [
         output.u ~ get_sampled_data(t, p)
     ]
-    return ODESystem(eqs, t, vars, [pars; p]; name, systems,
-        defaults = [
-            output.u => get_sampled_data(0.0, p)
-        ])
+    return ODESystem(eqs, t, vars, [pars; p]; name, systems)
 end
 
 """
@@ -701,8 +698,7 @@ data input component.
     eqs = [
         output.u ~ get_sampled_data(t, buffer)
     ]
-    return ODESystem(eqs, t, vars, pars; name, systems,
-        defaults = [output.u => get_sampled_data(0.0, buffer)])
+    return ODESystem(eqs, t, vars, pars; name, systems)
 end
 
 SampledData(x::SampledDataType.Option; kwargs...) = SampledData(Val(x); kwargs...)

test/Blocks/sources.jl

@@ -424,7 +424,7 @@ end
     x = @. time^2 + 1.0
 
     @testset "using Parameter type" begin
-        vars = @variables y(t)=1 dy(t)=0 ddy(t)=0
+        vars = @variables y(t) dy(t) ddy(t)
         @named src = SampledData(Float64)
         @named int = Integrator()
         @named iosys = ODESystem(
@@ -443,7 +443,7 @@ end
             tofloat = false)
         # prob = remake(prob; p = Parameter.(prob.p)) #<-- no longer needed with ModelingToolkit.jl PR #2231
 
-        sol = solve(prob, Rodas4(); initializealg = NoInit())
+        sol = solve(prob, Rodas4())
         @test sol.retcode == Success
         @test sol[src.output.u][1] == 1.0 #check correct initial condition
 
@@ -454,7 +454,7 @@ end
     end
 
     @testset "using Vector Based" begin
-        vars = @variables y(t)=1 dy(t)=0 ddy(t)=0
+        vars = @variables y(t) dy(t) ddy(t)
         @named src = SampledData(dt)
         @named int = Integrator()
         @named iosys = ODESystem(
@@ -472,7 +472,7 @@ end
             [s.src.buffer => x, s.src.sample_time => dt];
             tofloat = false)
 
-        sol = solve(prob, Rodas4(); initializealg = NoInit())
+        sol = solve(prob, Rodas4())
         @test sol.retcode == Success
         @test sol[src.output.u][1] == 1.0 #check correct initial condition