[cig-commits] r11291 - in long/3D/Gale/trunk: . src/Underworld/plugins/Output/VTKOutput
walter at geodynamics.org
walter at geodynamics.org
Thu Feb 28 15:57:34 PST 2008
Author: walter
Date: 2008-02-28 15:57:34 -0800 (Thu, 28 Feb 2008)
New Revision: 11291
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/src/Underworld/plugins/Output/VTKOutput/VTKOutput.c
Log:
r2016 at earth: boo | 2008-02-28 15:57:10 -0800
Make VTKOutput use Float64 instead of Float32, since the internal numbers are 8 byte floats
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2014
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2016
Modified: long/3D/Gale/trunk/src/Underworld/plugins/Output/VTKOutput/VTKOutput.c
===================================================================
--- long/3D/Gale/trunk/src/Underworld/plugins/Output/VTKOutput/VTKOutput.c 2008-02-28 20:53:10 UTC (rev 11290)
+++ long/3D/Gale/trunk/src/Underworld/plugins/Output/VTKOutput/VTKOutput.c 2008-02-28 23:57:34 UTC (rev 11291)
@@ -156,7 +156,7 @@
fprintf(fp," <Piece NumberOfPoints=\"%d\" NumberOfCells=\"1\">\n",
(num_particles-1)/stepping+1);
fprintf(fp," <Points>\n");
- fprintf(fp," <DataArray type=\"Float32\" NumberOfComponents=\"3\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" NumberOfComponents=\"3\" format=\"ascii\">\n");
if(myRank==0)
{
@@ -164,7 +164,7 @@
<VTKFile type=\"PUnstructuredGrid\" version=\"0.1\" byte_order=\"LittleEndian\">\n\
<PUnstructuredGrid GhostLevel=\"0\">\n\
<PPoints>\n\
- <DataArray type=\"Float32\" NumberOfComponents=\"3\" format=\"ascii\"/>\n\
+ <DataArray type=\"Float64\" NumberOfComponents=\"3\" format=\"ascii\"/>\n\
</PPoints>\n\
<PCellData></PCellData>\n");
}
@@ -293,42 +293,42 @@
fprintf(fp,"\n </DataArray>\n");
fprintf(fp," </Points>\n");
fprintf(fp," <PointData Scalars=\"Viscosity\" Tensors=\"Stress\">\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Viscosity\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Viscosity\" format=\"ascii\">\n");
if(myRank==0)
{
fprintf(pfp," <PPointData Scalars=\"Viscosity\" Tensors=\"Stress\">\n");
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Viscosity\" format=\"ascii\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Viscosity\" format=\"ascii\"/>\n");
}
break;
case 1:
fprintf(fp,"\n </DataArray>\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Yielding_fraction\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Yielding_fraction\" format=\"ascii\">\n");
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Yielding_fraction\" format=\"ascii\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Yielding_fraction\" format=\"ascii\"/>\n");
break;
case 2:
fprintf(fp,"\n </DataArray>\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Stress\" format=\"ascii\" NumberOfComponents=\"9\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Stress\" format=\"ascii\" NumberOfComponents=\"9\">\n");
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Stress\" format=\"ascii\" NumberOfComponents=\"9\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Stress\" format=\"ascii\" NumberOfComponents=\"9\"/>\n");
break;
case 3:
fprintf(fp,"\n </DataArray>\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Density\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Density\" format=\"ascii\">\n");
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Density\" format=\"ascii\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Density\" format=\"ascii\"/>\n");
break;
case 4:
fprintf(fp,"\n </DataArray>\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Alpha\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Alpha\" format=\"ascii\">\n");
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Alpha\" format=\"ascii\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Alpha\" format=\"ascii\"/>\n");
break;
case 5:
fprintf(fp,"\n </DataArray>\n");
- fprintf(fp," <DataArray type=\"Float32\" Name=\"Thermal_Diffusivity\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"Thermal_Diffusivity\" format=\"ascii\">\n");
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"Thermal_Diffusivity\" format=\"ascii\"/>\n");
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"Thermal_Diffusivity\" format=\"ascii\"/>\n");
break;
case 6:
fprintf(fp,"\n </DataArray>\n");
@@ -372,7 +372,7 @@
IJK ijk;
fprintf(fp," <Points>\n");
- fprintf(fp," <DataArray type=\"Float32\" NumberOfComponents=\"3\" format=\"ascii\">\n");
+ fprintf(fp," <DataArray type=\"Float64\" NumberOfComponents=\"3\" format=\"ascii\">\n");
for(ijk[2]=lower[2];ijk[2]<upper[2];++ijk[2])
for(ijk[1]=lower[1];ijk[1]<upper[1];++ijk[1])
@@ -539,7 +539,7 @@
<PStructuredGrid GhostLevel=\"0\" WholeExtent=\"0 %d 0 %d 0 %d\">\n\
<PCellData></PCellData>\n\
<PPoints>\n\
- <PDataArray type=\"Float32\" NumberOfComponents=\"3\"/>\n\
+ <PDataArray type=\"Float64\" NumberOfComponents=\"3\"/>\n\
</PPoints>\n\
<PPointData Scalars=\"StrainRateInvariantField\" Vectors=\"VelocityField\" Tensors=\"Stress\">\n",global[0]-1,global[1]-1,global[2]-1);
fprintf(ppressure_fp,"<?xml version=\"1.0\"?>\n\
@@ -547,7 +547,7 @@
<PStructuredGrid GhostLevel=\"0\" WholeExtent=\"0 %d 0 %d 0 %d\">\n\
<PCellData></PCellData>\n\
<PPoints>\n\
- <PDataArray type=\"Float32\" NumberOfComponents=\"3\"/>\n\
+ <PDataArray type=\"Float64\" NumberOfComponents=\"3\"/>\n\
</PPoints>\n\
<PPointData Scalars=\"PressureField\">\n",
global[0]-2,global[1]-2,p_global[2]-1);
@@ -583,25 +583,25 @@
/* Scalars */
case 2:
case 3:
- fprintf(fp," <DataArray type=\"Float32\" Name=\"%s\" format=\"ascii\">\n",feVar->name);
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"%s\" format=\"ascii\">\n",feVar->name);
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"%s\" format=\"ascii\"/>\n",feVar->name);
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"%s\" format=\"ascii\"/>\n",feVar->name);
break;
/* Vectors */
case 4:
case 9:
- fprintf(fp," <DataArray type=\"Float32\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"3\">\n",feVar->name);
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"3\">\n",feVar->name);
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"3\"/>\n",feVar->name);
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"3\"/>\n",feVar->name);
break;
/* Rank 2 Tensors */
case 6:
case 8:
case 18:
case 27:
- fprintf(fp," <DataArray type=\"Float32\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"9\">\n",feVar->name);
+ fprintf(fp," <DataArray type=\"Float64\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"9\">\n",feVar->name);
if(myRank==0)
- fprintf(pfp," <PDataArray type=\"Float32\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"9\"/>\n",feVar->name);
+ fprintf(pfp," <PDataArray type=\"Float64\" Name=\"%s\" format=\"ascii\" NumberOfComponents=\"9\"/>\n",feVar->name);
break;
/* Unknown */
default:
More information about the cig-commits
mailing list