void FluidSimulationSaveState::saveState(std::string filename, FluidSimulation *_fluidsim) {
std::ofstream erasefile;
erasefile.open(filename, std::ofstream::out | std::ofstream::trunc);
erasefile.close();
std::ofstream state(filename.c_str(), std::ios::out | std::ios::binary);
int i, j, k;
_fluidsim->getGridDimensions(&i, &j, &k);
double dx = _fluidsim->getCellSize();
_width = i;
_height = j;
_depth = k;
// i, j, k, dx
_writeInt(&i, &state);
_writeInt(&j, &state);
_writeInt(&k, &state);
_writeDouble(&dx, &state);
// number of marker particles
int n = _fluidsim->getNumMarkerParticles();
_writeInt(&n, &state);
// next frame to be processed
int currentFrame = _fluidsim->getCurrentFrame();
_writeInt(¤tFrame, &state);
// number of solid cell indices
int numIndices = _getNumSolidCells(_fluidsim);
_writeInt(&numIndices, &state);
// floats: marker particle positions in form [x1, y1, z1, x2, y2, z2, ...]
_writeBinaryMarkerParticlePositions(_fluidsim, &state);
// floats: marker particle velocities in form [x1, y1, z1, x2, y2, z2, ...]
_writeBinaryMarkerParticleVelocities(_fluidsim, &state);
// ints: solid cell indicies in form [i1, j1, k1, i2, j2, k2, ...]
_writeSolidCellIndices(_fluidsim, &state);
state.close();
}
void Writer::writeDouble(double value) {
_writeByte(trace::TYPE_DOUBLE);
_writeDouble(value);
}