diff --git a/src/OpenFOAM/matrices/lduMatrix/preconditioners/DILUPreconditioner/DILUPreconditioner.C b/src/OpenFOAM/matrices/lduMatrix/preconditioners/DILUPreconditioner/DILUPreconditioner.C
index 11f0da27c36267bc1fd9a4a0a0921412303f6c31..07cb690ce700c7c44d002a074a4f2e367abbf04c 100644
--- a/src/OpenFOAM/matrices/lduMatrix/preconditioners/DILUPreconditioner/DILUPreconditioner.C
+++ b/src/OpenFOAM/matrices/lduMatrix/preconditioners/DILUPreconditioner/DILUPreconditioner.C
@@ -107,8 +107,6 @@ void Foam::DILUPreconditioner::precondition
         solver_.matrix().lduAddr().upperAddr().begin();
     const label* const __restrict__ lPtr =
         solver_.matrix().lduAddr().lowerAddr().begin();
-    const label* const __restrict__ losortPtr =
-        solver_.matrix().lduAddr().losortAddr().begin();
 
     const scalar* const __restrict__ upperPtr =
         solver_.matrix().upper().begin();
@@ -126,9 +124,8 @@ void Foam::DILUPreconditioner::precondition
 
     for (label face=0; face<nFaces; face++)
     {
-        const label sface = losortPtr[face];
-        wAPtr[uPtr[sface]] -=
-            rDPtr[uPtr[sface]]*lowerPtr[sface]*wAPtr[lPtr[sface]];
+        wAPtr[uPtr[face]] -=
+            rDPtr[uPtr[face]]*lowerPtr[face]*wAPtr[lPtr[face]];
     }
 
     for (label face=nFacesM1; face>=0; face--)