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_dgbmv.c 2.3 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081
  1. /*
  2. *
  3. * cblas_dgbmv.c
  4. * This program is a C interface to dgbmv.
  5. * Written by Keita Teranishi
  6. * 4/6/1998
  7. *
  8. */
  9. #include "cblas.h"
  10. #include "cblas_f77.h"
  11. void cblas_dgbmv(const CBLAS_LAYOUT layout,
  12. const CBLAS_TRANSPOSE TransA, const int M, const int N,
  13. const int KL, const int KU,
  14. const double alpha, const double *A, const int lda,
  15. const double *X, const int incX, const double beta,
  16. double *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_dgbmv","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_dgbmv(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_dgbmv","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_dgbmv(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_dgbmv", "Illegal layout setting, %d\n", layout);
  78. CBLAS_CallFromC = 0;
  79. RowMajorStrg = 0;
  80. }