Re: Maple's linear algebra routine and higher precision arithmetic

On Mon, 29 May 2006 13:12:03 -0700, Ronald Bruck wrote:

In article <pan.2006.05.29.18.47.40.640507@xxxxxxxxxx>, Ronald Haynes
<ronald.haynes@xxxxxxxxxx> wrote:

Hi, many thanks... I just want to test the validity of some lower
precision results with results generated by increasing levels of precision.

Would defining the matrix entries as -1.0 4.0 1.0 do the trick to force
floats?

In Maple--yes. If you do this in Mathematica you're well-advised to
use a notation like N[-1,64], N[4, 64], N[1,64]. The "4.0" will be
read as a standard double-precision floating-point number.
(Mathematica doesn't have a "Digits" command--each number has its own
intrinsic precision.)

But you asked about Maple, not Mathematica, so this is just for your
general edification...

Hi, upon simple testing this does not seem to work. i.e. setting
Digits:=1 and computing the inverse does not give different results from
setting Digits:=100 and computing the inverse.

Sometime ago I recall hearing that Maple was changing that way it dealt
with linear algebra routines --- maybe by using NAG routines (or similar)
under the hood? If Maple is actually using compiled routines how could
Digits:= whatever actually have an effect, I would presume the routines
would be double precision by default and not be affected by a change of
Digits (except possibly by rounding the input matrix). In my case the
matrix is integer valued. In fact only a slight modification of the
tridiag(-1,4,-1) matrix.

Any further thoughts?
Thanks,
R Haynes

.

Relevant Pages

  • Re: Ambiguous interfaces
    ... > I have a question about ambiguous interfaces. ... > subroutine DoSomething_DefaultInteger ... if you combine two or more specific routines ... default real precision and at double precision. ...
    (comp.lang.fortran)
  • Re: qfloat extended precision in Cephes library
    ... for Mathematical Functions (Ellis Horwood Series in Mathematics and Its ... precision offering in particular. ... qfloat.h C++ wrapper makes it easier to write math code that looks like ... Some of Mr. Moshier's built-in routines, ...
    (sci.math.num-analysis)
  • Re: high-precision eigenvalue solver
    ... >> precision. ... Rather than hand-coding the routines, I would use a package like Gnu ... and recode the lapack routines to call the ... Since it's quadratically convergent, this isn't as ...
    (sci.math.num-analysis)
  • Re: Q: how to compute modulo Pi for large arguments?
    ... >>want to use it for periodics of large arguments in a C pgm. ... > But that's only got ~16 digits of precision at the best of times. ... much better (GSL can partially be interfaced with Maple), ...
    (sci.math.num-analysis)
  • Re: Maples linear algebra routine and higher precision arithmetic
    ... precision results with results generated by increasing levels of precision. ... with linear algebra routines --- maybe by using NAG routines ... Digits. ... meaning that hardware is used if Digits <= evalhf(i.e. ...
    (sci.math.symbolic)