void VTKWriter::initializeOutput(int numParticles) {
vtkFile = new VTKFile_t("UnstructuredGrid");
// per point, we add type, position, velocity and force
PointData pointData;
DataArray_t mass(type::Float32, "mass", 1);
DataArray_t velocity(type::Float32, "velocity", 3);
DataArray_t forces(type::Float32, "force", 3);
DataArray_t type(type::Int32, "type", 1);
pointData.DataArray().push_back(mass);
pointData.DataArray().push_back(velocity);
pointData.DataArray().push_back(forces);
pointData.DataArray().push_back(type);
CellData cellData; // we don't have cell data => leave it empty
// 3 coordinates
Points points;
DataArray_t pointCoordinates(type::Float32, "points", 3);
points.DataArray().push_back(pointCoordinates);
Cells cells; // we don't have cells, => leave it empty
// for some reasons, we have to add a dummy entry for paraview
DataArray_t cells_data(type::Float32, "types", 0);
cells.DataArray().push_back(cells_data);
PieceUnstructuredGrid_t piece(pointData, cellData, points, cells, numParticles, 0);
UnstructuredGrid_t unstructuredGrid(piece);
vtkFile->UnstructuredGrid(unstructuredGrid);
}
void VTKWriter::initializeOutput(int numParticles) {
vtkFile = new VTKFile_t("UnstructuredGrid");
// per point, we add mass, velocity, force, type and stress
// and in DEBUG mode the id of the molecules
PointData pointData;
DataArray_t mass(type::Float32, "mass", 1);
DataArray_t velocity(type::Float32, "velocity", 3);
DataArray_t forces(type::Float32, "force", 3);
DataArray_t type(type::Int32, "type", 1);
DataArray_t stress(type::Float32, "stress", 1);
DataArray_t flag(type::Int32, "flag", 1);
pointData.DataArray().push_back(mass);
pointData.DataArray().push_back(velocity);
pointData.DataArray().push_back(forces);
pointData.DataArray().push_back(type);
pointData.DataArray().push_back(stress);
pointData.DataArray().push_back(flag);
#ifdef DEBUG
DataArray_t id(type::Int32, "id", 1);
pointData.DataArray().push_back(id);
#endif
CellData cellData; // we don't have cell data => leave it empty
// 3 coordinates
Points points;
DataArray_t pointCoordinates(type::Float32, "points", 3);
points.DataArray().push_back(pointCoordinates);
Cells cells; // we don't have cells, => leave it empty
// for some reasons, we have to add a dummy entry for paraview
DataArray_t cells_data(type::Float32, "types", 0);
cells.DataArray().push_back(cells_data);
PieceUnstructuredGrid_t piece(pointData, cellData, points, cells, numParticles, 0);
UnstructuredGrid_t unstructuredGrid(piece);
vtkFile->UnstructuredGrid(unstructuredGrid);
}
