libflame revision_anchor
Functions
FLASH_Copy_other.c File Reference

(r)

Functions

FLA_Error FLASH_Copy_buffer_to_hier (dim_t m, dim_t n, void *buffer, dim_t rs, dim_t cs, dim_t i, dim_t j, FLA_Obj H)
 
FLA_Error FLASH_Copy_hier_to_buffer (dim_t i, dim_t j, FLA_Obj H, dim_t m, dim_t n, void *buffer, dim_t rs, dim_t cs)
 
FLA_Error FLASH_Copy_flat_to_hier (FLA_Obj F, dim_t i, dim_t j, FLA_Obj H)
 
FLA_Error FLASH_Copy_hier_to_flat (dim_t i, dim_t j, FLA_Obj H, FLA_Obj F)
 
FLA_Error FLASH_Copy_hierarchy (int direction, FLA_Obj F, FLA_Obj *H)
 

Function Documentation

◆ FLASH_Copy_buffer_to_hier()

FLA_Error FLASH_Copy_buffer_to_hier ( dim_t m,
dim_t n,
void * buffer,
dim_t rs,
dim_t cs,
dim_t i,
dim_t j,
FLA_Obj H )
14{
15 FLA_Obj flat_matrix;
16 FLA_Datatype datatype;
17 FLA_Error e_val;
18
19 if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
20 {
21 e_val = FLA_Check_matrix_strides( m, n, rs, cs );
22 FLA_Check_error_code( e_val );
23
24 e_val = FLA_Check_submatrix_dims_and_offset( m, n, i, j, H );
25 FLA_Check_error_code( e_val );
26 }
27
28 // Acquire the datatype from the hierarchical matrix object.
29 datatype = FLASH_Obj_datatype( H );
30
31 // Create a temporary conventional matrix object of the requested datatype
32 // and dimensions and attach the given buffer containing the incoming data.
33 FLA_Obj_create_without_buffer( datatype, m, n, &flat_matrix );
34 FLA_Obj_attach_buffer( buffer, rs, cs, &flat_matrix );
35
36 // Recurse through H, adding in the corresponding elements of flat_matrix,
37 // starting at the (i,j) element offset.
38 FLASH_Copy_flat_to_hier( flat_matrix, i, j, H );
39
40 // Free the object (but don't free the buffer!).
41 FLA_Obj_free_without_buffer( &flat_matrix );
42
43 return FLA_SUCCESS;
44}
FLASH_Copy_flat_to_hier
FLA_Error FLASH_Copy_flat_to_hier(FLA_Obj F, dim_t i, dim_t j, FLA_Obj H)
Definition FLASH_Copy_other.c:81
FLASH_Obj_datatype
FLA_Datatype FLASH_Obj_datatype(FLA_Obj H)
Definition FLASH_Obj.c:14
FLA_Check_matrix_strides
FLA_Error FLA_Check_matrix_strides(dim_t m, dim_t n, dim_t rs, dim_t cs)
Definition FLA_Check.c:1146
FLA_Obj_create_without_buffer
FLA_Error FLA_Obj_create_without_buffer(FLA_Datatype datatype, dim_t m, dim_t n, FLA_Obj *obj)
Definition FLA_Obj.c:362
FLA_Obj_free_without_buffer
FLA_Error FLA_Obj_free_without_buffer(FLA_Obj *obj)
Definition FLA_Obj.c:615
FLA_Obj_attach_buffer
FLA_Error FLA_Obj_attach_buffer(void *buffer, dim_t rs, dim_t cs, FLA_Obj *obj)
Definition FLA_Obj.c:522
FLA_Check_submatrix_dims_and_offset
FLA_Error FLA_Check_submatrix_dims_and_offset(dim_t m, dim_t n, dim_t i, dim_t j, FLA_Obj A)
Definition FLA_Check.c:834
FLA_Check_error_level
unsigned int FLA_Check_error_level(void)
Definition FLA_Check.c:18
FLA_Error
int FLA_Error
Definition FLA_type_defs.h:47
FLA_Datatype
int FLA_Datatype
Definition FLA_type_defs.h:49
i
int i
Definition bl1_axmyv2.c:145
FLA_Obj_view
Definition FLA_type_defs.h:159

