imstkSphState.h

/*
** This file is part of the Interactive Medical Simulation Toolkit (iMSTK)
** iMSTK is distributed under the Apache License, Version 2.0.
** See accompanying NOTICE for details.
*/

#pragma once

#include "imstkMath.h"

namespace imstk
{
template<typename T> class DataArray;
template<typename T, int N> class VecDataArray;

struct NeighborInfo
{
    Vec3d relativePos;      
    double density;         
};

class SphState
{
public:
    SphState() = default;
    SphState(const int numElements);
    virtual ~SphState() = default;

public:
    void setBoundaryParticlePositions(std::shared_ptr<VecDataArray<double, 3>> positions) { m_boundaryParticlePositions = positions; }

    size_t getNumParticles() const;

    std::shared_ptr<VecDataArray<double, 3>> getPositions() const { return m_positions; }
    void setPositions(std::shared_ptr<VecDataArray<double, 3>> positions) { m_positions = positions; }

    std::shared_ptr<VecDataArray<double, 3>> getFullStepVelocities() const { return m_fullStepVelocities; }

    std::shared_ptr<VecDataArray<double, 3>> getHalfStepVelocities() const { return m_halfStepVelocities; }

    std::shared_ptr<VecDataArray<double, 3>> getVelocities() const { return m_velocities; }
    void setVelocities(std::shared_ptr<VecDataArray<double, 3>> velocities) { m_velocities = velocities; }

    std::shared_ptr<VecDataArray<double, 3>> getBoundaryParticlePositions() const { return m_boundaryParticlePositions; }

    std::shared_ptr<VecDataArray<double, 3>> getNormals() const { return m_normals; }

    std::shared_ptr<DataArray<double>> getDensities() const { return m_densities; }

    std::shared_ptr<VecDataArray<double, 3>> getAccelerations() const { return m_acceleration; }

    std::shared_ptr<VecDataArray<double, 3>> getDiffuseVelocities() const { return m_diffuseVelocities; }

    std::vector<std::vector<size_t>>& getFluidNeighborLists() { return m_neighborLists; }
    const std::vector<std::vector<size_t>>& getFluidNeighborLists() const { return m_neighborLists; }

    std::vector<std::vector<size_t>>& getBoundaryNeighborLists() { return m_boundaryParticleNeighborLists; }
    const std::vector<std::vector<size_t>>& getBoundaryNeighborLists() const { return m_boundaryParticleNeighborLists; }

    std::vector<std::vector<NeighborInfo>>& getNeighborInfo() { return m_neighborInfo; }
    const std::vector<std::vector<NeighborInfo>>& getNeighborInfo() const { return m_neighborInfo; }

    void setState(std::shared_ptr<SphState> rhs);

private:
    std::shared_ptr<VecDataArray<double, 3>> m_positions;
    std::shared_ptr<VecDataArray<double, 3>> m_fullStepVelocities;
    std::shared_ptr<VecDataArray<double, 3>> m_halfStepVelocities;
    std::shared_ptr<VecDataArray<double, 3>> m_velocities;

    std::shared_ptr<VecDataArray<double, 3>> m_boundaryParticlePositions;   

    std::shared_ptr<DataArray<double>>       m_densities;                   
    std::shared_ptr<VecDataArray<double, 3>> m_normals;                     
    std::shared_ptr<VecDataArray<double, 3>> m_acceleration;                
    std::shared_ptr<VecDataArray<double, 3>> m_diffuseVelocities;           

    std::vector<std::vector<size_t>>       m_neighborLists;                 
    std::vector<std::vector<size_t>>       m_boundaryParticleNeighborLists; 
    std::vector<std::vector<NeighborInfo>> m_neighborInfo;                  
};
} // namespace imstk