From 618cc3da6255c8df25c1cd07a21e0aa356fad1de Mon Sep 17 00:00:00 2001 From: Alexis Montoison Date: Tue, 30 Apr 2024 02:32:06 -0400 Subject: [PATCH] Clean adjoint_forward_substitution! --- src/application.jl | 23 ----------------------- 1 file changed, 23 deletions(-) diff --git a/src/application.jl b/src/application.jl index 6b20cfa..017e890 100644 --- a/src/application.jl +++ b/src/application.jl @@ -181,34 +181,11 @@ end function adjoint_forward_substitution!(F::ILUFactorization, y::AbstractVector) ldiv!(LowerTriangular(F.U'), y) - # U = F.U - # @inbounds for col = 1 : U.n - # y[col] /= U.nzval[U.colptr[col]] - # end - # @inbounds for col = 1 : U.n - # for idx = U.colptr[col] + 1 : U.colptr[col + 1] - 1 - # y[U.rowval[idx]] -= (U.nzval[idx] / U.nzval[U.colptr[U.rowval[idx]]]) * y[col] - # end - # end - y end function adjoint_forward_substitution!(F::ILUFactorization, y::AbstractMatrix) ldiv!(LowerTriangular(F.U'), y) - # U = F.U - # p = size(y, 2) - # @inbounds for c = 1 : p - # @inbounds for col = 1 : U.n - # y[col,c] /= U.nzval[U.colptr[col]] - # end - # @inbounds for col = 1 : U.n - # for idx = U.colptr[col] + 1 : U.colptr[col + 1] - 1 - # y[U.rowval[idx],c] -= (U.nzval[idx] / U.nzval[U.colptr[U.rowval[idx]]]) * y[col,c] - # end - # end - # end - y end