semi_lagrangian3.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_LAGRANGIAN3_H_
#define INCLUDE_JET_SEMI_LAGRANGIAN3_H_
 
#include <jet/advection_solver3.h>
#include <limits>
 
namespace jet {
 
class SemiLagrangian3 : public AdvectionSolver3 {
 public:
    SemiLagrangian3();
 
    virtual ~SemiLagrangian3();
 
    void advect(const ScalarGrid3& input, const VectorField3& flow, double dt,
                ScalarGrid3* output,
                const ScalarField3& boundarySdf = ConstantScalarField3(
                    std::numeric_limits<double>::max())) final;
 
    void advect(const CollocatedVectorGrid3& input, const VectorField3& flow,
                double dt, CollocatedVectorGrid3* output,
                const ScalarField3& boundarySdf = ConstantScalarField3(
                    std::numeric_limits<double>::max())) final;
 
    void advect(const FaceCenteredGrid3& input, const VectorField3& flow,
                double dt, FaceCenteredGrid3* output,
                const ScalarField3& boundarySdf = ConstantScalarField3(
                    std::numeric_limits<double>::max())) final;
 
 protected:
    virtual std::function<double(const Vector3D&)> getScalarSamplerFunc(
        const ScalarGrid3& input) const;
 
    virtual std::function<Vector3D(const Vector3D&)> getVectorSamplerFunc(
        const CollocatedVectorGrid3& input) const;
 
    virtual std::function<Vector3D(const Vector3D&)> getVectorSamplerFunc(
        const FaceCenteredGrid3& input) const;
 
 private:
    Vector3D backTrace(const VectorField3& flow, double dt, double h,
                       const Vector3D& pt0, const ScalarField3& boundarySdf);
};
 
typedef std::shared_ptr<SemiLagrangian3> SemiLagrangian3Ptr;
 
}  // namespace jet
 
#endif  // INCLUDE_JET_SEMI_LAGRANGIAN3_H_