[Lapackpp-devel] LaBandFactDouble Question
Status: Beta
Brought to you by:
cstim
|
From: Christophe P. <pi...@tl...> - 2005-04-14 18:18:12
|
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
|