template<typename V >

void vtkfig::DataSet::SetRectilinearGrid	(	const V &	x,
		const V &	y,
		const V &	z
	)

inline

Enter data of a 3D rectilinear grid.

Template Parameters

V	Vector class counting from zero with member functions size() and operator[]. std::vector will work.

Parameters

x	x-coordinates
y	y-coordinates
z	z-coordinates

Definition at line 498 of file vtkfigDataSet.h.

   {
     vtkSmartPointer<vtkDoubleArray> xcoord;
     vtkSmartPointer<vtkDoubleArray> ycoord;
     vtkSmartPointer<vtkDoubleArray> zcoord;
     assert(spacedim!=2);
     spacedim=3;
     int Nx = x.size();
     int Ny = y.size();
     int Nz = z.size();
     
     if (this->data==NULL)
     {
       this->data=vtkSmartPointer<vtkRectilinearGrid>::New();
       auto rdata=vtkRectilinearGrid::SafeDownCast(this->data);
       assert(rdata);
       xcoord = vtkSmartPointer<vtkDoubleArray>::New();
       ycoord = vtkSmartPointer<vtkDoubleArray>::New();
       zcoord = vtkSmartPointer<vtkDoubleArray>::New();
       rdata->SetXCoordinates(xcoord);
       rdata->SetYCoordinates(ycoord);
       rdata->SetZCoordinates(zcoord);
       rdata->SetDimensions(Nx, Ny, Nz );
     }
     else
     {
       auto rdata=vtkRectilinearGrid::SafeDownCast(this->data);
       assert(rdata);
       xcoord=vtkDoubleArray::SafeDownCast(rdata->GetXCoordinates());
       ycoord=vtkDoubleArray::SafeDownCast(rdata->GetYCoordinates());
       zcoord=vtkDoubleArray::SafeDownCast(rdata->GetZCoordinates());
       xcoord->Initialize();
       ycoord->Initialize();
       zcoord->Initialize();
       rdata->SetDimensions(Nx, Ny, Nz );
     }
     
     xcoord->SetNumberOfComponents(1);
     xcoord->SetNumberOfTuples(Nx);
     
     ycoord->SetNumberOfComponents(1);
     ycoord->SetNumberOfTuples(Ny);
     
     zcoord->SetNumberOfComponents(1);
     zcoord->SetNumberOfTuples(Nz);
     
     for (int i=0; i<Nx; i++)
       xcoord->InsertComponent(i, 0, x[i]*coordinate_scale_factor*coordinate_scale_factor_xyz[0]);
     for (int i=0; i<Ny; i++)
       ycoord->InsertComponent(i, 0, y[i]*coordinate_scale_factor*coordinate_scale_factor_xyz[1]);
     for (int i=0; i<Nz; i++)
       zcoord->InsertComponent(i, 0, z[i]*coordinate_scale_factor*coordinate_scale_factor_xyz[2]);
     
     this->data->Modified();
   }