References FLA_Check_error_level(), FLA_Check_matrix_strides(), FLA_Check_submatrix_dims_and_offset(), FLA_Obj_attach_buffer(), FLA_Obj_create_without_buffer(), FLA_Obj_free_without_buffer(), FLASH_Copy_flat_to_hier(), FLASH_Obj_datatype(), and i.

◆ FLASH_Copy_flat_to_hier()

FLA_Error FLASH_Copy_flat_to_hier ( FLA_Obj F,
dim_t i,
dim_t j,
FLA_Obj H )
82{
83 FLA_Obj HTL, HTR,
84 HBL, HBR;
85 FLA_Obj HBR_tl, HBR_tr,
86 HBR_bl, HBR_br;
87 dim_t m, n;
88
89 m = FLA_Obj_length( F );
90 n = FLA_Obj_width( F );
91
92 FLASH_Part_create_2x2( H, &HTL, &HTR,
93 &HBL, &HBR, i, j, FLA_TL );
94
95 FLASH_Part_create_2x2( HBR, &HBR_tl, &HBR_tr,
96 &HBR_bl, &HBR_br, m, n, FLA_TL );
97
98 FLASH_Copy_hierarchy( FLA_FLAT_TO_HIER, F, &HBR_tl );
99
100 FLASH_Part_free_2x2( &HBR_tl, &HBR_tr,
101 &HBR_bl, &HBR_br );
102
103 FLASH_Part_free_2x2( &HTL, &HTR,
104 &HBL, &HBR );
105
106 return FLA_SUCCESS;
107}
FLASH_Copy_hierarchy
FLA_Error FLASH_Copy_hierarchy(int direction, FLA_Obj F, FLA_Obj *H)
Definition FLASH_Copy_other.c:139
FLASH_Part_create_2x2
FLA_Error FLASH_Part_create_2x2(FLA_Obj A, FLA_Obj *ATL, FLA_Obj *ATR, FLA_Obj *ABL, FLA_Obj *ABR, dim_t n_rows, dim_t n_cols, FLA_Side side)
Definition FLASH_View.c:177
FLASH_Part_free_2x2
FLA_Error FLASH_Part_free_2x2(FLA_Obj *ATL, FLA_Obj *ATR, FLA_Obj *ABL, FLA_Obj *ABR)
Definition FLASH_View.c:589
FLA_Obj_width
dim_t FLA_Obj_width(FLA_Obj obj)
Definition FLA_Query.c:123
FLA_Obj_length
dim_t FLA_Obj_length(FLA_Obj obj)
Definition FLA_Query.c:116
dim_t
unsigned long dim_t
Definition FLA_type_defs.h:71

References FLA_Obj_length(), FLA_Obj_width(), FLASH_Copy_hierarchy(), FLASH_Part_create_2x2(), FLASH_Part_free_2x2(), and i.

Referenced by FLA_LQ_UT_macro_task(), FLA_LU_piv_macro_task(), FLA_QR_UT_macro_task(), FLASH_Copy_buffer_to_hier(), FLASH_Obj_create_hier_copy_of_flat(), FLASH_Obj_create_hier_copy_of_flat_ext(), and FLASH_Obj_hierarchify().

◆ FLASH_Copy_hier_to_buffer()

FLA_Error FLASH_Copy_hier_to_buffer ( dim_t i,
dim_t j,
FLA_Obj H,
dim_t m,
dim_t n,
void * buffer,
dim_t rs,
dim_t cs )
48{
49 FLA_Obj flat_matrix;
50 FLA_Datatype datatype;
51 FLA_Error e_val;
52
53 if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
54 {
55 e_val = FLA_Check_matrix_strides( m, n, rs, cs );
56 FLA_Check_error_code( e_val );
57
58 e_val = FLA_Check_submatrix_dims_and_offset( m, n, i, j, H );
59 FLA_Check_error_code( e_val );
60 }
61
62 // Acquire the datatype from the hierarchical matrix object.
63 datatype = FLASH_Obj_datatype( H );
64
65 // Create a temporary conventional matrix object of the requested datatype
66 // and dimensions and attach the given buffer containing the incoming data.
67 FLA_Obj_create_without_buffer( datatype, m, n, &flat_matrix );
68 FLA_Obj_attach_buffer( buffer, rs, cs, &flat_matrix );
69
70 // Recurse through H, adding in the corresponding elements of flat_matrix,
71 // starting at the (i,j) element offset.
72 FLASH_Copy_hier_to_flat( i, j, H, flat_matrix );
73
74 // Free the object (but don't free the buffer!).
75 FLA_Obj_free_without_buffer( &flat_matrix );
76
77 return FLA_SUCCESS;
78}
FLASH_Copy_hier_to_flat
FLA_Error FLASH_Copy_hier_to_flat(dim_t i, dim_t j, FLA_Obj H, FLA_Obj F)
Definition FLASH_Copy_other.c:110

