i use the sparse matrix multiplication in Magma. Yesterday I implemented a case for that: A*B = C. I multiply 2 matrices A and B. A has 24 rows and columns. B has 24 rows and 3 columns. I wondered when i realized that C is 24 x 24. So i looked in the code from magma_scuspmm and found in line 68 and 69 the following:
Code: Select all
C.num_rows = A.num_rows;
C.num_cols = A.num_cols;
Code: Select all
C.num_rows = A.num_rows;
C.num_cols = B.num_cols;