36{
37 FLA_Obj ATL, ATR, A00, a01, A02,
38 ABL, ABR, a10t, alpha11, a12t,
39 A20, a21, A22;
40 FLA_Obj YTL, YTR, Y00, y01, Y02,
41 YBL, YBR, y10t, psi11, y12t,
42 Y20, y21, Y22;
43 FLA_Obj ZTL, ZTR, Z00, z01, Z02,
44 ZBL, ZBR, z10t, zeta11, z12t,
45 Z20, z21, Z22;
46 FLA_Obj TTL, TTR, T00, t01, T02,
47 TBL, TBR, t10t, tau11, t12t,
48 T20, t21, T22;
49 FLA_Obj STL, STR, S00, s01, S02,
50 SBL, SBR, s10t, sigma11, s12t,
51 S20, s21, S22;
53 uB, upsilon11,
54 u21;
56 vB, nu11,
57 v21;
59 dB, delta1,
60 d2;
62 eB, epsilon1,
63 e2;
65 fB, phi1,
66 f2;
68 gB, ghi1,
69 g2;
72
77
79 a01_b;
82 v21_b;
84
88
89
91
95
106
109
111 &ABL, &ABR, 0, 0, FLA_TL );
113 &YBL, &YBR, 0, 0, FLA_TL );
115 &ZBL, &ZBR, 0, 0, FLA_TL );
117 &TBL, &TBR, 0, 0, FLA_TL );
119 &SBL, &SBR, 0, 0, FLA_TL );
121 &uB, 0, FLA_TOP );
123 &vB, 0, FLA_TOP );
125 &dB, 0, FLA_TOP );
127 &eB, 0, FLA_TOP );
129 &fB, 0, FLA_TOP );
131 &gB, 0, FLA_TOP );
132
134 {
136
137 &a10t, &alpha11, &a12t,
138 ABL, ABR, &A20, &a21, &A22,
139 1, 1, FLA_BR );
141
142 &y10t, &psi11, &y12t,
143 YBL, YBR, &Y20, &y21, &Y22,
144 1, 1, FLA_BR );
146
147 &z10t, &zeta11, &z12t,
148 ZBL, ZBR, &Z20, &z21, &Z22,
149 1, 1, FLA_BR );
151
152 &t10t, &tau11, &t12t,
153 TBL, TBR, &T20, &t21, &T22,
154 1, 1, FLA_BR );
156
157 &s10t, &sigma11, &s12t,
158 SBL, SBR, &S20, &s21, &S22,
159 1, 1, FLA_BR );
161
162 &upsilon11,
163 uB, &u21, 1, FLA_BOTTOM );
165
166 &nu11,
167 vB, &v21, 1, FLA_BOTTOM );
169
170 &delta1,
171 dB, &d2, 1, FLA_BOTTOM );
173
174 &epsilon1,
175 eB, &e2, 1, FLA_BOTTOM );
177
178 &phi1,
179 fB, &f2, 1, FLA_BOTTOM );
181
182 &ghi1,
183 gB, &g2, 1, FLA_BOTTOM );
184
185
186
187
188
189
191 {
193 &a01_b, 1, FLA_BOTTOM );
196 }
197
199 a21, &a2 );
200
201
202
205
206
209
210
212 {
214 }
215
216
218 alpha11,
219 a21, tau11 );
221
223 {
224
225
226 FLA_Copyt( FLA_CONJ_TRANSPOSE, a12t, y21 );
228
229
232
233
234 FLA_Copyt( FLA_CONJ_TRANSPOSE, a10t, t01 );
236
239
240
242
243
245
248 &v21_b, 1, FLA_TOP );
249
250
252
253
254
256 FLA_Copyt( FLA_TRANSPOSE, a12t_r, v21_b );
257
258
259 FLA_Dotc( FLA_CONJUGATE, y21, v21, beta );
261
262
265
266
269
272
273
275
276
278 }
279 else
280 {
281
282 FLA_Copyt( FLA_CONJ_TRANSPOSE, a10t, t01 );
284 }
285
286
287
289 a10t, alpha11, a12t,
290
291 &ABL, &ABR, A20, a21, A22,
292 FLA_TL );
294 y10t, psi11, y12t,
295
296 &YBL, &YBR, Y20, y21, Y22,
297 FLA_TL );
299 z10t, zeta11, z12t,
300
301 &ZBL, &ZBR, Z20, z21, Z22,
302 FLA_TL );
304 t10t, tau11, t12t,
305
306 &TBL, &TBR, T20, t21, T22,
307 FLA_TL );
309 s10t, sigma11, s12t,
310
311 &SBL, &SBR, S20, s21, S22,
312 FLA_TL );
314 upsilon11,
315
316 &uB, u21, FLA_TOP );
318 nu11,
319
320 &vB, v21, FLA_TOP );
322 delta1,
323
324 &dB, d2, FLA_TOP );
326 epsilon1,
327
328 &eB, e2, FLA_TOP );
330 phi1,
331
332 &fB, f2, FLA_TOP );
334 ghi1,
335
336 &gB, g2, FLA_TOP );
337 }
338
349
350 return FLA_SUCCESS;
351}
FLA_Error FLA_Copy(FLA_Obj A, FLA_Obj B)
Definition FLA_Copy.c:15
FLA_Error FLA_Scal(FLA_Obj alpha, FLA_Obj A)
Definition FLA_Scal.c:15
FLA_Error FLA_Dotc(FLA_Conj conj, FLA_Obj x, FLA_Obj y, FLA_Obj rho)
Definition FLA_Dotc.c:13
FLA_Error FLA_Axpy(FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
Definition FLA_Axpy.c:15
FLA_Error FLA_Axpyt(FLA_Trans trans, FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
Definition FLA_Axpyt.c:15
FLA_Error FLA_Copyt(FLA_Trans trans, FLA_Obj A, FLA_Obj B)
Definition FLA_Copyt.c:15
FLA_Error FLA_Inv_scalc(FLA_Conj conjalpha, FLA_Obj alpha, FLA_Obj A)
Definition FLA_Inv_scalc.c:13
FLA_Error FLA_Gemvc(FLA_Trans transa, FLA_Conj conjx, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
Definition FLA_Gemvc.c:13
FLA_Error FLA_Gemv(FLA_Trans transa, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
Definition FLA_Gemv.c:15
FLA_Obj FLA_MINUS_ONE
Definition FLA_Init.c:22
FLA_Obj FLA_ZERO
Definition FLA_Init.c:20
FLA_Obj FLA_ONE
Definition FLA_Init.c:18
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
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
dim_t FLA_Obj_width(FLA_Obj obj)
Definition FLA_Query.c:123
FLA_Error FLA_Obj_create(FLA_Datatype datatype, dim_t m, dim_t n, dim_t rs, dim_t cs, FLA_Obj *obj)
Definition FLA_Obj.c:55
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_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
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
FLA_Error FLA_Merge_2x1(FLA_Obj AT, FLA_Obj AB, FLA_Obj *A)
Definition FLA_View.c:541
FLA_Error FLA_Obj_free(FLA_Obj *obj)
Definition FLA_Obj.c:588
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition FLA_Query.c:13
int FLA_Datatype
Definition FLA_type_defs.h:49
unsigned long dim_t
Definition FLA_type_defs.h:71
FLA_Error FLA_Set(FLA_Obj alpha, FLA_Obj A)
Definition FLA_Set.c:13
FLA_Error FLA_Househ2_UT(FLA_Side side, FLA_Obj chi_1, FLA_Obj x2, FLA_Obj tau)
Definition FLA_Househ2_UT.c:16
Definition FLA_type_defs.h:159