Merge pull request #588 from SciML/deprecations

Removed deprecated functionality

Author: ChrisRackauckas

src/LinearSolve.jl

@@ -138,7 +138,15 @@ include("solve_function.jl")
 include("default.jl")
 include("init.jl")
 include("adjoint.jl")
-include("deprecated.jl")
+
+## Deprecated, remove in July 2025
+
+@static if isdefined(SciMLBase, :DiffEqArrayOperator)
+    function defaultalg(A::SciMLBase.DiffEqArrayOperator, b,
+            assump::OperatorAssumptions{Bool})
+        defaultalg(A.A, b, assump)
+    end
+end
 
 @inline function _notsuccessful(F::LinearAlgebra.QRCompactWY)
     (m, n) = size(F)

src/default.jl

@@ -45,7 +45,7 @@ end
 # For SciML algorithms already using `defaultalg`, all assume square matrix.
 defaultalg(A, b) = defaultalg(A, b, OperatorAssumptions(true))
 
-function defaultalg(A::Union{DiffEqArrayOperator, MatrixOperator}, b,
+function defaultalg(A::MatrixOperator, b,
         assump::OperatorAssumptions{Bool})
     defaultalg(A.A, b, assump)
 end

src/deprecated.jl

@@ -54,17 +54,17 @@ prob = LinearProblem(copy(A), copy(b1))
 
 linsolve = init(prob, UMFPACKFactorization())
 sol11 = solve!(linsolve)
-linsolve = LinearSolve.set_b(sol11.cache, copy(b2))
+linsolve.b = copy(b2)
 sol12 = solve!(linsolve)
-linsolve = LinearSolve.set_A(sol12.cache, copy(A2))
+linsolve.A = copy(A2)
 sol13 = solve!(linsolve)
 
 for alg in algs
     linsolve = init(prob, alg)
     sol31 = solve!(linsolve)
-    linsolve = LinearSolve.set_b(sol31.cache, copy(b2))
+    linsolve.b = copy(b2)
     sol32 = solve!(linsolve)
-    linsolve = LinearSolve.set_A(sol32.cache, copy(A2))
+    linsolve.A = copy(A2)
     sol33 = solve!(linsolve)
     @test sol11.u ≈ sol31.u
     @test sol12.u ≈ sol32.u