point_kernel.cpp

/*---------------------------------------------------------------------------*\
 *
 *  bitpit
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of bitpit.
 *
 *  bitpit is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License v3 (LGPL)
 *  as published by the Free Software Foundation.
 *
 *  bitpit 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 Lesser General Public
 *  License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with bitpit. If not, see <http://www.gnu.org/licenses/>.
 *
\*---------------------------------------------------------------------------*/
 
#include "point_kernel.hpp"
 
namespace bitpit {
 
#if BITPIT_ENABLE_MPI==1
PointKernel::PointKernel(MPI_Comm communicator,
                         AdaptionMode adaptionMode, PartitioningMode partitioningMode)
    : PatchKernel(communicator, 0, adaptionMode, partitioningMode)
#else
PointKernel::PointKernel(AdaptionMode adaptionMode)
    : PatchKernel(adaptionMode)
#endif
{
    initialize();
}
 
#if BITPIT_ENABLE_MPI==1
PointKernel::PointKernel(int dimension, MPI_Comm communicator,
                         AdaptionMode adaptionMode, PartitioningMode partitioningMode)
    : PatchKernel(dimension, communicator, 0, adaptionMode, partitioningMode)
#else
PointKernel::PointKernel(int dimension, AdaptionMode adaptionMode)
    : PatchKernel(dimension, adaptionMode)
#endif
{
    initialize();
}
 
#if BITPIT_ENABLE_MPI==1
PointKernel::PointKernel(int id, int dimension, MPI_Comm communicator,
                         AdaptionMode adaptionMode, PartitioningMode partitioningMode)
    : PatchKernel(id, dimension, communicator, 0, adaptionMode, partitioningMode)
#else
PointKernel::PointKernel(int id, int dimension, AdaptionMode adaptionMode)
    : PatchKernel(id, dimension, adaptionMode)
#endif
{
}
 
void PointKernel::initialize()
{
    // Nothing to do
}
 
int PointKernel::getVolumeCodimension() const
{
    return 3;
}
 
int PointKernel::getSurfaceCodimension() const
{
    return 2;
}
 
int PointKernel::getLineCodimension() const
{
    return 1;
}
 
int PointKernel::getPointCodimension() const
{
    return 0;
}
 
double PointKernel::evalCellSize(long id) const
{
    BITPIT_UNUSED(id);
 
    return 0.;
}
 
double PointKernel::evalPointsDistance(long id1, long id2) const
{
    const Cell &cell1 = m_cells[id1];
    const Cell &cell2 = m_cells[id2];
 
    ConstProxyVector<long> cellVertexIds1 = cell1.getVertexIds();
    const Vertex &vertex1 = getVertex(cellVertexIds1[0]);
 
    ConstProxyVector<long> cellVertexIds2 = cell2.getVertexIds();
    const Vertex &vertex2 = getVertex(cellVertexIds2[0]);
 
    double distance = norm2(vertex2.getCoords() - vertex1.getCoords());
 
    return distance;
}
 
std::array<double, 3> PointKernel::evalPointsDirection(long id1, long id2) const
{
    const Cell &cell1 = m_cells[id1];
    const Cell &cell2 = m_cells[id2];
 
    ConstProxyVector<long> cellVertexIds1 = cell1.getVertexIds();
    const Vertex &vertex1 = getVertex(cellVertexIds1[0]);
 
    ConstProxyVector<long> cellVertexIds2 = cell2.getVertexIds();
    const Vertex &vertex2 = getVertex(cellVertexIds2[0]);
 
    std::array<double, 3> direction = vertex2.getCoords() - vertex1.getCoords();
 
    double distance = norm2(direction);
    if (!utils::DoubleFloatingEqual()(distance, 0., getTol())) {
        direction /= distance;
    }
 
    return direction;
}
 
}