UNBDB1: simultaneously bidiagonalizes the blocks of a tall and skinny matrix X with orthonomal columns: [ B11 ] [ X11 ] [ P1 | ] [ 0 ] [--—] = [------—] [--—] Q1**T . [ X21 ] [ | P2 ] [ B21 ] [ 0 ] X11 is P-by-Q, and X21 is (M-P)-by-Q. Q must be no larger than P, M-P, or M-Q. Routines CUNBDB2, CUNBDB3, and CUNBDB4 handle cases in which Q is not the minimum dimension. The unitary matrices P1, P2, and Q1 are P-by-P, (M-P)-by-(M-P), and (M-Q)-by-(M-Q), respectively. They are represented implicitly by Householder vectors. B11 and B12 are Q-by-Q bidiagonal matrices represented implicitly by angles THETA, PHI. More...

Public Member Functions
subroutine	cunbdb1 (m, p, q, x11, ldx11, x21, ldx21, theta, phi, taup1, taup2, tauq1, work, lwork, info)

	la_cunbdb1

	la_wunbdb1

subroutine	zunbdb1 (m, p, q, x11, ldx11, x21, ldx21, theta, phi, taup1, taup2, tauq1, work, lwork, info)

	la_zunbdb1

Detailed Description

UNBDB1: simultaneously bidiagonalizes the blocks of a tall and skinny matrix X with orthonomal columns: [ B11 ] [ X11 ] [ P1 | ] [ 0 ] [--—] = [------—] [--—] Q1**T . [ X21 ] [ | P2 ] [ B21 ] [ 0 ] X11 is P-by-Q, and X21 is (M-P)-by-Q. Q must be no larger than P, M-P, or M-Q. Routines CUNBDB2, CUNBDB3, and CUNBDB4 handle cases in which Q is not the minimum dimension. The unitary matrices P1, P2, and Q1 are P-by-P, (M-P)-by-(M-P), and (M-Q)-by-(M-Q), respectively. They are represented implicitly by Householder vectors. B11 and B12 are Q-by-Q bidiagonal matrices represented implicitly by angles THETA, PHI.

Member Function/Subroutine Documentation

◆ cunbdb1()

subroutine la_lapack::unbdb1::cunbdb1	(	integer(ilp), intent(in)	m,
		integer(ilp), intent(in)	p,
		integer(ilp), intent(in)	q,
		complex(sp), dimension(ldx11,*), intent(inout)	x11,
		integer(ilp), intent(in)	ldx11,
		complex(sp), dimension(ldx21,*), intent(inout)	x21,
		integer(ilp), intent(in)	ldx21,
		real(sp), dimension(*), intent(out)	theta,
		real(sp), dimension(*), intent(out)	phi,
		complex(sp), dimension(*), intent(out)	taup1,
		complex(sp), dimension(*), intent(out)	taup2,
		complex(sp), dimension(*), intent(out)	tauq1,
		complex(sp), dimension(*), intent(out)	work,
		integer(ilp), intent(in)	lwork,
		integer(ilp), intent(out)	info
	)

◆ la_cunbdb1()

la_lapack::unbdb1::la_cunbdb1

◆ la_wunbdb1()

la_lapack::unbdb1::la_wunbdb1

◆ la_zunbdb1()

la_lapack::unbdb1::la_zunbdb1

◆ zunbdb1()

subroutine la_lapack::unbdb1::zunbdb1	(	integer(ilp), intent(in)	m,
		integer(ilp), intent(in)	p,
		integer(ilp), intent(in)	q,
		complex(dp), dimension(ldx11,*), intent(inout)	x11,
		integer(ilp), intent(in)	ldx11,
		complex(dp), dimension(ldx21,*), intent(inout)	x21,
		integer(ilp), intent(in)	ldx21,
		real(dp), dimension(*), intent(out)	theta,
		real(dp), dimension(*), intent(out)	phi,
		complex(dp), dimension(*), intent(out)	taup1,
		complex(dp), dimension(*), intent(out)	taup2,
		complex(dp), dimension(*), intent(out)	tauq1,
		complex(dp), dimension(*), intent(out)	work,
		integer(ilp), intent(in)	lwork,
		integer(ilp), intent(out)	info
	)

The documentation for this interface was generated from the following file:

la_lapack.f90