F77_cgbmv()

void F77_cgbmv	(	CBLAS_INT *	layout,
		char *	transp,
		CBLAS_INT *	m,
		CBLAS_INT *	n,
		CBLAS_INT *	kl,
		CBLAS_INT *	ku,
		CBLAS_TEST_COMPLEX *	alpha,
		CBLAS_TEST_COMPLEX *	a,
		CBLAS_INT *	lda,
		CBLAS_TEST_COMPLEX *	x,
		CBLAS_INT *	incx,
		CBLAS_TEST_COMPLEX *	beta,
		CBLAS_TEST_COMPLEX *	y,
		CBLAS_INT *	incy )

Definition at line 45 of file c_cblas2.c.

  {
 
  CBLAS_TEST_COMPLEX *A;
  CBLAS_INT i,j,irow,jcol,LDA;
  CBLAS_TRANSPOSE trans;
 
  get_transpose_type(transp, &trans);
  if (*layout == TEST_ROW_MJR) {
     LDA = *ku+*kl+2;
     A=( CBLAS_TEST_COMPLEX* )malloc((*n+*kl)*LDA*sizeof(CBLAS_TEST_COMPLEX));
     for( i=0; i<*ku; i++ ){
        irow=*ku+*kl-i;
        jcol=(*ku)-i;
        for( j=jcol; j<*n; j++ ){
           A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].real;
           A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].imag;
        }
     }
     i=*ku;
     irow=*ku+*kl-i;
     for( j=0; j<*n; j++ ){
        A[ LDA*j+irow ].real=a[ (*lda)*j+i ].real;
        A[ LDA*j+irow ].imag=a[ (*lda)*j+i ].imag;
     }
     for( i=*ku+1; i<*ku+*kl+1; i++ ){
        irow=*ku+*kl-i;
        jcol=i-(*ku);
        for( j=jcol; j<(*n+*kl); j++ ){
           A[ LDA*j+irow ].real=a[ (*lda)*(j-jcol)+i ].real;
           A[ LDA*j+irow ].imag=a[ (*lda)*(j-jcol)+i ].imag;
        }
     }
     cblas_cgbmv( CblasRowMajor, trans, *m, *n, *kl, *ku, alpha, A, LDA, x,
                  *incx, beta, y, *incy );
     free(A);
  }
  else if (*layout == TEST_COL_MJR)
     cblas_cgbmv( CblasColMajor, trans, *m, *n, *kl, *ku, alpha, a, *lda, x,
                  *incx, beta, y, *incy );
  else
     cblas_cgbmv( UNDEFINED, trans, *m, *n, *kl, *ku, alpha, a, *lda, x,
                  *incx, beta, y, *incy );
}

Here is the call graph for this function: