F77_dgbmv()

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

Definition at line 187 of file c_dblas2.c.

                                                          {
 
  double *A;
  CBLAS_INT i,irow,j,jcol,LDA;
  CBLAS_TRANSPOSE trans;
 
  get_transpose_type(transp, &trans);
 
  if (*layout == TEST_ROW_MJR) {
     LDA = *ku+*kl+2;
     A   = ( double* )malloc( (*n+*kl)*LDA*sizeof( double ) );
     for( i=0; i<*ku; i++ ){
        irow=*ku+*kl-i;
        jcol=(*ku)-i;
        for( j=jcol; j<*n; j++ )
           A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
     }
     i=*ku;
     irow=*ku+*kl-i;
     for( j=0; j<*n; j++ )
        A[ LDA*j+irow ]=a[ (*lda)*j+i ];
     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 ]=a[ (*lda)*(j-jcol)+i ];
     }
     cblas_dgbmv( CblasRowMajor, trans, *m, *n, *kl, *ku, *alpha,
                  A, LDA, x, *incx, *beta, y, *incy );
     free(A);
  }
  else
     cblas_dgbmv( CblasColMajor, trans, *m, *n, *kl, *ku, *alpha,
                  a, *lda, x, *incx, *beta, y, *incy );
}

Here is the call graph for this function: