Re: Using DTRTRS (solve triangular system) from Octave

[c.]

On 28 Nov 2013, at 15:51, c. <address@hidden> wrote:

> 
> On 28 Nov 2013, at 15:43, José Luis García Pallero <address@hidden> wrote:
> 
>> Hello:
>> 
>> I need to solve a triangular system in Octave. Exists a way in Octave
>> to use the DTRTRS Lapack routine, which is the one to solve a
>> triangular system?
>> 
>> Cheers
> 
> IIRC, Yes. 
> 
> If you use mldivide to solve a triangular system Octave should automatically
> choose DTRTRS to solve it, i.e.:
> 
> A = tril (randn (10, 10));
> b = ones (10, 1);
> 
> x = A\b;
> 
> Should solve the system using DTRTRS.
> 
> If you want to check more precisely what is actually done, you can have a 
> look 
> at the source code of the file dMatrix.cc.
> 
> HTH,
> c.

Here's a more direct link to the relevant line in the source code:

http://hg.savannah.gnu.org/hgweb/octave/file/2ef2043a0601/liboctave/array/dMatrix.cc#l1727

c.