181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/* 281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_dgbmv.c 481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * This program is a C interface to dgbmv. 581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * Written by Keita Teranishi 681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 4/6/1998 781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray */ 981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas.h" 1081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas_f77.h" 1181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murrayvoid cblas_dgbmv(const enum CBLAS_ORDER order, 1281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const enum CBLAS_TRANSPOSE TransA, const int M, const int N, 1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const int KL, const int KU, 1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const double alpha, const double *A, const int lda, 1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const double *X, const int incX, const double beta, 1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray double *Y, const int incY) 1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{ 1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray char TA; 1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_CHAR 2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_CHAR F77_TA; 2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_TA &TA 2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 2481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT 2581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY; 2681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_KL=KL,F77_KU=KU; 2781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_M M 2981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_N N 3081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_lda lda 3181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_KL KL 3281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_KU KU 3381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incX incX 3481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incY incY 3581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 3681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int CBLAS_CallFromC; 3781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int RowMajorStrg; 3881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 3981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 4081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 1; 4181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (order == CblasColMajor) 4281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasNoTrans) TA = 'N'; 4481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasTrans) TA = 'T'; 4581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasConjTrans) TA = 'C'; 4681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 4781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_dgbmv","Illegal TransA setting, %d\n", TransA); 4981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 5081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 5181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 5281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 5481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 5581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 5681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_dgbmv(F77_TA, &F77_M, &F77_N, &F77_KL, &F77_KU, &alpha, 5781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray A, &F77_lda, X, &F77_incX, &beta, Y, &F77_incY); 5881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (order == CblasRowMajor) 6081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 6181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 1; 6281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasNoTrans) TA = 'T'; 6381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasTrans) TA = 'N'; 6481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasConjTrans) TA = 'N'; 6581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 6681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 6781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_dgbmv","Illegal TransA setting, %d\n", TransA); 6881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 6981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 7081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 7181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 7281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 7381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 7481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 7581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_dgbmv(F77_TA, &F77_N, &F77_M, &F77_KU, &F77_KL, &alpha, 7681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray A ,&F77_lda, X,&F77_incX, &beta, Y, &F77_incY); 7781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 7881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else cblas_xerbla(1, "cblas_dgbmv", "Illegal Order setting, %d\n", order); 7981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 8081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 8181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray} 82