1 /*
2 * Written by D.P. Manley, Digital Equipment Corporation.
3 * Prefixed "C_" to BLAS routines and their declarations.
4 *
5 * Modified by T. H. Do, 4/15/98, SGI/CRAY Research.
6 */
7 #include <stdlib.h>
8 #include "cblas.h"
9 #include "cblas_test.h"
10 #define TEST_COL_MJR 0
11 #define TEST_ROW_MJR 1
12 #define UNDEFINED -1
13
F77_zgemm(int * order,char * transpa,char * transpb,int * m,int * n,int * k,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb,CBLAS_TEST_ZOMPLEX * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)14 void F77_zgemm(int *order, char *transpa, char *transpb, int *m, int *n,
15 int *k, CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
16 CBLAS_TEST_ZOMPLEX *b, int *ldb, CBLAS_TEST_ZOMPLEX *beta,
17 CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
18
19 CBLAS_TEST_ZOMPLEX *A, *B, *C;
20 int i,j,LDA, LDB, LDC;
21 enum CBLAS_TRANSPOSE transa, transb;
22
23 get_transpose_type(transpa, &transa);
24 get_transpose_type(transpb, &transb);
25
26 if (*order == TEST_ROW_MJR) {
27 if (transa == CblasNoTrans) {
28 LDA = *k+1;
29 A=(CBLAS_TEST_ZOMPLEX*)malloc((*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
30 for( i=0; i<*m; i++ )
31 for( j=0; j<*k; j++ ) {
32 A[i*LDA+j].real=a[j*(*lda)+i].real;
33 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
34 }
35 }
36 else {
37 LDA = *m+1;
38 A=(CBLAS_TEST_ZOMPLEX* )malloc(LDA*(*k)*sizeof(CBLAS_TEST_ZOMPLEX));
39 for( i=0; i<*k; i++ )
40 for( j=0; j<*m; j++ ) {
41 A[i*LDA+j].real=a[j*(*lda)+i].real;
42 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
43 }
44 }
45
46 if (transb == CblasNoTrans) {
47 LDB = *n+1;
48 B=(CBLAS_TEST_ZOMPLEX* )malloc((*k)*LDB*sizeof(CBLAS_TEST_ZOMPLEX) );
49 for( i=0; i<*k; i++ )
50 for( j=0; j<*n; j++ ) {
51 B[i*LDB+j].real=b[j*(*ldb)+i].real;
52 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
53 }
54 }
55 else {
56 LDB = *k+1;
57 B=(CBLAS_TEST_ZOMPLEX* )malloc(LDB*(*n)*sizeof(CBLAS_TEST_ZOMPLEX));
58 for( i=0; i<*n; i++ )
59 for( j=0; j<*k; j++ ) {
60 B[i*LDB+j].real=b[j*(*ldb)+i].real;
61 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
62 }
63 }
64
65 LDC = *n+1;
66 C=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDC*sizeof(CBLAS_TEST_ZOMPLEX));
67 for( j=0; j<*n; j++ )
68 for( i=0; i<*m; i++ ) {
69 C[i*LDC+j].real=c[j*(*ldc)+i].real;
70 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
71 }
72 cblas_zgemm( CblasRowMajor, transa, transb, *m, *n, *k, alpha, A, LDA,
73 B, LDB, beta, C, LDC );
74 for( j=0; j<*n; j++ )
75 for( i=0; i<*m; i++ ) {
76 c[j*(*ldc)+i].real=C[i*LDC+j].real;
77 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
78 }
79 free(A);
80 free(B);
81 free(C);
82 }
83 else if (*order == TEST_COL_MJR)
84 cblas_zgemm( CblasColMajor, transa, transb, *m, *n, *k, alpha, a, *lda,
85 b, *ldb, beta, c, *ldc );
86 else
87 cblas_zgemm( UNDEFINED, transa, transb, *m, *n, *k, alpha, a, *lda,
88 b, *ldb, beta, c, *ldc );
89 }
F77_zhemm(int * order,char * rtlf,char * uplow,int * m,int * n,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb,CBLAS_TEST_ZOMPLEX * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)90 void F77_zhemm(int *order, char *rtlf, char *uplow, int *m, int *n,
91 CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
92 CBLAS_TEST_ZOMPLEX *b, int *ldb, CBLAS_TEST_ZOMPLEX *beta,
93 CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
94
95 CBLAS_TEST_ZOMPLEX *A, *B, *C;
96 int i,j,LDA, LDB, LDC;
97 enum CBLAS_UPLO uplo;
98 enum CBLAS_SIDE side;
99
100 get_uplo_type(uplow,&uplo);
101 get_side_type(rtlf,&side);
102
103 if (*order == TEST_ROW_MJR) {
104 if (side == CblasLeft) {
105 LDA = *m+1;
106 A= (CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
107 for( i=0; i<*m; i++ )
108 for( j=0; j<*m; j++ ) {
109 A[i*LDA+j].real=a[j*(*lda)+i].real;
110 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
111 }
112 }
113 else{
114 LDA = *n+1;
115 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
116 for( i=0; i<*n; i++ )
117 for( j=0; j<*n; j++ ) {
118 A[i*LDA+j].real=a[j*(*lda)+i].real;
119 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
120 }
121 }
122 LDB = *n+1;
123 B=(CBLAS_TEST_ZOMPLEX* )malloc( (*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX ) );
124 for( i=0; i<*m; i++ )
125 for( j=0; j<*n; j++ ) {
126 B[i*LDB+j].real=b[j*(*ldb)+i].real;
127 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
128 }
129 LDC = *n+1;
130 C=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDC*sizeof(CBLAS_TEST_ZOMPLEX ) );
131 for( j=0; j<*n; j++ )
132 for( i=0; i<*m; i++ ) {
133 C[i*LDC+j].real=c[j*(*ldc)+i].real;
134 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
135 }
136 cblas_zhemm( CblasRowMajor, side, uplo, *m, *n, alpha, A, LDA, B, LDB,
137 beta, C, LDC );
138 for( j=0; j<*n; j++ )
139 for( i=0; i<*m; i++ ) {
140 c[j*(*ldc)+i].real=C[i*LDC+j].real;
141 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
142 }
143 free(A);
144 free(B);
145 free(C);
146 }
147 else if (*order == TEST_COL_MJR)
148 cblas_zhemm( CblasColMajor, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
149 beta, c, *ldc );
150 else
151 cblas_zhemm( UNDEFINED, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
152 beta, c, *ldc );
153 }
F77_zsymm(int * order,char * rtlf,char * uplow,int * m,int * n,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb,CBLAS_TEST_ZOMPLEX * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)154 void F77_zsymm(int *order, char *rtlf, char *uplow, int *m, int *n,
155 CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
156 CBLAS_TEST_ZOMPLEX *b, int *ldb, CBLAS_TEST_ZOMPLEX *beta,
157 CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
158
159 CBLAS_TEST_ZOMPLEX *A, *B, *C;
160 int i,j,LDA, LDB, LDC;
161 enum CBLAS_UPLO uplo;
162 enum CBLAS_SIDE side;
163
164 get_uplo_type(uplow,&uplo);
165 get_side_type(rtlf,&side);
166
167 if (*order == TEST_ROW_MJR) {
168 if (side == CblasLeft) {
169 LDA = *m+1;
170 A=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
171 for( i=0; i<*m; i++ )
172 for( j=0; j<*m; j++ )
173 A[i*LDA+j]=a[j*(*lda)+i];
174 }
175 else{
176 LDA = *n+1;
177 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
178 for( i=0; i<*n; i++ )
179 for( j=0; j<*n; j++ )
180 A[i*LDA+j]=a[j*(*lda)+i];
181 }
182 LDB = *n+1;
183 B=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX ));
184 for( i=0; i<*m; i++ )
185 for( j=0; j<*n; j++ )
186 B[i*LDB+j]=b[j*(*ldb)+i];
187 LDC = *n+1;
188 C=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDC*sizeof(CBLAS_TEST_ZOMPLEX));
189 for( j=0; j<*n; j++ )
190 for( i=0; i<*m; i++ )
191 C[i*LDC+j]=c[j*(*ldc)+i];
192 cblas_zsymm( CblasRowMajor, side, uplo, *m, *n, alpha, A, LDA, B, LDB,
193 beta, C, LDC );
194 for( j=0; j<*n; j++ )
195 for( i=0; i<*m; i++ )
196 c[j*(*ldc)+i]=C[i*LDC+j];
197 free(A);
198 free(B);
199 free(C);
200 }
201 else if (*order == TEST_COL_MJR)
202 cblas_zsymm( CblasColMajor, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
203 beta, c, *ldc );
204 else
205 cblas_zsymm( UNDEFINED, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
206 beta, c, *ldc );
207 }
208
F77_zherk(int * order,char * uplow,char * transp,int * n,int * k,double * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,double * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)209 void F77_zherk(int *order, char *uplow, char *transp, int *n, int *k,
210 double *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
211 double *beta, CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
212
213 int i,j,LDA,LDC;
214 CBLAS_TEST_ZOMPLEX *A, *C;
215 enum CBLAS_UPLO uplo;
216 enum CBLAS_TRANSPOSE trans;
217
218 get_uplo_type(uplow,&uplo);
219 get_transpose_type(transp,&trans);
220
221 if (*order == TEST_ROW_MJR) {
222 if (trans == CblasNoTrans) {
223 LDA = *k+1;
224 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
225 for( i=0; i<*n; i++ )
226 for( j=0; j<*k; j++ ) {
227 A[i*LDA+j].real=a[j*(*lda)+i].real;
228 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
229 }
230 }
231 else{
232 LDA = *n+1;
233 A=(CBLAS_TEST_ZOMPLEX* )malloc((*k)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
234 for( i=0; i<*k; i++ )
235 for( j=0; j<*n; j++ ) {
236 A[i*LDA+j].real=a[j*(*lda)+i].real;
237 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
238 }
239 }
240 LDC = *n+1;
241 C=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDC*sizeof(CBLAS_TEST_ZOMPLEX ) );
242 for( i=0; i<*n; i++ )
243 for( j=0; j<*n; j++ ) {
244 C[i*LDC+j].real=c[j*(*ldc)+i].real;
245 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
246 }
247 cblas_zherk(CblasRowMajor, uplo, trans, *n, *k, *alpha, A, LDA, *beta,
248 C, LDC );
249 for( j=0; j<*n; j++ )
250 for( i=0; i<*n; i++ ) {
251 c[j*(*ldc)+i].real=C[i*LDC+j].real;
252 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
253 }
254 free(A);
255 free(C);
256 }
257 else if (*order == TEST_COL_MJR)
258 cblas_zherk(CblasColMajor, uplo, trans, *n, *k, *alpha, a, *lda, *beta,
259 c, *ldc );
260 else
261 cblas_zherk(UNDEFINED, uplo, trans, *n, *k, *alpha, a, *lda, *beta,
262 c, *ldc );
263 }
264
F77_zsyrk(int * order,char * uplow,char * transp,int * n,int * k,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)265 void F77_zsyrk(int *order, char *uplow, char *transp, int *n, int *k,
266 CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
267 CBLAS_TEST_ZOMPLEX *beta, CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
268
269 int i,j,LDA,LDC;
270 CBLAS_TEST_ZOMPLEX *A, *C;
271 enum CBLAS_UPLO uplo;
272 enum CBLAS_TRANSPOSE trans;
273
274 get_uplo_type(uplow,&uplo);
275 get_transpose_type(transp,&trans);
276
277 if (*order == TEST_ROW_MJR) {
278 if (trans == CblasNoTrans) {
279 LDA = *k+1;
280 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
281 for( i=0; i<*n; i++ )
282 for( j=0; j<*k; j++ ) {
283 A[i*LDA+j].real=a[j*(*lda)+i].real;
284 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
285 }
286 }
287 else{
288 LDA = *n+1;
289 A=(CBLAS_TEST_ZOMPLEX* )malloc((*k)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
290 for( i=0; i<*k; i++ )
291 for( j=0; j<*n; j++ ) {
292 A[i*LDA+j].real=a[j*(*lda)+i].real;
293 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
294 }
295 }
296 LDC = *n+1;
297 C=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDC*sizeof(CBLAS_TEST_ZOMPLEX ) );
298 for( i=0; i<*n; i++ )
299 for( j=0; j<*n; j++ ) {
300 C[i*LDC+j].real=c[j*(*ldc)+i].real;
301 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
302 }
303 cblas_zsyrk(CblasRowMajor, uplo, trans, *n, *k, alpha, A, LDA, beta,
304 C, LDC );
305 for( j=0; j<*n; j++ )
306 for( i=0; i<*n; i++ ) {
307 c[j*(*ldc)+i].real=C[i*LDC+j].real;
308 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
309 }
310 free(A);
311 free(C);
312 }
313 else if (*order == TEST_COL_MJR)
314 cblas_zsyrk(CblasColMajor, uplo, trans, *n, *k, alpha, a, *lda, beta,
315 c, *ldc );
316 else
317 cblas_zsyrk(UNDEFINED, uplo, trans, *n, *k, alpha, a, *lda, beta,
318 c, *ldc );
319 }
F77_zher2k(int * order,char * uplow,char * transp,int * n,int * k,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb,double * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)320 void F77_zher2k(int *order, char *uplow, char *transp, int *n, int *k,
321 CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
322 CBLAS_TEST_ZOMPLEX *b, int *ldb, double *beta,
323 CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
324 int i,j,LDA,LDB,LDC;
325 CBLAS_TEST_ZOMPLEX *A, *B, *C;
326 enum CBLAS_UPLO uplo;
327 enum CBLAS_TRANSPOSE trans;
328
329 get_uplo_type(uplow,&uplo);
330 get_transpose_type(transp,&trans);
331
332 if (*order == TEST_ROW_MJR) {
333 if (trans == CblasNoTrans) {
334 LDA = *k+1;
335 LDB = *k+1;
336 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ));
337 B=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDB*sizeof(CBLAS_TEST_ZOMPLEX ));
338 for( i=0; i<*n; i++ )
339 for( j=0; j<*k; j++ ) {
340 A[i*LDA+j].real=a[j*(*lda)+i].real;
341 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
342 B[i*LDB+j].real=b[j*(*ldb)+i].real;
343 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
344 }
345 }
346 else {
347 LDA = *n+1;
348 LDB = *n+1;
349 A=(CBLAS_TEST_ZOMPLEX* )malloc( LDA*(*k)*sizeof(CBLAS_TEST_ZOMPLEX ) );
350 B=(CBLAS_TEST_ZOMPLEX* )malloc( LDB*(*k)*sizeof(CBLAS_TEST_ZOMPLEX ) );
351 for( i=0; i<*k; i++ )
352 for( j=0; j<*n; j++ ){
353 A[i*LDA+j].real=a[j*(*lda)+i].real;
354 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
355 B[i*LDB+j].real=b[j*(*ldb)+i].real;
356 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
357 }
358 }
359 LDC = *n+1;
360 C=(CBLAS_TEST_ZOMPLEX* )malloc( (*n)*LDC*sizeof(CBLAS_TEST_ZOMPLEX ) );
361 for( i=0; i<*n; i++ )
362 for( j=0; j<*n; j++ ) {
363 C[i*LDC+j].real=c[j*(*ldc)+i].real;
364 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
365 }
366 cblas_zher2k(CblasRowMajor, uplo, trans, *n, *k, alpha, A, LDA,
367 B, LDB, *beta, C, LDC );
368 for( j=0; j<*n; j++ )
369 for( i=0; i<*n; i++ ) {
370 c[j*(*ldc)+i].real=C[i*LDC+j].real;
371 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
372 }
373 free(A);
374 free(B);
375 free(C);
376 }
377 else if (*order == TEST_COL_MJR)
378 cblas_zher2k(CblasColMajor, uplo, trans, *n, *k, alpha, a, *lda,
379 b, *ldb, *beta, c, *ldc );
380 else
381 cblas_zher2k(UNDEFINED, uplo, trans, *n, *k, alpha, a, *lda,
382 b, *ldb, *beta, c, *ldc );
383 }
F77_zsyr2k(int * order,char * uplow,char * transp,int * n,int * k,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb,CBLAS_TEST_ZOMPLEX * beta,CBLAS_TEST_ZOMPLEX * c,int * ldc)384 void F77_zsyr2k(int *order, char *uplow, char *transp, int *n, int *k,
385 CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a, int *lda,
386 CBLAS_TEST_ZOMPLEX *b, int *ldb, CBLAS_TEST_ZOMPLEX *beta,
387 CBLAS_TEST_ZOMPLEX *c, int *ldc ) {
388 int i,j,LDA,LDB,LDC;
389 CBLAS_TEST_ZOMPLEX *A, *B, *C;
390 enum CBLAS_UPLO uplo;
391 enum CBLAS_TRANSPOSE trans;
392
393 get_uplo_type(uplow,&uplo);
394 get_transpose_type(transp,&trans);
395
396 if (*order == TEST_ROW_MJR) {
397 if (trans == CblasNoTrans) {
398 LDA = *k+1;
399 LDB = *k+1;
400 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
401 B=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDB*sizeof(CBLAS_TEST_ZOMPLEX));
402 for( i=0; i<*n; i++ )
403 for( j=0; j<*k; j++ ) {
404 A[i*LDA+j].real=a[j*(*lda)+i].real;
405 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
406 B[i*LDB+j].real=b[j*(*ldb)+i].real;
407 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
408 }
409 }
410 else {
411 LDA = *n+1;
412 LDB = *n+1;
413 A=(CBLAS_TEST_ZOMPLEX* )malloc(LDA*(*k)*sizeof(CBLAS_TEST_ZOMPLEX));
414 B=(CBLAS_TEST_ZOMPLEX* )malloc(LDB*(*k)*sizeof(CBLAS_TEST_ZOMPLEX));
415 for( i=0; i<*k; i++ )
416 for( j=0; j<*n; j++ ){
417 A[i*LDA+j].real=a[j*(*lda)+i].real;
418 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
419 B[i*LDB+j].real=b[j*(*ldb)+i].real;
420 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
421 }
422 }
423 LDC = *n+1;
424 C=(CBLAS_TEST_ZOMPLEX* )malloc( (*n)*LDC*sizeof(CBLAS_TEST_ZOMPLEX));
425 for( i=0; i<*n; i++ )
426 for( j=0; j<*n; j++ ) {
427 C[i*LDC+j].real=c[j*(*ldc)+i].real;
428 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
429 }
430 cblas_zsyr2k(CblasRowMajor, uplo, trans, *n, *k, alpha, A, LDA,
431 B, LDB, beta, C, LDC );
432 for( j=0; j<*n; j++ )
433 for( i=0; i<*n; i++ ) {
434 c[j*(*ldc)+i].real=C[i*LDC+j].real;
435 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
436 }
437 free(A);
438 free(B);
439 free(C);
440 }
441 else if (*order == TEST_COL_MJR)
442 cblas_zsyr2k(CblasColMajor, uplo, trans, *n, *k, alpha, a, *lda,
443 b, *ldb, beta, c, *ldc );
444 else
445 cblas_zsyr2k(UNDEFINED, uplo, trans, *n, *k, alpha, a, *lda,
446 b, *ldb, beta, c, *ldc );
447 }
F77_ztrmm(int * order,char * rtlf,char * uplow,char * transp,char * diagn,int * m,int * n,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb)448 void F77_ztrmm(int *order, char *rtlf, char *uplow, char *transp, char *diagn,
449 int *m, int *n, CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a,
450 int *lda, CBLAS_TEST_ZOMPLEX *b, int *ldb) {
451 int i,j,LDA,LDB;
452 CBLAS_TEST_ZOMPLEX *A, *B;
453 enum CBLAS_SIDE side;
454 enum CBLAS_DIAG diag;
455 enum CBLAS_UPLO uplo;
456 enum CBLAS_TRANSPOSE trans;
457
458 get_uplo_type(uplow,&uplo);
459 get_transpose_type(transp,&trans);
460 get_diag_type(diagn,&diag);
461 get_side_type(rtlf,&side);
462
463 if (*order == TEST_ROW_MJR) {
464 if (side == CblasLeft) {
465 LDA = *m+1;
466 A=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
467 for( i=0; i<*m; i++ )
468 for( j=0; j<*m; j++ ) {
469 A[i*LDA+j].real=a[j*(*lda)+i].real;
470 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
471 }
472 }
473 else{
474 LDA = *n+1;
475 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
476 for( i=0; i<*n; i++ )
477 for( j=0; j<*n; j++ ) {
478 A[i*LDA+j].real=a[j*(*lda)+i].real;
479 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
480 }
481 }
482 LDB = *n+1;
483 B=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX));
484 for( i=0; i<*m; i++ )
485 for( j=0; j<*n; j++ ) {
486 B[i*LDB+j].real=b[j*(*ldb)+i].real;
487 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
488 }
489 cblas_ztrmm(CblasRowMajor, side, uplo, trans, diag, *m, *n, alpha,
490 A, LDA, B, LDB );
491 for( j=0; j<*n; j++ )
492 for( i=0; i<*m; i++ ) {
493 b[j*(*ldb)+i].real=B[i*LDB+j].real;
494 b[j*(*ldb)+i].imag=B[i*LDB+j].imag;
495 }
496 free(A);
497 free(B);
498 }
499 else if (*order == TEST_COL_MJR)
500 cblas_ztrmm(CblasColMajor, side, uplo, trans, diag, *m, *n, alpha,
501 a, *lda, b, *ldb);
502 else
503 cblas_ztrmm(UNDEFINED, side, uplo, trans, diag, *m, *n, alpha,
504 a, *lda, b, *ldb);
505 }
506
F77_ztrsm(int * order,char * rtlf,char * uplow,char * transp,char * diagn,int * m,int * n,CBLAS_TEST_ZOMPLEX * alpha,CBLAS_TEST_ZOMPLEX * a,int * lda,CBLAS_TEST_ZOMPLEX * b,int * ldb)507 void F77_ztrsm(int *order, char *rtlf, char *uplow, char *transp, char *diagn,
508 int *m, int *n, CBLAS_TEST_ZOMPLEX *alpha, CBLAS_TEST_ZOMPLEX *a,
509 int *lda, CBLAS_TEST_ZOMPLEX *b, int *ldb) {
510 int i,j,LDA,LDB;
511 CBLAS_TEST_ZOMPLEX *A, *B;
512 enum CBLAS_SIDE side;
513 enum CBLAS_DIAG diag;
514 enum CBLAS_UPLO uplo;
515 enum CBLAS_TRANSPOSE trans;
516
517 get_uplo_type(uplow,&uplo);
518 get_transpose_type(transp,&trans);
519 get_diag_type(diagn,&diag);
520 get_side_type(rtlf,&side);
521
522 if (*order == TEST_ROW_MJR) {
523 if (side == CblasLeft) {
524 LDA = *m+1;
525 A=(CBLAS_TEST_ZOMPLEX* )malloc( (*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
526 for( i=0; i<*m; i++ )
527 for( j=0; j<*m; j++ ) {
528 A[i*LDA+j].real=a[j*(*lda)+i].real;
529 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
530 }
531 }
532 else{
533 LDA = *n+1;
534 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
535 for( i=0; i<*n; i++ )
536 for( j=0; j<*n; j++ ) {
537 A[i*LDA+j].real=a[j*(*lda)+i].real;
538 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
539 }
540 }
541 LDB = *n+1;
542 B=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX));
543 for( i=0; i<*m; i++ )
544 for( j=0; j<*n; j++ ) {
545 B[i*LDB+j].real=b[j*(*ldb)+i].real;
546 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
547 }
548 cblas_ztrsm(CblasRowMajor, side, uplo, trans, diag, *m, *n, alpha,
549 A, LDA, B, LDB );
550 for( j=0; j<*n; j++ )
551 for( i=0; i<*m; i++ ) {
552 b[j*(*ldb)+i].real=B[i*LDB+j].real;
553 b[j*(*ldb)+i].imag=B[i*LDB+j].imag;
554 }
555 free(A);
556 free(B);
557 }
558 else if (*order == TEST_COL_MJR)
559 cblas_ztrsm(CblasColMajor, side, uplo, trans, diag, *m, *n, alpha,
560 a, *lda, b, *ldb);
561 else
562 cblas_ztrsm(UNDEFINED, side, uplo, trans, diag, *m, *n, alpha,
563 a, *lda, b, *ldb);
564 }
565