Re: [Lapackpp-devel] tridiagonal

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

Am Sonntag, 6. August 2006 13:25 schrieb Do bi:
> If I use LaTridiagMatDouble instead of LaGenMatDouble, will I save memory?

RTFM!!!!!!!!!!!! http://en.wikipedia.org/wiki/RTFM
Everything explained in the header include/trmd.h

> Moreover, if I use LaSymmTridiagMatDouble, will the new Lapackpp package
> solve the system?

There is no LaSymmTridiagMatDouble.

Christian

>
> Christian Stimming <sti...@tu...> wrote:  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 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.