unbdb5#

Functions

void cunbdb5(
    const INT           m1,
    const INT           m2,
    const INT           n,
          c64* restrict X1,
    const INT           incx1,
          c64* restrict X2,
    const INT           incx2,
          c64* restrict Q1,
    const INT           ldq1,
          c64* restrict Q2,
    const INT           ldq2,
          c64* restrict work,
    const INT           lwork,
          INT*          info
);
void cunbdb5(const INT m1, const INT m2, const INT n, c64 *restrict X1, const INT incx1, c64 *restrict X2, const INT incx2, c64 *restrict Q1, const INT ldq1, c64 *restrict Q2, const INT ldq2, c64 *restrict work, const INT lwork, INT *info)#

CUNBDB5 orthogonalizes the column vector X = [ X1 ] [ X2 ] with respect to the columns of Q = [ Q1 ] .

[ Q2 ] The columns of Q must be orthonormal.

If the projection is zero according to Kahan’s “twice is enough” criterion, then some other vector from the orthogonal complement is returned. This vector is chosen in an arbitrary but deterministic way.

Parameters

in
m1

The dimension of X1 and the number of rows in Q1. 0 <= M1.

in
m2

The dimension of X2 and the number of rows in Q2. 0 <= M2.

in
n

The number of columns in Q1 and Q2. 0 <= N.

inout
X1

Complex*16 array, dimension (M1). On entry, the top part of the vector to be orthogonalized. On exit, the top part of the projected vector.

in
incx1

Increment for entries of X1.

inout
X2

Complex*16 array, dimension (M2). On entry, the bottom part of the vector to be orthogonalized. On exit, the bottom part of the projected vector.

in
incx2

Increment for entries of X2.

in
Q1

Complex*16 array, dimension (LDQ1, N). The top part of the orthonormal basis matrix.

in
ldq1

The leading dimension of Q1. LDQ1 >= M1.

in
Q2

Complex*16 array, dimension (LDQ2, N). The bottom part of the orthonormal basis matrix.

in
ldq2

The leading dimension of Q2. LDQ2 >= M2.

out
work

Complex*16 array, dimension (LWORK).

in
lwork

The dimension of the array WORK. LWORK >= N.

out
info

= 0: successful exit. < 0: if INFO = -i, the i-th argument had an illegal value.

Functions

void zunbdb5(
    const INT            m1,
    const INT            m2,
    const INT            n,
          c128* restrict X1,
    const INT            incx1,
          c128* restrict X2,
    const INT            incx2,
          c128* restrict Q1,
    const INT            ldq1,
          c128* restrict Q2,
    const INT            ldq2,
          c128* restrict work,
    const INT            lwork,
          INT*           info
);
void zunbdb5(const INT m1, const INT m2, const INT n, c128 *restrict X1, const INT incx1, c128 *restrict X2, const INT incx2, c128 *restrict Q1, const INT ldq1, c128 *restrict Q2, const INT ldq2, c128 *restrict work, const INT lwork, INT *info)#

ZUNBDB5 orthogonalizes the column vector X = [ X1 ] [ X2 ] with respect to the columns of Q = [ Q1 ] .

[ Q2 ] The columns of Q must be orthonormal.

If the projection is zero according to Kahan’s “twice is enough” criterion, then some other vector from the orthogonal complement is returned. This vector is chosen in an arbitrary but deterministic way.

Parameters

in
m1

The dimension of X1 and the number of rows in Q1. 0 <= M1.

in
m2

The dimension of X2 and the number of rows in Q2. 0 <= M2.

in
n

The number of columns in Q1 and Q2. 0 <= N.

inout
X1

Complex*16 array, dimension (M1). On entry, the top part of the vector to be orthogonalized. On exit, the top part of the projected vector.

in
incx1

Increment for entries of X1.

inout
X2

Complex*16 array, dimension (M2). On entry, the bottom part of the vector to be orthogonalized. On exit, the bottom part of the projected vector.

in
incx2

Increment for entries of X2.

in
Q1

Complex*16 array, dimension (LDQ1, N). The top part of the orthonormal basis matrix.

in
ldq1

The leading dimension of Q1. LDQ1 >= M1.

in
Q2

Complex*16 array, dimension (LDQ2, N). The bottom part of the orthonormal basis matrix.

in
ldq2

The leading dimension of Q2. LDQ2 >= M2.

out
work

Complex*16 array, dimension (LWORK).

in
lwork

The dimension of the array WORK. LWORK >= N.

out
info

= 0: successful exit. < 0: if INFO = -i, the i-th argument had an illegal value.

On this page

This Page