add linearsolve

Author: ChrisRackauckas

Project.toml

@@ -43,6 +43,7 @@ julia = "1.6"
 
 [extras]
 FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
+LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 Parameters = "d96e819e-fc66-5662-9728-84c9c7592b0a"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
@@ -52,4 +53,4 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [targets]
-test = ["FillArrays", "OrdinaryDiffEq", "Parameters", "Random", "SafeTestsets", "SpecialFunctions", "Test", "Zygote"]
+test = ["FillArrays", "LinearSolve", "OrdinaryDiffEq", "Parameters", "Random", "SafeTestsets", "SpecialFunctions", "Test", "Zygote"]

test/Misc/jacvec_operators.jl

@@ -1,5 +1,5 @@
 using DiffEqBase,
-    DiffEqOperators, ForwardDiff, LinearAlgebra, SparseDiffTools, Test
+    DiffEqOperators, ForwardDiff, LinearAlgebra, SparseDiffTools, Test, LinearSolve
 const A = rand(300, 300)
 f(du, u) = mul!(du, A, u)
 f(u) = A * u