surface_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_SURFACE_KERNEL_HPP__
#define __BITPIT_SURFACE_KERNEL_HPP__
 
#include "patch_kernel.hpp"
#include "line_kernel.hpp"
 
#include "bitpit_common.hpp"
 
namespace bitpit {
 
class SurfaceKernel : public PatchKernel {
 
public:
    // Static constant
    BITPIT_PUBLIC_API static const unsigned short SELECT_TRIANGLE;
    BITPIT_PUBLIC_API static const unsigned short SELECT_QUAD;
    BITPIT_PUBLIC_API static const unsigned short SELECT_ALL;
 
    // Types definitions
    typedef double (SurfaceKernel::*eval_f_)(long, int &) const;
 
    void setSpaceDimension(int dimension);
    int getSpaceDimension(void) const;
 
    int getVolumeCodimension() const override;
    int getSurfaceCodimension() const override;
    int getLineCodimension() const override;
    int getPointCodimension() const override;
 
    void extractEnvelope(LineKernel &envelope) const;
 
    virtual double evalCellArea(long) const;
    virtual double evalEdgeLength(long, int) const;
    virtual double evalMinEdgeLength(long, int &) const;
    virtual double evalMaxEdgeLength(long, int &) const;
    virtual double evalAngleAtVertex(long, int) const;
    virtual double evalMinAngleAtVertex(long, int &) const;
    virtual double evalMaxAngleAtVertex(long, int &) const;
    virtual double evalAspectRatio(long, int &) const;
    virtual std::array<double, 3> evalFacetNormal(long, const std::array<double, 3> &orientation = {{0., 0., 1.}}) const;
    std::array<double, 3> evalEdgeNormal(long, int) const;
    std::array<double, 3> evalVertexNormal(long, int) const;
    std::array<double, 3> evalVertexNormal(long, int, std::size_t, const long *) const;
    std::array<double, 3> evalLimitedVertexNormal(long, int, double ) const;
    std::array<double, 3> evalLimitedVertexNormal(long, int, std::size_t, const long *, double ) const;
    virtual void evalVertexNormals(long id, int vertex, std::size_t nVertexNeighs, const long *vertexNeighs, double limit,
                                   std::array<double, 3> *unlimitedNormal, std::array<double, 3> *limitedNormal) const;
    virtual void evalEdgeNormals(long id, int edge, double limit, std::array<double, 3> *unlimitedNormal,
                                 std::array<double, 3> *limitedNormal) const;
    double evalCellSize(long id) const override;
    void evalBarycentricCoordinates(long id, const std::array<double, 3> &point, double *lambda) const;
 
    bool isCellOrientationConsistent() const;
    bool adjustCellOrientation();
    bool adjustCellOrientation(long id, bool invert = false);
    void flipCellOrientation(long id);
 
    ConstProxyVector<long> getFacetOrderedVertexIds(const Cell &facet) const;
    bool areFacetVerticesOrdered(const Cell &facet) const;
    int getFacetOrderedLocalVertex(const Cell &facet, std::size_t n) const;
 
    ConstProxyVector<long> getFacetOrderedEdgeIds(const Cell &facet) const;
    bool areFacetEdgesOrdered(const Cell &facet) const;
    int getFacetOrderedLocalEdge(const Cell &facet, std::size_t n) const;
 
    void displayQualityStats(std::ostream&, unsigned int padding = 0) const;
    std::vector<double> computeHistogram(eval_f_ funct_, std::vector<double> &bins, long &count, int n_intervals = 8, unsigned short mask = SELECT_ALL) const;
 
private:
    static const std::map<ElementType, unsigned short> m_selectionTypes;
 
    void initialize();
 
    bool compareSelectedTypes(unsigned short, ElementType) const;
    void displayHistogram(long, const std::vector<double>&, const std::vector<double>&, const std::string&, std::ostream&, unsigned int padding = 0) const;
 
    bool haveSameOrientation(const Cell &cell_A, int face_A, const Cell &cell_B, int face_B) const;
 
protected:
#if BITPIT_ENABLE_MPI==1
    SurfaceKernel(MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
    SurfaceKernel(int dimension, MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
    SurfaceKernel(int id, int dimension, MPI_Comm communicator, std::size_t haloSize, AdaptionMode adaptionMode, PartitioningMode partitioningMode);
#else
    SurfaceKernel(AdaptionMode adaptionMode);
    SurfaceKernel(int dimension, AdaptionMode adaptionMode);
    SurfaceKernel(int id, int dimension, AdaptionMode adaptionMode);
#endif
 
};
 
}
 
#endif