all computations to be done in extended precision?
Subject: all computations to be done in extended precision?
From: Richard Fateman
Date: Tue, 18 Dec 2007 12:00:57 -0800
how about
f(x,y):=block([fpprec:50], 0.0b0+x^2+y^2);
But this may not be what you want, since it computes 2.1^2 IN SINGLE
PRECISION, and then adds it up to 50 digits etc.
Probably you want
fpprec:50$
f(x,y):=x^2+y^2;
f(2.1b0, 3.2b0); which gives 1.465b1, not some bunch of garbage.
> -----Original Message-----
> From: maxima-bounces at math.utexas.edu
> [mailto:maxima-bounces at math.utexas.edu] On Behalf Of S. Newhouse
> Sent: Tuesday, December 18, 2007 11:53 AM
> To: maxima at math.utexas.edu
> Subject: all computations to be done in extended precision?
>
> Hello,
> Is there a declaration at the beginning of a program to have all
> calculations done in extended precision?
>
> It's a little cumbersome to do, e.g.
> %i39) fpprec: 50;
> Evaluation took 0.00 seconds (0.00 elapsed) using 0 bytes.
> (%o39) 50
> (%i40) f(x,y):= bfloat(x)^2 + bfloat(y)^2;
> Evaluation took 0.00 seconds (0.00 elapsed) using 0 bytes.
> 2 2
> (%o40) f(x, y) := bfloat (x) + bfloat (y)
> (%i41) f(2.1,3.2);
> Evaluation took 0.00 seconds (0.00 elapsed) using 0 bytes.
> (%o41) 1.4650000000000001509903313490212934419184329654106b1
>
> TIA,
> -sen
>
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>