F77_sgbmv()

void F77_sgbmv	(	int *	layout,
		char *	transp,
		int *	m,
		int *	n,
		int *	kl,
		int *	ku,
		float *	alpha,
		float *	a,
		int *	lda,
		float *	x,
		int *	incx,
		float *	beta,
		float *	y,
		int *	incy
	)

Definition at line 187 of file c_sblas2.c.

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

Here is the call graph for this function: