volume_kernel.hpp

/*---------------------------------------------------------------------------*\
 *
 *  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/>.
 *
\*---------------------------------------------------------------------------*/
 
#ifndef __BITPIT_VOLUME_KERNEL_HPP__
#define __BITPIT_VOLUME_KERNEL_HPP__
 
#include "patch_kernel.hpp"
#include "surface_kernel.hpp"
 
namespace bitpit {
 
class VolumeKernel : public PatchKernel {
 
public:
    int getVolumeCodimension() const override;
    int getSurfaceCodimension() const override;
    int getLineCodimension() const override;
    int getPointCodimension() const override;
 
    void extractEnvelope(SurfaceKernel &envelope) const;
 
    bool isPointInside(double x, double y, double z) const;
    virtual bool isPointInside(const std::array<double, 3> &point) const = 0;
    bool isPointInside(long id, double x, double y, double z) const;
    virtual bool isPointInside(long id, const std::array<double, 3> &point) const = 0;
 
    virtual double evalCellVolume(long id)const = 0;
 
    virtual double evalInterfaceArea(long id)const = 0;
        virtual std::array<double,3> evalInterfaceNormal(long id)const = 0;
 
    ConstProxyVector<long> getFaceOrderedVertexIds(const Cell &cell, int face) const;
    bool areFaceVerticesOrdered(const Cell &cell, int face) const;
    int getFaceOrderedLocalVertex(const Cell &cell, int face, std::size_t n) const;
 
protected:
#if BITPIT_ENABLE_MPI==1
    VolumeKernel(MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
    VolumeKernel(int dimension, MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
    VolumeKernel(int id, int dimension, MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
#else
    VolumeKernel(AdaptionMode adaptionMode);
    VolumeKernel(int dimension, AdaptionMode adaptionMode);
    VolumeKernel(int id, int dimension, AdaptionMode adaptionMode);
#endif
 
};
 
}
 
#endif