F77_dsbmv()

void F77_dsbmv	(	CBLAS_INT *	layout,
		char *	uplow,
		CBLAS_INT *	n,
		CBLAS_INT *	k,
		double *	alpha,
		double *	a,
		CBLAS_INT *	lda,
		double *	x,
		CBLAS_INT *	incx,
		double *	beta,
		double *	y,
		CBLAS_INT *	incy )

Definition at line 353 of file c_dblas2.c.

  {
  double *A;
  CBLAS_INT i,j,irow,jcol,LDA;
  CBLAS_UPLO uplo;
 
  get_uplo_type(uplow,&uplo);
 
  if (*layout == TEST_ROW_MJR) {
     LDA = *k+1;
     A   = ( double* )malloc( (*n+*k)*LDA*sizeof( double ) );
     if (uplo == CblasUpper) {
        for( i=0; i<*k; i++ ){
           irow=*k-i;
           jcol=(*k)-i;
           for( j=jcol; j<*n; j++ )
        A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
        }
        i=*k;
        irow=*k-i;
        for( j=0; j<*n; j++ )
           A[ LDA*j+irow ]=a[ (*lda)*j+i ];
     }
     else {
        i=0;
        irow=*k-i;
        for( j=0; j<*n; j++ )
           A[ LDA*j+irow ]=a[ (*lda)*j+i ];
        for( i=1; i<*k+1; i++ ){
           irow=*k-i;
           jcol=i;
           for( j=jcol; j<(*n+*k); j++ )
              A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
        }
     }
     cblas_dsbmv(CblasRowMajor, uplo, *n, *k, *alpha, A, LDA, x, *incx,
                 *beta, y, *incy );
     free(A);
   }
   else
     cblas_dsbmv(CblasColMajor, uplo, *n, *k, *alpha, a, *lda, x, *incx,
                 *beta, y, *incy );
}

Here is the call graph for this function: