Re: Cholesky factorization

Nicolas Neuss <[email protected]> 30 Oct 2002 10:47:00 +0100

Newsgroups	gmane.lisp.matlisp.user
Message-ID	<[email protected]>

"Michael A. Koerber" <[email protected]> writes:

> Rif,
> 
>   I should read and write more slowly.  WRT your first posting,
> the solution for multiple RHS would be to use (GETRS ...) or (GETRS! ...)
> 
> mike

Yes, something like

(multiple-value-bind (lu ipiv info)
    (getrf! (copy mat))
  (declare (ignore info))
  (getrs! lu ipiv rhs)) 

should work in recent CVS versions.

Nicolas.