On 07/22/2017 04:42 AM, Giovanni Di Ilio wrote:

Yes, I agree, the issue is definitely related to round-off error. But I still don't understand why the mapping fails, if the vertex coordinates read in from my .msh file are apparently stored as doubles... Anyway, yes, the problem is easily overcome by "manipulating" the values before mapping.

The numbers you consider differ in the last two or three bits when stored as double precision. Any computations you do with them, for example in the mapping, will then necessarily be completely off.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to