From a8591f65bf9bec05a6e7e5ec0a6ea8bb0f3c0a52 Mon Sep 17 00:00:00 2001 From: nightcityblade Date: Thu, 15 Oct 2026 12:56:38 +0800 Subject: [PATCH] fix: remove dead scalar boundary code --- julia/MOLE.jl/src/BCs/scalarBC.jl | 5 ----- 1 file changed, 5 deletions(-) diff --git a/julia/MOLE.jl/src/BCs/scalarBC.jl b/julia/MOLE.jl/src/BCs/scalarBC.jl index 797ddd8d..03ecd710 100644 --- a/julia/MOLE.jl/src/BCs/scalarBC.jl +++ b/julia/MOLE.jl/src/BCs/scalarBC.jl @@ -269,9 +269,6 @@ function addScalarBC!(A::SparseMatrixCSC, b::AbstractVector, bc::ScalarBC2D{T}, # remove rows of A associated to boundary Abc1 = Abcl .+ Abcr rowsbc1, _, _ = findnz(Abc1) - # For some reason, the below code doesn't work: - # rows1, cols1, s1 = findnz(A[rowsbc1, :]) - # A = A .- sparse(rows1, cols1, s1, size(A, 1), size(A, 2)) A[rowsbc1, :] = spzeros(size(A[rowsbc1, :])) # This does though # Update matrix A with boundary information A = A .+ Abc1 @@ -290,8 +287,6 @@ function addScalarBC!(A::SparseMatrixCSC, b::AbstractVector, bc::ScalarBC2D{T}, # Remove rows of A associated to boundary Abc2 = Abct .+ Abcb rowsbc2, _, _ = findnz(Abc2) - # rows2, cols2, s2 = findnz(A[rowsbc2, :]) - # A = A .- sparse(rows2, cols2, s2, size(A, 1), size(A, 2)) A[rowsbc2, :] = spzeros(size(A[rowsbc2, :])) # Update matriz A with boundary information A = A .+ Abc2