void VTKWriter::plotMDCell(double x,double y, double z, double rho, utils::Vector<double, 3> v) {
if (vtkFile->UnstructuredGrid().present()) {
//cout << "UnstructuredGrid is present" << endl;
} else {
cout << "ERROR: No UnstructuredGrid present" << endl;
}
PointData::DataArray_sequence& pointDataSequence = vtkFile->UnstructuredGrid()->Piece().PointData().DataArray();
PointData::DataArray_iterator dataIterator = pointDataSequence.begin();
if(rho>0.65)
std::cout <<"rho" <<rho << " x" << x<<" y"<<y<<" z"<<z<<" rrho"<<(rho*6.69/39.304*pow(10,4))<<std::endl;
dataIterator->push_back(rho*6.69/39.304*pow(10,4));
//cout << "Appended mass data in: " << dataIterator->Name();
dataIterator++;
//dataIterator->push_back(v[0]*157);
//dataIterator->push_back(v[1]*157);
//dataIterator->push_back(v[2]*157);
dataIterator->push_back(v[0]);
dataIterator->push_back(v[1]);
dataIterator->push_back(v[2]);
//cout << "Appended velocity data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(0.0);
dataIterator->push_back(0.0);
dataIterator->push_back(0.0);
//cout << "Appended force data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(199);
dataIterator++;
dataIterator->push_back(0);
dataIterator++;
dataIterator->push_back( 0);
#ifdef DEBUG
dataIterator++;
dataIterator->push_back(199);
#endif
Points::DataArray_sequence& pointsSequence = vtkFile->UnstructuredGrid()->Piece().Points().DataArray();
Points::DataArray_iterator pointsIterator = pointsSequence.begin();
pointsIterator->push_back(x);
pointsIterator->push_back(y);
pointsIterator->push_back(z);
}
void VTKWriter::plotParticle(Particle& p) {
if (vtkFile->UnstructuredGrid().present()) {
LOG4CXX_TRACE(iolog, "UnstructuredGrid is present");
} else {
LOG4CXX_ERROR(iolog, "No UnstructuredGrid present");
}
PointData::DataArray_sequence& pointDataSequence = vtkFile->UnstructuredGrid()->Piece().PointData().DataArray();
PointData::DataArray_iterator dataIterator = pointDataSequence.begin();
dataIterator->push_back(p.getM());
//cout << "Appended mass data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getV()[0]);
dataIterator->push_back(p.getV()[1]);
dataIterator->push_back(p.getV()[2]);
//cout << "Appended velocity data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getOldF()[0]);
dataIterator->push_back(p.getOldF()[1]);
dataIterator->push_back(p.getOldF()[2]);
//cout << "Appended force data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getType());
Points::DataArray_sequence& pointsSequence = vtkFile->UnstructuredGrid()->Piece().Points().DataArray();
Points::DataArray_iterator pointsIterator = pointsSequence.begin();
pointsIterator->push_back(p.getX()[0]);
pointsIterator->push_back(p.getX()[1]);
pointsIterator->push_back(p.getX()[2]);
}
void VTKWriter::plotCell(OffsetArray<double> *_p, OffsetArray<utils::Vector<double,2> > *_v, int i, int j, double cellSize, double xOffset, double yOffset, double lb2md) {
if (vtkFile->UnstructuredGrid().present()) {
//cout << "UnstructuredGrid is present" << endl;
} else {
std::cout << "ERROR: No UnstructuredGrid present" << std::endl;
}
PointData::DataArray_sequence& pointDataSequence = vtkFile->UnstructuredGrid()->Piece().PointData().DataArray();
PointData::DataArray_iterator dataIterator = pointDataSequence.begin();
dataIterator->push_back(_p->get(i,j));
//cout << "Appended mass data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(_v->get(i,j)[0]*901*lb2md);
dataIterator->push_back(_v->get(i,j)[1]*901*lb2md);
dataIterator->push_back(0.0);
//cout << "Appended velocity data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(0.0);
dataIterator->push_back(0.0);
dataIterator->push_back(0.0);
//cout << "Appended force data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(99);
dataIterator++;
dataIterator->push_back(0);
dataIterator++;
dataIterator->push_back( 0);
#ifdef DEBUG
dataIterator++;
dataIterator->push_back(99);
#endif
Points::DataArray_sequence& pointsSequence = vtkFile->UnstructuredGrid()->Piece().Points().DataArray();
Points::DataArray_iterator pointsIterator = pointsSequence.begin();
pointsIterator->push_back(((double)i*cellSize+xOffset*cellSize+cellSize/2.0));
pointsIterator->push_back(((double)j*cellSize+yOffset*cellSize+cellSize/2.0));
pointsIterator->push_back(cellSize/2.0);
}
void VTKWriter::plotParticle(Particle& p) {
if (vtkFile->UnstructuredGrid().present()) {
//cout << "UnstructuredGrid is present" << endl;
} else {
cout << "ERROR: No UnstructuredGrid present" << endl;
}
PointData::DataArray_sequence& pointDataSequence = vtkFile->UnstructuredGrid()->Piece().PointData().DataArray();
PointData::DataArray_iterator dataIterator = pointDataSequence.begin();
dataIterator->push_back(p.getM());
//cout << "Appended mass data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getV()[0]*157);
dataIterator->push_back(p.getV()[1]*157);
dataIterator->push_back(p.getV()[2]*157);
//cout << "Appended velocity data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getF()[0]);
dataIterator->push_back(p.getF()[1]);
dataIterator->push_back(p.getF()[2]);
//cout << "Appended force data in: " << dataIterator->Name();
dataIterator++;
dataIterator->push_back(p.getType());
dataIterator++;
dataIterator->push_back(p.getStress());
dataIterator++;
dataIterator->push_back(p.getFlag() == true? 1 : 0);
#ifdef DEBUG
dataIterator++;
dataIterator->push_back(p.getID());
#endif
Points::DataArray_sequence& pointsSequence = vtkFile->UnstructuredGrid()->Piece().Points().DataArray();
Points::DataArray_iterator pointsIterator = pointsSequence.begin();
pointsIterator->push_back(p.getX()[0]);
pointsIterator->push_back(p.getX()[1]);
pointsIterator->push_back(p.getX()[2]);
}