dgebal
dgebal - balance a general real matrix A
SUBROUTINE DGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO)
CHARACTER * 1 JOB
INTEGER N, LDA, ILO, IHI, INFO
DOUBLE PRECISION A(LDA,*), SCALE(*)
SUBROUTINE DGEBAL_64( JOB, N, A, LDA, ILO, IHI, SCALE, INFO)
CHARACTER * 1 JOB
INTEGER*8 N, LDA, ILO, IHI, INFO
DOUBLE PRECISION A(LDA,*), SCALE(*)
SUBROUTINE GEBAL( JOB, [N], A, [LDA], ILO, IHI, SCALE, [INFO])
CHARACTER(LEN=1) :: JOB
INTEGER :: N, LDA, ILO, IHI, INFO
REAL(8), DIMENSION(:) :: SCALE
REAL(8), DIMENSION(:,:) :: A
SUBROUTINE GEBAL_64( JOB, [N], A, [LDA], ILO, IHI, SCALE, [INFO])
CHARACTER(LEN=1) :: JOB
INTEGER(8) :: N, LDA, ILO, IHI, INFO
REAL(8), DIMENSION(:) :: SCALE
REAL(8), DIMENSION(:,:) :: A
#include <sunperf.h>
void dgebal(char job, int n, double *a, int lda, int *ilo, int *ihi, double *scale, int *info);
void dgebal_64(char job, long n, double *a, long lda, long *ilo, long *ihi, double *scale, long *info);
dgebal balances a general real matrix A. This involves, first,
permuting A by a similarity transformation to isolate eigenvalues
in the first 1 to ILO-1 and last IHI+1 to N elements on the
diagonal; and second, applying a diagonal similarity transformation
to rows and columns ILO to IHI to make the rows and columns as
close in norm as possible. Both steps are optional.
Balancing may reduce the 1-norm of the matrix, and improve the
accuracy of the computed eigenvalues and/or eigenvectors.
-
* JOB (input)
-
Specifies the operations to be performed on A:
-
* N (input)
-
The order of the matrix A. N >= 0.
-
* A (input/output)
-
On entry, the input matrix A.
On exit, A is overwritten by the balanced matrix.
If JOB = 'N', A is not referenced.
See Further Details.
-
* LDA (input)
-
The leading dimension of the array A. LDA >= max(1,N).
-
* ILO (output)
-
ILO and IHI are set to integers such that on exit
A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N.
If JOB = 'N' or 'S', ILO = 1 and IHI = N.
-
* IHI (output)
-
See the description for ILO.
-
* SCALE (output)
-
Details of the permutations and scaling factors applied to
A. If P(j) is the index of the row and column interchanged
with row and column j and D(j) is the scaling factor
applied to row and column j, then
SCALE(j) = P(j) for j = 1,...,ILO-1
-
* INFO (output)
-