<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7652.24">
<TITLE>Problem: can't write down double precision data to .vtp file using vtkXMLPolyDataWriter ! </TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>
<BR>
<P><FONT SIZE=2>Hi All,<BR>
<BR>
I am trying to write double precison data to a .vtp file using vtkXMLPolyDataWriter but it writes down only 5 digits after the point decimal.During debugging, I found cout behaving the same way and only a few digits after the decimal point were displayed on the screen which was sorted out by setting the precision of cout explicitly using cout.setprecision() function.This suggests, we probably need to set the precision of the vtkPolyDataWriter too.<BR>
<BR>
Any suggestion plese, I am really stuck badly just at the end of my job !<BR>
<BR>
The code which I am using is as follows:<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
#include <vld.h><BR>
#include "vtkXMLPolyDataWriter.h"<BR>
<BR>
#include "vtkXMLWriter.h"<BR>
<BR>
#include "vtkPolyData.h"<BR>
<BR>
#include "vtkPointData.h"<BR>
<BR>
#include "vtkIdType.h"<BR>
<BR>
#include <iostream><BR>
<BR>
#include <cmath><BR>
<BR>
using namespace std;<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
int main( int argc, char *argv[] )<BR>
<BR>
{<BR>
<BR>
<BR>
<BR>
vtkXMLPolyDataReader *reader1 = vtkXMLPolyDataReader::New();<BR>
<BR>
reader1->SetFileName("C:\\VtkProjects\\normal\\conemodel.vtp");<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
vtkPolyData *data1 =reader1->GetOutput();<BR>
<BR>
data1->BuildLinks();<BR>
<BR>
data1->Update();<BR>
<BR>
<BR>
<BR>
vtkPoints * pointSet = vtkPoints::New();<BR>
<BR>
pointSet->SetDataTypeToDouble();<BR>
<BR>
<BR>
<BR>
double point1 [3];<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
//Load pointSet from polydata <BR>
<BR>
for (int i =0; i< (data1->GetNumberOfPoints());i++)<BR>
<BR>
{<BR>
<BR>
data1->GetPoint(i,point1);<BR>
<BR>
pointSet->InsertPoint(i,point1);<BR>
<BR>
}<BR>
<BR>
<BR>
<BR>
double * buffer1 = new double[3];<BR>
<BR>
vtkIdType pointId = 0;<BR>
<BR>
buffer1[0]=1.809017843;<BR>
<BR>
buffer1[1]=1.809017490;<BR>
<BR>
buffer1[2]=1.809017840;<BR>
<BR>
<BR>
//To confirm the data type of pointSet<BR>
<BR>
cout.setprecision(10);<BR>
<BR>
cout<<"data type of pointSet"<<pointSet->GetDataType()<<endl;<BR>
<BR>
<BR>
// change the first element of pointSet <BR>
<BR>
<BR>
<BR>
double* aikPoint1 = new double[3];<BR>
<BR>
pointSet->SetPoint(pointId, buffer1);<BR>
<BR>
aikPoint1= pointSet->GetPoint(0);<BR>
<BR>
<BR>
// Display the contents of changed value<BR>
cout<<aikPoint1[0]<<endl;<BR>
<BR>
cout<<aikPoint1[1]<<endl;<BR>
<BR>
cout<<aikPoint1[2]<<endl<<endl;<BR>
<BR>
<BR>
<BR>
//linkNew pointSet with the data<BR>
<BR>
data1->SetPoints(pointSet);<BR>
<BR>
<BR>
// check if we see the updted value directly within the polydata<BR>
data1->GetPoint(0,aikPoint);<BR>
<BR>
<BR>
<BR>
cout<<aikPoint[0]<<endl;<BR>
<BR>
cout<<aikPoint[1]<<endl;<BR>
<BR>
cout<<aikPoint[2]<<endl;<BR>
<BR>
<BR>
//**************************************** SO FAR SO GOOD *************************************************<BR>
vtkPolyData * test = vtkPolyData::New();<BR>
<BR>
test->Allocate(data1,1000,1000);<BR>
<BR>
test->SetPoints(pointSet);<BR>
<BR>
test->SetPolys(data1->GetPolys());<BR>
<BR>
test->Update();<BR>
<BR>
<BR>
<BR>
vtkXMLPolyDataWriter *writer1 = vtkXMLPolyDataWriter::New();<BR>
<BR>
writer1->SetInput(test);<BR>
<BR>
<BR>
<BR>
writer1->SetDataModeToAscii();<BR>
<BR>
writer1->SetFileName("testDoubleFloatCone.vtp");<BR>
<BR>
writer1->Write();<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
return(0);<BR>
<BR>
<BR>
<BR>
}<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
Best Regards,<BR>
<BR>
Salman<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>