Constructor for LU: why dgesv?

[Thomas Hoffmann]

Excuse my ignorance, but when looking into the ctor for LU (e.g. dbleLU.cc),
I noticed that octave calls the Fortran subroutine dgesv, with a second 
parameter
(NRHS) of 0. 

dgesv calls in turn DGETRF (factorization) and DGETRS (solution). 

But if one uses NRHS=0, DGETRS only checks some parameters and returns.

Why does the ctor not call dgetrf directly, when only the factorization is 
necessary?

Regards,
Thomas Hoffmann.


=============================================================================
Thomas Hoffmann, Institut fuer Halbleiter- und Mikrosystemtechnik, TU Dresden 
E-mail: address@hidden



-----------------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.che.wisc.edu/octave/octave.html
How to fund new projects:  http://www.che.wisc.edu/octave/funding.html
Subscription information:  http://www.che.wisc.edu/octave/archive.html
-----------------------------------------------------------------------