latps#

slatps

Functions

void slatps(
    const char*          uplo,
    const char*          trans,
    const char*          diag,
    const char*          normin,
    const INT            n,
    const f32*  restrict AP,
          f32*  restrict X,
          f32*           scale,
          f32*  restrict cnorm,
          INT*           info
);

void slatps(const char *uplo, const char *trans, const char *diag, const char *normin, const INT n, const f32 *restrict AP, f32 *restrict X, f32 *scale, f32 *restrict cnorm, INT *info)#

SLATPS solves one of the triangular systems A * x = s*b or A**T * x = s*b with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form.

Here A**T denotes the transpose of A, x and b are n-element vectors, and s is a scaling factor, usually less than or equal to 1, chosen so that the components of x will be less than the overflow threshold. If the unscaled problem will not cause overflow, the Level 2 BLAS routine DTPSV is called. If the matrix A is singular (A(j,j) = 0 for some j), then s is set to 0 and a non-trivial solution to A*x = 0 is returned.

Parameters

in

uplo

‘U’: A is upper triangular; ‘L’: A is lower triangular.

in

trans

‘N’: Solve A * x = s*b; ‘T’/’C’: Solve A**T * x = s*b.

in

diag

‘N’: A is non-unit triangular; ‘U’: A is unit triangular.

in

normin

‘Y’: CNORM contains column norms on entry; ‘N’: CNORM is not set, will be computed.

in

n

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

in

AP

The triangular matrix A, packed columnwise. Array of dimension (n*(n+1)/2).

inout

X

On entry, the right hand side b. On exit, overwritten by the solution x. Array of dimension n.

out

scale

The scaling factor s for the triangular system.

inout

cnorm

If normin=’Y’, cnorm contains column norms on entry. If normin=’N’, cnorm returns the 1-norm of offdiagonal part of column j. Array of dimension n.

out

info

Exit status:

= 0: successful exit
< 0: if info = -k, the k-th argument had an illegal value

dlatps

Functions

void dlatps(
    const char*          uplo,
    const char*          trans,
    const char*          diag,
    const char*          normin,
    const INT            n,
    const f64*  restrict AP,
          f64*  restrict X,
          f64*           scale,
          f64*  restrict cnorm,
          INT*           info
);

void dlatps(const char *uplo, const char *trans, const char *diag, const char *normin, const INT n, const f64 *restrict AP, f64 *restrict X, f64 *scale, f64 *restrict cnorm, INT *info)#

DLATPS solves one of the triangular systems A * x = s*b or A**T * x = s*b with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form.

Here A**T denotes the transpose of A, x and b are n-element vectors, and s is a scaling factor, usually less than or equal to 1, chosen so that the components of x will be less than the overflow threshold. If the unscaled problem will not cause overflow, the Level 2 BLAS routine DTPSV is called. If the matrix A is singular (A(j,j) = 0 for some j), then s is set to 0 and a non-trivial solution to A*x = 0 is returned.

Parameters

in

uplo

‘U’: A is upper triangular; ‘L’: A is lower triangular.

in

trans

‘N’: Solve A * x = s*b; ‘T’/’C’: Solve A**T * x = s*b.

in

diag

‘N’: A is non-unit triangular; ‘U’: A is unit triangular.

in

normin

‘Y’: CNORM contains column norms on entry; ‘N’: CNORM is not set, will be computed.

in

n

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

in

AP

The triangular matrix A, packed columnwise. Array of dimension (n*(n+1)/2).

inout

X

On entry, the right hand side b. On exit, overwritten by the solution x. Array of dimension n.

out

scale

The scaling factor s for the triangular system.

inout

cnorm

If normin=’Y’, cnorm contains column norms on entry. If normin=’N’, cnorm returns the 1-norm of offdiagonal part of column j. Array of dimension n.

out

info

Exit status:

= 0: successful exit
< 0: if info = -k, the k-th argument had an illegal value

clatps

Functions

