libflame  revision_anchor
Functions
FLASH_QR_UT_inc_solve.c File Reference

(r)

Functions

FLA_Error FLASH_QR_UT_inc_solve (FLA_Obj A, FLA_Obj TW, FLA_Obj B, FLA_Obj X)
 

Function Documentation

◆ FLASH_QR_UT_inc_solve()

FLA_Error FLASH_QR_UT_inc_solve ( FLA_Obj  A,
FLA_Obj  TW,
FLA_Obj  B,
FLA_Obj  X 
)
14 {
15  FLA_Obj W, Y;
16  FLA_Obj AT, AB;
17  FLA_Obj YT, YB;
18 
19  // Check parameters.
20  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
21  FLA_QR_UT_inc_solve_check( A, TW, B, X );
22 
24 
25  FLASH_Obj_create_copy_of( FLA_NO_TRANSPOSE, B, &Y );
26 
27  FLASH_Apply_Q_UT_inc( FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_FORWARD, FLA_COLUMNWISE,
28  A, TW, W, Y );
29 
30  // Create a temporary hierarchical view of only the top n-by-n part of A in
31  // case m > n so that AT captures the upper triangular factor R. We do the
32  // same for Y to ensure conformality.
33  FLASH_Part_create_2x1( A, &AT,
34  &AB, FLASH_Obj_scalar_width( A ), FLA_TOP );
35  FLASH_Part_create_2x1( Y, &YT,
36  &YB, FLASH_Obj_scalar_width( A ), FLA_TOP );
37 
38  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_NONUNIT_DIAG,
39  FLA_ONE, AT, YT );
40 
41  FLASH_Copy( YT, X );
42 
43  // Free the temporary hierarchical views.
45  &AB );
47  &YB );
48 
49  FLASH_Obj_free( &Y );
50  FLASH_Obj_free( &W );
51 
52  return FLA_SUCCESS;
53 }
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
FLA_Error FLASH_Apply_Q_UT_inc(FLA_Side side, FLA_Trans trans, FLA_Direct direct, FLA_Store storev, FLA_Obj A, FLA_Obj TW, FLA_Obj W1, FLA_Obj B)
Definition: FLASH_Apply_Q_UT_inc.c:18
void FLASH_Obj_free(FLA_Obj *H)
Definition: FLASH_Obj.c:638
FLA_Error FLASH_Obj_create_copy_of(FLA_Trans trans, FLA_Obj H_cur, FLA_Obj *H_new)
Definition: FLASH_Obj.c:561
dim_t FLASH_Obj_scalar_width(FLA_Obj H)
Definition: FLASH_View.c:641
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
FLA_Error FLASH_Part_free_2x1(FLA_Obj *AT, FLA_Obj *AB)
Definition: FLASH_View.c:572
FLA_Error FLASH_Apply_Q_UT_inc_create_workspace(FLA_Obj TW, FLA_Obj B, FLA_Obj *W)
Definition: FLASH_Apply_Q_UT_inc_create_workspace.c:13
FLA_Error FLA_QR_UT_inc_solve_check(FLA_Obj A, FLA_Obj TW, FLA_Obj B, FLA_Obj X)
Definition: FLA_QR_UT_inc_solve_check.c:13
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
Definition: FLA_type_defs.h:159

References FLA_Check_error_level(), FLA_ONE, FLA_QR_UT_inc_solve_check(), FLASH_Apply_Q_UT_inc(), FLASH_Apply_Q_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().