libflame  revision_anchor
Functions
FLA_CAQR_UT_inc_factorize_panels.c File Reference

(r)

Functions

FLA_Error FLA_CAQR_UT_inc_factorize_panels (dim_t nb_part, FLA_Obj A, FLA_Obj TW)
 

Function Documentation

◆ FLA_CAQR_UT_inc_factorize_panels()

FLA_Error FLA_CAQR_UT_inc_factorize_panels ( dim_t  nb_part,
FLA_Obj  A,
FLA_Obj  TW 
)
14 {
15  FLA_Obj AT, A0,
16  AB, A1,
17  A2;
18 
19  FLA_Obj TWT, TW0,
20  TWB, TW1,
21  TW2;
22 
23  dim_t b;
24 
25  FLA_Part_2x1( A, &AT,
26  &AB, 0, FLA_TOP );
27 
28  FLA_Part_2x1( TW, &TWT,
29  &TWB, 0, FLA_TOP );
30 
31  while ( FLA_Obj_length( AB ) > 0 ){
32 
33  b = min( nb_part, FLA_Obj_length( AB ) );
34 
35  FLA_Repart_2x1_to_3x1( AT, &A0,
36  /* ** */ /* ** */
37  &A1,
38  AB, &A2, b, FLA_BOTTOM );
39 
40  FLA_Repart_2x1_to_3x1( TWT, &TW0,
41  /* ** */ /* ** */
42  &TW1,
43  TWB, &TW2, b, FLA_BOTTOM );
44 
45  /*------------------------------------------------------------*/
46 
47  // Perform an incremental QR factorization on A1, writing triangular
48  // block Householder factors to T in TW1.
49  FLASH_QR_UT_inc( A1, TW1 );
50 
51  /*------------------------------------------------------------*/
52 
53  FLA_Cont_with_3x1_to_2x1( &AT, A0,
54  A1,
55  /* ** */ /* ** */
56  &AB, A2, FLA_TOP );
57 
58  FLA_Cont_with_3x1_to_2x1( &TWT, TW0,
59  TW1,
60  /* ** */ /* ** */
61  &TWB, TW2, FLA_TOP );
62  }
63 
64  return FLA_SUCCESS;
65 }
FLA_Error FLASH_QR_UT_inc(FLA_Obj A, FLA_Obj TW)
Definition: FLASH_QR_UT_inc.c:13
FLA_Error FLA_Cont_with_3x1_to_2x1(FLA_Obj *AT, FLA_Obj A0, FLA_Obj A1, FLA_Obj *AB, FLA_Obj A2, FLA_Side side)
Definition: FLA_View.c:428
FLA_Error FLA_Repart_2x1_to_3x1(FLA_Obj AT, FLA_Obj *A0, FLA_Obj *A1, FLA_Obj AB, FLA_Obj *A2, dim_t mb, FLA_Side side)
Definition: FLA_View.c:226
FLA_Error FLA_Part_2x1(FLA_Obj A, FLA_Obj *A1, FLA_Obj *A2, dim_t mb, FLA_Side side)
Definition: FLA_View.c:76
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
unsigned long dim_t
Definition: FLA_type_defs.h:71
Definition: FLA_type_defs.h:159

References FLA_Cont_with_3x1_to_2x1(), FLA_Obj_length(), FLA_Part_2x1(), FLA_Repart_2x1_to_3x1(), and FLASH_QR_UT_inc().

Referenced by FLASH_CAQR_UT_inc_noopt().