void clatps(
    const char*          uplo,
    const char*          trans,
    const char*          diag,
    const char*          normin,
    const INT            n,
    const c64*  restrict AP,
          c64*  restrict X,
          f32*           scale,
          f32*  restrict cnorm,
          INT*           info
);

void clatps(const char *uplo, const char *trans, const char *diag, const char *normin, const INT n, const c64 *restrict AP, c64 *restrict X, f32 *scale, f32 *restrict cnorm, INT *info)#

CLATPS solves one of the triangular systems A * x = s*b, A**T * x = s*b, or A**H * x = s*b with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form.

Here A**T denotes the transpose of A, A**H denotes the conjugate transpose of A, x and b are n-element vectors, and s is a scaling factor, usually less than or equal to 1, chosen so that the components of x will be less than the overflow threshold. If the unscaled problem will not cause overflow, the Level 2 BLAS routine ZTPSV is called. If the matrix A is singular (A(j,j) = 0 for some j), then s is set to 0 and a non-trivial solution to A*x = 0 is returned.

Parameters

in

uplo

‘U’: A is upper triangular; ‘L’: A is lower triangular.

in

trans

‘N’: Solve A * x = s*b; ‘T’: Solve A**T * x = s*b; ‘C’: Solve A**H * x = s*b.

in

diag

‘N’: A is non-unit triangular; ‘U’: A is unit triangular.

in

normin

‘Y’: CNORM contains column norms on entry; ‘N’: CNORM is not set, will be computed.

in

n

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

in

AP

The triangular matrix A, packed columnwise. Array of dimension (n*(n+1)/2).

inout

X

On entry, the right hand side b. On exit, overwritten by the solution x. Array of dimension n.

out

scale

The scaling factor s for the triangular system.

inout

cnorm

If normin=’Y’, cnorm contains column norms on entry. If normin=’N’, cnorm returns the 1-norm of offdiagonal part of column j. Array of dimension n.

out

info

= 0: successful exit < 0: if info = -k, the k-th argument had an illegal value

zlatps

Functions

void zlatps(
    const char*          uplo,
    const char*          trans,
    const char*          diag,
    const char*          normin,
    const INT            n,
    const c128* restrict AP,
          c128* restrict X,
          f64*           scale,
          f64*  restrict cnorm,
          INT*           info
);

void zlatps(const char *uplo, const char *trans, const char *diag, const char *normin, const INT n, const c128 *restrict AP, c128 *restrict X, f64 *scale, f64 *restrict cnorm, INT *info)#

ZLATPS solves one of the triangular systems A * x = s*b, A**T * x = s*b, or A**H * x = s*b with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form.

Here A**T denotes the transpose of A, A**H denotes the conjugate transpose of A, x and b are n-element vectors, and s is a scaling factor, usually less than or equal to 1, chosen so that the components of x will be less than the overflow threshold. If the unscaled problem will not cause overflow, the Level 2 BLAS routine ZTPSV is called. If the matrix A is singular (A(j,j) = 0 for some j), then s is set to 0 and a non-trivial solution to A*x = 0 is returned.

Parameters

in

uplo

‘U’: A is upper triangular; ‘L’: A is lower triangular.

in

trans

‘N’: Solve A * x = s*b; ‘T’: Solve A**T * x = s*b; ‘C’: Solve A**H * x = s*b.

in

diag

‘N’: A is non-unit triangular; ‘U’: A is unit triangular.

in

normin

‘Y’: CNORM contains column norms on entry; ‘N’: CNORM is not set, will be computed.

in

n

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

in

AP

The triangular matrix A, packed columnwise. Array of dimension (n*(n+1)/2).

inout

X

On entry, the right hand side b. On exit, overwritten by the solution x. Array of dimension n.

out

scale

The scaling factor s for the triangular system.

inout

cnorm

If normin=’Y’, cnorm contains column norms on entry. If normin=’N’, cnorm returns the 1-norm of offdiagonal part of column j. Array of dimension n.

out

info

= 0: successful exit < 0: if info = -k, the k-th argument had an illegal value

latps#

This Page