Feature: regionCoupling with AMG solvers

src/OpenFOAM/Make/files

@@ -276,6 +276,7 @@ $(GAMGInterfaces)/processorGAMGInterface/processorGAMGInterface.C
 $(GAMGInterfaces)/cyclicGAMGInterface/cyclicGAMGInterface.C
 $(GAMGInterfaces)/ggiGAMGInterface/ggiGAMGInterface.C
 $(GAMGInterfaces)/cyclicGgiGAMGInterface/cyclicGgiGAMGInterface.C
+$(GAMGInterfaces)/regionCoupleGAMGInterface/regionCoupleGAMGInterface.C
 /*$(GAMGInterfaces)/mixingPlaneGAMGInterface/mixingPlaneGAMGInterface.C*/
 
 GAMGInterfaceFields = $(GAMG)/interfaceFields
@@ -285,6 +286,7 @@ $(GAMGInterfaceFields)/processorGAMGInterfaceField/processorGAMGInterfaceField.C
 $(GAMGInterfaceFields)/cyclicGAMGInterfaceField/cyclicGAMGInterfaceField.C
 $(GAMGInterfaceFields)/ggiGAMGInterfaceField/ggiGAMGInterfaceField.C
 $(GAMGInterfaceFields)/cyclicGgiGAMGInterfaceField/cyclicGgiGAMGInterfaceField.C
+$(GAMGInterfaceFields)/regionCoupleGAMGInterfaceField/regionCoupleGAMGInterfaceField.C
 /*$(GAMGInterfaceFields)/mixingPlaneGAMGInterfaceField/mixingPlaneGAMGInterfaceField.C*/
 
 GAMGAgglomerations = $(GAMG)/GAMGAgglomerations

src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterface/lduInterface.H

@@ -84,6 +84,9 @@ public:
 
         // Access
 
+            //- Return true if interface is coupled
+            virtual bool coupled() const = 0;
+
             //- Return faceCell addressing
             virtual const unallocLabelList& faceCells() const = 0;
 

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomerateLduAddressing.C

@@ -217,11 +217,14 @@ void Foam::GAMGAgglomeration::agglomerateLduAddressing
     {
         if (fineInterfaces.set(inti))
         {
-            fineInterfaces[inti].initInternalFieldTransfer
+            if (fineInterfaces[inti].coupled())
-            (
+            {
-                Pstream::blocking,
+                fineInterfaces[inti].initInternalFieldTransfer
-                restrictMap
+                (
-            );
+                    Pstream::blocking,
+                    restrictMap
+                );
+            }
         }
     }
 
@@ -233,18 +236,21 @@ void Foam::GAMGAgglomeration::agglomerateLduAddressing
     {
         if (fineInterfaces.set(inti))
         {
-            fineInterfaceAddr.set
+            if (fineInterfaces[inti].coupled())
-            (
+            {
-                inti,
+                fineInterfaceAddr.set
-                new labelField
                 (
-                    fineInterfaces[inti].internalFieldTransfer
+                    inti,
+                    new labelField
                     (
-                        Pstream::blocking,
+                        fineInterfaces[inti].internalFieldTransfer
-                        restrictMap
+                        (
+                            Pstream::blocking,
+                            restrictMap
+                        )
                     )
-                )
+                );
-            );
+            }
         }
     }
 
@@ -279,19 +285,25 @@ void Foam::GAMGAgglomeration::agglomerateLduAddressing
     {
         if (fineInterfaces.set(inti))
         {
-            coarseInterfaces.set
+            if (fineInterfaces[inti].coupled())
-            (
+            {
-                inti,
+                coarseInterfaces.set
-                GAMGInterface::New
                 (
-                    meshLevels_[fineLevelIndex],
+                    inti,
-                    fineInterfaces[inti],
+                    GAMGInterface::New
-                    fineInterfaces[inti].interfaceInternalField(restrictMap),
+                    (
-                    fineInterfaceAddr[inti]
+                        meshLevels_[fineLevelIndex],
-                ).ptr()
+                        fineInterfaces[inti],
-            );
+                        fineInterfaces[inti].interfaceInternalField
+                        (
+                            restrictMap
+                        ),
+                        fineInterfaceAddr[inti]
+                    ).ptr()
+                );
 
-            coarseInterfaceAddr[inti] = coarseInterfaces[inti].faceCells();
+                coarseInterfaceAddr[inti] = coarseInterfaces[inti].faceCells();
+            }
         }
     }
 

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaceFields/regionCoupleGAMGInterfaceField/regionCoupleGAMGInterfaceField.C

