Hi, Using the C++ VTK library for writing an .vtr file, I'm able to create the mesh, but not assign any skalars/vectors to the grid points. This is more or less given in the example file RGrid.cxx.
But how do I assign values at the mesh-points? Write the mesh: #include <stdio.h> #include <vtkRectilinearGrid.h> #include <vtkXMLRectilinearGridWriter.h> #include <vtkDoubleArray.h> #include <vtkIntArray.h> #include <vtkSmartPointer.h> int main(){ int nx = 3, ny = 3, nz = 3; // coordinates vtkDoubleArray *xCoords = vtkDoubleArray::New(); for (int i=0; i<nx; i++) xCoords->InsertNextValue(i); vtkDoubleArray *yCoords = vtkDoubleArray::New(); for (int i=0; i<ny; i++) yCoords->InsertNextValue(i); vtkDoubleArray *zCoords = vtkDoubleArray::New(); for (int i=0; i<nz; i++) zCoords->InsertNextValue(i); // The coordinates are assigned to the rectilinear grid. vtkRectilinearGrid *rgrid = vtkRectilinearGrid::New(); rgrid->SetDimensions(nx,ny,nz); rgrid->SetXCoordinates(xCoords); rgrid->SetYCoordinates(yCoords); rgrid->SetZCoordinates(zCoords); /* Write to file */ vtkSmartPointer<vtkXMLRectilinearGridWriter> writer = vtkSmartPointer<vtkXMLRectilinearGridWriter>::New(); writer->SetFileName("test.vtr"); #if VTK_MAJOR_VERSION <= 5 writer->SetInput(rgrid); #else writer->SetInputData(rgrid); #endif writer->Write(); /* clean up */ xCoords->Delete(); yCoords->Delete(); zCoords->Delete(); rgrid->Delete(); printf("DONE printing\n"); } For setting values at the mesh-points, I assumed something like: vtkIntArray* temperature = vtkIntArray::New(); temperature->SetName("Temperature"); temperature->SetNumberOfComponents(1); temperature->SetNumberOfTuples(nx*ny*nz); for(int i=0;i<mesh.nn;i++){ temperature->SetValue(i,10); // set everything to 10 } rgrid->GetPointData()->SetScalars(temperature); But that doesn't work. As expected, when looking at http://www.vtk.org/doc/nightly/html/classvtkRectilinearGrid.html So how do I set values/vectors for all the points in the mesh. Regards, Paw
_______________________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView Follow this link to subscribe/unsubscribe: http://www.paraview.org/mailman/listinfo/paraview