PCISPHSolver3.h

Go to the documentation of this file.

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

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

namespace CubbyFlow
{
    class PCISPHSolver3 : public SPHSolver3
    {
    public:
        class Builder;

        PCISPHSolver3();

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

        virtual ~PCISPHSolver3();

        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;

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

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

    using PCISPHSolver3Ptr = std::shared_ptr<PCISPHSolver3>;

    class PCISPHSolver3::Builder final : public SPHSolverBuilderBase3<PCISPHSolver3::Builder>
    {
    public:
        PCISPHSolver3 Build() const;

        PCISPHSolver3Ptr MakeShared() const;
    };
}

#endif