References FLA_Check_error_level(), FLA_Check_matrix_strides(), FLA_Check_submatrix_dims_and_offset(), FLA_Obj_attach_buffer(), FLA_Obj_create_without_buffer(), FLA_Obj_free_without_buffer(), FLASH_Copy_hier_to_flat(), FLASH_Obj_datatype(), and i.

◆ FLASH_Copy_hier_to_flat()

FLA_Error FLASH_Copy_hier_to_flat ( dim_t i,
dim_t j,
FLA_Obj H,
FLA_Obj F )
111{
112 FLA_Obj HTL, HTR,
113 HBL, HBR;
114 FLA_Obj HBR_tl, HBR_tr,
115 HBR_bl, HBR_br;
116 dim_t m, n;
117
118 m = FLA_Obj_length( F );
119 n = FLA_Obj_width( F );
120
121 FLASH_Part_create_2x2( H, &HTL, &HTR,
122 &HBL, &HBR, i, j, FLA_TL );
123
124 FLASH_Part_create_2x2( HBR, &HBR_tl, &HBR_tr,
125 &HBR_bl, &HBR_br, m, n, FLA_TL );
126
127 FLASH_Copy_hierarchy( FLA_HIER_TO_FLAT, F, &HBR_tl );
128
129 FLASH_Part_free_2x2( &HBR_tl, &HBR_tr,
130 &HBR_bl, &HBR_br );
131
132 FLASH_Part_free_2x2( &HTL, &HTR,
133 &HBL, &HBR );
134
135 return FLA_SUCCESS;
136}

References FLA_Obj_length(), FLA_Obj_width(), FLASH_Copy_hierarchy(), FLASH_Part_create_2x2(), FLASH_Part_free_2x2(), and i.

Referenced by FLASH_Copy_hier_to_buffer(), FLASH_Obj_create_flat_copy_of_hier(), and FLASH_Obj_flatten().

◆ FLASH_Copy_hierarchy()

