Karl Stiefvater - 2011-06-21

can be fixed by adding "+ C_EPSILON" to computation above.