@@ -96,6 +96,7 @@ function xp_linsolve_rhs(
96
96
b_list
97
97
end
98
98
99
+ #=
99
100
function SciMLBase.solve(prob::DualAbstractLinearProblem, args...; kwargs...)
100
101
return solve(prob, nothing, args...; kwargs...)
101
102
end
@@ -109,6 +110,7 @@ function SciMLBase.solve(prob::DualAbstractLinearProblem,
109
110
alg::LinearSolve.SciMLLinearSolveAlgorithm, args...; kwargs...)
110
111
solve!(init(prob, alg, args...; kwargs...))
111
112
end
113
+ =#
112
114
113
115
function linearsolve_dual_solution (
114
116
u:: Number , partials, dual_type)
@@ -122,6 +124,7 @@ function linearsolve_dual_solution(
122
124
zip (u, partials_list[i, :] for i in 1 : length (partials_list[1 ])))
123
125
end
124
126
127
+ #=
125
128
function SciMLBase.init(
126
129
prob::DualAbstractLinearProblem, alg::LinearSolve.SciMLLinearSolveAlgorithm,
127
130
args...;
@@ -170,6 +173,7 @@ function SciMLBase.solve!(cache::DualLinearCache, args...; kwargs...)
170
173
cache.alg, dual_sol, sol.resid, cache; sol.retcode, sol.iters, sol.stats
171
174
)
172
175
end
176
+ =#
173
177
174
178
# If setting A or b for DualLinearCache, put the Dual-stripped versions in the LinearCache
175
179
# Also "forwards" setproperty so that
0 commit comments