------- Comment #4 from macracan at yahoo dot com 2007-10-01 17:43 ------- (In reply to comment #3) > For once a real floating point bug. In this code > double p = 0.422244 * f[a.e]; > if (f[0] < p) > > with a.e=1, f[1]=0.285433, we should calculate p=0.285433. Since > f[0]=0.0461109, we shouldn't enter the if-clause, but we do when
Sorry, but this isn't quite accurate; the math is off, and the code _should_ enter the if branch the first time around, just not the second time around. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33611