@@ -5,34 +5,6 @@ using LinearSolve.LinearAlgebra, LinearSolve.ArrayInterface, RecursiveFactorizat
5
5
6
6
LinearSolve. userecursivefactorization (A:: Union{Nothing,AbstractMatrix} ) = true
7
7
8
- function LinearSolve. init_cacheval (alg:: RFLUFactorization , A, b, u, Pl, Pr, maxiters:: Int ,
9
- abstol, reltol, verbose:: Bool , assumptions:: OperatorAssumptions )
10
- ipiv = Vector {LinearAlgebra.BlasInt} (undef, min (size (A)... ))
11
- ArrayInterface. lu_instance (convert (AbstractMatrix, A)), ipiv
12
- end
13
-
14
- function LinearSolve. init_cacheval (
15
- alg:: RFLUFactorization , A:: Matrix{Float64} , b, u, Pl, Pr,
16
- maxiters:: Int ,
17
- abstol, reltol, verbose:: Bool , assumptions:: OperatorAssumptions )
18
- ipiv = Vector {LinearAlgebra.BlasInt} (undef, 0 )
19
- LinearSolve. PREALLOCATED_LU, ipiv
20
- end
21
-
22
- function LinearSolve. init_cacheval (alg:: RFLUFactorization ,
23
- A:: Union{LinearSolve.SparseArrays.AbstractSparseArray, LinearSolve.SciMLOperators.AbstractSciMLOperator} , b, u, Pl, Pr,
24
- maxiters:: Int ,
25
- abstol, reltol, verbose:: Bool , assumptions:: OperatorAssumptions )
26
- nothing , nothing
27
- end
28
-
29
- function LinearSolve. init_cacheval (alg:: RFLUFactorization ,
30
- A:: Union{Diagonal, SymTridiagonal, Tridiagonal} , b, u, Pl, Pr,
31
- maxiters:: Int ,
32
- abstol, reltol, verbose:: Bool , assumptions:: OperatorAssumptions )
33
- nothing , nothing
34
- end
35
-
36
8
function SciMLBase. solve! (cache:: LinearSolve.LinearCache , alg:: RFLUFactorization{P, T} ;
37
9
kwargs... ) where {P, T}
38
10
A = cache. A
0 commit comments