@@ -0,0 +1,121 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright held by original author
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 2 of the License, or (at your
+    option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM; if not, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+\*---------------------------------------------------------------------------*/
+
+#include "regionCoupleGAMGInterfaceField.H"
+#include "addToRunTimeSelectionTable.H"
+#include "lduMatrix.H"
+
+#include "OSspecific.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    defineTypeNameAndDebug(regionCoupleGAMGInterfaceField, 0);
+    addToRunTimeSelectionTable
+    (
+        GAMGInterfaceField,
+        regionCoupleGAMGInterfaceField,
+        lduInterface
+    );
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::regionCoupleGAMGInterfaceField::regionCoupleGAMGInterfaceField
+(
+    const GAMGInterface& GAMGCp,
+    const lduInterfaceField& fineInterface
+)
+:
+    ggiGAMGInterfaceField(GAMGCp, fineInterface),
+    regionCoupleInterface_(refCast<const regionCoupleGAMGInterface>(GAMGCp))
+{}
+
+
+// * * * * * * * * * * * * * * * * Desstructor * * * * * * * * * * * * * * * //
+
+Foam::regionCoupleGAMGInterfaceField::~regionCoupleGAMGInterfaceField()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+void Foam::regionCoupleGAMGInterfaceField::initInterfaceMatrixUpdate
+(
+    const scalarField& psiInternal,
+    scalarField& result,
+    const lduMatrix& m,
+    const scalarField& coeffs,
+    const direction cmpt,
+    const Pstream::commsTypes commsType
+) const
+{
+    // This must have a reduce in it.  HJ, 15/May/2009
+    if (regionCoupleInterface_.coupled())
+    {
+        ggiGAMGInterfaceField::initInterfaceMatrixUpdate
+        (
+            psiInternal,
+            result,
+            m,
+            coeffs,
+            cmpt,
+            commsType
+        );
+    }
+}
+
+
+void Foam::regionCoupleGAMGInterfaceField::updateInterfaceMatrix
+(
+    const scalarField& psiInternal,
+    scalarField& result,
+    const lduMatrix& m,
+    const scalarField& coeffs,
+    const direction cmpt,
+    const Pstream::commsTypes commsType
+) const
+{
+    // Get expanded data to zone size.  No global reduce allowed
+    // HJ, 15/May/2009
+    if (regionCoupleInterface_.coupled())
+    {
+        ggiGAMGInterfaceField::updateInterfaceMatrix
+        (
+            psiInternal,
+            result,
+            m,
+            coeffs,
+            cmpt,
+            commsType
+        );
+    }
+}
+
+
+// ************************************************************************* //

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaceFields/regionCoupleGAMGInterfaceField/regionCoupleGAMGInterfaceField.H

@@ -0,0 +1,120 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright held by original author
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 2 of the License, or (at your
+    option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM; if not, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+Class
+    Foam::regionCoupleGAMGInterfaceField
+
+Description
+    GAMG agglomerated regionCouple interface field.
+
+Author
+    Hrvoje Jasak, Wikki Ltd.  All rights reserved.
+
+SourceFiles
+    regionCoupleGAMGInterfaceField.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef regionCoupleGAMGInterfaceField_H
+#define regionCoupleGAMGInterfaceField_H
+
+#include "regionCoupleGAMGInterface.H"
+#include "ggiGAMGInterfaceField.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                  Class regionCoupleGAMGInterfaceField Declaration
+\*---------------------------------------------------------------------------*/
+
+class regionCoupleGAMGInterfaceField
+:
+    public ggiGAMGInterfaceField
+{
+    // Private data
+
+        //- Local reference cast into the regionCouple interface
+        const regionCoupleGAMGInterface& regionCoupleInterface_;
+
+public:
+
+    //- Runtime type information
+    TypeName("regionCouple");
+
+
+    // Constructors
+
+        //- Construct from GAMG interface and fine level interface field
+        regionCoupleGAMGInterfaceField
+        (
+            const GAMGInterface& GAMGCp,
+            const lduInterfaceField& fineInterfaceField
+        );
+
+
+    // Destructor
+
+        virtual ~regionCoupleGAMGInterfaceField();
+
+
+    // Member Functions
+
+        // Interface matrix update
+
+            //- Initialise neighbour matrix update
+            virtual void initInterfaceMatrixUpdate
+            (
+                const scalarField& psiInternal,
+                scalarField& result,
+                const lduMatrix& m,
+                const scalarField& coeffs,
+                const direction cmpt,
+                const Pstream::commsTypes commsType
+            ) const;
+
+            //- Update result field based on interface functionality
+            virtual void updateInterfaceMatrix
+            (
+                const scalarField& psiInternal,
+                scalarField& result,
+                const lduMatrix&,
+                const scalarField& coeffs,
+                const direction cmpt,
+                const Pstream::commsTypes commsType
+            ) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaces/cyclicGAMGInterface/cyclicGAMGInterface.H

@@ -95,6 +95,15 @@ public:
 
     // Member Functions
 
+        // Access
+
+            //- Return true if interface is coupled
+            virtual bool coupled() const
+            {
+                return true;
+            }
+
+
         // Interface transfer functions
 
             //- Transfer and return neighbour field

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaces/ggiGAMGInterface/ggiGAMGInterface.H

@@ -138,6 +138,15 @@ public:
 
     // Member Functions
 
+        // Access
+
+            //- Return true if interface is coupled
+            virtual bool coupled() const
+            {
+                return true;
+            }
+
+
         // Agglomeration
 
             //- Agglomerating the given fine-level coefficients and return

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaces/processorGAMGInterface/processorGAMGInterface.H

@@ -95,6 +95,15 @@ public:
 
     // Member Functions
 
+        // Access
+
+            //- Return true if interface is coupled
+            virtual bool coupled() const
+            {
+                return true;
+            }
+
+
         // Interface transfer functions
 
             //- Initialise interface data transfer

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaces/regionCoupleGAMGInterface/regionCoupleGAMGInterface.C

@@ -0,0 +1,71 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright held by original author
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 2 of the License, or (at your
+    option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM; if not, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+\*---------------------------------------------------------------------------*/
+
+#include "regionCoupleGAMGInterface.H"
+#include "addToRunTimeSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    defineTypeNameAndDebug(regionCoupleGAMGInterface, 0);
+    addToRunTimeSelectionTable
+    (
+        GAMGInterface,
+        regionCoupleGAMGInterface,
+        lduInterface
+    );
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::regionCoupleGAMGInterface::regionCoupleGAMGInterface
+(
+    const lduPrimitiveMesh& lduMesh,
+    const lduInterface& fineInterface,
+    const labelField& localRestrictAddressing,
+    const labelField& neighbourRestrictAddressing
+)
+:
+    ggiGAMGInterface
+    (
+        lduMesh,
+        fineInterface,
+        localRestrictAddressing,
+        neighbourRestrictAddressing
+    ),
+    coupled_(fineInterface.coupled())
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
+
+Foam::regionCoupleGAMGInterface::~regionCoupleGAMGInterface()
+{}
+
+
+// ************************************************************************* //

src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/interfaces/regionCoupleGAMGInterface/regionCoupleGAMGInterface.H

@@ -0,0 +1,108 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright held by original author
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 2 of the License, or (at your
+    option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM; if not, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+Class
+    Foam::regionCoupleGAMGInterface
+
+Description
+    GAMG agglomerated regionCouple interface.
+
+Author
+    Hrvoje Jasak, Wikki Ltd.  All rights reserved.
+
+SourceFiles
+    regionCoupleGAMGInterface.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef regionCoupleGAMGInterface_H
+#define regionCoupleGAMGInterface_H
+
+#include "ggiGAMGInterface.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                     Class regionCoupleGAMGInterface Declaration
+\*---------------------------------------------------------------------------*/
+
+class regionCoupleGAMGInterface
+:
+    public ggiGAMGInterface
+{
+
+    // Private data
+
+        //- Is the interface coupled?
+        const bool coupled_;
+
+
+public:
+
+    //- Runtime type information
+    TypeName("regionCouple");
+
+
+    // Constructors
+
+        //- Construct from fine level interface,
+        //  local and neighbour restrict addressing
+        regionCoupleGAMGInterface
+        (
+            const lduPrimitiveMesh& lduMesh,
+            const lduInterface& fineInterface,
+            const labelField& localRestrictAddressing,
+            const labelField& neighbourRestrictAddressing
+        );
+
+
+    // Destructor
+
+        virtual ~regionCoupleGAMGInterface();
+
+
+    // Member Functions
+
+        // Access
+
+            //- Return true if interface is coupled
+            virtual bool coupled() const
+            {Info<< "BANANA" << endl;
+                return coupled_;
+            }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/OpenFOAM/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H

@@ -156,7 +156,7 @@ public:
 
         // Access
 
-            //- Return true because this patch is coupled
+            //- Return true if patch is coupled
             virtual bool coupled() const
             {
                 return true;

src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C

@@ -85,8 +85,12 @@ void Foam::ggiPolyPatch::calcZoneAddressing() const
     // Check zone addressing
     if (zAddr.size() > 0 && min(zAddr) < 0)
     {
+        Info<< "myZone: " << myZone << nl
+            << "my start and size: " << start() << " and " << size() << nl
+            << "zAddr: " << zAddr << endl;
+
         FatalErrorIn("void ggiPolyPatch::calcZoneAddressing() const")
-            << "Problem with patch-to zone addressing: some patch faces "
+            << "Problem with patch-to-zone addressing: some patch faces "
             << "not found in interpolation zone"
             << abort(FatalError);
     }
@@ -439,6 +443,9 @@ void Foam::ggiPolyPatch::clearOut()
 {
     clearGeom();
 
+    shadowIndex_ = -1;
+    zoneIndex_ = -1;
+
     deleteDemandDrivenData(zoneAddressingPtr_);
     deleteDemandDrivenData(patchToPatchPtr_);
     deleteDemandDrivenData(localParallelPtr_);

src/finiteVolume/fvMesh/fvPatches/basic/coupled/coupledFvPatch.H

@@ -102,7 +102,7 @@ public:
 
         // Access
 
-            //- Return true because this patch is coupled
+            //- Return true if patch is coupled
             virtual bool coupled() const
             {
                 return coupledPolyPatch_.coupled();