1/*
2 *
3 * cblas_sgbmv.c
4 * This program is a C interface to sgbmv.
5 * Written by Keita Teranishi
6 * 4/6/1998
7 *
8 */
9
10#include "cblas.h"
11#include "cblas_f77.h"
12void cblas_sgbmv(const enum CBLAS_ORDER order,
13                 const enum CBLAS_TRANSPOSE TransA, const int M, const int N,
14                 const int KL, const int KU,
15                 const float alpha, const float *A, const int lda,
16                 const float  *X, const int incX, const float beta,
17                 float  *Y, const int incY)
18{
19   char TA;
20#ifdef F77_CHAR
21   F77_CHAR F77_TA;
22#else
23   #define F77_TA &TA
24#endif
25#ifdef F77_INT
26   F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY;
27   F77_INT F77_KL=KL,F77_KU=KU;
28#else
29   #define F77_M M
30   #define F77_N N
31   #define F77_lda lda
32   #define F77_KL KL
33   #define F77_KU KU
34   #define F77_incX incX
35   #define F77_incY incY
36#endif
37   extern int CBLAS_CallFromC;
38   extern int RowMajorStrg;
39   RowMajorStrg = 0;
40
41   CBLAS_CallFromC = 1;
42   if (order == CblasColMajor)
43   {
44      if (TransA == CblasNoTrans) TA = 'N';
45      else if (TransA == CblasTrans) TA = 'T';
46      else if (TransA == CblasConjTrans) TA = 'C';
47      else
48      {
49         cblas_xerbla(2, "cblas_sgbmv","Illegal TransA setting, %d\n", TransA);
50         CBLAS_CallFromC = 0;
51         RowMajorStrg = 0;
52         return;
53      }
54      #ifdef F77_CHAR
55         F77_TA = C2F_CHAR(&TA);
56      #endif
57      F77_sgbmv(F77_TA, &F77_M, &F77_N, &F77_KL, &F77_KU, &alpha,
58                     A, &F77_lda, X, &F77_incX, &beta, Y, &F77_incY);
59   }
60   else if (order == CblasRowMajor)
61   {
62      RowMajorStrg = 1;
63      if (TransA == CblasNoTrans) TA = 'T';
64      else if (TransA == CblasTrans) TA = 'N';
65      else if (TransA == CblasConjTrans) TA = 'N';
66      else
67      {
68         cblas_xerbla(2, "cblas_sgbmv","Illegal TransA setting, %d\n", TransA);
69         CBLAS_CallFromC = 0;
70         RowMajorStrg = 0;
71         return;
72      }
73      #ifdef F77_CHAR
74         F77_TA = C2F_CHAR(&TA);
75      #endif
76      F77_sgbmv(F77_TA, &F77_N, &F77_M, &F77_KU, &F77_KL, &alpha,
77                     A ,&F77_lda, X, &F77_incX, &beta, Y, &F77_incY);
78   }
79   else cblas_xerbla(1, "cblas_sgbmv", "Illegal Order setting, %d\n", order);
80   CBLAS_CallFromC = 0;
81   RowMajorStrg = 0;
82   return;
83}
84