Re: LAPACK, et al., vs. Numerical Recipes



James Amundson <amundson@xxxxxxxxxxxxxxxxxxxxx> wrote:

Hi,

I am looking for comparisons (i.e., with regards to robustness,
accuracy and speed) between the linear algebra routines from Numerical
Recipes and LAPACK, preferably including ATLAS and/or Goto BLAS.

If NR uses the reference blas (just writing out the loops) then anything
using Atlas or Goto Blas will beat the pants off it in terms of speed.

In terms of accuracy, the solver routines in Lapack are probably better
in error checking, but not necessarily much better. On the other hand,
the eigenvalue routines in Lapack are way more up to date.

Sorry, I know of no comparative studies.

Victor.
--
Victor Eijkhout -- eijkhout at tacc utexas edu
.



Relevant Pages

  • Re: Matrix Multiplication in BLAS or Lapack
    ... > so you need not to link the lapack library. ... > LAPACK library includes the blas routines already due to ... > BLAS offers the generic and simple routines, e.g. daxpy, dgemv, dgemm, ... > playing around with the optimization options of the compiler - ...
    (sci.math.num-analysis)
  • Re: sin(10^300), ignorance, collusion. was Re: TILE64 embedded multicore processors - Guy Macon
    ... And LAPACK gives reliable answers so long as the arithmetic isn't ... count routines for bisection, ... slower fall-backs for crap arithmetics. ... And there may be some optimizations ...
    (comp.arch)
  • Re: Matrix Multiplication in BLAS or Lapack
    ... BLAS-library and not of LAPACK; ... LAPACK library includes the blas routines already due to ... BLAS offers the generic and simple routines, e.g. daxpy, dgemv, dgemm, ... playing around with the optimization options of the compiler - ...
    (sci.math.num-analysis)
  • Re: Lapack ?
    ... Martin Raspaud writes: ... > Does anyone know if there is a 100% Common lisp implementation of Lapack ... > routines for matrix inversion? ...
    (comp.lang.lisp)