libflame  revision_anchor
Functions
FLA_CAQR_UT_inc.h File Reference

(r)

Go to the source code of this file.

Functions

FLA_Error FLASH_CAQR_UT_inc (dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW)
 
FLA_Error FLASH_CAQR_UT_inc_noopt (dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW)
 
FLA_Error FLASH_CAQR_UT_inc_create_hier_matrices (dim_t p, FLA_Obj A_flat, dim_t depth, dim_t *b_flash, dim_t b_alg, FLA_Obj *A, FLA_Obj *ATW, FLA_Obj *R, FLA_Obj *RTW)
 
dim_t FLASH_CAQR_UT_inc_determine_alg_blocksize (FLA_Obj A)
 
FLA_Error FLASH_CAQR_UT_inc_adjust_views (FLA_Obj A, FLA_Obj TW)
 
void FLA_CAQR_UT_inc_init_structure (dim_t p, dim_t nb_part, FLA_Obj R)
 
dim_t FLA_CAQR_UT_inc_compute_blocks_per_part (dim_t p, FLA_Obj A)
 
FLA_Error FLA_CAQR_UT_inc_factorize_panels (dim_t nb_part, FLA_Obj A, FLA_Obj ATW)
 
FLA_Error FLA_CAQR_UT_inc_copy_triangles (dim_t nb_part, FLA_Obj A, FLA_Obj R)
 
FLA_Error FLA_CAQR_UT_inc_blk_var1 (FLA_Obj R, FLA_Obj TW, fla_caqrutinc_t *cntl)
 
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

◆ FLA_CAQR_UT_inc_blk_var1()

FLA_Error FLA_CAQR_UT_inc_blk_var1 ( FLA_Obj  R,
FLA_Obj  TW,
fla_caqrutinc_t cntl 
)
14 {
15  FLA_Obj ATL, ATR, A00, A01, A02,
16  ABL, ABR, A10, A11, A12,
17  A20, A21, A22;
18 
19  FLA_Obj TTL, WTR, T00, W01, W02,
20  TBL, TBR, T10, T11, W12,
21  T20, T21, T22;
22 
23  dim_t b;
24 
25  FLA_Part_2x2( A, &ATL, &ATR,
26  &ABL, &ABR, 0, 0, FLA_TL );
27 
28  FLA_Part_2x2( TW, &TTL, &WTR,
29  &TBL, &TBR, 0, 0, FLA_TL );
30 
31  while ( FLA_Obj_min_dim( ABR ) > 0 ){
32 
33  b = FLA_Determine_blocksize( ABR, FLA_BR, FLA_Cntl_blocksize( cntl ) );
34 
35  FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &A01, &A02,
36  /* ************* */ /* ******************** */
37  &A10, /**/ &A11, &A12,
38  ABL, /**/ ABR, &A20, /**/ &A21, &A22,
39  b, b, FLA_BR );
40 
41  FLA_Repart_2x2_to_3x3( TTL, /**/ WTR, &T00, /**/ &W01, &W02,
42  /* ************* */ /* ******************** */
43  &T10, /**/ &T11, &W12,
44  TBL, /**/ TBR, &T20, /**/ &T21, &T22,
45  b, b, FLA_BR );
46 
47  /*------------------------------------------------------------*/
48 
50  A21, T21,
51  FLA_Cntl_sub_caqr2ut( cntl ) );
52 
53 
54  if ( FLA_Obj_width( A12 ) > 0 )
55  {
56  FLA_Apply_CAQ2_UT_internal( FLA_LEFT, FLA_CONJ_TRANSPOSE, FLA_FORWARD, FLA_COLUMNWISE,
57  A21, T21, W12, A12,
58  A22,
59  FLA_Cntl_sub_apcaq2ut( cntl ) );
60  }
61 
62  /*------------------------------------------------------------*/
63 
64  FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, A01, /**/ A02,
65  A10, A11, /**/ A12,
66  /* ************** */ /* ****************** */
67  &ABL, /**/ &ABR, A20, A21, /**/ A22,
68  FLA_TL );
69 
70  FLA_Cont_with_3x3_to_2x2( &TTL, /**/ &WTR, T00, W01, /**/ W02,
71  T10, T11, /**/ W12,
72  /* ************** */ /* ****************** */
73  &TBL, /**/ &TBR, T20, T21, /**/ T22,
74  FLA_TL );
75 
76  }
77 
78  return FLA_SUCCESS;
79 }
FLA_Error FLA_Apply_CAQ2_UT_internal(FLA_Side side, FLA_Trans trans, FLA_Direct direct, FLA_Store storev, FLA_Obj D, FLA_Obj T, FLA_Obj W, FLA_Obj C, FLA_Obj E, fla_apcaq2ut_t *cntl)
Definition: FLA_Apply_CAQ2_UT_internal.c:17
FLA_Error FLA_CAQR2_UT_internal(FLA_Obj U, FLA_Obj D, FLA_Obj T, fla_caqr2ut_t *cntl)
Definition: FLA_CAQR2_UT_internal.c:16
FLA_Error FLA_Cont_with_3x3_to_2x2(FLA_Obj *ATL, FLA_Obj *ATR, FLA_Obj A00, FLA_Obj A01, FLA_Obj A02, FLA_Obj A10, FLA_Obj A11, FLA_Obj A12, FLA_Obj *ABL, FLA_Obj *ABR, FLA_Obj A20, FLA_Obj A21, FLA_Obj A22, FLA_Quadrant quadrant)
Definition: FLA_View.c:304
FLA_Error FLA_Part_2x2(FLA_Obj A, FLA_Obj *A11, FLA_Obj *A12, FLA_Obj *A21, FLA_Obj *A22, dim_t mb, dim_t nb, FLA_Quadrant quadrant)
Definition: FLA_View.c:17
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
FLA_Error FLA_Repart_2x2_to_3x3(FLA_Obj ATL, FLA_Obj ATR, FLA_Obj *A00, FLA_Obj *A01, FLA_Obj *A02, FLA_Obj *A10, FLA_Obj *A11, FLA_Obj *A12, FLA_Obj ABL, FLA_Obj ABR, FLA_Obj *A20, FLA_Obj *A21, FLA_Obj *A22, dim_t mb, dim_t nb, FLA_Quadrant quadrant)
Definition: FLA_View.c:142
dim_t FLA_Determine_blocksize(FLA_Obj A_unproc, FLA_Quadrant to_dir, fla_blocksize_t *cntl_blocksizes)
Definition: FLA_Blocksize.c:234
dim_t FLA_Obj_min_dim(FLA_Obj obj)
Definition: FLA_Query.c:153
unsigned long dim_t
Definition: FLA_type_defs.h:71
Definition: FLA_type_defs.h:159

