orbdb5#

Functions

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

SORBDB5 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.

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

Double precision 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

Double precision 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

Double precision array, dimension (ldq1, n). The top part of the orthonormal basis matrix.

in
ldq1

The leading dimension of Q1. ldq1 >= m1.

in
Q2

Double precision array, dimension (ldq2, n). The bottom part of the orthonormal basis matrix.

in
ldq2

The leading dimension of Q2. ldq2 >= m2.

out
work

Double precision 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 dorbdb5(
    const INT           m1,
    const INT           m2,
    const INT           n,
          f64* restrict X1,
    const INT           incx1,
          f64* restrict X2,
    const INT           incx2,
    const f64* restrict Q1,
    const INT           ldq1,
    const f64* restrict Q2,
    const INT           ldq2,
          f64* restrict work,
    const INT           lwork,
          INT*          info
);
void dorbdb5(const INT m1, const INT m2, const INT n, f64 *restrict X1, const INT incx1, f64 *restrict X2, const INT incx2, const f64 *restrict Q1, const INT ldq1, const f64 *restrict Q2, const INT ldq2, f64 *restrict work, const INT lwork, INT *info)#

DORBDB5 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.

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

Double precision 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

Double precision 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

Double precision array, dimension (ldq1, n). The top part of the orthonormal basis matrix.

in
ldq1

The leading dimension of Q1. ldq1 >= m1.

in
Q2

Double precision array, dimension (ldq2, n). The bottom part of the orthonormal basis matrix.

in
ldq2

The leading dimension of Q2. ldq2 >= m2.

out
work

Double precision 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