libflame
revision_anchor
|
Functions | |
FLA_Error | FLA_CAQR_UT_inc_factorize_panels (dim_t nb_part, FLA_Obj A, FLA_Obj TW) |
FLA_Error FLA_CAQR_UT_inc_factorize_panels | ( | dim_t | nb_part, |
FLA_Obj | A, | ||
FLA_Obj | TW | ||
) |
References FLA_Cont_with_3x1_to_2x1(), FLA_Obj_length(), FLA_Part_2x1(), FLA_Repart_2x1_to_3x1(), and FLASH_QR_UT_inc().
Referenced by FLASH_CAQR_UT_inc_noopt().
{ FLA_Obj AT, A0, AB, A1, A2; FLA_Obj TWT, TW0, TWB, TW1, TW2; dim_t b; FLA_Part_2x1( A, &AT, &AB, 0, FLA_TOP ); FLA_Part_2x1( TW, &TWT, &TWB, 0, FLA_TOP ); while ( FLA_Obj_length( AB ) > 0 ){ b = min( nb_part, FLA_Obj_length( AB ) ); FLA_Repart_2x1_to_3x1( AT, &A0, /* ** */ /* ** */ &A1, AB, &A2, b, FLA_BOTTOM ); FLA_Repart_2x1_to_3x1( TWT, &TW0, /* ** */ /* ** */ &TW1, TWB, &TW2, b, FLA_BOTTOM ); /*------------------------------------------------------------*/ // Perform an incremental QR factorization on A1, writing triangular // block Householder factors to T in TW1. FLASH_QR_UT_inc( A1, TW1 ); /*------------------------------------------------------------*/ FLA_Cont_with_3x1_to_2x1( &AT, A0, A1, /* ** */ /* ** */ &AB, A2, FLA_TOP ); FLA_Cont_with_3x1_to_2x1( &TWT, TW0, TW1, /* ** */ /* ** */ &TWB, TW2, FLA_TOP ); } return FLA_SUCCESS; }