Bug in magma_scuspmm

[Noran]

Hello @ all,
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;

Shouldn't there stand the following?

Code: Select all

C.num_rows = A.num_rows;
C.num_cols = B.num_cols;

Greetings Noran

[hartwig anzt]

Yes - great thanks for pointing that out!

Hartwig