bodies.h

Go to the documentation of this file.

#pragma once

#include "domain.h"
#include "parameterDB.h"
#include "body.h"


template <typename memoryType>
class bodies
{
public:
    int  numBodies,   
         totalPoints; 

    bool bodiesMove;  

    cusp::array1d<int, memoryType>
      numPoints,    
      offsets,      
      I,            
      J;            
    
    cusp::array1d<int, memoryType>
      startI,       
      startJ,       
      numCellsX,    
      numCellsY;    
        
    cusp::array1d<real, memoryType>
      xmin,  
      xmax,  
      ymin,  
      ymax;  
    
    cusp::array1d<real, memoryType>
      forceX,   
      forceY;   

    cusp::array1d<real, memoryType>
      X,     
      Y,     
      ds,    
      ones,  
      x,     
      y,     
      uB,    
      vB;    

    // set initial position and velocity of each body
    void initialise(parameterDB &db, domain &D);
    
    // store index of each cell that contains a boundary point
    void calculateCellIndices(domain &D);
    
    // store indices of the bounding box of each body
    void calculateBoundingBoxes(parameterDB &db, domain &D);
    
    // update position, velocity and neighbors of each body
    void update(parameterDB &db, domain &D, real Time);

    // write body coordinates into a file
    void writeToFile(std::string &caseFolder, int timeStep);

    // write body coordinates into a file called \a bodies
    void writeToFile(real *bx, real *by, std::string &caseFolder, int timeStep);
}; // bodies