sph_system_data3.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_SPH_SYSTEM_DATA3_H_
#define INCLUDE_JET_SPH_SYSTEM_DATA3_H_
 
#include <jet/constants.h>
#include <jet/particle_system_data3.h>
#include <vector>
 
namespace jet {
 
class SphSystemData3 : public ParticleSystemData3 {
 public:
    SphSystemData3();
 
    explicit SphSystemData3(size_t numberOfParticles);
 
    SphSystemData3(const SphSystemData3& other);
 
    virtual ~SphSystemData3();
 
    void setRadius(double newRadius) override;
 
    void setMass(double newMass) override;
 
    ConstArrayAccessor1<double> densities() const;
 
    ArrayAccessor1<double> densities();
 
    ConstArrayAccessor1<double> pressures() const;
 
    ArrayAccessor1<double> pressures();
 
    void updateDensities();
 
    void setTargetDensity(double targetDensity);
 
    double targetDensity() const;
 
    void setTargetSpacing(double spacing);
 
    double targetSpacing() const;
 
    void setRelativeKernelRadius(double relativeRadius);
 
    void setKernelRadius(double kernelRadius);
 
    double relativeKernelRadius() const;
 
    double kernelRadius() const;
 
    double sumOfKernelNearby(const Vector3D& position) const;
 
    double interpolate(const Vector3D& origin,
                       const ConstArrayAccessor1<double>& values) const;
 
    Vector3D interpolate(const Vector3D& origin,
                         const ConstArrayAccessor1<Vector3D>& values) const;
 
    Vector3D gradientAt(size_t i,
                        const ConstArrayAccessor1<double>& values) const;
 
    double laplacianAt(size_t i,
                       const ConstArrayAccessor1<double>& values) const;
 
    Vector3D laplacianAt(size_t i,
                         const ConstArrayAccessor1<Vector3D>& values) const;
 
    void buildNeighborSearcher();
 
    void buildNeighborLists();
 
    void serialize(std::vector<uint8_t>* buffer) const override;
 
    void deserialize(const std::vector<uint8_t>& buffer) override;
 
    void set(const SphSystemData3& other);
 
    SphSystemData3& operator=(const SphSystemData3& other);
 
 private:
    double _targetDensity = kWaterDensity;
 
    double _targetSpacing = 0.1;
 
    double _kernelRadiusOverTargetSpacing = 1.8;
 
    double _kernelRadius;
 
    size_t _pressureIdx;
 
    size_t _densityIdx;
 
    void computeMass();
};
 
typedef std::shared_ptr<SphSystemData3> SphSystemData3Ptr;
 
}  // namespace jet
 
#endif  // INCLUDE_JET_SPH_SYSTEM_DATA3_H_