GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-01 01:32:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix the embarrassingly incomplete fix to the embarrassing bug in blocked transpose.

Browse Source
This commit is contained in:
Rasmus Munk Larsen 2020-04-29 17:27:36 +00:00
parent a32923a439
commit 923ee9aba3
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/Core/Transpose.h
View File

@ -285,8 +285,8 @@ void BlockedInPlaceTranspose(MatrixType& m) {
internal::ptranspose(A); internal::ptranspose(A);
internal::ptranspose(B); internal::ptranspose(B);
for (Index i=0; i<PacketSize; ++i) { for (Index i=0; i<PacketSize; ++i) {
m.template writePacket<Alignment>(m.rowIndexByOuterInner(row_start + i, col_start), m.colIndexByOuterInner(row_start + i,col_start), A.packet[i]); m.template writePacket<Alignment>(m.rowIndexByOuterInner(row_start + i, col_start), m.colIndexByOuterInner(row_start + i,col_start), B.packet[i]);
m.template writePacket<Alignment>(m.rowIndexByOuterInner(col_start + i, row_start), m.colIndexByOuterInner(col_start + i,row_start), B.packet[i]); m.template writePacket<Alignment>(m.rowIndexByOuterInner(col_start + i, row_start), m.colIndexByOuterInner(col_start + i,row_start), A.packet[i]);
} }
} }
} }