Shape.h

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// Copyright (c) 2021-2025 MPM-Geomechanics Development Team
 
#ifndef SHAPE_H_
#define SHAPE_H_
 
#include "Eigen/Core"
using Eigen::Vector3d;
 
class Shape {
 
public:
    
    Shape();
    
    virtual ~Shape();
    
    virtual void update(const Vector3d& particle_position, const Vector3d& nodal_position, const Vector3d& cell_dimension, const Vector3d& particle_size) = 0;
 
    inline const Vector3d& getShape() const { return this->shape; }
    
    inline const Vector3d& getDerivate() const { return this->derivate; }
    
    inline void setShape(double sx, double sy, double sz) { this->shape=Vector3d(sx,sy,sz); }
    
    inline void setDerivate(double gx, double gy, double gz) { this->derivate=Vector3d(gx,gy,gz); }
    
private:
    
    virtual double computeGradient(double pI_position, double cell_dimension, double lp) = 0;
    
    virtual double computeShape(double pI_position, double cell_dimension, double lp) = 0;
 
    Vector3d shape; 
    
    Vector3d derivate; 
};
 
inline Shape::Shape() {
    
    this->shape.setZero();
    this->derivate.setZero();
}
 
inline Shape::~Shape() {
    
}
 
#endif /* SHAPE_H_ */