semi_lagrangian2.h

Go to the documentation of this file.

// Copyright (c) 2018 Doyub Kim
//
// I am making my contributions/submissions to this project solely in my
// personal capacity and am not conveying any rights to any intellectual
// property of any third parties.
 
#ifndef INCLUDE_JET_SEMI_LAGRANGIAN2_H_
#define INCLUDE_JET_SEMI_LAGRANGIAN2_H_
 
#include <jet/advection_solver2.h>
#include <limits>
 
namespace jet {
 
class SemiLagrangian2 : public AdvectionSolver2 {
 public:
    SemiLagrangian2();
 
    virtual ~SemiLagrangian2();
 
    void advect(const ScalarGrid2& input, const VectorField2& flow, double dt,
                ScalarGrid2* output, const ScalarField2& boundarySdf =
                                         ConstantScalarField2(kMaxD)) final;
 
    void advect(
        const CollocatedVectorGrid2& input, const VectorField2& flow, double dt,
        CollocatedVectorGrid2* output,
        const ScalarField2& boundarySdf = ConstantScalarField2(kMaxD)) final;
 
    void advect(
        const FaceCenteredGrid2& input, const VectorField2& flow, double dt,
        FaceCenteredGrid2* output,
        const ScalarField2& boundarySdf = ConstantScalarField2(kMaxD)) final;
 
 protected:
    virtual std::function<double(const Vector2D&)> getScalarSamplerFunc(
        const ScalarGrid2& input) const;
 
    virtual std::function<Vector2D(const Vector2D&)> getVectorSamplerFunc(
        const CollocatedVectorGrid2& input) const;
 
    virtual std::function<Vector2D(const Vector2D&)> getVectorSamplerFunc(
        const FaceCenteredGrid2& input) const;
 
 private:
    Vector2D backTrace(const VectorField2& flow, double dt, double h,
                       const Vector2D& pt0, const ScalarField2& boundarySdf);
};
 
typedef std::shared_ptr<SemiLagrangian2> SemiLagrangian2Ptr;
 
}  // namespace jet
 
#endif  // INCLUDE_JET_SEMI_LAGRANGIAN2_H_