Question regarding PDGBTRF

Open discussion regarding features, bugs, issues, vendors, etc.

Question regarding PDGBTRF

Postby ashtonaut » Thu Jan 27, 2005 6:59 pm

I have a banded matrix distributed across 6 processes in the block-column fashion required. When I call PDGBTRF, I get an INFO=1 error code on all 6 processes.

From Scalapack guide:
If INFO = K<=NPROCS, the submatrix stored on processor INFO and factored locally was not nonsingular, and the factorization was not completed


Is it a requirement of this routine that all local submatrices can be factorised individually? If it is, this may explain the error, as I think the first column block is not diagonally dominant, even though the global matrix as a whole will be.

Any help would be most appreciated.

Ashton Peters
ashtonaut
 
Posts: 13
Joined: Thu Jan 27, 2005 6:53 pm

Postby Stan Tomov » Mon Jan 31, 2005 3:37 pm

Yes, it is a requirement that the blocks can be factored individually.
For example, the matrix
/ 1 1 | 0 0 \
| 1 1 | 1 0 |
| 0 1 | 2 0 |
\ 0 0 | 0 2 /
is non-singular, but if solved on 2 processors with
splitting as indicated, the method will fail since it
will try to factor (as a first step of the algorithm)
the singular block A(1:2,1:2).
The method implemented is described in
http://www.netlib.org/lapack/lawns/lawn125.ps
It should work for diagonally dominant nonsymmetric banded
matrices.
Stan Tomov
 
Posts: 13
Joined: Thu Dec 09, 2004 1:28 pm


Return to User Discussion

Who is online

Users browsing this forum: Yahoo [Bot] and 4 guests