pci_sph_solver3.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_PCI_SPH_SOLVER3_H_
#define INCLUDE_JET_PCI_SPH_SOLVER3_H_
 
#include <jet/sph_solver3.h>
 
namespace jet {
 
class PciSphSolver3 : public SphSolver3 {
 public:
    class Builder;
 
    PciSphSolver3();
 
    PciSphSolver3(
        double targetDensity,
        double targetSpacing,
        double relativeKernelRadius);
 
    virtual ~PciSphSolver3();
 
    double maxDensityErrorRatio() const;
 
    void setMaxDensityErrorRatio(double ratio);
 
    unsigned int maxNumberOfIterations() const;
 
    void setMaxNumberOfIterations(unsigned int n);
 
    static Builder builder();
 
 protected:
    void accumulatePressureForce(double timeIntervalInSeconds) override;
 
    void onBeginAdvanceTimeStep(double timeStepInSeconds) override;
 
 private:
    double _maxDensityErrorRatio = 0.01;
    unsigned int _maxNumberOfIterations = 5;
 
    ParticleSystemData3::VectorData _tempPositions;
    ParticleSystemData3::VectorData _tempVelocities;
    ParticleSystemData3::VectorData _pressureForces;
    ParticleSystemData3::ScalarData _densityErrors;
 
    double computeDelta(double timeStepInSeconds);
    double computeBeta(double timeStepInSeconds);
};
 
typedef std::shared_ptr<PciSphSolver3> PciSphSolver3Ptr;
 
class PciSphSolver3::Builder final
    : public SphSolverBuilderBase3<PciSphSolver3::Builder> {
 public:
    PciSphSolver3 build() const;
 
    PciSphSolver3Ptr makeShared() const;
};
 
}  // namespace jet
 
#endif  // INCLUDE_JET_PCI_SPH_SOLVER3_H_