blockDescriptorEdges.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) 2019-2021 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 "blockDescriptor.H"
#include "lineEdge.H"
#include "lineDivide.H"
#include "hexCell.H"

// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //

int Foam::blockDescriptor::calcEdgePointsWeights
(
    pointField& edgePoints,
    scalarList& edgeWeights,
    const Foam::edge& cellModelEdge,
    const label nDiv,
    const gradingDescriptors& expand
) const
{
    // The topological edge on the block
    const Foam::edge thisEdge(blockShape_, cellModelEdge);

    if (blockEdge::debug && !thisEdge.good())
    {
        Info<< "Collapsed edge:" << thisEdge;
        if (index_ >= 0)
        {
            Info << " block:" << index_;
        }
        Info<< " model edge:" << cellModelEdge << nl;
    }

    // FUTURE: skip point generation for collapsed edge


    // Set the edge points/weights
    // The edge is a straight-line if it is not in the list of blockEdges

    for (const blockEdge& cedge : blockEdges_)
    {
        const int cmp = cedge.compare(thisEdge);

        if (cmp > 0)
        {
            // Curve has the same orientation

            // Divide the line
            const lineDivide divEdge(cedge, nDiv, expand);

            edgePoints  = divEdge.points();
            edgeWeights = divEdge.lambdaDivisions();

            return 1;  // Found curved-edge: done
        }
        else if (cmp < 0)
        {
            // Curve has the opposite orientation

            // Divide the line
            const lineDivide divEdge(cedge, nDiv, expand.inv());

            const pointField& p = divEdge.points();
            const scalarList& d = divEdge.lambdaDivisions();

            edgePoints.resize(p.size());
            edgeWeights.resize(d.size());

            // Copy in reverse order
            const label pn = (p.size() - 1);
            forAll(p, pi)
            {
                edgePoints[pi]  = p[pn - pi];
                edgeWeights[pi] = 1 - d[pn - pi];
            }

            return 1;  // Found curved-edge: done
        }
    }

    // Not curved-edge: divide the edge as a straight line

    // Get list of points for this block
    const pointField blockPoints(blockShape_.points(vertices_));

    lineDivide divEdge
    (
        blockEdges::lineEdge(blockPoints, cellModelEdge),
        nDiv,
        expand
    );

    edgePoints = divEdge.points();
    edgeWeights = divEdge.lambdaDivisions();

    return 0;
}


// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

int Foam::blockDescriptor::edgesPointsWeights
(
    pointField (&edgesPoints)[12],
    scalarList (&edgesWeights)[12]
) const
{
    int nCurved = 0;

    for (label edgei = 0; edgei < 12; ++edgei)  //< hexCell::nEdges()
    {
        nCurved += calcEdgePointsWeights
        (
            edgesPoints[edgei],
            edgesWeights[edgei],
            hexCell::modelEdges()[edgei],

            sizes()[edgei/4],   // 12 edges -> 3 components (x,y,z)
            expand_[edgei]
        );
    }

    return nCurved;
}


bool Foam::blockDescriptor::edgePointsWeights
(
    const label edgei,
    pointField& edgePoints,
    scalarList& edgeWeights,
    const label nDiv,
    const gradingDescriptors& gd
) const
{
    if (edgei < 0 || edgei >= 12)  //< hexCell::nEdges()
    {
        FatalErrorInFunction
            << "Edge label " << edgei
            << " out of range 0..11"
            << exit(FatalError);
    }

    const int nCurved = calcEdgePointsWeights
    (
        edgePoints,
        edgeWeights,
        hexCell::modelEdges()[edgei],

        nDiv,
        gd
    );

    return nCurved;
}


bool Foam::blockDescriptor::edgePointsWeights
(
    const label edgei,
    pointField& edgePoints,
    scalarList& edgeWeights
) const
{
    if (edgei < 0 || edgei >= 12)  //< hexCell::nEdges()
    {
        FatalErrorInFunction
            << "Edge label " << edgei
            << " out of range 0..11"
            << exit(FatalError);
    }

    const int nCurved = calcEdgePointsWeights
    (
        edgePoints,
        edgeWeights,
        hexCell::modelEdges()[edgei],

        sizes()[edgei/4],   // 12 edges -> 3 components (x,y,z)
        expand_[edgei]
    );

    return nCurved;
}


// ************************************************************************* //