void ExplicitTriang( AbstractDistMatrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_ONLY(CallStackEntry cse("qr::ExplicitTriang"))
DistMatrix<F,MD,STAR> t(A.Grid());
DistMatrix<Base<F>,MD,STAR> d(A.Grid());
if( ctrl.colPiv )
{
DistMatrix<Int,VC,STAR> p(A.Grid());
BusingerGolub( A, t, d, p, ctrl );
}
else
Householder( A, t, d );
A.Resize( t.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}
void ExplicitTriang( Matrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_ONLY(CallStackEntry cse("qr::ExplicitTriang"))
Matrix<F> t;
Matrix<Base<F>> d;
if( ctrl.colPiv )
{
Matrix<Int> p;
BusingerGolub( A, t, d, p, ctrl );
}
else
Householder( A, t, d );
A.Resize( t.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}
void ExplicitTriang( ElementalMatrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_CSE
DistMatrix<F,MD,STAR> householderScalars(A.Grid());
DistMatrix<Base<F>,MD,STAR> signature(A.Grid());
if( ctrl.colPiv )
{
DistPermutation Omega(A.Grid());
BusingerGolub( A, householderScalars, signature, Omega, ctrl );
}
else
Householder( A, householderScalars, signature );
A.Resize( householderScalars.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}
void ExplicitTriang( Matrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_CSE
Matrix<F> householderScalars;
Matrix<Base<F>> signature;
if( ctrl.colPiv )
{
Permutation Omega;
BusingerGolub( A, householderScalars, signature, Omega, ctrl );
}
else
Householder( A, householderScalars, signature );
A.Resize( householderScalars.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}
void ExplicitTriang( ElementalMatrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_ONLY(CSE cse("qr::ExplicitTriang"))
DistMatrix<F,MD,STAR> t(A.Grid());
DistMatrix<Base<F>,MD,STAR> d(A.Grid());
if( ctrl.colPiv )
{
DistPermutation Omega(A.Grid());
BusingerGolub( A, t, d, Omega, ctrl );
}
else
Householder( A, t, d );
A.Resize( t.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}
void ExplicitTriang( Matrix<F>& A, const QRCtrl<Base<F>>& ctrl )
{
DEBUG_ONLY(CSE cse("qr::ExplicitTriang"))
Matrix<F> t;
Matrix<Base<F>> d;
if( ctrl.colPiv )
{
Permutation Omega;
BusingerGolub( A, t, d, Omega, ctrl );
}
else
Householder( A, t, d );
A.Resize( t.Height(), A.Width() );
MakeTrapezoidal( UPPER, A );
}