[Lapackpp-devel] LaBandFactDouble Question

[Christophe P. <pi...@tl...>]

Hello,

I am trying to use the banded matrices to solve pde, and I notice that the
return arguments in the function LaBandMatFactorize are weird. I mean that
the pivot and upper/lower matrix should be save in the member B and not in
the original matrix:

So my question is in bfd.h :

00138 inline void LaBandMatFactorize(LaBandMatDouble &A, LaBandFactDouble
&AF)
00139 {
00140     integer n = A.size(1), m = n, LDA = A.gdim(0);
00141     integer KL = A.subdiags(), KU = A.superdiags(), info=0;
00142
00143     F77NAME(dgbtrf)(&m, &n, &KL, &KU, &A(0,0), &LDA,
&(AF.pivot()(0)), &info);

should it be replace with something more along this (even maybe overload
the function)

00138 inline void LaBandMatFactorize(LaBandMatDouble &A, LaBandFactDouble
&AF)
00139 {
00140     integer n = A.size(1), m = n, LDA = A.gdim(0);
00141     integer KL = A.subdiags(), KU = A.superdiags(), info=0;
00142
00143     F77NAME(dgbtrf)(&m, &n, &KL, &KU, &AF.B(0,0), &LDA,
&(AF.pivot()(0)), &info);

otherwise the factorization is not saved into AF but in A, making AF.B
useless, and the function LaLinearSolve unusable.
This is just a question/remark.
Thanks
Christophe