baseMaterial.h

//  Copyright (c) 2019 Prashant K. Jha
//  Copyright (c) 2019 Patrick Diehl
//
//  Distributed under the Boost Software License, Version 1.0. (See accompanying
//  file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
 
#ifndef MATERIAL_PD_BASEMATERIAL_H
#define MATERIAL_PD_BASEMATERIAL_H
 
#include <cstring>
#include <vector>
 
#include "util/point.h"
 
namespace material {
 
namespace pd {
 
class BaseMaterial {
 
public:
  BaseMaterial(const size_t &dim, const double &horizon)
      : d_horizon(horizon), d_dimension(dim), d_stateActive(false), d_density
      (0.), d_name (""){};
 
  bool isStateActive() const { return d_stateActive; };
 
  std::string name() const {return d_name;};
 
  virtual std::pair<util::Point3,double> getBondEF(size_t i , size_t j){
 
      return std::make_pair(util::Point3(), 0.);
 
  };
 
  virtual double getS(const util::Point3 &dx, const util::Point3 &du) {
    return
        0.;
  };
 
  virtual double getS(size_t i, size_t j) { return 0.; };
 
  virtual double getSc(const double &r) { return 0.; };
 
  virtual double getSc(size_t i , size_t j) { return 0.; };
 
  virtual util::Point3 getBondForceDirection(const util::Point3 &dx,
                                     const util::Point3 &du) const {return
                                     util::Point3(); }
 
  virtual util::Matrix33 getStrain(size_t i) {
 
    return util::Matrix33();
  }
 
  virtual util::Matrix33 getStress(size_t i) {
 
    return util::Matrix33();
  }
 
  virtual void update(){}
 
  virtual double getInfFn(const double &r) const {return 0.; };
 
  double getHorizon() const {
    return d_horizon;
  }
 
  double getDensity() const {return d_density; };
 
  virtual std::string printStr(int nt = 0, int lvl = 0) const {
    // TODO implement
    return "BaseMaterial";
  };
 
  virtual void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt,
      lvl); };
 
protected:
 
  double d_horizon;
 
  size_t d_dimension;
 
  bool d_stateActive;
 
  double d_density;
 
  std::string d_name;
};
 
} // namespace pd
 
} // namespace material
 
#endif // MATERIAL_PD_BASEMATERIAL_H