From 887e07a3943b6833a9da4b16fbeb1ba8e87b484b Mon Sep 17 00:00:00 2001
From: Thorsten Liebig <Thorsten.Liebig@gmx.de>
Date: Mon, 4 Oct 2010 09:30:24 +0200
Subject: [PATCH] define vtk dump data type

Using double as vtk data type will fix some errors/warnings in Paraview about a datasize mismatch.
---
 FDTD/processfields.cpp | 12 ++++++------
 FDTD/processfields.h   |  2 ++
 2 files changed, 8 insertions(+), 6 deletions(-)

diff --git a/FDTD/processfields.cpp b/FDTD/processfields.cpp
index 16c8b38..be2c7a9 100644
--- a/FDTD/processfields.cpp
+++ b/FDTD/processfields.cpp
@@ -305,15 +305,15 @@ void ProcessFields::WriteVTKCartesianGridHeader(ofstream &file, double const* co
 	file << "ASCII" << endl;
 	file << "DATASET RECTILINEAR_GRID " << endl;
 	file << "DIMENSIONS " << numLines[0] << " " << numLines[1] << " " << numLines[2] << endl;
-	file << "X_COORDINATES " << numLines[0] << " float" << endl;
+	file << "X_COORDINATES " << numLines[0] << " " << __VTK_DATA_TYPE__ << endl;
 	for (unsigned int i=0;i<numLines[0];++i)
 		file << setprecision(precision) << discLines[0][i] * discLines_scaling << " ";
 	file << endl;
-	file << "Y_COORDINATES " << numLines[1] << " float" << endl;
+	file << "Y_COORDINATES " << numLines[1] << " " << __VTK_DATA_TYPE__ << endl;
 	for (unsigned int i=0;i<numLines[1];++i)
 		file << setprecision(precision) << discLines[1][i] * discLines_scaling << " ";
 	file << endl;
-	file << "Z_COORDINATES " << numLines[2] << " float" << endl;
+	file << "Z_COORDINATES " << numLines[2] << " " << __VTK_DATA_TYPE__ << endl;
 	for (unsigned int i=0;i<numLines[2];++i)
 		file << setprecision(precision) << discLines[2][i] * discLines_scaling << " ";
 	file << endl << endl;
@@ -330,7 +330,7 @@ void ProcessFields::WriteVTKCylindricalGridHeader(ofstream &file, double const*
 	file << "ASCII" << endl;
 	file << "DATASET STRUCTURED_GRID " << endl;
 	file << "DIMENSIONS " << numLines[0] << " " << numLines[1] << " " << numLines[2] << endl;
-	file << "POINTS " << numLines[0]*numLines[1]*numLines[2] << " float" << endl;
+	file << "POINTS " << numLines[0]*numLines[1]*numLines[2] << " " << __VTK_DATA_TYPE__ << endl;
 	for (unsigned int k=0;k<numLines[2];++k)
 		for (unsigned int j=0;j<numLines[1];++j)
 			for (unsigned int i=0;i<numLines[0];++i)
@@ -347,7 +347,7 @@ void ProcessFields::WriteVTKCylindricalGridHeader(ofstream &file, double const*
 
 void ProcessFields::WriteVTKVectorArray(ofstream &file, string name, FDTD_FLOAT const* const* const* const* array, double const* const* discLines, unsigned int const* numLines, unsigned int precision, MeshType meshT)
 {
-	file << "VECTORS " << name << " float " << endl;
+	file << "VECTORS " << name << " " << __VTK_DATA_TYPE__ << endl;
 
 	if (g_settings.NativeFieldDumps())
 		meshT = CARTESIAN_MESH; //dump field components as they are...
@@ -407,7 +407,7 @@ bool ProcessFields::DumpMultiVectorArray2VTK(ofstream &file, string names[], FDT
 
 void ProcessFields::WriteVTKScalarArray(ofstream &file, string name, FDTD_FLOAT const* const* const* array, unsigned int const* numLines, unsigned int precision)
 {
-	file << "SCALARS " << name << " float " << 1 << endl;
+	file << "SCALARS " << name << " " << __VTK_DATA_TYPE__ << 1 << endl;
 	file << "LOOKUP_TABLE default" << endl;
 	unsigned int pos[3];
 	int count=0;
diff --git a/FDTD/processfields.h b/FDTD/processfields.h
index 3ac523e..58b6dac 100644
--- a/FDTD/processfields.h
+++ b/FDTD/processfields.h
@@ -21,6 +21,8 @@
 #include "processing.h"
 #include "tools/array_ops.h"
 
+#define __VTK_DATA_TYPE__ "double"
+
 class ProcessFields : public Processing
 {
 public: