IterativeLevelSetSolver2.h

Go to the documentation of this file.

/*************************************************************************
> File Name: IterativeLevelSetSolver2.h
> Project Name: CubbyFlow
> Author: Chan-Ho Chris Ohk
> Purpose: Abstract base class for 2-D PDE-based iterative level set solver.
> Created Time: 2017/08/30
> Copyright (c) 2018, Chan-Ho Chris Ohk
*************************************************************************/
#ifndef CUBBYFLOW_ITERATIVE_LEVEL_SET_SOLVER2_H
#define CUBBYFLOW_ITERATIVE_LEVEL_SET_SOLVER2_H

#include <Core/Solver/LevelSet/LevelSetSolver2.h>

namespace CubbyFlow
{
    class IterativeLevelSetSolver2 : public LevelSetSolver2
    {
    public:
        IterativeLevelSetSolver2();

        virtual ~IterativeLevelSetSolver2();

        void Reinitialize(
            const ScalarGrid2& inputSDF,
            double maxDistance,
            ScalarGrid2* outputSDF) override;

        void Extrapolate(
            const ScalarGrid2& input,
            const ScalarField2& sdf,
            double maxDistance,
            ScalarGrid2* output) override;

        void Extrapolate(
            const CollocatedVectorGrid2& input,
            const ScalarField2& sdf,
            double maxDistance,
            CollocatedVectorGrid2* output) override;

        void Extrapolate(
            const FaceCenteredGrid2& input,
            const ScalarField2& sdf,
            double maxDistance,
            FaceCenteredGrid2* output) override;

        double GetMaxCFL() const;

        void SetMaxCFL(double newMaxCfl);

    protected:
        virtual void GetDerivatives(
            ConstArrayAccessor2<double> grid,
            const Vector2D& gridSpacing,
            size_t i, size_t j,
            std::array<double, 2>* dx, std::array<double, 2>* dy) const = 0;

    private:
        double m_maxCFL = 0.5;

        void Extrapolate(
            const ConstArrayAccessor2<double>& input,
            const ConstArrayAccessor2<double>& sdf,
            const Vector2D& gridSpacing,
            double maxDistance,
            ArrayAccessor2<double> output);

        static unsigned int DistanceToNumberOfIterations(double distance, double dtau);

        static double Sign(
            const ConstArrayAccessor2<double>& sdf,
            const Vector2D& gridSpacing,
            size_t i, size_t j);

        double PseudoTimeStep(
            ConstArrayAccessor2<double> sdf,
            const Vector2D& gridSpacing) const;
    };

    using IterativeLevelSetSolver2Ptr = std::shared_ptr<IterativeLevelSetSolver2>;
}

#endif