Page 1 of 1

Bug in magma_scuspmm

Posted: Thu Nov 06, 2014 3:54 am
by 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

Re: Bug in magma_scuspmm

Posted: Thu Nov 06, 2014 6:30 am
by hartwig anzt
Yes - great thanks for pointing that out!

Hartwig