|
libflame
12600
|
Go to the source code of this file.
Functions | |
| FLA_Error | FLA_Copyr_l_blk_var1 (FLA_Obj A, FLA_Obj B, fla_copyr_t *cntl) |
| FLA_Error | FLA_Copyr_l_blk_var2 (FLA_Obj A, FLA_Obj B, fla_copyr_t *cntl) |
| FLA_Error | FLA_Copyr_l_blk_var3 (FLA_Obj A, FLA_Obj B, fla_copyr_t *cntl) |
| FLA_Error | FLA_Copyr_l_blk_var4 (FLA_Obj A, FLA_Obj B, fla_copyr_t *cntl) |
| FLA_Error FLA_Copyr_l_blk_var1 | ( | FLA_Obj | A, |
| FLA_Obj | B, | ||
| fla_copyr_t * | cntl | ||
| ) |
References FLA_Cont_with_3x3_to_2x2(), FLA_Copy_internal(), FLA_Copyr_internal(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Part_2x2(), and FLA_Repart_2x2_to_3x3().
Referenced by FLA_Copyr_l().
{
FLA_Obj ATL, ATR, A00, A01, A02,
ABL, ABR, A10, A11, A12,
A20, A21, A22;
FLA_Obj BTL, BTR, B00, B01, B02,
BBL, BBR, B10, B11, B12,
B20, B21, B22;
dim_t b;
FLA_Part_2x2( A, &ATL, &ATR,
&ABL, &ABR, 0, 0, FLA_TL );
FLA_Part_2x2( B, &BTL, &BTR,
&BBL, &BBR, 0, 0, FLA_TL );
while ( FLA_Obj_min_dim( ABR ) > 0 ){
b = FLA_Determine_blocksize( ABR, FLA_BR, FLA_Cntl_blocksize( cntl ) );
FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &A01, &A02,
/* ************* */ /* ******************** */
&A10, /**/ &A11, &A12,
ABL, /**/ ABR, &A20, /**/ &A21, &A22,
b, b, FLA_BR );
FLA_Repart_2x2_to_3x3( BTL, /**/ BTR, &B00, /**/ &B01, &B02,
/* ************* */ /* ******************** */
&B10, /**/ &B11, &B12,
BBL, /**/ BBR, &B20, /**/ &B21, &B22,
b, b, FLA_BR );
/*------------------------------------------------------------*/
// B11 = tril( A11 );
FLA_Copyr_internal( FLA_LOWER_TRIANGULAR, A11, B11,
FLA_Cntl_sub_copyr( cntl ) );
// B10 = A10;
FLA_Copy_internal( A10, B10,
FLA_Cntl_sub_copy( cntl ) );
/*------------------------------------------------------------*/
FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, A01, /**/ A02,
A10, A11, /**/ A12,
/* ************** */ /* ****************** */
&ABL, /**/ &ABR, A20, A21, /**/ A22,
FLA_TL );
FLA_Cont_with_3x3_to_2x2( &BTL, /**/ &BTR, B00, B01, /**/ B02,
B10, B11, /**/ B12,
/* ************** */ /* ****************** */
&BBL, /**/ &BBR, B20, B21, /**/ B22,
FLA_TL );
}
return FLA_SUCCESS;
}
| FLA_Error FLA_Copyr_l_blk_var2 | ( | FLA_Obj | A, |
| FLA_Obj | B, | ||
| fla_copyr_t * | cntl | ||
| ) |
References FLA_Cont_with_3x3_to_2x2(), FLA_Copy_internal(), FLA_Copyr_internal(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Part_2x2(), and FLA_Repart_2x2_to_3x3().
Referenced by FLA_Copyr_l().
{
FLA_Obj ATL, ATR, A00, A01, A02,
ABL, ABR, A10, A11, A12,
A20, A21, A22;
FLA_Obj BTL, BTR, B00, B01, B02,
BBL, BBR, B10, B11, B12,
B20, B21, B22;
dim_t b;
FLA_Part_2x2( A, &ATL, &ATR,
&ABL, &ABR, 0, 0, FLA_BR );
FLA_Part_2x2( B, &BTL, &BTR,
&BBL, &BBR, 0, 0, FLA_BR );
while ( FLA_Obj_min_dim( ATL ) > 0 ){
b = FLA_Determine_blocksize( ATL, FLA_TL, FLA_Cntl_blocksize( cntl ) );
FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, &A01, /**/ &A02,
&A10, &A11, /**/ &A12,
/* ************* */ /* ******************** */
ABL, /**/ ABR, &A20, &A21, /**/ &A22,
b, b, FLA_TL );
FLA_Repart_2x2_to_3x3( BTL, /**/ BTR, &B00, &B01, /**/ &B02,
&B10, &B11, /**/ &B12,
/* ************* */ /* ******************** */
BBL, /**/ BBR, &B20, &B21, /**/ &B22,
b, b, FLA_TL );
/*------------------------------------------------------------*/
// B11 = tril( A11 );
FLA_Copyr_internal( FLA_LOWER_TRIANGULAR, A11, B11,
FLA_Cntl_sub_copyr( cntl ) );
// B10 = A10;
FLA_Copy_internal( A10, B10,
FLA_Cntl_sub_copy( cntl ) );
/*------------------------------------------------------------*/
FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, /**/ A01, A02,
/* ************** */ /* ****************** */
A10, /**/ A11, A12,
&ABL, /**/ &ABR, A20, /**/ A21, A22,
FLA_BR );
FLA_Cont_with_3x3_to_2x2( &BTL, /**/ &BTR, B00, /**/ B01, B02,
/* ************** */ /* ****************** */
B10, /**/ B11, B12,
&BBL, /**/ &BBR, B20, /**/ B21, B22,
FLA_BR );
}
return FLA_SUCCESS;
}
| FLA_Error FLA_Copyr_l_blk_var3 | ( | FLA_Obj | A, |
| FLA_Obj | B, | ||
| fla_copyr_t * | cntl | ||
| ) |
References FLA_Cont_with_3x3_to_2x2(), FLA_Copy_internal(), FLA_Copyr_internal(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Part_2x2(), and FLA_Repart_2x2_to_3x3().
Referenced by FLA_Copyr_l().
{
FLA_Obj ATL, ATR, A00, A01, A02,
ABL, ABR, A10, A11, A12,
A20, A21, A22;
FLA_Obj BTL, BTR, B00, B01, B02,
BBL, BBR, B10, B11, B12,
B20, B21, B22;
dim_t b;
FLA_Part_2x2( A, &ATL, &ATR,
&ABL, &ABR, 0, 0, FLA_TL );
FLA_Part_2x2( B, &BTL, &BTR,
&BBL, &BBR, 0, 0, FLA_TL );
while ( FLA_Obj_min_dim( ABR ) > 0 ){
b = FLA_Determine_blocksize( ABR, FLA_BR, FLA_Cntl_blocksize( cntl ) );
FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &A01, &A02,
/* ************* */ /* ******************** */
&A10, /**/ &A11, &A12,
ABL, /**/ ABR, &A20, /**/ &A21, &A22,
b, b, FLA_BR );
FLA_Repart_2x2_to_3x3( BTL, /**/ BTR, &B00, /**/ &B01, &B02,
/* ************* */ /* ******************** */
&B10, /**/ &B11, &B12,
BBL, /**/ BBR, &B20, /**/ &B21, &B22,
b, b, FLA_BR );
/*------------------------------------------------------------*/
// B11 = tril( A11 );
FLA_Copyr_internal( FLA_LOWER_TRIANGULAR, A11, B11,
FLA_Cntl_sub_copyr( cntl ) );
// B21 = A21;
FLA_Copy_internal( A21, B21,
FLA_Cntl_sub_copy( cntl ) );
/*------------------------------------------------------------*/
FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, A01, /**/ A02,
A10, A11, /**/ A12,
/* ************** */ /* ****************** */
&ABL, /**/ &ABR, A20, A21, /**/ A22,
FLA_TL );
FLA_Cont_with_3x3_to_2x2( &BTL, /**/ &BTR, B00, B01, /**/ B02,
B10, B11, /**/ B12,
/* ************** */ /* ****************** */
&BBL, /**/ &BBR, B20, B21, /**/ B22,
FLA_TL );
}
return FLA_SUCCESS;
}
| FLA_Error FLA_Copyr_l_blk_var4 | ( | FLA_Obj | A, |
| FLA_Obj | B, | ||
| fla_copyr_t * | cntl | ||
| ) |
References FLA_Cont_with_3x3_to_2x2(), FLA_Copy_internal(), FLA_Copyr_internal(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Part_2x2(), and FLA_Repart_2x2_to_3x3().
Referenced by FLA_Copyr_l().
{
FLA_Obj ATL, ATR, A00, A01, A02,
ABL, ABR, A10, A11, A12,
A20, A21, A22;
FLA_Obj BTL, BTR, B00, B01, B02,
BBL, BBR, B10, B11, B12,
B20, B21, B22;
dim_t b;
FLA_Part_2x2( A, &ATL, &ATR,
&ABL, &ABR, 0, 0, FLA_BR );
FLA_Part_2x2( B, &BTL, &BTR,
&BBL, &BBR, 0, 0, FLA_BR );
while ( FLA_Obj_min_dim( ATL ) > 0 ){
b = FLA_Determine_blocksize( ATL, FLA_TL, FLA_Cntl_blocksize( cntl ) );
FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, &A01, /**/ &A02,
&A10, &A11, /**/ &A12,
/* ************* */ /* ******************** */
ABL, /**/ ABR, &A20, &A21, /**/ &A22,
b, b, FLA_TL );
FLA_Repart_2x2_to_3x3( BTL, /**/ BTR, &B00, &B01, /**/ &B02,
&B10, &B11, /**/ &B12,
/* ************* */ /* ******************** */
BBL, /**/ BBR, &B20, &B21, /**/ &B22,
b, b, FLA_TL );
/*------------------------------------------------------------*/
// B11 = tril( A11 );
FLA_Copyr_internal( FLA_LOWER_TRIANGULAR, A11, B11,
FLA_Cntl_sub_copyr( cntl ) );
// B21 = A21;
FLA_Copy_internal( A21, B21,
FLA_Cntl_sub_copy( cntl ) );
/*------------------------------------------------------------*/
FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, /**/ A01, A02,
/* ************** */ /* ****************** */
A10, /**/ A11, A12,
&ABL, /**/ &ABR, A20, /**/ A21, A22,
FLA_BR );
FLA_Cont_with_3x3_to_2x2( &BTL, /**/ &BTR, B00, /**/ B01, B02,
/* ************** */ /* ****************** */
B10, /**/ B11, B12,
&BBL, /**/ &BBR, B20, /**/ B21, B22,
FLA_BR );
}
return FLA_SUCCESS;
}
1.7.6.1