level_set_shape.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
 *                              SPHinXsys                                   *
 * --------------------------------------------------------------------------*
 * SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle    *
 * Hydrodynamics for industrial compleX systems. It provides C++ APIs for   *
 * physical accurate simulation and aims to model coupled industrial dynamic *
 * systems including fluid, solid, multi-body dynamics and beyond with SPH  *
 * (smoothed particle hydrodynamics), a meshless computational method using *
 * particle discretization.                                                 *
 *                                                                          *
 * SPHinXsys is partially funded by German Research Foundation              *
 * (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1                *
 * and HU1527/12-1.                                                         *
 *                                                                           *
 * Portions copyright (c) 2017-2020 Technical University of Munich and      *
 * the authors' affiliations.                                               *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may   *
 * not use this file except in compliance with the License. You may obtain a *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
 *                                                                           *
 * --------------------------------------------------------------------------*/
#ifndef LEVEL_SET_SHAPE_H
#define LEVEL_SET_SHAPE_H

#include "base_geometry.h"
#include "level_set.h"

#include <string>

namespace SPH
{

    class SPHBody;

    class LevelSetShape : public Shape
    {
    private:
        UniquePtrKeeper<BaseLevelSet> level_set_keeper_;

    public:
        LevelSetShape(SPHBody *sph_body, Shape &shape, Real refinement_ratio = 1.0);
        virtual ~LevelSetShape(){};

        virtual bool checkContain(const Vecd &input_pnt, bool BOUNDARY_INCLUDED = true) override;
        virtual Vecd findClosestPoint(const Vecd &input_pnt) override;

        Vecd findLevelSetGradient(const Vecd &input_pnt);
        Real computeKernelIntegral(const Vecd &input_pnt, Real h_ratio = 1.0);
        Vecd computeKernelGradientIntegral(const Vecd &input_pnt, Real h_ratio = 1.0);
        LevelSetShape *cleanLevelSet(Real small_shift_factor = 1.0);
        void writeLevelSet(SPHBody &sph_body);

    protected:
        BaseLevelSet *level_set_; 
        virtual BoundingBox findBounds() override;
    };
}
#endif // LEVEL_SET_SHAPE_H