bug fixing

Author: anastasia21112

benchmarks/LinearSolve/MatrixDepot.jmd

@@ -51,16 +51,24 @@ for z in 1:length(allmatrices_md.content[1].rows)
         matrix_size[z] = n
 
         mtx_copy = copy(A)
-    
+
+        @info "$n × $n"
+        n > 500 && error("Skipping too large matrices")
+
+        rows, cols = size(mtx_copy)
+        new_rows = div(rows, 2)
+        new_cols = div(cols, 2)
+        condensed = zeros(Int, new_rows, new_cols) 
+
         while size(condensed, 1) > 32 || size(condensed, 2) > 32
             rows, cols = size(mtx_copy)
             new_rows = div(rows, 2)
             new_cols = div(cols, 2)
-            condensed = zeros(Int, new_rows, new_cols) 
+            condensed = sparse(zeros(Int, new_rows, new_cols))
 
             for r in 1:2:rows-1
                 for c in 1:2:cols-1
-                    block = mtx_copy[r:r+1, c:c+1]
+                    block = sparse(mtx_copy[r:r+1, c:c+1])
                     condensed[div(r-1, 2) + 1, div(c-1, 2) + 1] = (length(nonzeros(block)) >= 3) ? 1 : 0
                 end
             end 
@@ -70,11 +78,10 @@ for z in 1:length(allmatrices_md.content[1].rows)
 
 
         percentage_sparsity[z] = length(nonzeros(A)) / n^2
-        @info "$n × $n"
-
+        
         spaced_out_sparsity[z] = length(nonzeros(mtx_copy)) * percentage_sparsity[z]
 
-        n > 500 && error("Skipping too large matrices")
+        
 
         b = rand(rng, n)
         u0 = rand(rng, n)
@@ -99,6 +106,8 @@ for z in 1:length(allmatrices_md.content[1].rows)
         println(e)
     end
 end
+
+print(percentage_sparsity)
 ```
 
 ```julia