[Getfem-users] Glitch with VTK export

[Roman Putanowicz]

Hi All,

I have encountered problem with VTK export. 
In getfem_export.h  the type of exported data
is set to "float" (e.g. line 275).

Apparently that causes problems when the data
is saved with with more than 6 digits after
the decimal dot (this is just guess - I went to
this conclusion by manually truncating the digits)

Well, the fix is simple -- declare data type as double.

Now there is a question -- should we just change
"float" to "double" in vtk_export::write_dataset_
or shall we go for more sophisticated solution
of checking what is the actual type of elements of
written vector and then adjust the VTK data type ?

The fist solution is trivial the second might be a bit complex.
There is also another solution to add new property
to vtk_export to explicitly choose the VTK data type.

What is your opinion on this?

In the attachment there is the VTK file exported form GetFEM
that causes the problems. I am using ParaView 3.6.2.
and the SVN getfem version.

Regards,

Roman
-- 
Roman Putanowicz, PhD  < address@hidden  >
Institute for Computational Civil Engng (L-5)
Dept. of Civil Engng, Cracow Univ. of Technology
www.l5.pk.edu.pl, tel. +48 12 628 2569, fax 2034

annulus.vtk
Description: Text document