From cf492d42a11df04db713ee25f4446d3fc63a10c4 Mon Sep 17 00:00:00 2001 From: Kutalmis Bercin Date: Mon, 31 Jan 2022 12:00:06 +0000 Subject: [PATCH] ENH: Matrix: improve row access in inner/outer products --- src/OpenFOAM/matrices/Matrix/Matrix.C | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/OpenFOAM/matrices/Matrix/Matrix.C b/src/OpenFOAM/matrices/Matrix/Matrix.C index fc262b707b..0ba6168888 100644 --- a/src/OpenFOAM/matrices/Matrix/Matrix.C +++ b/src/OpenFOAM/matrices/Matrix/Matrix.C @@ -1020,9 +1020,9 @@ operator& 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) { @@ -1065,9 +1065,9 @@ operator^ 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)); }