Re: Cholesky factorization

>>>>> "rif" == rif  <[email protected]> writes:

    rif> I assume it would be a lot of work to expose the ability to do work in
    rif> single-precision floating point rather than double?

I suppose with some hacking of the macros we could make it work,
assuming that single-precision routines always want single-precision
arrays and numbers.

This would be a bit tedious, and I'm reluctant to do this, however.  I
decided long ago that double-precision almost always allowed me to
ignore round-off issues for the things I do.  The extra time and
memory were not an issue.  Are you on an x86?  Then single-precision
only buys you memory.

Ray