On Fri, Oct 25, 2002 at 05:12:56PM -0700, John Hawkes wrote:
> From: "Siva Koti Reddy" <siva.kotireddy@wipro.com>
> > 39 new_raph Unable to solve equation in 100 tries. P =
> 1.5708, P0
> > = 1.5708, delta = 6.12574e-17
> > new_raph: Success
> > *** Failed to execute new_raph ***
>
> The AIM7/AIM9 new_raph is broken code. The convergence loop termination
> conditional looks something like:
> if (delta == 0) break;
> for a type "double" delta. You ought to change that to be something
> like:
> if (delta <= 0.00000001L) break;
I usually specify the compiler flag -ffloat-store and that fixes the
issue for me.
-- Nate Straz nstraz@sgi.com sgi, inc http://www.sgi.com/ Linux Test Project http://ltp.sf.net/ - To unsubscribe from this list: send the line "unsubscribe linux-kernel" in the body of a message to majordomo@vger.kernel.org More majordomo info at http://vger.kernel.org/majordomo-info.html Please read the FAQ at http://www.tux.org/lkml/
This archive was generated by hypermail 2b29 : Thu Oct 31 2002 - 22:00:39 EST