GridForwardEulerDiffusionSolver2.h

Go to the documentation of this file.

/*************************************************************************
> File Name: GridForwardEulerDiffusionSolver2.h
> Project Name: CubbyFlow
> Author: Chan-Ho Chris Ohk
> Purpose: 2-D grid-based forward Euler diffusion solver.
> Created Time: 2017/08/11
> Copyright (c) 2018, Chan-Ho Chris Ohk
*************************************************************************/
#ifndef CUBBYFLOW_GRID_FORWARD_EULER_DIFFUSION_SOLVER2_H
#define CUBBYFLOW_GRID_FORWARD_EULER_DIFFUSION_SOLVER2_H

#include <Core/Solver/Grid/GridDiffusionSolver2.h>

namespace CubbyFlow
{
    class GridForwardEulerDiffusionSolver2 final : public GridDiffusionSolver2
    {
    public:
        GridForwardEulerDiffusionSolver2();

        void Solve(
            const ScalarGrid2& source,
            double diffusionCoefficient,
            double timeIntervalInSeconds,
            ScalarGrid2* dest,
            const ScalarField2& boundarySDF = ConstantScalarField2(std::numeric_limits<double>::max()),
            const ScalarField2& fluidSDF = ConstantScalarField2(-std::numeric_limits<double>::max())) override;

        void Solve(
            const CollocatedVectorGrid2& source,
            double diffusionCoefficient,
            double timeIntervalInSeconds,
            CollocatedVectorGrid2* dest,
            const ScalarField2& boundarySDF = ConstantScalarField2(std::numeric_limits<double>::max()),
            const ScalarField2& fluidSDF = ConstantScalarField2(-std::numeric_limits<double>::max())) override;

        void Solve(
            const FaceCenteredGrid2& source,
            double diffusionCoefficient,
            double timeIntervalInSeconds,
            FaceCenteredGrid2* dest,
            const ScalarField2& boundarySDF = ConstantScalarField2(std::numeric_limits<double>::max()),
            const ScalarField2& fluidSDF = ConstantScalarField2(-std::numeric_limits<double>::max())) override;

    private:
        Array2<char> m_markers;

        void BuildMarkers(
            const Size2& size,
            const std::function<Vector2D(size_t, size_t)>& pos,
            const ScalarField2& boundarySDF,
            const ScalarField2& fluidSDF);
    };

    using GridForwardEulerDiffusionSolver2Ptr = std::shared_ptr<GridForwardEulerDiffusionSolver2>;
}

#endif