@@ -728,15 +728,12 @@ struct JoinCompressed{R<:Real,M1<:AbstractMatrix{R},M2<:AbstractMatrix{R}} <:
728
728
AbstractMatrix{R}
729
729
m:: Int
730
730
n:: Int
731
- c:: Int
732
731
Br:: M1
733
732
Bc:: M2
734
733
row_color_ind:: Dict{Int,Int}
735
734
col_color_ind:: Dict{Int,Int}
736
735
end
737
736
738
- Base. size (B:: JoinCompressed ) = (B. m + B. n, B. c)
739
-
740
737
function Base. getindex (B:: JoinCompressed , i:: Integer , j:: Integer )
741
738
(; n, Br, Bc, row_color_ind, col_color_ind) = B
742
739
if i ≤ n
@@ -758,8 +755,7 @@ function decompress!(
758
755
)
759
756
(; row_color_ind, col_color_ind, symmetric_result) = result
760
757
m, n = size (A)
761
- c = ncolors (result)
762
- Br_and_Bc = JoinCompressed (m, n, c, Br, Bc, row_color_ind, col_color_ind)
758
+ Br_and_Bc = JoinCompressed (m, n, Br, Bc, row_color_ind, col_color_ind)
763
759
A_and_Aᵀ = decompress (Br_and_Bc, symmetric_result)
764
760
copyto! (A, A_and_Aᵀ[(n + 1 ): (n + m), 1 : n]) # original matrix in bottom left corner
765
761
return A
@@ -770,8 +766,7 @@ function decompress!(
770
766
)
771
767
(; row_color_ind, col_color_ind, symmetric_result, large_colptr, large_rowval) = result
772
768
m, n = size (A)
773
- c = ncolors (result)
774
- Br_and_Bc = JoinCompressed (m, n, c, Br, Bc, row_color_ind, col_color_ind)
769
+ Br_and_Bc = JoinCompressed (m, n, Br, Bc, row_color_ind, col_color_ind)
775
770
# pretend A is larger
776
771
A_and_noAᵀ = SparseMatrixCSC (m + n, m + n, large_colptr, large_rowval, A. nzval)
777
772
# decompress lower triangle only
0 commit comments