lapll#

Functions

void slapll(
    const INT           n,
          f32* restrict x,
    const INT           incx,
          f32* restrict y,
    const INT           incy,
          f32*          ssmin
);
void slapll(const INT n, f32 *restrict x, const INT incx, f32 *restrict y, const INT incy, f32 *ssmin)#

Given two column vectors X and Y, let.

                 A = ( X Y ).
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.

Parameters

in
n

The length of the vectors X and Y.

inout
x

On entry, X contains the N-vector X. On exit, X is overwritten.

in
incx

The increment between successive elements of X. INCX > 0.

inout
y

On entry, Y contains the N-vector Y. On exit, Y is overwritten.

in
incy

The increment between successive elements of Y. INCY > 0.

out
ssmin

The smallest singular value of the N-by-2 matrix A = ( X Y ).

Functions

void dlapll(
    const INT           n,
          f64* restrict x,
    const INT           incx,
          f64* restrict y,
    const INT           incy,
          f64*          ssmin
);
void dlapll(const INT n, f64 *restrict x, const INT incx, f64 *restrict y, const INT incy, f64 *ssmin)#

Given two column vectors X and Y, let.

                 A = ( X Y ).
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.

Parameters

in
n

The length of the vectors X and Y.

inout
x

On entry, X contains the N-vector X. On exit, X is overwritten.

in
incx

The increment between successive elements of X. INCX > 0.

inout
y

On entry, Y contains the N-vector Y. On exit, Y is overwritten.

in
incy

The increment between successive elements of Y. INCY > 0.

out
ssmin

The smallest singular value of the N-by-2 matrix A = ( X Y ).

Functions

void clapll(
    const INT           n,
          c64* restrict x,
    const INT           incx,
          c64* restrict y,
    const INT           incy,
          f32*          ssmin
);
void clapll(const INT n, c64 *restrict x, const INT incx, c64 *restrict y, const INT incy, f32 *ssmin)#

Given two column vectors X and Y, let.

                 A = ( X Y ).
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.

Parameters

in
n

The length of the vectors X and Y.

inout
x

On entry, X contains the N-vector X. On exit, X is overwritten.

in
incx

The increment between successive elements of X. INCX > 0.

inout
y

On entry, Y contains the N-vector Y. On exit, Y is overwritten.

in
incy

The increment between successive elements of Y. INCY > 0.

out
ssmin

The smallest singular value of the N-by-2 matrix A = ( X Y ).

Functions

void zlapll(
    const INT            n,
          c128* restrict x,
    const INT            incx,
          c128* restrict y,
    const INT            incy,
          f64*           ssmin
);
void zlapll(const INT n, c128 *restrict x, const INT incx, c128 *restrict y, const INT incy, f64 *ssmin)#

Given two column vectors X and Y, let.

                 A = ( X Y ).
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.

Parameters

in
n

The length of the vectors X and Y.

inout
x

On entry, X contains the N-vector X. On exit, X is overwritten.

in
incx

The increment between successive elements of X. INCX > 0.

inout
y

On entry, Y contains the N-vector Y. On exit, Y is overwritten.

in
incy

The increment between successive elements of Y. INCY > 0.

out
ssmin

The smallest singular value of the N-by-2 matrix A = ( X Y ).