Merge be73e8f into 1bb2a3a

SciML · Apr 18, 2022 · 7313c17 · 7313c17
2 parents 1bb2a3a + be73e8f
commit 7313c17
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/common.jl b/src/common.jl
@@ -8,6 +8,7 @@ struct LinearCache{TA,Tb,Tu,Tp,Talg,Tc,Tl,Tr,Ttol}
     isfresh::Bool # false => cacheval is set wrt A, true => update cacheval wrt A
     Pl::Tl        # preconditioners
     Pr::Tr
+    solve_adjoint::Bool # solve adjoint problem
     abstol::Ttol
     reltol::Ttol
     maxiters::Int
@@ -77,6 +78,7 @@ default_tol(::Type{Any}) = 0
 
 function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorithm,Nothing}, args...;
                         alias_A = false, alias_b = false,
+                        solve_adjoint = false,
                         abstol=default_tol(eltype(prob.A)),
                         reltol=default_tol(eltype(prob.A)),
                         maxiters=length(prob.b),
@@ -94,6 +96,11 @@ function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorith
         fill!(u0,false)
     end
 
+    # TODO
+    if solve_adjoint === true
+        A = A'
+    end
+
     cacheval = init_cacheval(alg, A, b, u0, Pl, Pr, maxiters, abstol, reltol, verbose)
     isfresh = true
     Tc = typeof(cacheval)
@@ -121,6 +128,7 @@ function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorith
         isfresh,
         Pl,
         Pr,
+        solve_adjoint,
         abstol,
         reltol,
         maxiters,