libflame
revision_anchor
|
Functions | |
FLA_Error | FLA_Apply_pivots_ln_blk_var2 (FLA_Obj p, FLA_Obj A, fla_appiv_t *cntl) |
FLA_Error FLA_Apply_pivots_ln_blk_var2 | ( | FLA_Obj | p, |
FLA_Obj | A, | ||
fla_appiv_t * | cntl | ||
) |
References FLA_Apply_pivots_internal(), FLA_Cont_with_3x1_to_2x1(), FLA_Determine_blocksize(), FLA_Obj_length(), FLA_Part_2x1(), and FLA_Repart_2x1_to_3x1().
Referenced by FLA_Apply_pivots_ln().
{ FLA_Obj AT, A0, AB, A1, A2; FLA_Obj pT, p0, pB, pi1, p2; dim_t b; FLA_Part_2x1( A, &AT, &AB, 0, FLA_TOP ); FLA_Part_2x1( p, &pT, &pB, 0, FLA_TOP ); while ( FLA_Obj_length( AT ) < FLA_Obj_length( A ) ) { b = FLA_Determine_blocksize( AB, FLA_BOTTOM, FLA_Cntl_blocksize( cntl ) ); FLA_Repart_2x1_to_3x1( AT, &A0, /* ** */ /* ** */ &A1, AB, &A2, b, FLA_BOTTOM ); FLA_Repart_2x1_to_3x1( pT, &p0, /* ** */ /* ** */ &pi1, pB, &p2, b, FLA_BOTTOM ); /*------------------------------------------------------------*/ /* Apply pivots to a block and matrix */ FLA_Apply_pivots_internal( FLA_LEFT, FLA_NO_TRANSPOSE, pi1, AB, FLA_Cntl_sub_appiv( cntl ) ); /*------------------------------------------------------------*/ FLA_Cont_with_3x1_to_2x1( &AT, A0, A1, /* ** */ /* ** */ &AB, A2, FLA_TOP ); FLA_Cont_with_3x1_to_2x1( &pT, p0, pi1, /* ** */ /* ** */ &pB, p2, FLA_TOP ); } return FLA_SUCCESS; }