26 A = (
double* )malloc( (*m)*LDA*
sizeof( double ) );
29 A[ LDA*i+j ]=a[ (*lda)*j+i ];
31 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy );
36 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
39 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
50 A = (
double* )malloc( (*m)*LDA*
sizeof( double ) );
52 for( i=0; i<*m; i++ ) {
54 A[ LDA*i+j ]=a[ (*lda)*j+i ];
60 a[ (*lda)*j+i ]=A[ LDA*i+j ];
85 A = (
double* )malloc( (*n)*LDA*
sizeof( double ) );
88 A[ LDA*i+j ]=a[ (*lda)*j+i ];
117 A = (
double* )malloc( (*n)*LDA*
sizeof( double ) );
118 for( i=0; i<*n; i++ )
119 for( j=0; j<*n; j++ )
120 A[ LDA*i+j ]=a[ (*lda)*j+i ];
142 A = (
double* )malloc( (*n)*LDA*
sizeof( double ) );
143 for( i=0; i<*n; i++ )
144 for( j=0; j<*n; j++ )
145 A[ LDA*i+j ]=a[ (*lda)*j+i ];
169 A = (
double* )malloc( (*n)*LDA*
sizeof( double ) );
170 for( i=0; i<*n; i++ )
171 for( j=0; j<*n; j++ )
172 A[ LDA*i+j ]=a[ (*lda)*j+i ];
174 for( i=0; i<*n; i++ )
175 for( j=0; j<*n; j++ )
176 a[ (*lda)*j+i ]=A[ LDA*i+j ];
197 A = (
double* )malloc( (*n)*LDA*
sizeof( double ) );
198 for( i=0; i<*n; i++ )
199 for( j=0; j<*n; j++ )
200 A[ LDA*i+j ]=a[ (*lda)*j+i ];
202 for( i=0; i<*n; i++ )
203 for( j=0; j<*n; j++ )
204 a[ (*lda)*j+i ]=A[ LDA*i+j ];
227 A = (
double* )malloc( (*n+*kl)*LDA*
sizeof( double ) );
228 for( i=0; i<*ku; i++ ){
231 for( j=jcol; j<*n; j++ )
232 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
236 for( j=0; j<*n; j++ )
237 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
238 for( i=*ku+1; i<*ku+*kl+1; i++ ){
241 for( j=jcol; j<(*n+*kl); j++ )
242 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
245 A, LDA, x, *incx, *beta, y, *incy );
250 a, *lda, x, *incx, *beta, y, *incy );
271 A = (
double* )malloc( (*n+*k)*LDA*
sizeof( double ) );
273 for( i=0; i<*k; i++ ){
276 for( j=jcol; j<*n; j++ )
277 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
281 for( j=0; j<*n; j++ )
282 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
287 for( j=0; j<*n; j++ )
288 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
289 for( i=1; i<*k+1; i++ ){
292 for( j=jcol; j<(*n+*k); j++ )
293 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
321 A = (
double* )malloc( (*n+*k)*LDA*
sizeof( double ) );
323 for( i=0; i<*k; i++ ){
326 for( j=jcol; j<*n; j++ )
327 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
331 for( j=0; j<*n; j++ )
332 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
337 for( j=0; j<*n; j++ )
338 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
339 for( i=1; i<*k+1; i++ ){
342 for( j=jcol; j<(*n+*k); j++ )
343 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
368 A = (
double* )malloc( (*n+*k)*LDA*
sizeof( double ) );
370 for( i=0; i<*k; i++ ){
373 for( j=jcol; j<*n; j++ )
374 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
378 for( j=0; j<*n; j++ )
379 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
384 for( j=0; j<*n; j++ )
385 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
386 for( i=1; i<*k+1; i++ ){
389 for( j=jcol; j<(*n+*k); j++ )
390 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
416 A = (
double* )malloc( LDA*LDA*
sizeof(
double ) );
417 AP = (
double* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
double ) );
419 for( j=0, k=0; j<*n; j++ )
420 for( i=0; i<j+1; i++, k++ )
421 A[ LDA*i+j ]=ap[ k ];
422 for( i=0, k=0; i<*n; i++ )
423 for( j=i; j<*n; j++, k++ )
424 AP[ k ]=A[ LDA*i+j ];
427 for( j=0, k=0; j<*n; j++ )
428 for( i=j; i<*n; i++, k++ )
429 A[ LDA*i+j ]=ap[ k ];
430 for( i=0, k=0; i<*n; i++ )
431 for( j=0; j<i+1; j++, k++ )
432 AP[ k ]=A[ LDA*i+j ];
462 A = (
double* )malloc( LDA*LDA*
sizeof(
double ) );
463 AP = (
double* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
double ) );
465 for( j=0, k=0; j<*n; j++ )
466 for( i=0; i<j+1; i++, k++ )
467 A[ LDA*i+j ]=ap[ k ];
468 for( i=0, k=0; i<*n; i++ )
469 for( j=i; j<*n; j++, k++ )
470 AP[ k ]=A[ LDA*i+j ];
473 for( j=0, k=0; j<*n; j++ )
474 for( i=j; i<*n; i++, k++ )
475 A[ LDA*i+j ]=ap[ k ];
476 for( i=0, k=0; i<*n; i++ )
477 for( j=0; j<i+1; j++, k++ )
478 AP[ k ]=A[ LDA*i+j ];
506 A = (
double* )malloc( LDA*LDA*
sizeof(
double ) );
507 AP = (
double* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
double ) );
509 for( j=0, k=0; j<*n; j++ )
510 for( i=0; i<j+1; i++, k++ )
511 A[ LDA*i+j ]=ap[ k ];
512 for( i=0, k=0; i<*n; i++ )
513 for( j=i; j<*n; j++, k++ )
514 AP[ k ]=A[ LDA*i+j ];
518 for( j=0, k=0; j<*n; j++ )
519 for( i=j; i<*n; i++, k++ )
520 A[ LDA*i+j ]=ap[ k ];
521 for( i=0, k=0; i<*n; i++ )
522 for( j=0; j<i+1; j++, k++ )
523 AP[ k ]=A[ LDA*i+j ];
547 A = (
double* )malloc( LDA*LDA*
sizeof(
double ) );
548 AP = (
double* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
double ) );
550 for( j=0, k=0; j<*n; j++ )
551 for( i=0; i<j+1; i++, k++ )
552 A[ LDA*i+j ]=ap[ k ];
553 for( i=0, k=0; i<*n; i++ )
554 for( j=i; j<*n; j++, k++ )
555 AP[ k ]=A[ LDA*i+j ];
558 for( j=0, k=0; j<*n; j++ )
559 for( i=j; i<*n; i++, k++ )
560 A[ LDA*i+j ]=ap[ k ];
561 for( i=0, k=0; i<*n; i++ )
562 for( j=0; j<i+1; j++, k++ )
563 AP[ k ]=A[ LDA*i+j ];
567 for( i=0, k=0; i<*n; i++ )
568 for( j=i; j<*n; j++, k++ )
569 A[ LDA*i+j ]=AP[ k ];
570 for( j=0, k=0; j<*n; j++ )
571 for( i=0; i<j+1; i++, k++ )
572 ap[ k ]=A[ LDA*i+j ];
575 for( i=0, k=0; i<*n; i++ )
576 for( j=0; j<i+1; j++, k++ )
577 A[ LDA*i+j ]=AP[ k ];
578 for( j=0, k=0; j<*n; j++ )
579 for( i=j; i<*n; i++, k++ )
580 ap[ k ]=A[ LDA*i+j ];
603 A = (
double* )malloc( LDA*LDA*
sizeof(
double ) );
604 AP = (
double* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
double ) );
606 for( j=0, k=0; j<*n; j++ )
607 for( i=0; i<j+1; i++, k++ )
608 A[ LDA*i+j ]=ap[ k ];
609 for( i=0, k=0; i<*n; i++ )
610 for( j=i; j<*n; j++, k++ )
611 AP[ k ]=A[ LDA*i+j ];
614 for( j=0, k=0; j<*n; j++ )
615 for( i=j; i<*n; i++, k++ )
616 A[ LDA*i+j ]=ap[ k ];
617 for( i=0, k=0; i<*n; i++ )
618 for( j=0; j<i+1; j++, k++ )
619 AP[ k ]=A[ LDA*i+j ];
623 for( i=0, k=0; i<*n; i++ )
624 for( j=i; j<*n; j++, k++ )
625 A[ LDA*i+j ]=AP[ k ];
626 for( j=0, k=0; j<*n; j++ )
627 for( i=0; i<j+1; i++, k++ )
628 ap[ k ]=A[ LDA*i+j ];
631 for( i=0, k=0; i<*n; i++ )
632 for( j=0; j<i+1; j++, k++ )
633 A[ LDA*i+j ]=AP[ k ];
634 for( j=0, k=0; j<*n; j++ )
635 for( i=j; i<*n; i++, k++ )
636 ap[ k ]=A[ LDA*i+j ];
void cblas_dgemv(CBLAS_LAYOUT layout, CBLAS_TRANSPOSE TransA, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT lda, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void cblas_dtbmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const CBLAS_INT K, const double *A, const CBLAS_INT lda, double *X, const CBLAS_INT incX)
void cblas_dsyr(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, double *A, const CBLAS_INT lda)
void cblas_dspr2(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, const double *Y, const CBLAS_INT incY, double *A)
void cblas_dspr(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, double *Ap)
void cblas_dtbsv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const CBLAS_INT K, const double *A, const CBLAS_INT lda, double *X, const CBLAS_INT incX)
void cblas_dsbmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const CBLAS_INT K, const double alpha, const double *A, const CBLAS_INT lda, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void cblas_dtrsv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const double *A, const CBLAS_INT lda, double *X, const CBLAS_INT incX)
void cblas_dger(CBLAS_LAYOUT layout, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, const double *Y, const CBLAS_INT incY, double *A, const CBLAS_INT lda)
void cblas_dspmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *Ap, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void cblas_dgbmv(CBLAS_LAYOUT layout, CBLAS_TRANSPOSE TransA, const CBLAS_INT M, const CBLAS_INT N, const CBLAS_INT KL, const CBLAS_INT KU, const double alpha, const double *A, const CBLAS_INT lda, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void cblas_dtrmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const double *A, const CBLAS_INT lda, double *X, const CBLAS_INT incX)
void cblas_dsymv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT lda, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void cblas_dsyr2(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, const double *Y, const CBLAS_INT incY, double *A, const CBLAS_INT lda)
void cblas_dtpsv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const double *Ap, double *X, const CBLAS_INT incX)
void cblas_dtpmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const double *Ap, double *X, const CBLAS_INT incX)
#define BLAS_FORTRAN_STRLEN_END
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)