FLA_Error FLASH_Copy_hierarchy ( int direction,
FLA_Obj F,
FLA_Obj * H )
140{
141 // Once we get down to a submatrix whose elements are scalars, we are down
142 // to our base case.
143 if ( FLA_Obj_elemtype( *H ) == FLA_SCALAR )
144 {
145 // Depending on which top-level function invoked us, we either copy
146 // the source data in the flat matrix to the leaf-level submatrix of
147 // the hierarchical matrix, or copy the data in the hierarchical
148 // submatrix to the flat matrix.
149 if ( direction == FLA_FLAT_TO_HIER )
150 {
151#ifdef FLA_ENABLE_SCC
152 if ( FLA_is_owner() )
153#endif
154 FLA_Copy_external( F, *H );
155 }
156 else if ( direction == FLA_HIER_TO_FLAT )
157 {
158#ifdef FLA_ENABLE_SCC
159 if ( FLA_is_owner() )
160#endif
161 FLA_Copy_external( *H, F );
162 }
163 }
164 else
165 {
166 FLA_Obj HL, HR, H0, H1, H2;
167 FLA_Obj FL, FR, F0, F1, F2;
168
169 FLA_Obj H1T, H01,
170 H1B, H11,
171 H21;
172 FLA_Obj F1T, F01,
173 F1B, F11,
174 F21;
175
176 dim_t b_m;
177 dim_t b_n;
178
179 FLA_Part_1x2( *H, &HL, &HR, 0, FLA_LEFT );
180 FLA_Part_1x2( F, &FL, &FR, 0, FLA_LEFT );
181
182 while ( FLA_Obj_width( HL ) < FLA_Obj_width( *H ) )
183 {
184 FLA_Repart_1x2_to_1x3( HL, /**/ HR, &H0, /**/ &H1, &H2,
185 1, FLA_RIGHT );
186
187 // Get the scalar width of H1 and use that to determine the
188 // width of F1.
189 b_n = FLASH_Obj_scalar_width( H1 );
190
191 FLA_Repart_1x2_to_1x3( FL, /**/ FR, &F0, /**/ &F1, &F2,
192 b_n, FLA_RIGHT );
193
194 // -------------------------------------------------------------
195
196 FLA_Part_2x1( H1, &H1T,
197 &H1B, 0, FLA_TOP );
198 FLA_Part_2x1( F1, &F1T,
199 &F1B, 0, FLA_TOP );
200
201 while ( FLA_Obj_length( H1T ) < FLA_Obj_length( H1 ) )
202 {
203 FLA_Repart_2x1_to_3x1( H1T, &H01,
204 /* ** */ /* *** */
205 &H11,
206 H1B, &H21, 1, FLA_BOTTOM );
207
208 // Get the scalar length of H11 and use that to determine the
209 // length of F11.
210 b_m = FLASH_Obj_scalar_length( H11 );
211
212 FLA_Repart_2x1_to_3x1( F1T, &F01,
213 /* ** */ /* *** */
214 &F11,
215 F1B, &F21, b_m, FLA_BOTTOM );
216 // -------------------------------------------------------------
217
218 // Recursively copy between F11 and H11.
219 FLASH_Copy_hierarchy( direction, F11,
220 FLASH_OBJ_PTR_AT( H11 ) );
221
222 // -------------------------------------------------------------
223
224 FLA_Cont_with_3x1_to_2x1( &H1T, H01,
225 H11,
226 /* ** */ /* *** */
227 &H1B, H21, FLA_TOP );
228 FLA_Cont_with_3x1_to_2x1( &F1T, F01,
229 F11,
230 /* ** */ /* *** */
231 &F1B, F21, FLA_TOP );
232 }
233
234 // -------------------------------------------------------------
235
236 FLA_Cont_with_1x3_to_1x2( &HL, /**/ &HR, H0, H1, /**/ H2,
237 FLA_LEFT );
238 FLA_Cont_with_1x3_to_1x2( &FL, /**/ &FR, F0, F1, /**/ F2,
239 FLA_LEFT );
240 }
241 }
242
243 return FLA_SUCCESS;
244}
FLASH_Obj_scalar_length
dim_t FLASH_Obj_scalar_length(FLA_Obj H)
Definition FLASH_View.c:600
FLASH_Obj_scalar_width
dim_t FLASH_Obj_scalar_width(FLA_Obj H)
Definition FLASH_View.c:641
FLA_Copy_external
FLA_Error FLA_Copy_external(FLA_Obj A, FLA_Obj B)
Definition FLA_Copy_external.c:13
FLA_Cont_with_3x1_to_2x1
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_Repart_2x1_to_3x1
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_Cont_with_1x3_to_1x2
FLA_Error FLA_Cont_with_1x3_to_1x2(FLA_Obj *AL, FLA_Obj *AR, FLA_Obj A0, FLA_Obj A1, FLA_Obj A2, FLA_Side side)
Definition FLA_View.c:475
FLA_Part_1x2
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
FLA_Part_2x1
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
FLA_Repart_1x2_to_1x3
FLA_Error FLA_Repart_1x2_to_1x3(FLA_Obj AL, FLA_Obj AR, FLA_Obj *A0, FLA_Obj *A1, FLA_Obj *A2, dim_t nb, FLA_Side side)
Definition FLA_View.c:267
FLA_is_owner
FLA_Bool FLA_is_owner(void)
Definition FLA_Obj.c:33
FLA_Obj_elemtype
FLA_Elemtype FLA_Obj_elemtype(FLA_Obj obj)
Definition FLA_Query.c:51

References FLA_Cont_with_1x3_to_1x2(), FLA_Cont_with_3x1_to_2x1(), FLA_Copy_external(), FLA_is_owner(), FLA_Obj_elemtype(), FLA_Obj_length(), FLA_Obj_width(), FLA_Part_1x2(), FLA_Part_2x1(), FLA_Repart_1x2_to_1x3(), FLA_Repart_2x1_to_3x1(), FLASH_Copy_hierarchy(), FLASH_Obj_scalar_length(), and FLASH_Obj_scalar_width().

Referenced by FLASH_Copy_flat_to_hier(), FLASH_Copy_hier_to_flat(), and FLASH_Copy_hierarchy().

Generated on Thu Sep 23 2021 17:46:17 for libflame by doxygen 1.12.0