Lines Matching refs:LDA
12 double *a, int *lda, double *x, int *incx, double *beta,
16 int i,j,LDA;
21 LDA = *n+1;
22 A = ( double* )malloc( (*m)*LDA*sizeof( double ) );
25 A[ LDA*i+j ]=a[ (*lda)*j+i ];
27 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy );
32 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
35 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
39 double *y, int *incy, double *a, int *lda ) {
42 int i,j,LDA;
45 LDA = *n+1;
46 A = ( double* )malloc( (*m)*LDA*sizeof( double ) );
50 A[ LDA*i+j ]=a[ (*lda)*j+i ];
53 cblas_dger(CblasRowMajor, *m, *n, *alpha, x, *incx, y, *incy, A, LDA );
56 a[ (*lda)*j+i ]=A[ LDA*i+j ];
60 cblas_dger( CblasColMajor, *m, *n, *alpha, x, *incx, y, *incy, a, *lda );
64 int *n, double *a, int *lda, double *x, int *incx) {
66 int i,j,LDA;
76 LDA = *n+1;
77 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
80 A[ LDA*i+j ]=a[ (*lda)*j+i ];
81 cblas_dtrmv(CblasRowMajor, uplo, trans, diag, *n, A, LDA, x, *incx);
85 cblas_dtrmv(CblasColMajor, uplo, trans, diag, *n, a, *lda, x, *incx);
87 cblas_dtrmv(UNDEFINED, uplo, trans, diag, *n, a, *lda, x, *incx);
92 int *n, double *a, int *lda, double *x, int *incx ) {
94 int i,j,LDA;
104 LDA = *n+1;
105 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
108 A[ LDA*i+j ]=a[ (*lda)*j+i ];
109 cblas_dtrsv(CblasRowMajor, uplo, trans, diag, *n, A, LDA, x, *incx );
113 cblas_dtrsv(CblasColMajor, uplo, trans, diag, *n, a, *lda, x, *incx );
116 int *lda, double *x, int *incx, double *beta, double *y,
119 int i,j,LDA;
125 LDA = *n+1;
126 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
129 A[ LDA*i+j ]=a[ (*lda)*j+i ];
130 cblas_dsymv(CblasRowMajor, uplo, *n, *alpha, A, LDA, x, *incx,
135 cblas_dsymv(CblasColMajor, uplo, *n, *alpha, a, *lda, x, *incx,
140 int *incx, double *a, int *lda) {
142 int i,j,LDA;
148 LDA = *n+1;
149 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
152 A[ LDA*i+j ]=a[ (*lda)*j+i ];
153 cblas_dsyr(CblasRowMajor, uplo, *n, *alpha, x, *incx, A, LDA);
156 a[ (*lda)*j+i ]=A[ LDA*i+j ];
160 cblas_dsyr(CblasColMajor, uplo, *n, *alpha, x, *incx, a, *lda);
164 int *incx, double *y, int *incy, double *a, int *lda) {
166 int i,j,LDA;
172 LDA = *n+1;
173 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
176 A[ LDA*i+j ]=a[ (*lda)*j+i ];
177 cblas_dsyr2(CblasRowMajor, uplo, *n, *alpha, x, *incx, y, *incy, A, LDA);
180 a[ (*lda)*j+i ]=A[ LDA*i+j ];
184 cblas_dsyr2(CblasColMajor, uplo, *n, *alpha, x, *incx, y, *incy, a, *lda);
188 double *alpha, double *a, int *lda, double *x, int *incx,
192 int i,irow,j,jcol,LDA;
198 LDA = *ku+*kl+2;
199 A = ( double* )malloc( (*n+*kl)*LDA*sizeof( double ) );
204 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
209 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
214 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
217 A, LDA, x, *incx, *beta, y, *incy );
222 a, *lda, x, *incx, *beta, y, *incy );
226 int *n, int *k, double *a, int *lda, double *x, int *incx) {
228 int irow, jcol, i, j, LDA;
238 LDA = *k+1;
239 A = ( double* )malloc( (*n+*k)*LDA*sizeof( double ) );
245 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
250 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
256 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
261 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
264 cblas_dtbmv(CblasRowMajor, uplo, trans, diag, *n, *k, A, LDA, x, *incx);
268 cblas_dtbmv(CblasColMajor, uplo, trans, diag, *n, *k, a, *lda, x, *incx);
272 int *n, int *k, double *a, int *lda, double *x, int *incx) {
274 int irow, jcol, i, j, LDA;
284 LDA = *k+1;
285 A = ( double* )malloc( (*n+*k)*LDA*sizeof( double ) );
291 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
296 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
302 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
307 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
310 cblas_dtbsv(CblasRowMajor, uplo, trans, diag, *n, *k, A, LDA, x, *incx);
314 cblas_dtbsv(CblasColMajor, uplo, trans, diag, *n, *k, a, *lda, x, *incx);
318 double *a, int *lda, double *x, int *incx, double *beta,
321 int i,j,irow,jcol,LDA;
327 LDA = *k+1;
328 A = ( double* )malloc( (*n+*k)*LDA*sizeof( double ) );
334 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
339 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
345 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
350 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
353 cblas_dsbmv(CblasRowMajor, uplo, *n, *k, *alpha, A, LDA, x, *incx,
358 cblas_dsbmv(CblasColMajor, uplo, *n, *k, *alpha, a, *lda, x, *incx,
365 int i,j,k,LDA;
371 LDA = *n;
372 A = ( double* )malloc( LDA*LDA*sizeof( double ) );
373 AP = ( double* )malloc( (((LDA+1)*LDA)/2)*sizeof( double ) );
377 A[ LDA*i+j ]=ap[ k ];
380 AP[ k ]=A[ LDA*i+j ];
385 A[ LDA*i+j ]=ap[ k ];
388 AP[ k ]=A[ LDA*i+j ];
403 int i, j, k, LDA;
413 LDA = *n;
414 A = ( double* )malloc( LDA*LDA*sizeof( double ) );
415 AP = ( double* )malloc( (((LDA+1)*LDA)/2)*sizeof( double ) );
419 A[ LDA*i+j ]=ap[ k ];
422 AP[ k ]=A[ LDA*i+j ];
427 A[ LDA*i+j ]=ap[ k ];
430 AP[ k ]=A[ LDA*i+j ];
443 int i, j, k, LDA;
453 LDA = *n;
454 A = ( double* )malloc( LDA*LDA*sizeof( double ) );
455 AP = ( double* )malloc( (((LDA+1)*LDA)/2)*sizeof( double ) );
459 A[ LDA*i+j ]=ap[ k ];
462 AP[ k ]=A[ LDA*i+j ];
468 A[ LDA*i+j ]=ap[ k ];
471 AP[ k ]=A[ LDA*i+j ];
484 int i,j,k,LDA;
490 LDA = *n;
491 A = ( double* )malloc( LDA*LDA*sizeof( double ) );
492 AP = ( double* )malloc( (((LDA+1)*LDA)/2)*sizeof( double ) );
496 A[ LDA*i+j ]=ap[ k ];
499 AP[ k ]=A[ LDA*i+j ];
504 A[ LDA*i+j ]=ap[ k ];
507 AP[ k ]=A[ LDA*i+j ];
513 A[ LDA*i+j ]=AP[ k ];
516 ap[ k ]=A[ LDA*i+j ];
521 A[ LDA*i+j ]=AP[ k ];
524 ap[ k ]=A[ LDA*i+j ];
536 int i,j,k,LDA;
542 LDA = *n;
543 A = ( double* )malloc( LDA*LDA*sizeof( double ) );
544 AP = ( double* )malloc( (((LDA+1)*LDA)/2)*sizeof( double ) );
548 A[ LDA*i+j ]=ap[ k ];
551 AP[ k ]=A[ LDA*i+j ];
556 A[ LDA*i+j ]=ap[ k ];
559 AP[ k ]=A[ LDA*i+j ];
565 A[ LDA*i+j ]=AP[ k ];
568 ap[ k ]=A[ LDA*i+j ];
573 A[ LDA*i+j ]=AP[ k ];
576 ap[ k ]=A[ LDA*i+j ];