remove special problems

Author: ChrisRackauckas

src/solve.jl

@@ -49,7 +49,11 @@ function init{uType,tType,isinplace,algType<:Union{AbstractRODEAlgorithm,Abstrac
     save_everystep = save_timeseries
   end
 
-  if prob.mass_matrix != I
+  if typeof(prob.f)<:Tuple
+    if min((mm != I for mm in prob.mass_matrix)...)
+      error("This solver is not able to use mass matrices.")
+    end
+  elseif prob.mass_matrix != I
     error("This solver is not able to use mass matrices.")
   end
 

test/iif_methods.jl

@@ -11,8 +11,8 @@ no_noise = (t,u) -> 0.0
 f1_no_noise = (t,u) -> μ
 (p::typeof(f1_no_noise))(::Type{Val{:analytic}},t,u0,W) = u0.*exp.(2μ*t)
 
-prob = SplitSDEProblem((f1,f2),σ,1/2,(0.0,1.0))
-no_noise_prob = SplitSDEProblem((f1_no_noise,f2),no_noise,1/2,(0.0,1.0))
+prob = SDEProblem((f1,f2),σ,1/2,(0.0,1.0))
+no_noise_prob = SDEProblem((f1_no_noise,f2),no_noise,1/2,(0.0,1.0))
 
 sol = solve(prob,IIF1M(),dt=1/10)
 
@@ -69,7 +69,7 @@ function (p::typeof(f1))(::Type{Val{:analytic}},t,u0,W)
 end
 f2 = (t,u,du) -> du .= μ .* u
 
-prob = SplitSDEProblem((f1,f2),σ,u0,(0.0,1.0),noise_rate_prototype=rand(2,2))
+prob = SDEProblem((f1,f2),σ,u0,(0.0,1.0),noise_rate_prototype=rand(2,2))
 
 f1_no_noise = (t,u,du) -> A
 f2 = (t,u,du) -> du .= μ .* u
@@ -80,7 +80,7 @@ function (p::typeof(f1_no_noise))(::Type{Val{:analytic}},t,u0,W)
  tmp = (A+1.01I)*t
  expm(tmp)*u0
 end
-prob_no_noise = SplitSDEProblem((f1_no_noise,f2),σ,u0,(0.0,1.0),noise_rate_prototype=rand(2,2))
+prob_no_noise = SDEProblem((f1_no_noise,f2),σ,u0,(0.0,1.0),noise_rate_prototype=rand(2,2))
 
 
 sol = solve(prob,IIF1M(),dt=1/10)

test/split_tests.jl

@@ -6,7 +6,7 @@ f2 = (t,u) -> (1.01)/2 * u
 σ = (t,u) -> 0.87u
 #(p::typeof(f))(::Type{Val{:analytic}},t,u0,W) = u0.*exp.(0.63155t+0.87W)
 
-prob = SplitSDEProblem((f1,f2),σ,1/2,(0.0,1.0))
+prob = SDEProblem((f1,f2),σ,1/2,(0.0,1.0))
 
 sol = solve(prob,SplitEM(),dt=1/10)
 
@@ -17,7 +17,7 @@ sol2 = solve(prob,EM(),dt=1/10)
 @test sol[:] ≈ sol2[:]
 
 u0 = rand(4)
-prob = SplitSDEProblem((f1,f2),σ,u0,(0.0,1.0))
+prob = SDEProblem((f1,f2),σ,u0,(0.0,1.0))
 
 sol = solve(prob,SplitEM(),dt=1/10)