References FLA_Apply_CAQ2_UT_internal(), FLA_CAQR2_UT_internal(), FLA_Cont_with_3x3_to_2x2(), FLA_Determine_blocksize(), FLA_Obj_min_dim(), FLA_Obj_width(), FLA_Part_2x2(), and FLA_Repart_2x2_to_3x3().

Referenced by FLASH_CAQR_UT_inc_noopt().

◆ FLA_CAQR_UT_inc_compute_blocks_per_part()

dim_t FLA_CAQR_UT_inc_compute_blocks_per_part ( dim_t  p,
FLA_Obj  A 
)
14 {
15  dim_t nb_part;
16  dim_t nb_left;
17  dim_t num_blocks;
18 
19  // Query the element (not scalar) length of A.
20  num_blocks = FLA_Obj_length( A );
21 
22  // Compute the number of blocks per partitions.
23  nb_part = num_blocks / p;
24  nb_left = num_blocks % p;
25 
26  // If there are leftover blocks, increase nb_part by one.
27  if ( nb_left > 0 ) nb_part += 1;
28 
29  return nb_part;
30 }
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116

References FLA_Obj_length().

Referenced by FLASH_Apply_CAQ_UT_inc(), FLASH_CAQR_UT_inc_create_hier_matrices(), and FLASH_CAQR_UT_inc_noopt().

◆ FLA_CAQR_UT_inc_copy_triangles()

