boundaryShapes3D.h

Go to the documentation of this file.

/* This file is part of the Palabos library.
 *
 * Copyright (C) 2011 FlowKit Sarl
 * Avenue de Chailly 23
 * 1012 Lausanne, Switzerland
 * E-mail contact: contact@flowkit.com
 *
 * The most recent release of Palabos can be downloaded at 
 * <http://www.palabos.org/>
 *
 * The library Palabos is free software: you can redistribute it and/or
 * modify it under the terms of the GNU Affero General Public License as
 * published by the Free Software Foundation, either version 3 of the
 * License, or (at your option) any later version.
 *
 * The library 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 Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef BOUNDARY_SHAPES_3D_H
#define BOUNDARY_SHAPES_3D_H

#include "core/globalDefs.h"
#include "atomicBlock/atomicBlock3D.h"
#include "core/geometry3D.h"
#include "core/array.h"

namespace plb {

struct OffBoundary {
    enum Type {dirichlet, neumann, freeSlip, constRhoInlet, densityNeumann, flux, isolation};
};

template<typename T, class SurfaceData>
struct BoundaryShape3D {
    virtual ~BoundaryShape3D() { }

    virtual bool isInside(Dot3D const& location) const =0;
    virtual bool pointOnSurface (
            Array<T,3> const& fromPoint, Array<T,3> const& direction,
            Array<T,3>& locatedPoint, T& distance,
            Array<T,3>& wallNormal, SurfaceData& surfaceData,
            OffBoundary::Type& bdType, plint& id ) const =0;
    virtual bool gridPointOnSurface (
            Dot3D const& fromPoint, Dot3D const& direction,
            Array<T,3>& locatedPoint, T& distance,
            Array<T,3>& wallNormal, SurfaceData& surfaceData,
            OffBoundary::Type& bdType, plint& id ) const
    {
        return pointOnSurface(Array<T,3>((T)fromPoint.x, (T)fromPoint.y, (T)fromPoint.z),
                              Array<T,3>((T)direction.x, (T)direction.y, (T)direction.z),
                              locatedPoint, distance, wallNormal, surfaceData, bdType, id);
    }
    virtual Array<T,3> computeContinuousNormal (
            Array<T,3> const& p, plint id, bool isAreaWeighted = false ) const =0;
    virtual bool intersectsSurface (
            Array<T,3> const& p1, Array<T,3> const& p2, plint& id ) const =0;
    virtual bool gridIntersectsSurface (
            Dot3D const& p1, Dot3D const& p2, plint& id ) const
    {
        return intersectsSurface( Array<T,3>((T)p1.x, (T)p1.y, (T)p1.z),
                                  Array<T,3>((T)p2.x, (T)p2.y, (T)p2.z), id );
    }
    virtual bool distanceToSurface( Array<T,3> const& point,
                                    T& distance, bool& isBehind ) const =0;
    virtual bool gridDistanceToSurface( Dot3D const& point,
                                        T& distance, bool& isBehind ) const
    {
        return distanceToSurface( Array<T,3>((T)point.x, (T)point.y, (T)point.z),
                                  distance, isBehind );
    }
    virtual plint getTag(plint id) const =0;
    virtual BoundaryShape3D<T,SurfaceData>* clone() const =0;
    virtual BoundaryShape3D<T,SurfaceData>* clone(std::vector<AtomicBlock3D*> args) const {
        return clone();
    }
};

template<typename T, class SurfaceData>
class BoundaryShapeIntersection3D : public BoundaryShape3D<T,SurfaceData> {
public:
    BoundaryShapeIntersection3D();
    virtual ~BoundaryShapeIntersection3D();
    BoundaryShapeIntersection3D(BoundaryShapeIntersection3D<T,SurfaceData> const& rhs);
    BoundaryShapeIntersection3D<T,SurfaceData>& operator=(BoundaryShapeIntersection3D<T,SurfaceData> const& rhs);
    void swap(BoundaryShapeIntersection3D<T,SurfaceData>& rhs);
    void addComponent(BoundaryShape3D<T,SurfaceData>* component);

    virtual bool isInside(Array<T,3> const& location) const;
    virtual bool isInside(Array<T,3> const& location, T epsilon) const;
    virtual bool pointOnSurface (
            Array<T,3> const& fromPoint, Array<T,3> const& direction,
            Array<T,3>& locatedPoint, T& distance,
            Array<T,3>& wallNormal, SurfaceData& surfaceData,
            OffBoundary::Type& bdType, plint& id ) const;
    virtual bool intersectsSurface (
            Array<T,3> const& p1, Array<T,3> const& p2, plint& id ) const;
    virtual bool distanceToSurface( Array<T,3> const& point,
                                    T& distance, bool& isBehind ) const;
    virtual plint getTag(plint id) const;
    virtual BoundaryShapeIntersection3D<T,SurfaceData>* clone() const;
private:
    std::vector<BoundaryShape3D<T,SurfaceData>*> components;
};

}  // namespace plb

#endif  // BOUNDARY_SHAPES_3D_H