libflame  revision_anchor
Functions
FLASH_Apply_Q_UT_inc_create_workspace.c File Reference

(r)

Functions

FLA_Error FLASH_Apply_Q_UT_inc_create_workspace (FLA_Obj TW, FLA_Obj B, FLA_Obj *W)
 

Function Documentation

◆ FLASH_Apply_Q_UT_inc_create_workspace()

FLA_Error FLASH_Apply_Q_UT_inc_create_workspace ( FLA_Obj  TW,
FLA_Obj  B,
FLA_Obj W 
)

References FLA_Abort(), FLA_Obj_datatype(), FLA_Obj_width(), FLA_Print_message(), FLASH_Obj_create_ext(), FLASH_Obj_depth(), FLASH_Obj_scalar_length_tl(), and FLASH_Obj_scalar_width_tl().

Referenced by FLASH_QR_UT_inc_solve().

14 {
15  FLA_Datatype datatype;
16  dim_t depth;
17  dim_t b_alg;
18  dim_t b_flash;
19  dim_t m, n;
20 
21  // Query the depth.
22  depth = FLASH_Obj_depth( TW );
23 
24  // *** The current Apply_Q_UT_inc algorithm implemented assumes that
25  // the matrix has a hierarchical depth of 1. We check for that here
26  // because we anticipate that we'll use a more general algorithm in the
27  // future, and we don't want to forget to remove the constraint. ***
28  if ( depth != 1 )
29  {
30  FLA_Print_message( "FLASH_Apply_Q_UT_inc() currently only supports matrices of depth 1",
31  __FILE__, __LINE__ );
32  FLA_Abort();
33  }
34 
35  // Query the datatype of matrix TW.
36  datatype = FLA_Obj_datatype( TW );
37 
38  // Inspect the length of a the top-left element of TW to get the
39  // algorithmic blocksize we'll use throughout the Apply_Q_UT_inc
40  // algorithm.
41  b_alg = FLASH_Obj_scalar_length_tl( TW );
42 
43  // The width of the top-left element gives us the storage blocksize.
44  b_flash = FLASH_Obj_scalar_width_tl( TW );
45 
46  // The element length of W is 1.
47  m = 1;
48 
49  // Query the element (not scalar) width of the right-hand side
50  // matrix B. This is done so we can create W with full blocks for the
51  // right "edge cases" of B.
52  n = FLA_Obj_width( B );
53 
54  // Create hierarchical matrix W.
55  FLASH_Obj_create_ext( datatype, m * b_alg, n * b_flash,
56  depth, &b_alg, &b_flash,
57  W );
58 
59  return FLA_SUCCESS;
60 }
unsigned long dim_t
Definition: FLA_type_defs.h:71
dim_t FLASH_Obj_depth(FLA_Obj H)
Definition: FLASH_Obj.c:20
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_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
dim_t FLASH_Obj_scalar_width_tl(FLA_Obj H)
Definition: FLASH_View.c:737
void FLA_Abort(void)
Definition: FLA_Error.c:248
void FLA_Print_message(char *str, char *file, int line)
Definition: FLA_Error.c:234
int FLA_Datatype
Definition: FLA_type_defs.h:49
dim_t FLASH_Obj_scalar_length_tl(FLA_Obj H)
Definition: FLASH_View.c:723