Mentions légales du service

Skip to content
Snippets Groups Projects
Commit f5a7f449 authored by hhakim's avatar hhakim
Browse files

Fix MatBSR<FPP, Cpu>::transpose for complex case (it was adjoint instead of transpose).

parent 2ad805e4
No related branches found
No related tags found
No related merge requests found
...@@ -1011,7 +1011,7 @@ BSRMat<T, BlockStorageOrder> BSRMat<T, BlockStorageOrder>::apply_op(const char o ...@@ -1011,7 +1011,7 @@ BSRMat<T, BlockStorageOrder> BSRMat<T, BlockStorageOrder>::apply_op(const char o
if(op == 'H') if(op == 'H')
block = DenseMatMap<T>(this->data+block_offset*this->bm*this->bn, this->bm, this->bn).adjoint().eval(); block = DenseMatMap<T>(this->data+block_offset*this->bm*this->bn, this->bm, this->bn).adjoint().eval();
else if(op == 'T') else if(op == 'T')
block = DenseMatMap<T>(this->data+block_offset*this->bm*this->bn, this->bm, this->bn).adjoint().eval(); block = DenseMatMap<T>(this->data+block_offset*this->bm*this->bn, this->bm, this->bn).transpose().eval();
bcolinds[t_block_offset++] = mat_row_id/this->bm; bcolinds[t_block_offset++] = mat_row_id/this->bm;
i++; i++;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment