test for correctness from problem expr builders

Author: ChrisRackauckas

test/sdesystem.jl

@@ -23,6 +23,11 @@ f = SDEFunction(de)
 prob = SDEProblem(SDEFunction(de),f.g,[1.0,0.0,0.0],(0.0,100.0),(10.0,26.0,2.33))
 sol = solve(prob,SRIW1())
 
+probexpr = SDEProblem(SDEFunction(de),f.g,[1.0,0.0,0.0],(0.0,100.0),(10.0,26.0,2.33))
+solexpr = solve(eval(probexpr),SRIW1())
+
+@test all(x->x==0,Array(sol-solexpr))
+
 noiseeqs_nd = [0.01*x 0.01*x*y 0.02*x*z
                σ      0.01*y   0.02*x*z
                ρ      β        0.01*z  ]

test/steadystatesystems.jl

@@ -14,4 +14,7 @@ for factor in [1e-1, 1e0, 1e10], u0_p in [(2.34,2.676),(22.34,1.632),(.3,15.676)
     ss_prob = SteadyStateProblem(de,u0,p)
     sol = solve(ss_prob,SSRootfind()).u[1]
     @test abs(sol^2 - factor*u0_p[2]) < 1e-8
+    ss_prob = SteadyStateProblemExpr(de,u0,p)
+    sol_expr = solve(eval(ss_prob),SSRootfind()).u[1]
+    @test all(x->x==0,Array(sol-sol_expr))
 end