[Scilab-Dev] In intdgesv3.f : there is something I don't understand

[bernard.hugueney at scilab.org]

Hi,

On Fri, 5 Jun 2009 23:27:07 +0200, Vincent Lejeune
<vincent.lejeune at institutoptique.fr> wrote:

> of the previous LU decomposition. The matter is, that DGETRS is called on
> the 
> lB pointer : according to the documentation of the DGETRS function, the
lB 
> matrix should be modified, that means that the B matrix is not the same
> after 
> the computation.
> 
> A quick run into Scilab prooves me that this assertion is wrong, B stay 
> unchanged before and after using A\B.
> 
> Which step miss I ?
> 

Sad, but true : getrhsvar() does a *copy* of the argument matrix. So dgetrs
modifies the copy, not the original data.

Regards,

Bernard