 }
 
 /**
-* \brief Update the location of the body points.
-*
-* This is done using the formulae:
-*
-* \f$ x_{i,m} = X^c_m + (X_{i,m} - X^0_m) \cos\theta - (Y_{i,m} - Y^0_m) \sin\theta \f$
-*
-* and
-*
-* \f$ y_{i,m} = Y^c_m + (X_{i,m} - X^0_m) \sin\theta + (Y_{i,m} - Y^0_m) \cos\theta \f$
-*
-* \param db database that contains all the simulation parameters
-* \param D information about the computational grid
-* \param Time the time
-*/
+ * \brief Updates position and velocity of each body point and its neighbor.
+ *
+ * This is done using the formulae:
+ *
+ * \f$ x_{i,m} = X^c_m + (X_{i,m} - X^0_m) \cos\theta - (Y_{i,m} - Y^0_m) \sin\theta \f$
+ *
+ * and
+ *
+ * \f$ y_{i,m} = Y^c_m + (X_{i,m} - X^0_m) \sin\theta + (Y_{i,m} - Y^0_m) \cos\theta \f$
+ *
+ * \param db database that contains all the simulation parameters
+ * \param D information about the computational grid
+ * \param Time the time
+ */
 template <typename memoryType>
 void bodies<memoryType>::update(parameterDB &db, domain &D, real Time)
 {

/***************************************************************************//**
* \file bodies.cu
* \author Krishnan, A. (anush@bu.edu)
* \brief Definition of the class \c bodies
*/


#include "bodies.h"


/**
* \brief Initialize the arrays in the class with information from \c body instances.
*
* Information regarding the coordinates of the body points and motion of the bodies
* is stored on the host as an array of instances of the class \c body.
* This function transfers that information to arrays on the device,
* where they are stored as a structure of arrays.
* This makes computation more efficient.
*
* \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)
{

}

/**
 * \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)
{
