Functions | |
FLA_Error | FLA_Apply_Q_UT_inc_lhc_blk_var1 (FLA_Obj A, FLA_Obj TW, FLA_Obj W1, FLA_Obj B, fla_apqutinc_t *cntl) |
FLA_Error FLA_Apply_Q_UT_inc_lhc_blk_var1 | ( | FLA_Obj | A, | |
FLA_Obj | TW, | |||
FLA_Obj | W1, | |||
FLA_Obj | B, | |||
fla_apqutinc_t * | cntl | |||
) |
References FLA_Apply_Q_UT_internal(), FLA_Apply_Q_UT_UD_internal(), FLA_Cont_with_3x1_to_2x1(), FLA_Cont_with_3x3_to_2x2(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Part_2x1(), FLA_Part_2x2(), FLA_Repart_2x1_to_3x1(), and FLA_Repart_2x2_to_3x3().
Referenced by FLA_Apply_Q_UT_inc_lhc().
00036 { 00037 FLA_Obj ATL, ATR, A00, A01, A02, 00038 ABL, ABR, A10, A11, A12, 00039 A20, A21, A22; 00040 00041 FLA_Obj TTL, WTR, T00, W01, W02, 00042 TBL, TBR, T10, T11, W12, 00043 T20, T21, T22; 00044 00045 FLA_Obj BT, B0, 00046 BB, B1, 00047 B2; 00048 00049 dim_t b; 00050 00051 FLA_Part_2x2( A, &ATL, &ATR, 00052 &ABL, &ABR, 0, 0, FLA_TL ); 00053 00054 FLA_Part_2x2( TW, &TTL, &WTR, 00055 &TBL, &TBR, 0, 0, FLA_TL ); 00056 00057 FLA_Part_2x1( B, &BT, 00058 &BB, 0, FLA_TOP ); 00059 00060 while ( FLA_Obj_min_dim( ABR ) > 0 ){ 00061 00062 b = FLA_Determine_blocksize( ABR, FLA_BR, FLA_Cntl_blocksize( cntl ) ); 00063 00064 FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &A01, &A02, 00065 /* ************* */ /* ******************** */ 00066 &A10, /**/ &A11, &A12, 00067 ABL, /**/ ABR, &A20, /**/ &A21, &A22, 00068 b, b, FLA_BR ); 00069 00070 FLA_Repart_2x2_to_3x3( TTL, /**/ WTR, &T00, /**/ &W01, &W02, 00071 /* ************* */ /* ******************** */ 00072 &T10, /**/ &T11, &W12, 00073 TBL, /**/ TBR, &T20, /**/ &T21, &T22, 00074 b, b, FLA_BR ); 00075 00076 FLA_Repart_2x1_to_3x1( BT, &B0, 00077 /* ** */ /* ** */ 00078 &B1, 00079 BB, &B2, b, FLA_BOTTOM ); 00080 00081 /*------------------------------------------------------------*/ 00082 00083 /* 00084 Apply Q^H to B1 from the left: 00085 00086 B1 = Q^H * B1 00087 00088 where Q is formed from A11 and T11. Note that W1 refers 00089 to a row-panel of blocks where each block refers to an 00090 nb_alg-by-b row-panel of workspace. 00091 */ 00092 00093 FLA_Apply_Q_UT_internal( FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_COLUMNWISE, 00094 A11, T11, W1, B1, 00095 FLA_Cntl_sub_apqut( cntl ) ); 00096 00097 /* 00098 Apply Q^H to B1 and B2 from the left: 00099 00100 / B1 \ = Q^H * / B1 \ 00101 \ B2 / \ B2 / 00102 00103 where Q is formed from A21 and T21. 00104 */ 00105 00106 FLA_Apply_Q_UT_UD_internal( FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_COLUMNWISE, 00107 A21, T21, W1, B1, 00108 B2, 00109 FLA_Cntl_sub_apqutud( cntl ) ); 00110 00111 /*------------------------------------------------------------*/ 00112 00113 FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, A01, /**/ A02, 00114 A10, A11, /**/ A12, 00115 /* ************** */ /* ****************** */ 00116 &ABL, /**/ &ABR, A20, A21, /**/ A22, 00117 FLA_TL ); 00118 00119 FLA_Cont_with_3x3_to_2x2( &TTL, /**/ &WTR, T00, W01, /**/ W02, 00120 T10, T11, /**/ W12, 00121 /* ************** */ /* ****************** */ 00122 &TBL, /**/ &TBR, T20, T21, /**/ T22, 00123 FLA_TL ); 00124 00125 FLA_Cont_with_3x1_to_2x1( &BT, B0, 00126 B1, 00127 /* ** */ /* ** */ 00128 &BB, B2, FLA_TOP ); 00129 00130 } 00131 00132 return FLA_SUCCESS; 00133 }