Merge pull request #89 from SciML/myb/init

More robust inits

Author: YingboMa

test/Mechanical/rotational.jl

@@ -20,7 +20,7 @@ D = Differential(t)
 
     @named model = ODESystem(connections, t, systems=[fixed, inertia1, inertia2, spring, damper])
     sys = structural_simplify(model)
-    prob = DAEProblem(sys, D.(states(sys)) .=> 0.0, [], (0, 10.0))
+    prob = DAEProblem(sys, D.(states(sys)) .=> 0, D.(states(model)) .=> 0.0, (0, 10.0))
     sol = solve(prob, DFBDF())
 
     # Plots.plot(sol; vars=[inertia1.w, inertia2.w])
@@ -53,8 +53,7 @@ end
 
     @named model = ODESystem(connections, t, systems=[fixed, torque, inertia1, inertia2, spring, damper, sine])
     sys = structural_simplify(model)
-    #prob = ODAEProblem(sys, Pair[], (0, 1.0))
-    prob = DAEProblem(sys, D.(states(sys)) .=> 0.0, [D(D(inertia2.phi)) => 1.0], (0, 10.0))
+    prob = DAEProblem(sys, D.(states(sys)) .=> 0.0, [D(D(inertia2.phi)) => 1.0; D.(states(model)) .=> 0.0], (0, 10.0))
     sol = solve(prob, DFBDF())
 
     # Plots.plot(sol; vars=[inertia1.w, -inertia2.w*2])