libflame  revision_anchor
Functions
FLASH_CAQR_UT_inc_solve.c File Reference

(r)

Functions

FLA_Error FLASH_CAQR_UT_inc_solve (dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW, FLA_Obj B, FLA_Obj X)
 

Function Documentation

◆ FLASH_CAQR_UT_inc_solve()

FLA_Error FLASH_CAQR_UT_inc_solve ( dim_t  p,
FLA_Obj  A,
FLA_Obj  ATW,
FLA_Obj  R,
FLA_Obj  RTW,
FLA_Obj  B,
FLA_Obj  X 
)

References FLA_CAQR_UT_inc_solve_check(), FLA_Check_error_level(), FLA_ONE, FLASH_Apply_CAQ_UT_inc(), FLASH_Apply_CAQ_UT_inc_create_workspace(), FLASH_Copy(), FLASH_Obj_create_copy_of(), FLASH_Obj_free(), FLASH_Obj_scalar_width(), FLASH_Part_create_2x1(), FLASH_Part_free_2x1(), and FLASH_Trsm().

14 {
15  FLA_Obj W, Y;
16  FLA_Obj RT, RB;
17  FLA_Obj YT, YB;
18 
19  // Check parameters.
20  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
21  FLA_CAQR_UT_inc_solve_check( p, A, ATW, R, RTW, B, X );
22 
24 
25  FLASH_Obj_create_copy_of( FLA_NO_TRANSPOSE, B, &Y );
26 
27  // Create a temporary hierarchical view of only the top n-by-n part of R in
28  // case m > n so that RT captures the upper triangle. We do the same for Y
29  // to ensure conformality.
30  FLASH_Part_create_2x1( R, &RT,
31  &RB, FLASH_Obj_scalar_width( R ), FLA_TOP );
32  FLASH_Part_create_2x1( Y, &YT,
33  &YB, FLASH_Obj_scalar_width( R ), FLA_TOP );
34 
36  FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_FORWARD, FLA_COLUMNWISE,
37  A, ATW, R, RTW, W, Y );
38 
39  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_NONUNIT_DIAG,
40  FLA_ONE, RT, YT );
41 
42  FLASH_Copy( YT, X );
43 
44  // Free the temporary hierarchical views.
46  &RB );
48  &YB );
49 
50  FLASH_Obj_free( &Y );
51  FLASH_Obj_free( &W );
52 
53  return FLA_SUCCESS;
54 }
FLA_Error FLASH_Apply_CAQ_UT_inc(dim_t p, FLA_Side side, FLA_Trans trans, FLA_Direct direct, FLA_Store storev, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW, FLA_Obj W, FLA_Obj B)
Definition: FLASH_Apply_CAQ_UT_inc.c:15
FLA_Error FLASH_Obj_create_copy_of(FLA_Trans trans, FLA_Obj H_cur, FLA_Obj *H_new)
Definition: FLASH_Obj.c:561
FLA_Error FLASH_Apply_CAQ_UT_inc_create_workspace(dim_t p, FLA_Obj TW, FLA_Obj B, FLA_Obj *W)
Definition: FLASH_Apply_CAQ_UT_inc_create_workspace.c:13
FLA_Error FLA_CAQR_UT_inc_solve_check(dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW, FLA_Obj B, FLA_Obj X)
Definition: FLA_CAQR_UT_inc_solve_check.c:13
void FLASH_Obj_free(FLA_Obj *H)
Definition: FLASH_Obj.c:638
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
Definition: FLA_type_defs.h:158
FLA_Error FLASH_Part_create_2x1(FLA_Obj A, FLA_Obj *AT, FLA_Obj *AB, dim_t n_rows, FLA_Side side)
Definition: FLASH_View.c:13
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
FLA_Error FLASH_Part_free_2x1(FLA_Obj *AT, FLA_Obj *AB)
Definition: FLASH_View.c:572
dim_t FLASH_Obj_scalar_width(FLA_Obj H)
Definition: FLASH_View.c:641
FLA_Error FLASH_Copy(FLA_Obj A, FLA_Obj B)
Definition: FLASH_Copy.c:15
FLA_Error FLASH_Trsm(FLA_Side side, FLA_Uplo uplo, FLA_Trans trans, FLA_Diag diag, FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
Definition: FLASH_Trsm.c:15