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
Branches
Tags
No related merge requests found
......@@ -1011,7 +1011,7 @@ BSRMat<T, BlockStorageOrder> BSRMat<T, BlockStorageOrder>::apply_op(const char o
if(op == 'H')
block = DenseMatMap<T>(this->data+block_offset*this->bm*this->bn, this->bm, this->bn).adjoint().eval();
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;
i++;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment