element_reference.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_ELEMENT_REFERENCE_HPP__
#define __BITPIT_ELEMENT_REFERENCE_HPP__
 
#include <assert.h>
#include <vector>
 
#include "bitpit_common.hpp"
 
#include "element_type.hpp"
 
namespace bitpit {
 
class ReferenceElementInfo {
 
public:
    BITPIT_PUBLIC_API static const int MAX_ELEM_VERTICES = 8;
    BITPIT_PUBLIC_API static const int MAX_ELEM_FACES    = 6;
    BITPIT_PUBLIC_API static const int MAX_ELEM_EDGES    = 12;
 
    static bool hasInfo(ElementType type);
    BITPIT_PUBLIC_API static const ReferenceElementInfo & getInfo(ElementType type);
 
    int dimension;
    ElementType type;
 
    int nVertices;
    int nFaces;
    int nEdges;
 
    std::array<ElementType, MAX_ELEM_FACES> faceTypeStorage;
    std::array<std::array<int, MAX_ELEM_VERTICES>, MAX_ELEM_FACES> faceConnectStorage;
    std::array<std::array<int, MAX_ELEM_FACES>, MAX_ELEM_FACES> faceEdgeStorage;
 
    std::array<ElementType, MAX_ELEM_EDGES> edgeTypeStorage;
    std::array<std::array<int, MAX_ELEM_VERTICES>, MAX_ELEM_EDGES> edgeConnectStorage;
 
    virtual ~ReferenceElementInfo() = default;
 
    virtual double evalSize(const std::array<double, 3> *vertexCoords) const = 0;
 
    virtual void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const = 0;
    virtual double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const = 0;
 
protected:
    ReferenceElementInfo(int _dimension, ElementType _type, int _nVertices, int _nFaces, int _nEdges);
 
    ReferenceElementInfo(ReferenceElementInfo const &) = delete;
    ReferenceElementInfo & operator=(ReferenceElementInfo const &) = delete;
 
    void initializeFaceEdges(const std::vector<const ReferenceElementInfo *> &facesInfo, const std::vector<const ReferenceElementInfo *> &edgesInfo);
 
};
 
class Reference3DElementInfo : public ReferenceElementInfo {
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    virtual double evalVolume(const std::array<double, 3> *vertexCoords) const = 0;
    double evalSurfaceArea(const std::array<double, 3> *vertexCoords) const;
    double evalEdgePerimeter(const std::array<double, 3> *vertexCoords) const;
 
    double evalFaceArea(int face, const std::array<double, 3> *vertexCoords) const;
    double evalEdgeLength(int edge, const std::array<double, 3> *vertexCoords) const;
 
    void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const override;
    double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const override;
 
protected:
    Reference3DElementInfo(ElementType type, int nVertices, int nFaces);
 
};
 
class ReferenceTetraInfo : public Reference3DElementInfo {
 
friend class ReferenceElementInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    double evalVolume(const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceTetraInfo info;
 
protected:
    ReferenceTetraInfo();
 
    ReferenceTetraInfo(ReferenceTetraInfo const &) = delete;
    ReferenceTetraInfo & operator=(ReferenceTetraInfo const &) = delete;
 
};
 
class ReferenceVoxelInfo : public Reference3DElementInfo {
 
friend class ReferenceElementInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    double evalVolume(const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceVoxelInfo info;
 
protected:
    ReferenceVoxelInfo();
 
    ReferenceVoxelInfo(ReferenceVoxelInfo const &) = delete;
    ReferenceVoxelInfo & operator=(ReferenceVoxelInfo const &) = delete;
 
};
 
class ReferenceHexahedronInfo : public Reference3DElementInfo {
 
friend class ReferenceElementInfo;
 
public:
    double evalVolume(const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceHexahedronInfo info;
 
protected:
    ReferenceHexahedronInfo();
 
    ReferenceHexahedronInfo(ReferenceHexahedronInfo const &) = delete;
    ReferenceHexahedronInfo & operator=(ReferenceHexahedronInfo const &) = delete;
 
};
 
class ReferencePyramidInfo : public Reference3DElementInfo {
 
friend class ReferenceElementInfo;
 
public:
    double evalVolume(const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferencePyramidInfo info;
 
protected:
    ReferencePyramidInfo();
 
    ReferencePyramidInfo(ReferencePyramidInfo const &) = delete;
    ReferencePyramidInfo & operator=(ReferencePyramidInfo const &) = delete;
 
};
 
class ReferenceWedgeInfo : public Reference3DElementInfo {
 
friend class ReferenceElementInfo;
 
public:
    double evalVolume(const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceWedgeInfo info;
 
protected:
    ReferenceWedgeInfo();
 
    ReferenceWedgeInfo(ReferenceWedgeInfo const &) = delete;
    ReferenceWedgeInfo & operator=(ReferenceWedgeInfo const &) = delete;
 
};
 
class Reference2DElementInfo : public ReferenceElementInfo {
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    virtual double evalArea(const std::array<double, 3> *vertexCoords) const = 0;
    double evalPerimeter(const std::array<double, 3> *vertexCoords) const;
 
    double evalFaceLength(int face, const std::array<double, 3> *vertexCoords) const;
 
    virtual std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const = 0;
 
    void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const override;
    double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const override;
 
    virtual bool areVerticesCCWOrdered() const;
    virtual int getCCWOrderedVertex(int n) const;
 
    virtual bool areFacesCCWOrdered() const;
    virtual int getCCWOrderedFace(int n) const;
 
protected:
    Reference2DElementInfo(ElementType type, int nVertices);
 
    void getCCWVertexCoords(const std::array<double, 3> *vertexCoords, const std::array<double, 3> **ccwVertexCoords, std::array<double, 3> *ccwVertexCoordsStorage) const;
 
};
 
class ReferenceTriangleInfo : public Reference2DElementInfo {
 
friend class ReferenceElementInfo;
friend class ReferenceTetraInfo;
friend class ReferencePyramidInfo;
friend class ReferenceWedgeInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    double evalArea(const std::array<double, 3> *vertexCoords) const override;
 
    std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const override;
 
    void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const override;
    double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceTriangleInfo info;
 
protected:
    ReferenceTriangleInfo();
 
    ReferenceTriangleInfo(ReferenceTriangleInfo const &) = delete;
    ReferenceTriangleInfo & operator=(ReferenceTriangleInfo const &) = delete;
 
};
 
class ReferencePixelInfo : public Reference2DElementInfo {
 
friend class ReferenceElementInfo;
friend class ReferenceVoxelInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    double evalArea(const std::array<double, 3> *vertexCoords) const override;
 
    std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const override;
 
    bool areVerticesCCWOrdered() const override;
    int getCCWOrderedVertex(int n) const override;
 
    bool areFacesCCWOrdered() const override;
    int getCCWOrderedFace(int n) const override;
 
    BITPIT_PUBLIC_API static const ReferencePixelInfo info;
 
protected:
    ReferencePixelInfo();
 
    ReferencePixelInfo(ReferencePixelInfo const &) = delete;
    ReferencePixelInfo & operator=(ReferencePixelInfo const &) = delete;
 
};
 
class ReferenceQuadInfo : public Reference2DElementInfo {
 
friend class ReferenceElementInfo;
friend class ReferenceHexahedronInfo;
friend class ReferencePyramidInfo;
friend class ReferenceWedgeInfo;
 
public:
    double evalArea(const std::array<double, 3> *vertexCoords) const override;
 
    std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const override;
 
    BITPIT_PUBLIC_API static const ReferenceQuadInfo info;
 
protected:
    ReferenceQuadInfo();
 
    ReferenceQuadInfo(ReferenceQuadInfo const &) = delete;
    ReferenceQuadInfo & operator=(ReferenceQuadInfo const &) = delete;
 
};
 
class Reference1DElementInfo : public ReferenceElementInfo {
 
public:
    virtual double evalLength(const std::array<double, 3> *vertexCoords) const = 0;
 
    virtual std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &orientation = {{0., 0., 1.}}, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const = 0;
 
protected:
    Reference1DElementInfo(ElementType type);
 
};
 
class ReferenceLineInfo : public Reference1DElementInfo {
 
friend class ReferenceElementInfo;
friend class ReferenceTriangleInfo;
friend class ReferencePixelInfo;
friend class ReferenceQuadInfo;
friend class ReferenceTetraInfo;
friend class ReferenceVoxelInfo;
friend class ReferenceHexahedronInfo;
friend class ReferencePyramidInfo;
friend class ReferenceWedgeInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    double evalLength(const std::array<double, 3> *vertexCoords) const override;
 
    std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &orientation = {{0., 0., 1.}}, const std::array<double, 3> &point = {{0.5, 0.5, 0.5}}) const override;
 
    void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const override;
    double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceLineInfo info;
 
protected:
    ReferenceLineInfo();
 
    ReferenceLineInfo(ReferenceLineInfo const &) = delete;
    ReferenceLineInfo & operator=(ReferenceLineInfo const &) = delete;
 
};
 
class Reference0DElementInfo : public ReferenceElementInfo {
 
public:
    virtual std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &orientation = {{1., 0., 0.}}) const = 0;
 
protected:
    Reference0DElementInfo(ElementType type);
 
};
 
class ReferenceVertexInfo : public Reference0DElementInfo {
 
friend class ReferenceElementInfo;
friend class ReferenceLineInfo;
friend class ReferenceTriangleInfo;
friend class ReferencePixelInfo;
friend class ReferenceQuadInfo;
 
public:
    double evalSize(const std::array<double, 3> *vertexCoords) const override;
 
    std::array<double, 3> evalNormal(const std::array<double, 3> *vertexCoords, const std::array<double, 3> &orientation = {{1., 0., 0.}}) const override;
 
    void evalPointProjection(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords, std::array<double, 3> *projection, double *distance) const override;
    double evalPointDistance(const std::array<double, 3> &point, const std::array<double, 3> *vertexCoords) const override;
 
    BITPIT_PUBLIC_API static const ReferenceVertexInfo info;
 
protected:
    ReferenceVertexInfo();
 
    ReferenceVertexInfo(ReferenceVertexInfo const &) = delete;
    ReferenceVertexInfo & operator=(ReferenceVertexInfo const &) = delete;
 
};
 
}
 
#endif