ENH: Matrix: improve row access in inner/outer products

This commit is contained in:
Kutalmis Bercin
2022-01-31 12:00:06 +00:00
committed by Mark Olesen
parent 98b4779793
commit cf492d42a1

View File

@ -1020,9 +1020,9 @@ operator&
Zero Zero
); );
for (label i = 0; i < AB.m(); ++i) for (label k = 0; k < B.m(); ++k)
{ {
for (label k = 0; k < B.m(); ++k) for (label i = 0; i < AB.m(); ++i)
{ {
for (label j = 0; j < AB.n(); ++j) for (label j = 0; j < AB.n(); ++j)
{ {
@ -1065,9 +1065,9 @@ operator^
for (label i = 0; i < AB.m(); ++i) for (label i = 0; i < AB.m(); ++i)
{ {
for (label k = 0; k < BT.n(); ++k) for (label j = 0; j < AB.n(); ++j)
{ {
for (label j = 0; j < AB.n(); ++j) for (label k = 0; k < BT.n(); ++k)
{ {
AB(i, j) += A(i, k)*Detail::conj(BT(j, k)); AB(i, j) += A(i, k)*Detail::conj(BT(j, k));
} }