FECCellToFaceStencil.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2011-2016 OpenFOAM Foundation
    Copyright (C) 2018-2020 OpenCFD Ltd.
-------------------------------------------------------------------------------
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 3 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, see <http://www.gnu.org/licenses/>.
 
\*---------------------------------------------------------------------------*/
 
#include "FECCellToFaceStencil.H"
#include "syncTools.H"
#include "emptyPolyPatch.H"
#include "dummyTransform.H"
 
// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
void Foam::FECCellToFaceStencil::calcEdgeBoundaryData
(
    const boolList& isValidBFace,
    const labelList& boundaryEdges,
    EdgeMap<labelList>& neiGlobal
) const
{
    neiGlobal.resize(2*boundaryEdges.size());
 
    labelHashSet edgeGlobals;
 
    for (const label edgeI : boundaryEdges)
    {
        neiGlobal.insert
        (
            mesh().edges()[edgeI],
            calcFaceCells
            (
                isValidBFace,
                mesh().edgeFaces(edgeI),
                edgeGlobals
            )
        );
    }
 
    syncTools::syncEdgeMap
    (
        mesh(),
        neiGlobal,
        ListOps::unionEqOp(),
        dummyTransform()
    );
}
 
 
void Foam::FECCellToFaceStencil::calcFaceStencil
(
    labelListList& faceStencil
) const
{
    const polyBoundaryMesh& patches = mesh().boundaryMesh();
    const label nBnd = mesh().nBoundaryFaces();
    const labelList& own = mesh().faceOwner();
    const labelList& nei = mesh().faceNeighbour();
 
 
 
    // Determine neighbouring global cell
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
    labelList neiGlobalCell(nBnd);
    forAll(patches, patchi)
    {
        const polyPatch& pp = patches[patchi];
 
        if (pp.coupled())
        {
            label facei = pp.start();
 
            forAll(pp, i)
            {
                neiGlobalCell[facei-mesh().nInternalFaces()] =
                    globalNumbering().toGlobal(own[facei]);
                facei++;
            }
        }
    }
    syncTools::swapBoundaryFaceList(mesh(), neiGlobalCell);
 
 
 
    // Determine on coupled edges the edge cells on the other side
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
    // Calculate edges on coupled patches
    labelList boundaryEdges
    (
        allCoupledFacesPatch()().meshEdges
        (
            mesh().edges(),
            mesh().pointEdges()
        )
    );
 
    // Mark boundary faces to be included in stencil (i.e. not coupled or empty)
    boolList isValidBFace;
    validBoundaryFaces(isValidBFace);
 
    // Swap edgeCells for coupled edges. Note: use EdgeMap for now since we've
    // got syncTools::syncEdgeMap for those. Should be replaced with Map and
    // syncTools functionality to handle those.
    EdgeMap<labelList> neiGlobal;
    calcEdgeBoundaryData
    (
        isValidBFace,
        boundaryEdges,
        neiGlobal
    );
 
 
 
    // Construct stencil in global numbering
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
    faceStencil.setSize(mesh().nFaces());
 
    // Do coupled edges first
 
    forAll(boundaryEdges, i)
    {
        label edgeI = boundaryEdges[i];
 
        const labelList& eGlobals = neiGlobal[mesh().edges()[edgeI]];
 
        // Distribute to all edgeFaces
        const labelList& eFaces = mesh().edgeFaces(edgeI);
 
        forAll(eFaces, j)
        {
            label facei = eFaces[j];
 
            // Insert eGlobals into faceStencil.
            merge(-1, -1, eGlobals, faceStencil[facei]);
        }
    }
    neiGlobal.clear();
 
 
    // Do remaining edges by looping over all faces
 
    // Work arrays
    DynamicList<label> fEdgesSet;
    DynamicList<label> eFacesSet;
    labelHashSet faceStencilSet;
 
    for (label facei = 0; facei < mesh().nInternalFaces(); facei++)
    {
        label globalOwn = globalNumbering().toGlobal(own[facei]);
        label globalNei = globalNumbering().toGlobal(nei[facei]);
 
        // Convert any existing faceStencil (from coupled edges) into
        // set and operate on this.
 
        faceStencilSet.clear();
 
        // Insert all but global owner and neighbour
        forAll(faceStencil[facei], i)
        {
            label globalI = faceStencil[facei][i];
            if (globalI != globalOwn && globalI != globalNei)
            {
                faceStencilSet.insert(globalI);
            }
        }
        faceStencil[facei].clear();
 
        // Collect all edge connected (internal) cells
        const labelList& fEdges = mesh().faceEdges(facei, fEdgesSet);
 
        forAll(fEdges, i)
        {
            label edgeI = fEdges[i];
 
            insertFaceCells
            (
                globalOwn,
                globalNei,
                isValidBFace,
                mesh().edgeFaces(edgeI, eFacesSet),
                faceStencilSet
            );
        }
 
        // Extract, guarantee owner first, neighbour second.
        faceStencil[facei].setSize(faceStencilSet.size()+2);
        label n = 0;
        faceStencil[facei][n++] = globalOwn;
        faceStencil[facei][n++] = globalNei;
        for (const label stencili : faceStencilSet)
        {
            if (stencili == globalOwn || stencili == globalNei)
            {
                FatalErrorInFunction
                    << "problem:" << faceStencilSet
                    << abort(FatalError);
            }
            faceStencil[facei][n++] = stencili;
        }
    }
    forAll(patches, patchi)
    {
        const polyPatch& pp = patches[patchi];
        label facei = pp.start();
 
        if (pp.coupled())
        {
            forAll(pp, i)
            {
                label globalOwn = globalNumbering().toGlobal(own[facei]);
                label globalNei = neiGlobalCell[facei-mesh().nInternalFaces()];
 
                // Convert any existing faceStencil (from coupled edges) into
                // set and operate on this.
 
                faceStencilSet.clear();
 
                // Insert all but global owner and neighbour
                forAll(faceStencil[facei], i)
                {
                    label globalI = faceStencil[facei][i];
                    if (globalI != globalOwn && globalI != globalNei)
                    {
                        faceStencilSet.insert(globalI);
                    }
                }
                faceStencil[facei].clear();
 
                // Collect all edge connected (internal) cells
                const labelList& fEdges = mesh().faceEdges(facei, fEdgesSet);
 
                forAll(fEdges, i)
                {
                    label edgeI = fEdges[i];
 
                    insertFaceCells
                    (
                        globalOwn,
                        globalNei,
                        isValidBFace,
                        mesh().edgeFaces(edgeI, eFacesSet),
                        faceStencilSet
                    );
                }
 
                // Extract, guarantee owner first, neighbour second.
                faceStencil[facei].setSize(faceStencilSet.size()+2);
                label n = 0;
                faceStencil[facei][n++] = globalOwn;
                faceStencil[facei][n++] = globalNei;
                for (const label stencili : faceStencilSet)
                {
                    if (stencili == globalOwn || stencili == globalNei)
                    {
                        FatalErrorInFunction
                            << "problem:" << faceStencilSet
                            << abort(FatalError);
                    }
                    faceStencil[facei][n++] = stencili;
                }
 
                if (n != faceStencil[facei].size())
                {
                    FatalErrorInFunction
                        << " size:" << faceStencil[facei].size()
                        << abort(FatalError);
                }
 
                facei++;
            }
        }
        else if (!isA<emptyPolyPatch>(pp))
        {
            forAll(pp, i)
            {
                label globalOwn = globalNumbering().toGlobal(own[facei]);
 
                // Convert any existing faceStencil (from coupled edges) into
                // set and operate on this.
 
                faceStencilSet.clear();
 
                // Insert all but global owner and neighbour
                forAll(faceStencil[facei], i)
                {
                    label globalI = faceStencil[facei][i];
                    if (globalI != globalOwn)
                    {
                        faceStencilSet.insert(globalI);
                    }
                }
                faceStencil[facei].clear();
 
                // Collect all edge connected (internal) cells
                const labelList& fEdges = mesh().faceEdges(facei, fEdgesSet);
 
                forAll(fEdges, i)
                {
                    label edgeI = fEdges[i];
 
                    insertFaceCells
                    (
                        globalOwn,
                        -1,
                        isValidBFace,
                        mesh().edgeFaces(edgeI, eFacesSet),
                        faceStencilSet
                    );
                }
 
                // Extract, guarantee owner first, neighbour second.
                faceStencil[facei].setSize(faceStencilSet.size()+1);
                label n = 0;
                faceStencil[facei][n++] = globalOwn;
                for (const label stencili : faceStencilSet)
                {
                    if (stencili == globalOwn)
                    {
                        FatalErrorInFunction
                            << "problem:" << faceStencilSet
                            << abort(FatalError);
                    }
                    faceStencil[facei][n++] = stencili;
                }
 
                facei++;
            }
        }
    }
 
 
    for (label facei = 0; facei < mesh().nInternalFaces(); facei++)
    {
        label globalOwn = globalNumbering().toGlobal(own[facei]);
        if (faceStencil[facei][0] != globalOwn)
        {
            FatalErrorInFunction
                << "problem:" << faceStencil[facei]
                << " globalOwn:" << globalOwn
                << abort(FatalError);
        }
        label globalNei = globalNumbering().toGlobal(nei[facei]);
        if (faceStencil[facei][1] != globalNei)
        {
            FatalErrorInFunction
                << "problem:" << faceStencil[facei]
                << " globalNei:" << globalNei
                << abort(FatalError);
        }
    }
 
 
    forAll(patches, patchi)
    {
        const polyPatch& pp = patches[patchi];
 
        if (pp.coupled())
        {
            forAll(pp, i)
            {
                label facei = pp.start()+i;
 
                label globalOwn = globalNumbering().toGlobal(own[facei]);
                if (faceStencil[facei][0] != globalOwn)
                {
                    FatalErrorInFunction
                        << "problem:" << faceStencil[facei]
                        << " globalOwn:" << globalOwn
                        << abort(FatalError);
                }
                label globalNei = neiGlobalCell[facei-mesh().nInternalFaces()];
                if (faceStencil[facei][1] != globalNei)
                {
                    FatalErrorInFunction
                        << "problem:" << faceStencil[facei]
                        << " globalNei:" << globalNei
                        << abort(FatalError);
                }
            }
        }
        else if (!isA<emptyPolyPatch>(pp))
        {
            forAll(pp, i)
            {
                label facei = pp.start()+i;
 
                label globalOwn = globalNumbering().toGlobal(own[facei]);
                if (faceStencil[facei][0] != globalOwn)
                {
                    FatalErrorInFunction
                        << "problem:" << faceStencil[facei]
                        << " globalOwn:" << globalOwn
                        << abort(FatalError);
                }
            }
        }
    }
}
 
 
// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
Foam::FECCellToFaceStencil::FECCellToFaceStencil(const polyMesh& mesh)
:
    cellToFaceStencil(mesh)
{
    // Calculate per face the (edge) connected cells (in global numbering)
    labelListList faceStencil;
    calcFaceStencil(faceStencil);
 
    // Transfer to *this
    transfer(faceStencil);
}
 
 
// ************************************************************************* //