PCISPHSolver2.h

Go to the documentation of this file.

/*************************************************************************
> File Name: PCISPHSolver2.h
> Project Name: CubbyFlow
> Author: Chan-Ho Chris Ohk
> Purpose: 2-D PCISPH solver.
> Created Time: 2017/06/09
> Copyright (c) 2018, Chan-Ho Chris Ohk
*************************************************************************/
#ifndef CUBBYFLOW_PCISPH_SOLVER2_H
#define CUBBYFLOW_PCISPH_SOLVER2_H

#include <Core/Solver/Particle/SPH/SPHSolver2.h>

namespace CubbyFlow
{
    class PCISPHSolver2 : public SPHSolver2
    {
    public:
        class Builder;

        PCISPHSolver2();

        PCISPHSolver2(double targetDensity, double targetSpacing, double relativeKernelRadius);

        virtual ~PCISPHSolver2();

        double GetMaxDensityErrorRatio() const;

        void SetMaxDensityErrorRatio(double ratio);

        unsigned int GetMaxNumberOfIterations() const;

        void SetMaxNumberOfIterations(unsigned int n);

        static Builder GetBuilder();

    protected:
        void AccumulatePressureForce(double timeIntervalInSeconds) override;

        void OnBeginAdvanceTimeStep(double timeStepInSeconds) override;

    private:
        double m_maxDensityErrorRatio = 0.01;
        unsigned int m_maxNumberOfIterations = 5;

        ParticleSystemData2::VectorData m_tempPositions;
        ParticleSystemData2::VectorData m_tempVelocities;
        ParticleSystemData2::VectorData m_pressureForces;
        ParticleSystemData2::ScalarData m_densityErrors;

        double ComputeDelta(double timeStepInSeconds) const;
        double ComputeBeta(double timeStepInSeconds) const;
    };

    using PCISPHSolver2Ptr = std::shared_ptr<PCISPHSolver2>;

    class PCISPHSolver2::Builder final : public SPHSolverBuilderBase2<PCISPHSolver2::Builder>
    {
    public:
        PCISPHSolver2 Build() const;

        PCISPHSolver2Ptr MakeShared() const;
    };
}

#endif