Tiny memory save

JuliaLinearAlgebra · Oct 5, 2017 · 3a914bf · 3a914bf
1 parent b3e9bdb
commit 3a914bf
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 5 deletions.
diff --git a/src/classical.jl b/src/classical.jl
@@ -42,12 +42,11 @@ function direct_interpolation(A, T, splitting)
     fill!(T.nzval, 1.)
     T = A .* T
     Pp = rs_direct_interpolation_pass1(T, A, splitting)
-    Pp = Pp .+ 1
+    Pp .= Pp .+ 1
 
     Px, Pj, Pp = rs_direct_interpolation_pass2(A, T, splitting, Pp)
 
-    # Px .= abs.(Px)
-    Pj = Pj .+ 1
+    Pj .= Pj .+ 1
 
     R = SparseMatrixCSC(maximum(Pj), size(A, 1), Pp, Pj, Px)
     P = R'

diff --git a/src/splitting.jl b/src/splitting.jl
@@ -116,8 +116,6 @@ function RS_CF_splitting(S::SparseMatrixCSC, T::SparseMatrixCSC)
 
 					lambda[row] == 0 && continue
 
-					# assert(lambda[j] > 0);//this would cause problems!
-
 					# move j to the beginning of its current interval
 					lambda_j = lambda[row] + 1
 					old_pos  = node_to_index[row]