Skip to content

Peano: Geometrical comparisons must always use relative tolerance

The tarch::la:: routines use usually an absolute tolerance (default: machine precision) for comparing floating point numbers. This is problematic for numbers with large absolute values.

Fix: These numbers must be normalised and scaled to range [-1,1] before comparing with respect to the tolerance.