FLA_Error FLA_CAQR_UT_inc_copy_triangles ( dim_t  nb_part,
FLA_Obj  A,
FLA_Obj  R 
)
14 {
15  FLA_Obj AT, A0,
16  AB, A1,
17  A2;
18 
19  FLA_Obj RT, R0,
20  RB, R1,
21  R2;
22 
23  dim_t b;
24 
25  FLA_Part_2x1( A, &AT,
26  &AB, 0, FLA_TOP );
27 
28  FLA_Part_2x1( R, &RT,
29  &RB, 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( RT, &R0,
41  /* ** */ /* ** */
42  &R1,
43  RB, &R2, b, FLA_BOTTOM );
44 
45  /*------------------------------------------------------------*/
46 
47  // Copy the individual upper triangles in A into R.
48  FLASH_Copyr( FLA_UPPER_TRIANGULAR, A1, R1 );
49 
50  /*------------------------------------------------------------*/
51 
52  FLA_Cont_with_3x1_to_2x1( &AT, A0,
53  A1,
54  /* ** */ /* ** */
55  &AB, A2, FLA_TOP );
56 
57  FLA_Cont_with_3x1_to_2x1( &RT, R0,
58  R1,
59  /* ** */ /* ** */
60  &RB, R2, FLA_TOP );
61  }
62 
63  return FLA_SUCCESS;
64 }
FLA_Error FLASH_Copyr(FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLASH_Copyr.c:15
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

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

Referenced by FLASH_CAQR_UT_inc_noopt().

◆ FLA_CAQR_UT_inc_factorize_panels()

FLA_Error FLA_CAQR_UT_inc_factorize_panels ( dim_t  nb_part,
FLA_Obj  A,
FLA_Obj  ATW 
)
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

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().

◆ FLA_CAQR_UT_inc_init_structure()

void FLA_CAQR_UT_inc_init_structure ( dim_t  p,
dim_t  nb_part,
FLA_Obj  R 
)
14 {
15  dim_t m, n;
16  dim_t rs, cs;
17  dim_t i, j, ip;
18  FLA_Obj* buff_R;
19 
20  m = FLA_Obj_length( R );
21  n = FLA_Obj_width( R );
22  rs = FLA_Obj_row_stride( R );
23  cs = FLA_Obj_col_stride( R );
24  buff_R = FLA_Obj_buffer_at_view( R );
25 
26  // Fill in R by row panels.
27  for ( ip = 0; ip < p; ++ip )
28  {
29  FLA_Obj* buff_R1 = buff_R + (ip*nb_part)*rs;
30 
31  int m_behind = ip*nb_part;
32  int m_ahead = m - m_behind;
33 
34  int m_cur = min( nb_part, m_ahead );
35  int n_cur = n;
36 
37  // Iterate across columns for the current panel.
38  for ( j = 0; j < n_cur; ++j )
39  {
40  FLA_Obj* rho = buff_R1 + j*cs;
41 
42  // Mark the above-diagonal blocks as full.
43  for ( i = 0; i < j; ++i )
44  {
45  rho->base->uplo = FLA_FULL_MATRIX;
46  rho += rs;
47  }
48 
49  // Mark the diagonal block as triangular.
50  rho->base->uplo = FLA_UPPER_TRIANGULAR;
51  rho += rs;
52 
53  // Mark the below-diagonal blocks as zero.
54  for ( i = j + 1; i < m_cur; ++i )
55  {
56  rho->base->uplo = FLA_ZERO_MATRIX;
57  rho += rs;
58  }
59  }
60  }
61 }
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
void * FLA_Obj_buffer_at_view(FLA_Obj obj)
Definition: FLA_Query.c:215
int i
Definition: bl1_axmyv2.c:145
* rho
Definition: bl1_axpyv2bdotaxpy.c:322

References FLA_Obj_buffer_at_view(), FLA_Obj_col_stride(), FLA_Obj_length(), FLA_Obj_row_stride(), FLA_Obj_width(), i, and rho.

Referenced by FLASH_CAQR_UT_inc_create_hier_matrices().

◆ FLASH_CAQR_UT_inc()

FLA_Error FLASH_CAQR_UT_inc ( dim_t  p,
FLA_Obj  A,
FLA_Obj  ATW,
FLA_Obj  R,
FLA_Obj  RTW 
)
14 {
15  FLA_Error r_val;
16 
17  //if ( FLASH_Queue_stack_depth() == 0 )
18  // r_val = FLASH_CAQR_UT_inc_opt1( A, ATW, R, RTW );
19  //else
20  r_val = FLASH_CAQR_UT_inc_noopt( p, A, ATW, R, RTW );
21 
22  return r_val;
23 }
FLA_Error FLASH_CAQR_UT_inc_noopt(dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW)
Definition: FLASH_CAQR_UT_inc_noopt.c:15
int FLA_Error
Definition: FLA_type_defs.h:47

References FLASH_CAQR_UT_inc_noopt().

◆ FLASH_CAQR_UT_inc_adjust_views()

FLA_Error FLASH_CAQR_UT_inc_adjust_views ( FLA_Obj  A,
FLA_Obj  TW 
)
76 {
77  dim_t b_flash;
78  dim_t n, n_last;
79 
80  // We can query b_flash as the width of the top-left element of TW.
81  b_flash = FLASH_Obj_scalar_width_tl( TW );
82 
83  // Query the element (not scalar) n dimension of A.
84  n = FLA_Obj_width( A );
85 
86  // If the bottom-right-most block along the diagonal is a partial block,
87  // adjust the view of the corresponding T block.
88  n_last = FLASH_Obj_scalar_width( A ) % b_flash;
89 
90  if ( n_last > 0 )
91  {
92  FLA_Obj TWTL, TWTR,
93  TWBL, TWBR;
94  FLA_Obj TWL, TWR;
95  FLA_Obj TWT, TW0,
96  TWB, TW1,
97  TW2;
98  FLA_Obj* TW1p;
99 
100  FLA_Part_2x2( TW, &TWTL, &TWTR,
101  &TWBL, &TWBR, n-1, n-1, FLA_TL );
102 
103  FLA_Part_2x1( TWBR, &TWT,
104  &TWB, 0, FLA_TOP );
105 
106  while ( FLA_Obj_length( TWB ) > 0 )
107  {
108  FLA_Repart_2x1_to_3x1( TWT, &TW0,
109  /* *** */ /* *** */
110  &TW1,
111  TWB, &TW2, 1, FLA_BOTTOM );
112 
113  // -----------------------------------------------------------
114 
115  TW1p = FLASH_OBJ_PTR_AT( TW1 );
116 
117  FLA_Part_1x2( *TW1p, &TWL, &TWR, n_last, FLA_LEFT );
118 
119  *TW1p = TWL;
120  TW1p->m_inner = TW1p->m;
121  TW1p->n_inner = TW1p->n;
122 
123  // -----------------------------------------------------------
124 
125  FLA_Cont_with_3x1_to_2x1( &TWT, TW0,
126  TW1,
127  /* *** */ /* *** */
128  &TWB, TW2, FLA_TOP );
129  }
130  }
131 
132  return FLA_SUCCESS;
133 }
dim_t FLASH_Obj_scalar_width(FLA_Obj H)
Definition: FLASH_View.c:641
dim_t FLASH_Obj_scalar_width_tl(FLA_Obj H)
Definition: FLASH_View.c:737
FLA_Error FLA_Part_1x2(FLA_Obj A, FLA_Obj *A1, FLA_Obj *A2, dim_t nb, FLA_Side side)
Definition: FLA_View.c:110
dim_t n_inner
Definition: FLA_type_defs.h:166
dim_t m
Definition: FLA_type_defs.h:163
dim_t m_inner
Definition: FLA_type_defs.h:165
dim_t n
Definition: FLA_type_defs.h:164

References FLA_Cont_with_3x1_to_2x1(), FLA_Obj_length(), FLA_Obj_width(), FLA_Part_1x2(), FLA_Part_2x1(), FLA_Part_2x2(), FLA_Repart_2x1_to_3x1(), FLASH_Obj_scalar_width(), FLASH_Obj_scalar_width_tl(), FLA_Obj_view::m, FLA_Obj_view::m_inner, FLA_Obj_view::n, and FLA_Obj_view::n_inner.

Referenced by FLASH_CAQR_UT_inc_create_hier_matrices().

◆ FLASH_CAQR_UT_inc_create_hier_matrices()

FLA_Error FLASH_CAQR_UT_inc_create_hier_matrices ( dim_t  p,
FLA_Obj  A_flat,
dim_t  depth,
dim_t b_flash,
dim_t  b_alg,
FLA_Obj A,
FLA_Obj ATW,
FLA_Obj R,
FLA_Obj RTW 
)
14 {
15  FLA_Datatype datatype;
16  dim_t m, n;
17  dim_t nb_part;
18 
19  // *** The current CAQR_UT_inc algorithm implemented assumes that
20  // the matrix has a hierarchical depth of 1.
21  if ( depth != 1 )
22  {
23  FLA_Print_message( "FLASH_CAQR_UT_inc() currently only supports matrices of depth 1",
24  __FILE__, __LINE__ );
25  FLA_Abort();
26  }
27 
28  // Create hierarchical copy of matrix A_flat.
29  FLASH_Obj_create_hier_copy_of_flat( A_flat, depth, b_flash, A );
30 
31  // Create hierarchical copy of matrix A_flat.
32  FLASH_Obj_create_conf_to( FLA_NO_TRANSPOSE, *A, R );
33 
34  // Query the datatype of matrix A_flat.
35  datatype = FLA_Obj_datatype( A_flat );
36 
37  // If the user passed in zero for b_alg, then we need to set the
38  // algorithmic (inner) blocksize to a reasonable default value.
39  if ( b_alg == 0 )
40  {
42  }
43 
44  // Query the element (not scalar) dimensions of the new hierarchical
45  // matrix. This is done so we can create T with full blocks for the
46  // bottom and right "edge cases" of A.
47  m = FLA_Obj_length( *A );
48  n = FLA_Obj_width( *A );
49 
50  // Create hierarchical matrices T and W for both A and R. T is lower
51  // triangular where each block is b_alg-by-b_flash and W is strictly
52  // upper triangular where each block is b_alg-by-b_flash. So we can
53  // create them simultaneously as part of the same hierarchical matrix.
54  FLASH_Obj_create_ext( datatype, m * b_alg, n * b_flash[0],
55  depth, &b_alg, b_flash,
56  ATW );
57  FLASH_Obj_create_ext( datatype, m * b_alg, n * b_flash[0],
58  depth, &b_alg, b_flash,
59  RTW );
60 
61  // If the bottom-right-most block along the diagonal is a partial block,
62  // adjust the view of the corresponding T block.
65 
66  // Compute the partition length from the number of partitions.
67  nb_part = FLA_CAQR_UT_inc_compute_blocks_per_part( p, *A );
68 
69  // Encode block structure (upper tri, full, or zero) into blocks of R.
70  FLA_CAQR_UT_inc_init_structure( p, nb_part, *R );
71 
72  return FLA_SUCCESS;
73 }
dim_t FLASH_CAQR_UT_inc_determine_alg_blocksize(FLA_Obj A)
Definition: FLASH_CAQR_UT_inc_create_hier_matrices.c:136
FLA_Error FLASH_CAQR_UT_inc_adjust_views(FLA_Obj A, FLA_Obj TW)
Definition: FLASH_CAQR_UT_inc_create_hier_matrices.c:75
FLA_Error FLASH_Obj_create_ext(FLA_Datatype datatype, dim_t m, dim_t n, dim_t depth, dim_t *b_m, dim_t *b_n, FLA_Obj *H)
Definition: FLASH_Obj.c:151
FLA_Error FLASH_Obj_create_conf_to(FLA_Trans trans, FLA_Obj H_cur, FLA_Obj *H_new)
Definition: FLASH_Obj.c:406
FLA_Error FLASH_Obj_create_hier_copy_of_flat(FLA_Obj F, dim_t depth, dim_t *b_mn, FLA_Obj *H)
Definition: FLASH_Obj.c:591
void FLA_CAQR_UT_inc_init_structure(dim_t p, dim_t nb_part, FLA_Obj R)
Definition: FLA_CAQR_UT_inc_init_structure.c:13
dim_t FLA_CAQR_UT_inc_compute_blocks_per_part(dim_t p, FLA_Obj A)
Definition: FLA_CAQR_UT_inc_compute_p_length.c:13
void FLA_Abort(void)
Definition: FLA_Error.c:248
void FLA_Print_message(char *str, char *file, int line)
Definition: FLA_Error.c:234
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
int FLA_Datatype
Definition: FLA_type_defs.h:49

References FLA_Abort(), FLA_CAQR_UT_inc_compute_blocks_per_part(), FLA_CAQR_UT_inc_init_structure(), FLA_Obj_datatype(), FLA_Obj_length(), FLA_Obj_width(), FLA_Print_message(), FLASH_CAQR_UT_inc_adjust_views(), FLASH_CAQR_UT_inc_determine_alg_blocksize(), FLASH_Obj_create_conf_to(), FLASH_Obj_create_ext(), and FLASH_Obj_create_hier_copy_of_flat().

◆ FLASH_CAQR_UT_inc_determine_alg_blocksize()

dim_t FLASH_CAQR_UT_inc_determine_alg_blocksize ( FLA_Obj  A)
137 {
138  dim_t b_alg;
139  dim_t b_flash;
140 
141  // Acquire the storage blocksize.
142  b_flash = FLA_Obj_length( *FLASH_OBJ_PTR_AT( A ) );
143 
144  // Scale the storage blocksize by a pre-defined scalar to arrive at a
145  // reasonable algorithmic blocksize, but make sure it's at least 1.
146  b_alg = ( dim_t ) max( ( double ) b_flash * FLA_CAQR_INNER_TO_OUTER_B_RATIO, 1 );
147 
148  return b_alg;
149 }

References FLA_Obj_length().

Referenced by FLASH_CAQR_UT_inc_create_hier_matrices().

◆ FLASH_CAQR_UT_inc_noopt()

FLA_Error FLASH_CAQR_UT_inc_noopt ( dim_t  p,
FLA_Obj  A,
FLA_Obj  ATW,
FLA_Obj  R,
FLA_Obj  RTW 
)
16 {
17  FLA_Error r_val = FLA_SUCCESS;
18  dim_t nb_part;
19 
20  // Check parameters.
21  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
22  FLA_CAQR_UT_inc_check( p, A, ATW, R, RTW );
23 
24  // Compute the partition length from the number of partitions.
26 
27  // Begin a parallel region.
29 
30  // Perform incremental QR's on each of the p partitions.
31  FLA_CAQR_UT_inc_factorize_panels( nb_part, A, ATW );
32 
33  // Copy the triangles of A into R.
34  FLA_CAQR_UT_inc_copy_triangles( nb_part, A, R );
35 
36  // Perform an incremental CAQR on the resulting upper triangular R's in A.
38 
39  // End the parallel region.
41 
42  return r_val;
43 }
fla_caqrutinc_t * flash_caqrutinc_cntl
Definition: FLASH_CAQR_UT_inc_cntl_init.c:16
void FLASH_Queue_begin(void)
Definition: FLASH_Queue.c:59
void FLASH_Queue_end(void)
Definition: FLASH_Queue.c:81
FLA_Error FLA_CAQR_UT_inc_blk_var1(FLA_Obj R, FLA_Obj TW, fla_caqrutinc_t *cntl)
Definition: FLA_CAQR_UT_inc_blk_var1.c:13
FLA_Error FLA_CAQR_UT_inc_copy_triangles(dim_t nb_part, FLA_Obj A, FLA_Obj R)
Definition: FLA_CAQR_UT_inc_copy_triangles.c:13
FLA_Error FLA_CAQR_UT_inc_factorize_panels(dim_t nb_part, FLA_Obj A, FLA_Obj ATW)
Definition: FLA_CAQR_UT_inc_factorize_panels.c:13
FLA_Error FLA_CAQR_UT_inc_check(dim_t p, FLA_Obj A, FLA_Obj ATW, FLA_Obj R, FLA_Obj RTW)
Definition: FLA_CAQR_UT_inc_check.c:13
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18

References FLA_CAQR_UT_inc_blk_var1(), FLA_CAQR_UT_inc_check(), FLA_CAQR_UT_inc_compute_blocks_per_part(), FLA_CAQR_UT_inc_copy_triangles(), FLA_CAQR_UT_inc_factorize_panels(), FLA_Check_error_level(), flash_caqrutinc_cntl, FLASH_Queue_begin(), and FLASH_Queue_end().

Referenced by FLASH_CAQR_UT_inc().

◆ 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 
)
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_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
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
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_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_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
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18

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().