libflame
revision_anchor
|
Functions | |
FLA_Error | FLA_UDdate_UT_inc_blk_var1 (FLA_Obj R, FLA_Obj C, FLA_Obj D, FLA_Obj T, FLA_Obj W, fla_uddateutinc_t *cntl) |
FLA_Error FLA_UDdate_UT_inc_blk_var1 | ( | FLA_Obj | R, |
FLA_Obj | C, | ||
FLA_Obj | D, | ||
FLA_Obj | T, | ||
FLA_Obj | W, | ||
fla_uddateutinc_t * | cntl | ||
) |
References FLA_Apply_QUD_UT_internal(), FLA_Cont_with_1x3_to_1x2(), FLA_Cont_with_3x3_to_2x2(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Obj_width(), FLA_Part_1x2(), FLA_Part_2x2(), FLA_Repart_1x2_to_1x3(), FLA_Repart_2x2_to_3x3(), and FLA_UDdate_UT_internal().
Referenced by FLASH_UDdate_UT_inc().
{ FLA_Obj RTL, RTR, R00, R01, R02, RBL, RBR, R10, R11, R12, R20, R21, R22; FLA_Obj CL, CR, C0, C1, C2; FLA_Obj DL, DR, D0, D1, D2; FLA_Obj TL, TR, T0, T1, T2; FLA_Obj WTL, WTR, W00, W01, W02, WBL, WBR, W10, W11, W12, W20, W21, W22; dim_t b; FLA_Part_2x2( R, &RTL, &RTR, &RBL, &RBR, 0, 0, FLA_TL ); FLA_Part_1x2( C, &CL, &CR, 0, FLA_LEFT ); FLA_Part_1x2( D, &DL, &DR, 0, FLA_LEFT ); FLA_Part_1x2( T, &TL, &TR, 0, FLA_LEFT ); FLA_Part_2x2( W, &WTL, &WTR, &WBL, &WBR, 0, 0, FLA_TL ); while ( FLA_Obj_min_dim( RBR ) > 0 ){ b = FLA_Determine_blocksize( RBR, FLA_BR, FLA_Cntl_blocksize( cntl ) ); FLA_Repart_2x2_to_3x3( RTL, /**/ RTR, &R00, /**/ &R01, &R02, /* ************* */ /* ******************** */ &R10, /**/ &R11, &R12, RBL, /**/ RBR, &R20, /**/ &R21, &R22, b, b, FLA_BR ); FLA_Repart_1x2_to_1x3( CL, /**/ CR, &C0, /**/ &C1, &C2, b, FLA_RIGHT ); FLA_Repart_1x2_to_1x3( DL, /**/ DR, &D0, /**/ &D1, &D2, b, FLA_RIGHT ); FLA_Repart_1x2_to_1x3( TL, /**/ TR, &T0, /**/ &T1, &T2, b, FLA_RIGHT ); FLA_Repart_2x2_to_3x3( WTL, /**/ WTR, &W00, /**/ &W01, &W02, /* ************* */ /* ******************** */ &W10, /**/ &W11, &W12, WBL, /**/ WBR, &W20, /**/ &W21, &W22, b, b, FLA_BR ); /*------------------------------------------------------------*/ /* Perform an up/downdate of the upper triangular factor R11 via up/downdating UT Householder transformations: [ R11, ... C1, ... D1, T1 ] = FLA_UDdate_UT( R11, ... C1, ... D1, T1 ); by updating R11 in such a way that removes the contributions of the rows in D1 while simultaneously adding new contributions to the factorization from the rows of C1. Note that C1 and D1 are also updated in the process. */ FLA_UDdate_UT_internal( R11, C1, D1, T1, FLA_Cntl_sub_uddateut( cntl ) ); if ( FLA_Obj_width( R12 ) > 0 ) { /* Apply Q' to R12, C2, and D2 from the left: / R12 \ / R12 \ | C2 | = Q' * | C2 | \ D2 / \ D2 / where Q is formed from C1, D1, and T1. */ FLA_Apply_QUD_UT_internal( FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_FORWARD, FLA_COLUMNWISE, T1, W12, R12, C1, C2, D1, D2, FLA_Cntl_sub_apqudut( cntl ) ); } /*------------------------------------------------------------*/ FLA_Cont_with_3x3_to_2x2( &RTL, /**/ &RTR, R00, R01, /**/ R02, R10, R11, /**/ R12, /* ************** */ /* ****************** */ &RBL, /**/ &RBR, R20, R21, /**/ R22, FLA_TL ); FLA_Cont_with_1x3_to_1x2( &CL, /**/ &CR, C0, C1, /**/ C2, FLA_LEFT ); FLA_Cont_with_1x3_to_1x2( &DL, /**/ &DR, D0, D1, /**/ D2, FLA_LEFT ); FLA_Cont_with_1x3_to_1x2( &TL, /**/ &TR, T0, T1, /**/ T2, FLA_LEFT ); FLA_Cont_with_3x3_to_2x2( &WTL, /**/ &WTR, W00, W01, /**/ W02, W10, W11, /**/ W12, /* ************** */ /* ****************** */ &WBL, /**/ &WBR, W20, W21, /**/ W22, FLA_TL ); } return FLA_SUCCESS; }