grid.h

/* -*- c++ -*- ----------------------------------------------------------
 *
 *                    ***       Karamelo       ***
 *               Parallel Material Point Method Simulator
 * 
 * Copyright (2019) Alban de Vaucorbeil, alban.devaucorbeil@monash.edu
 * Materials Science and Engineering, Monash University
 * Clayton VIC 3800, Australia

 * This software is distributed under the GNU General Public License.
 *
 * ----------------------------------------------------------------------- */

#ifndef MPM_GRID_H
#define MPM_GRID_H

#include "pointers.h"
#include "grid.h"
#include "grid.h"
#include <vector>
#include <Eigen/Eigen>
#include <map>
#include <array>

using namespace Eigen;

struct Point {
  int owner;    
  tagint tag;   
  double x[3];  
  int ntype[3]; 
};


class Grid : protected Pointers {
 public:
  int ncells;            
  bigint nnodes;         
  bigint nnodes_local;   
  bigint nnodes_ghost;   
  vector<tagint> ntag;   
  map<int, int> map_ntag;

  int nx;                
  int ny;                
  int nz;                

  int nx_global;         
  int ny_global;         
  int nz_global;         

  int nshared;           
  vector<tagint> shared; 
  map<int, vector<tagint>> dest_nshared;   
  map<int, vector<tagint>> origin_nshared; 

  vector<int> nowner;    

  double cellsize;       

  vector<Eigen::Vector3d> x;            
  vector<Eigen::Vector3d> x0;           
  vector<Eigen::Vector3d> v;            
  vector<Eigen::Vector3d> v_update;     
  vector<Eigen::Vector3d> mb;           
  vector<Eigen::Vector3d> f;            

  vector<double> mass;              
  vector<int> mask;                 
  vector<bool> rigid;               
  vector<array<int, 3>> ntype;      

  vector<double> pH;

  MPI_Datatype Pointtype;    

  Grid(class MPM *);
  virtual ~Grid();
  void grow(int);              
  void grow_ghosts(int);       
  void setup(string);
  void init(double*, double*); 

  void reduce_mass_ghost_nodes();                  
  void reduce_mass_ghost_nodes_old();
  void reduce_ghost_nodes(bool only_v = false);    
  void reduce_ghost_nodes_old(bool only_v = false);
  void update_grid_velocities();                   
  void update_grid_positions();                    
  void reduce_regularized_variables();             
};

#endif