Re: [Lapackpp-devel] tridiagonal

[Christian S. <sti...@tu...>]

Am Freitag, 4. August 2006 15:24 schrieb Do bi:
> I want to solve the tridiagonal Ax=b. If this is not possible, then I
> proceed with LaMatDouble because I have already produced a tridiagonal with
> my code in the form of a LaMatDouble. Thanks.
>
> Christian Stimming <sti...@tu...> wrote:
> (...) Apart from this I actually don't quite understand what you
> want to do with the LaTridiagMat. There aren't any useful functions
> available for these, except for LaTridiagMatFactorize in trfd.h ...

If you want to solve tridiagonal Ax=b then the two functions at the end of 
include/trfd.h should be exactly what you need. You first create the 
LaTridiagMatDouble A with your data, then create its factorization by  
LaTridiagMatFactorize, then use that factorization to solve Ax=b by the 
LaLinearSolve in trfd.h. In numerical maths it is quite usual to solve Ax=b 
in these two separate steps.

Christian