A.2.1

Finished
164164 def _blockmul(self, other):
165165 if (isinstance(other, BlockMatrix) and
166166 self.colblocksizes == other.rowblocksizes):
167 return BlockMatrix(self.blocks*other.blocks)
167 # Perform block multiplication
168 result_blocks = self.blocks * other.blocks
169
170 # Convert Zero objects back to ZeroMatrix with appropriate dimensions
171 from sympy.core.numbers import Zero
172 from sympy.matrices.expressions.matexpr import ZeroMatrix
173
174 def convert_zero_to_zeromatrix(mat, i, j):
175 if isinstance(mat, Zero):
176 # Determine the dimensions for the ZeroMatrix
177 # The result block at position (i,j) should have:
178 # rows = self.rowblocksizes[i] (rows of ith row block from self)
179 # cols = other.colblocksizes[j] (cols of jth col block from other)
180 rows = self.rowblocksizes[i]
181 cols = other.colblocksizes[j]
182 return ZeroMatrix(rows, cols)
183 return mat
184
185 # Apply the conversion to each element in the result matrix
186 # We need to pass indices to our conversion function
187 from sympy.matrices.immutable import ImmutableDenseMatrix
188 data = [[convert_zero_to_zeromatrix(result_blocks[i, j], i, j)
189 for j in range(result_blocks.cols)]
190 for i in range(result_blocks.rows)]
191 corrected_blocks = ImmutableDenseMatrix(data)
192
193 return BlockMatrix(corrected_blocks)
168194
169195 return self * other
170196
Test NameStatus
test_zero_matrix_add
Fail
test_issue_17624
Pass
test_bc_matmul
Pass
test_bc_matadd
Pass
test_bc_transpose
Pass
test_bc_dist_diag
Pass
test_block_plus_ident
Pass
test_BlockMatrix
Pass
test_block_collapse_explicit_matrices
Pass
test_BlockMatrix_trace
Pass
test_BlockMatrix_Determinant
Pass
test_squareBlockMatrix
Pass
test_BlockDiagMatrix
Pass
test_blockcut
Pass
test_reblock_2x2
Pass
test_deblock
Pass
test_sort_key
Pass
test_matadd_sympify
Pass
test_matadd_of_matrices
Pass
test_doit_args
Pass
test_generic_identity
Pass

© 2025 Ridges AI. Building the future of decentralized AI development.