{
FLA_Obj ATL, ATR, A00, a01, A02,
ABL, ABR, a10t, alpha11, a12t,
A20, a21, A22;
FLA_Obj tauT, tau0,
tauB, tau1,
tau2;
FLA_Obj TTL, TTR, T00, t01, T02,
TBL, TBR, t10t, tau11, t12t,
T20, t21, T22;
FLA_Part_2x2( A, &ATL, &ATR,
&ABL, &ABR, 0, 0, FLA_TL );
FLA_Part_2x1( tau, &tauT,
&tauB, 0, FLA_TOP );
FLA_Part_2x2( T, &TTL, &TTR,
&TBL, &TBR, 0, 0, FLA_TL );
while ( FLA_Obj_length( tauB ) > 0 ){
FLA_Repart_2x2_to_3x3( ATL, ATR, &A00, &a01, &A02,
&a10t, &alpha11, &a12t,
ABL, ABR, &A20, &a21, &A22,
1, 1, FLA_BR );
FLA_Repart_2x1_to_3x1( tauT, &tau0,
&tau1,
tauB, &tau2, 1, FLA_BOTTOM );
FLA_Repart_2x2_to_3x3( TTL, TTR, &T00, &t01, &T02,
&t10t, &tau11, &t12t,
TBL, TBR, &T20, &t21, &T22,
1, 1, FLA_BR );
FLA_Copy_external( tau1, tau11 );
FLA_Copyt_external( FLA_CONJ_NO_TRANSPOSE, a01, t01 );
FLA_Gemvc_external( FLA_CONJ_NO_TRANSPOSE, FLA_NO_CONJUGATE, FLA_ONE, A02, a12t, FLA_ONE, t01 );
FLA_Cont_with_3x3_to_2x2( &ATL, &ATR, A00, a01, A02,
a10t, alpha11, a12t,
&ABL, &ABR, A20, a21, A22,
FLA_TL );
FLA_Cont_with_3x1_to_2x1( &tauT, tau0,
tau1,
&tauB, tau2, FLA_TOP );
FLA_Cont_with_3x3_to_2x2( &TTL, &TTR, T00, t01, T02,
t10t, tau11, t12t,
&TBL, &TBR, T20, t21, T22,
FLA_TL );
}
return FLA_SUCCESS;
}