 }
 
 /**
-* \brief Write the body coordinates into a file called \a bodies
-*
-* \param bx x-coordinates of the bodies
-* \param by y-coordinates of the bodies
-* \param caseFolder folder of the simulation
-* \param timeStep time-step of the simulation
-*/
+ * \brief Writes body coordinates into a file called \a bodies.
+ *
+ * \param bx x-coordinate of all points of all bodies
+ * \param by y-coordinate of all points of all bodies
+ * \param caseFolder directory of the simulation
+ * \param timeStep time-step of the simulation
+ */
 template <typename memoryType>
 void bodies<memoryType>::writeToFile(real *bx, real *by, std::string &caseFolder, int timeStep)
 {

}

/**
* \brief Write the body coordinates into a file (on the device)
*
* \param caseFolder folder of the simulation
* \param timeStep time-step of the simulation
*/
template <>
void bodies<device_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)
{
