laqsy#

Functions

void slaqsy(
    const char*          uplo,
    const INT            n,
          f32*  restrict A,
    const INT            lda,
    const f32*  restrict S,
    const f32            scond,
    const f32            amax,
          char*          equed
);
void slaqsy(const char *uplo, const INT n, f32 *restrict A, const INT lda, const f32 *restrict S, const f32 scond, const f32 amax, char *equed)#

SLAQSY equilibrates a symmetric matrix A using the scaling factors in the vector S.

Parameters

in
uplo

  • 'U': Upper triangular

  • 'L': Lower triangular

in
n

The order of the matrix A. n>=0.

inout
A

Array of dimension (lda, n). On entry, the symmetric matrix A. If uplo='U', the leading n by n upper triangular part of A contains the upper triangular part of the matrix A, and the strictly lower triangular part of A is not referenced. If uplo='L', the leading n by n lower triangular part of A contains the lower triangular part of the matrix A, and the strictly upper triangular part of A is not referenced. On exit, if equed='Y', the equilibrated matrix: diag(S) * A * diag(S).

in
lda

The leading dimension of the array A. lda>=max(n,1).

in
S

Array of dimension n. The scale factors for A.

in
scond

Ratio of the smallest S(i) to the largest S(i).

in
amax

Absolute value of largest matrix entry.

out
equed

  • 'N': No equilibration

  • 'Y': Equilibration was done, i.e., A has been replaced by diag(S) * A * diag(S)

Functions

void dlaqsy(
    const char*          uplo,
    const INT            n,
          f64*  restrict A,
    const INT            lda,
    const f64*  restrict S,
    const f64            scond,
    const f64            amax,
          char*          equed
);
void dlaqsy(const char *uplo, const INT n, f64 *restrict A, const INT lda, const f64 *restrict S, const f64 scond, const f64 amax, char *equed)#

DLAQSY equilibrates a symmetric matrix A using the scaling factors in the vector S.

Parameters

in
uplo

  • 'U': Upper triangular

  • 'L': Lower triangular

in
n

The order of the matrix A. n>=0.

inout
A

Array of dimension (lda, n). On entry, the symmetric matrix A. If uplo='U', the leading n by n upper triangular part of A contains the upper triangular part of the matrix A, and the strictly lower triangular part of A is not referenced. If uplo='L', the leading n by n lower triangular part of A contains the lower triangular part of the matrix A, and the strictly upper triangular part of A is not referenced. On exit, if equed='Y', the equilibrated matrix: diag(S) * A * diag(S).

in
lda

The leading dimension of the array A. lda>=max(n,1).

in
S

Array of dimension n. The scale factors for A.

in
scond

Ratio of the smallest S(i) to the largest S(i).

in
amax

Absolute value of largest matrix entry.

out
equed

  • 'N': No equilibration

  • 'Y': Equilibration was done, i.e., A has been replaced by diag(S) * A * diag(S)

Functions

void claqsy(
    const char*          uplo,
    const INT            n,
          c64*  restrict A,
    const INT            lda,
    const f32*  restrict S,
    const f32            scond,
    const f32            amax,
          char*          equed
);
void claqsy(const char *uplo, const INT n, c64 *restrict A, const INT lda, const f32 *restrict S, const f32 scond, const f32 amax, char *equed)#

CLAQSY equilibrates a symmetric matrix A using the scaling factors in the vector S.

Parameters

in
uplo

  • 'U': Upper triangular

  • 'L': Lower triangular

in
n

The order of the matrix A. n>=0.

inout
A

Array of dimension (lda, n). On entry, the symmetric matrix A. If uplo='U', the leading n by n upper triangular part of A contains the upper triangular part of the matrix A, and the strictly lower triangular part of A is not referenced. If uplo='L', the leading n by n lower triangular part of A contains the lower triangular part of the matrix A, and the strictly upper triangular part of A is not referenced. On exit, if equed='Y', the equilibrated matrix: diag(S) * A * diag(S).

in
lda

The leading dimension of the array A. lda>=max(n,1).

in
S

Array of dimension n. The scale factors for A.

in
scond

Ratio of the smallest S(i) to the largest S(i).

in
amax

Absolute value of largest matrix entry.

out
equed

  • 'N': No equilibration

  • 'Y': Equilibration was done, i.e., A has been replaced by diag(S) * A * diag(S)

Functions

void zlaqsy(
    const char*          uplo,
    const INT            n,
          c128* restrict A,
    const INT            lda,
    const f64*  restrict S,
    const f64            scond,
    const f64            amax,
          char*          equed
);
void zlaqsy(const char *uplo, const INT n, c128 *restrict A, const INT lda, const f64 *restrict S, const f64 scond, const f64 amax, char *equed)#

ZLAQSY equilibrates a symmetric matrix A using the scaling factors in the vector S.

Parameters

in
uplo

  • 'U': Upper triangular

  • 'L': Lower triangular

in
n

The order of the matrix A. n>=0.

inout
A

Array of dimension (lda, n). On entry, the symmetric matrix A. If uplo='U', the leading n by n upper triangular part of A contains the upper triangular part of the matrix A, and the strictly lower triangular part of A is not referenced. If uplo='L', the leading n by n lower triangular part of A contains the lower triangular part of the matrix A, and the strictly upper triangular part of A is not referenced. On exit, if equed='Y', the equilibrated matrix: diag(S) * A * diag(S).

in
lda

The leading dimension of the array A. lda>=max(n,1).

in
S

Array of dimension n. The scale factors for A.

in
scond

Ratio of the smallest S(i) to the largest S(i).

in
amax

Absolute value of largest matrix entry.

out
equed

  • 'N': No equilibration

  • 'Y': Equilibration was done, i.e., A has been replaced by diag(S) * A * diag(S)