libflame
revision_anchor
|
Go to the source code of this file.
Functions | |
FLA_Error | FLA_LU_piv_internal (FLA_Obj A, FLA_Obj p, fla_lu_t *cntl) |
FLA_Error | FLA_LU_piv_solve (FLA_Obj A, FLA_Obj p, FLA_Obj B, FLA_Obj X) |
FLA_Error | FLASH_LU_piv_solve (FLA_Obj A, FLA_Obj p, FLA_Obj B, FLA_Obj X) |
FLA_Error FLA_LU_piv_internal | ( | FLA_Obj | A, |
FLA_Obj | p, | ||
fla_lu_t * | cntl | ||
) |
References FLA_Check_error_level(), FLA_Check_null_pointer(), FLA_LU_piv_blk_ext(), FLA_LU_piv_blk_var3(), FLA_LU_piv_blk_var4(), FLA_LU_piv_blk_var5(), FLA_LU_piv_macro_task(), FLA_LU_piv_opt_var3(), FLA_LU_piv_opt_var4(), FLA_LU_piv_opt_var5(), FLA_LU_piv_unb_ext(), FLA_LU_piv_unb_var3(), FLA_LU_piv_unb_var4(), FLA_LU_piv_unb_var5(), and FLASH_Queue_get_enabled().
Referenced by FLA_LU_piv(), FLA_LU_piv_blk_var3(), FLA_LU_piv_blk_var4(), FLA_LU_piv_blk_var5(), FLA_LU_piv_task(), and FLASH_LU_piv().
{ FLA_Error r_val = FLA_SUCCESS; if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING ) { FLA_Error e_val = FLA_Check_null_pointer( ( void* ) cntl ); FLA_Check_error_code( e_val ); } if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER && FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM && FLASH_Queue_get_enabled( ) ) { // Enqueue ENQUEUE_FLASH_LU_piv_macro( A, *FLASH_OBJ_PTR_AT( p ), cntl ); } else { /* if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER && FLA_Obj_elemtype( A ) == FLA_SCALAR && !FLASH_Queue_get_enabled( ) ) { // Execute leaf cntl = fla_lu_piv_cntl_leaf; } */ // Choose implementation. if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_EXTERN || FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM ) { if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER ) { // When performing an unblocked operation on a hierarhical // object, we assume it's because we need to perform the // the operation on a macro block. r_val = FLA_LU_piv_macro_task( A, *FLASH_OBJ_PTR_AT( p ), cntl ); } else { r_val = FLA_LU_piv_blk_ext( A, p ); } } else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_EXTERN ) { if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER ) { // When performing an unblocked operation on a hierarhical // object, we assume it's because we need to perform the // the operation on a macro block. r_val = FLA_LU_piv_macro_task( A, *FLASH_OBJ_PTR_AT( p ), cntl ); } else { r_val = FLA_LU_piv_unb_ext( A, p ); } } #ifdef FLA_ENABLE_NON_CRITICAL_CODE else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT3 ) { r_val = FLA_LU_piv_unb_var3( A, p ); } else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT4 ) { r_val = FLA_LU_piv_unb_var4( A, p ); } #endif else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT5 ) { r_val = FLA_LU_piv_unb_var5( A, p ); } #ifdef FLA_ENABLE_NON_CRITICAL_CODE else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT3 ) { r_val = FLA_LU_piv_opt_var3( A, p ); } else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT4 ) { r_val = FLA_LU_piv_opt_var4( A, p ); } #endif else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT5 ) { r_val = FLA_LU_piv_opt_var5( A, p ); } #ifdef FLA_ENABLE_NON_CRITICAL_CODE else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 ) { r_val = FLA_LU_piv_blk_var3( A, p, cntl ); } else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT4 ) { r_val = FLA_LU_piv_blk_var4( A, p, cntl ); } #endif else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT5 ) { r_val = FLA_LU_piv_blk_var5( A, p, cntl ); } else { FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED ); } } return r_val; }
References FLA_Apply_pivots(), FLA_Check_error_level(), FLA_Copy_external(), FLA_LU_piv_solve_check(), FLA_ONE, and FLA_Trsm_external().
{ // Check parameters. if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING ) FLA_LU_piv_solve_check( A, p, B, X ); FLA_Copy_external( B, X ); FLA_Apply_pivots( FLA_LEFT, FLA_NO_TRANSPOSE, p, X ); FLA_Trsm_external( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_UNIT_DIAG, FLA_ONE, A, X ); FLA_Trsm_external( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_NONUNIT_DIAG, FLA_ONE, A, X ); return FLA_SUCCESS; }
References FLA_Check_error_level(), FLA_LU_piv_solve_check(), FLA_ONE, FLASH_Apply_pivots(), FLASH_Copy(), and FLASH_Trsm().
{ // Check parameters. if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING ) FLA_LU_piv_solve_check( A, p, B, X ); FLASH_Copy( B, X ); FLASH_Apply_pivots( FLA_LEFT, FLA_NO_TRANSPOSE, p, X ); FLASH_Trsm( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_UNIT_DIAG, FLA_ONE, A, X ); FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_NONUNIT_DIAG, FLA_ONE, A, X ); return FLA_SUCCESS; }