I&#39;m trying to visualize some scalar field data using VTK (Paraview, and some hand-rolled scripts), but I&#39;m having trouble outputting my data in VTK format.  I tried using the VTK Legacy ASCII format with the code below.  It works correctly when my spacing is 0.1 or greater, but when I decrease the spacing (I&#39;ve tried .08, .05, .01 and lower), The isosurfaces VTK generate are incorrect.  It looks as if the dimensions of the data are miscalculated, because I get bumpy, off-axis sheets offset from the true isosurface as shown here: <a href="http://www.cs.utah.edu/~benjones/isoFail.png">http://www.cs.utah.edu/~benjones/isoFail.png</a> .  The true isosurface is a triangel mesh created by contouring the same data I&#39;m dumping below.  Why might the code/sample data below be giving me these errors?<div>

<br></div><div>Output code:</div><div><br></div><div><div>   std::ofstream out(filename);</div><div>  </div><div>  out &lt;&lt; &quot;# vtk DataFile Version 3.0&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;Signed distance/biharmonic visualizer&quot; &lt;&lt; std::endl;</div>

<div>  out &lt;&lt; &quot;ASCII&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;DATASET STRUCTURED_POINTS&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;DIMENSIONS &quot; &lt;&lt; nx &lt;&lt; &quot; &quot; &lt;&lt; ny &lt;&lt; &quot; &quot; &lt;&lt; nz &lt;&lt; std::endl;</div>

<div>  out &lt;&lt; &quot;ORIGIN 0.0 0.0 0.0&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;SPACING &quot; &lt;&lt; h &lt;&lt; &quot; &quot; &lt;&lt; h &lt;&lt; &quot; &quot; &lt;&lt; h &lt;&lt; std::endl;</div>

<div>  out &lt;&lt; &quot;POINT_DATA &quot; &lt;&lt; nx*ny*nz &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;SCALARS signedDistance double 1&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;LOOKUP_TABLE default&quot; &lt;&lt; std::endl;</div>

<div>  for(size_t i = 0; i &lt; nx; ++i)</div><div>    for(size_t j = 0; j &lt; ny; ++j)</div><div>      for(size_t k = 0; k &lt; nz; ++k)</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>out &lt;&lt; tempPhi(i,j,k) &lt;&lt; std::endl;</div>

<div>  //out &lt;&lt; &quot;POINT_DATA &quot; &lt;&lt; nx*ny*nz &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;SCALARS biharmonic double 1&quot; &lt;&lt; std::endl;</div><div>  out &lt;&lt; &quot;LOOKUP_TABLE default&quot; &lt;&lt; std::endl;</div>

<div>  for(size_t i = 0; i &lt; nx; ++i)</div><div>    for(size_t j = 0; j &lt; ny; ++j)</div><div>      for(size_t k = 0; k &lt; nz; ++k)</div><div>      out &lt;&lt; biharmonic(i,j,k) &lt;&lt; std::endl;</div><div>  out.close();</div>

</div><div><br></div><div><br></div><div><br></div><div>Sample file:</div><div><br></div><div><br></div><div><div># vtk DataFile Version 3.0</div><div>Signed distance/biharmonic visualizer</div><div>ASCII</div><div>DATASET STRUCTURED_POINTS</div>

<div>DIMENSIONS 54 58 52</div><div>ORIGIN 0.0 0.0 0.0</div><div>SPACING 0.05 0.05 0.05</div><div>POINT_DATA 162864</div><div>SCALARS signedDistance double</div><div>LOOKUP_TABLE default</div><div>1.79769e+308</div><div>1.79769e+308</div>

<div>1.79769e+308</div></div><div>...</div><div>...</div><div><div>1.79769e+308</div><div>1.79769e+308</div><div>SCALARS biharmonic double</div><div>LOOKUP_TABLE default</div><div>0</div><div>0</div></div><div>...</div><div>

...</div><div>0</div><div>0</div><div><br></div>