[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## Re: Using DTRTRS (solve triangular system) from Octave

**From**: |
c. |

**Subject**: |
Re: Using DTRTRS (solve triangular system) from Octave |

**Date**: |
Thu, 28 Nov 2013 15:59:19 +0100 |

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.