You can not select more than 25 topics Topics must start with a chinese character,a letter or number, can include dashes ('-') and can be up to 35 characters long.

cblas_sgbmv.c 2.3 kB

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