From e3b53f2727f1f2e199820b9a7a95bc181ecf0c71 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Wed, 16 Aug 2023 16:03:32 +0100
Subject: [PATCH] multiphaseEuler: Fix indexing bug in partial elimination

---
 .../MomentumTransferPhaseSystem.C                           | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
index a732ec4be3..8450bd912b 100644
--- a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
+++ b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
@@ -1442,8 +1442,8 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialElimination
         {
             for (label j = i + 1; j < phases.size(); j++)
             {
-                KdByAs[i][j] /= KdByAs[i][i];
-                KdByAfs[i][j] /= KdByAfs[i][i];
+                KdByAs[j][i] /= KdByAs[i][i];
+                KdByAfs[j][i] /= KdByAfs[i][i];
                 for (label k = i + 1; k < phases.size(); ++ k)
                 {
                     KdByAs[j][k] -= KdByAs[j][i]*KdByAs[i][k];
@@ -1584,7 +1584,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialEliminationf
         {
             for (label j = i + 1; j < phases.size(); j++)
             {
-                phiKdfs[i][j] /= phiKdfs[i][i];
+                phiKdfs[j][i] /= phiKdfs[i][i];
                 for (label k = i + 1; k < phases.size(); ++ k)
                 {
                     phiKdfs[j][k] -= phiKdfs[j][i]*phiKdfs[i][k];