<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd"><html><head><meta name="qrichtext" content="1" /><style type="text/css">p, li { white-space: pre-wrap; }</style></head><body style=" font-family:'DejaVu Sans'; font-size:9pt; font-weight:400; font-style:normal;">Hi,<br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><br></p>I'm trying to reproduce the behaviour of the intdgesv3 function (that is, the "backslash" function between 2 double matrix).<br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><br></p>After several check of input argument, if the, let's say A matrix in the A\B expression, is a square matrix, the function copy the A matrix in another matrix AF, call DGETRF that compute the LU decomposition of the AL matrix inplace, and then call DGETRS, which resolve the A*X=B system, with the help 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.<br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><br></p>A quick run into Scilab prooves me that this assertion is wrong, B stay unchanged before and after using A\B.<br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><br></p>Which step miss I ?<br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><br></p>Thanks, Vincent.</p></body></html>