MeshDomain.h

Go to the documentation of this file.

#pragma once
 
 
#include "objects/containers/Volume.h"
#include "objects/finders/Bvh.h"
#include "objects/geometry/Domain.h"
#include "objects/geometry/Triangle.h"
 
NAMESPACE_SPH_BEGIN
 
class IScheduler;
 
struct MeshParams {
    AffineMatrix matrix = AffineMatrix::identity();
 
    bool precomputeInside = true;
 
    Size volumeResolution = 128;
};
 
 
class MeshDomain : public IDomain {
private:
    Bvh<BvhTriangle> bvh;
    Volume<char> mask;
 
    struct {
        Array<Vector> points;
        Array<Vector> normals;
        Box box;
        Float volume;
        Float area;
    } cached;
 
public:
    MeshDomain(IScheduler& scheduler, Array<Triangle>&& triangles, const MeshParams& params = MeshParams{});
 
    virtual Vector getCenter() const override {
        return cached.box.center();
    }
 
    virtual Box getBoundingBox() const override {
        return cached.box;
    }
 
    virtual Float getVolume() const override {
        return cached.volume;
    }
 
    virtual Float getSurfaceArea() const override {
        return cached.area;
    }
 
    virtual bool contains(const Vector& v) const override;
 
    virtual void getSubset(ArrayView<const Vector> vs,
        Array<Size>& output,
        const SubsetType type) const override;
 
    virtual void getDistanceToBoundary(ArrayView<const Vector>, Array<Float>&) const override {
        NOT_IMPLEMENTED;
    }
 
    virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices) const override;
 
    virtual void addGhosts(ArrayView<const Vector> vs,
        Array<Ghost>& ghosts,
        const Float eta,
        const Float eps) const override;
 
private:
    bool containImpl(const Vector& v) const;
};
 
NAMESPACE_SPH_END