slapll (l) - Linux Manuals

slapll: two column vectors X and Y, let A = ( X Y )

NAME

SLAPLL - two column vectors X and Y, let A = ( X Y )

SYNOPSIS

SUBROUTINE SLAPLL(
N, X, INCX, Y, INCY, SSMIN )

    
INTEGER INCX, INCY, N

    
REAL SSMIN

    
REAL X( * ), Y( * )

PURPOSE

Given two column vectors X and Y, let The subroutine first computes the QR factorization of A = Q*R, and then computes the SVD of the 2-by-2 upper triangular matrix R. The smaller singular value of R is returned in SSMIN, which is used as the measurement of the linear dependency of the vectors X and Y.

ARGUMENTS

N (input) INTEGER
The length of the vectors X and Y.
X (input/output) REAL array,
dimension (1+(N-1)*INCX) On entry, X contains the N-vector X. On exit, X is overwritten.
INCX (input) INTEGER
The increment between successive elements of X. INCX > 0.
Y (input/output) REAL array,
dimension (1+(N-1)*INCY) On entry, Y contains the N-vector Y. On exit, Y is overwritten.
INCY (input) INTEGER
The increment between successive elements of Y. INCY > 0.
SSMIN (output) REAL
The smallest singular value of the N-by-2 matrix A = ( X Y ).

Pages related to slapll

  • slapll (3)
  • slapmt (l) - rearranges the columns of the M by N matrix X as specified by the permutation K(1),K(2),...,K(N) of the integers 1,...,N
  • sla_gbamv (l) - performs one of the matrix-vector operations y := alpha*abs(A)*abs(x) + beta*abs(y),
  • sla_gbrcond (l) - SLA_GERCOND Estimate the Skeel condition number of op(A) * op2(C) where op2 is determined by CMODE as follows CMODE = 1 op2(C) = C CMODE = 0 op2(C) = I CMODE = -1 op2(C) = inv(C) The Skeel condition number cond(A) = norminf( |inv(A)||A| ) is computed by computing scaling factors R such that diag(R)*A*op2(C) is row equilibrated and computing the standard infinity-norm condition number
  • sla_gbrfsx_extended (l) - computes ..
  • sla_geamv (l) - performs one of the matrix-vector operations y := alpha*abs(A)*abs(x) + beta*abs(y),
  • sla_gercond (l) - SLA_GERCOND estimate the Skeel condition number of op(A) * op2(C) where op2 is determined by CMODE as follows CMODE = 1 op2(C) = C CMODE = 0 op2(C) = I CMODE = -1 op2(C) = inv(C) The Skeel condition number cond(A) = norminf( |inv(A)||A| ) is computed by computing scaling factors R such that diag(R)*A*op2(C) is row equilibrated and computing the standard infinity-norm condition number
  • sla_gerfsx_extended (l) - computes ..