More fixes!

Author: amontoison

src/decompression.jl

@@ -724,7 +724,7 @@ end
 #     Bc::AbstractMatrix,
 #     result::StarSetBicoloringResult,
 # )
-#     (; ag, S, A_indices, compressed_indices, pos_Bc) = result
+#     (; S, A_indices, compressed_indices, pos_Bc) = result
 #     fill!(A, zero(eltype(A)))
 
 #     ind_Br = 0
@@ -754,7 +754,7 @@ function decompress!(
     Bc::AbstractMatrix,
     result::StarSetBicoloringResult,
 )
-    (; ag, A_indices, compressed_indices, pos_Bc) = result
+    (; A_indices, compressed_indices, pos_Bc) = result
     nzA = nonzeros(A)
     for k in 1:pos_Bc
         nzA[A_indices[k]] = Bc[compressed_indices[k]]
@@ -774,7 +774,6 @@ function decompress!(
     result::TreeSetBicoloringResult,
 ) where {R<:Real}
     (;
-        ag,
         symmetric_color,
         symmetric_to_row,
         symmetric_to_column,
@@ -783,7 +782,6 @@ function decompress!(
         nt,
         buffer,
     ) = result
-    (; S) = ag
 
     m, n = size(A)
     fill!(A, zero(R))
@@ -801,6 +799,7 @@ function decompress!(
 
         # Reset the buffer to zero for all vertices in the tree (except the root)
         for pos in first:last
+            (vertex, _) = reverse_bfs_orders[pos]
             buffer_right_type[vertex] = zero(R)
         end
         # Reset the buffer to zero for the root vertex
@@ -831,7 +830,6 @@ function decompress!(
     result::TreeSetBicoloringResult,
 ) where {R<:Real}
     (;
-        ag,
         symmetric_color,
         symmetric_to_column,
         symmetric_to_row,
@@ -841,7 +839,6 @@ function decompress!(
         A_indices,
         buffer,
     ) = result
-    (; S) = ag
 
     m, n = size(A)
     A_colptr = A.colptr
@@ -864,6 +861,7 @@ function decompress!(
 
         # Reset the buffer to zero for all vertices in the tree (except the root)
         for pos in first:last
+            (vertex, _) = reverse_bfs_orders[pos]
             buffer_right_type[vertex] = zero(R)
         end
         # Reset the buffer to zero for the root vertex

src/result.jl

@@ -754,6 +754,11 @@ function TreeSetBicoloringResult(
             (leaf, neighbor) = reverse_bfs_orders[pos]
             index += 1
 
+            i = leaf
+            j = neighbor
+            col_i = view(rv, nzrange(S, i))
+            col_j = view(rv, nzrange(S, j))
+
             #! format: off
             # S[i,j] is in the lower triangular part of S
             if in_triangle(i, j, :L)