fracture.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 GEOM_FRACTURE_H
#define GEOM_FRACTURE_H
 
#include "util/point.h" // definition of Point3
#include <inp/decks/fractureDeck.h>
#include <stdint.h> // uint8_t type
#include <string.h> // size_t type
#include <vector>
 
// forward declaration of fracture deck
namespace inp {
struct EdgeCrack;
struct FractureDeck;
} // namespace inp
 
namespace geometry {
 
class Fracture {
 
public:
  Fracture(inp::FractureDeck *deck, const std::vector<util::Point3> *nodes,
           const std::vector<std::vector<size_t>> *neighbor_list);
 
  explicit Fracture(inp::FractureDeck *deck);
 
  bool addCrack(const double &time, const std::vector<util::Point3> *nodes,
            const std::vector<std::vector<size_t>> *neighbor_list);
 
  void setBondState(const size_t &i, const size_t &j, const bool &state);
 
  bool getBondState(const size_t &i, const size_t &j) const;
 
  std::vector<uint8_t> getBonds(const size_t &i) const;
 
  std::string printStr(int nt = 0, int lvl = 0) const;
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); };
 
private:
  void computeFracturedBondFd(const size_t &i, inp::EdgeCrack *crack,
                              const std::vector<util::Point3> *nodes,
                              const std::vector<size_t> *neighbors);
 
  inp::FractureDeck *d_fractureDeck_p;
 
  std::vector<std::vector<uint8_t>> d_fracture;
};
 
} // namespace geometry
 
#endif // GEOM_FRACTURE_H