 }
 
 /**
-* \brief Calculate the indices of the cells in which the boundary points are located.
-*
-* It calculates the index of the x-coordinate and the index of the y-coordinate
-* of the bottom-left node of the containing cell.
-* This information is useful when transferring data between the boundary points
-* and the computational grid.
-*
-* \param D information about the computational grid
-*/
+ * \brief Stores index of each cell that contains a boundary point.
+ *
+ * It calculates the index of the x-coordinate and the index of the y-coordinate
+ * of the bottom-left node of each cell that contains a boundary point.
+ * This information is useful when transferring data between the boundary points
+ * and the computational grid.
+ *
+ * \param D information about the computational grid
+ */
 template <typename memoryType>
 void bodies<memoryType>::calculateCellIndices(domain &D)
 {

}

/**
* \brief Write the body coordinates into a file (on the host)
*
* \param caseFolder folder of the simulation
* \param timeStep time-step of the simulation
*/
template <>
void bodies<host_memory>::writeToFile(std::string &caseFolder, int timeStep)
{

/***************************************************************************//**
 * \file bodies.cu
 * \author Anush Krishnan (anush@bu.edu)
 * \brief Implementation of the methods of the class \c bodies.
 */


#include "bodies.h"


/**
 * \brief Stores initial position and velocity of each body in the flow.
 *
 * \param db database that contains all the simulation parameters
 * \param D information about the computational grid
 */
template <typename memoryType>
void bodies<memoryType>::initialise(parameterDB &db, domain &D)
{
