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_ssymm.c 2.7 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108
  1. /*
  2. *
  3. * cblas_ssymm.c
  4. * This program is a C interface to ssymm.
  5. * Written by Keita Teranishi
  6. * 4/8/1998
  7. *
  8. */
  9. #include "cblas.h"
  10. #include "cblas_f77.h"
  11. void cblas_ssymm(const CBLAS_LAYOUT layout, const CBLAS_SIDE Side,
  12. const CBLAS_UPLO Uplo, const int M, const int N,
  13. const float alpha, const float *A, const int lda,
  14. const float *B, const int ldb, const float beta,
  15. float *C, const int ldc)
  16. {
  17. char SD, UL;
  18. #ifdef F77_CHAR
  19. F77_CHAR F77_SD, F77_UL;
  20. #else
  21. #define F77_SD &SD
  22. #define F77_UL &UL
  23. #endif
  24. #ifdef F77_INT
  25. F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_ldb=ldb;
  26. F77_INT F77_ldc=ldc;
  27. #else
  28. #define F77_M M
  29. #define F77_N N
  30. #define F77_lda lda
  31. #define F77_ldb ldb
  32. #define F77_ldc ldc
  33. #endif
  34. extern int CBLAS_CallFromC;
  35. extern int RowMajorStrg;
  36. RowMajorStrg = 0;
  37. CBLAS_CallFromC = 1;
  38. if( layout == CblasColMajor )
  39. {
  40. if( Side == CblasRight) SD='R';
  41. else if ( Side == CblasLeft ) SD='L';
  42. else
  43. {
  44. cblas_xerbla(2, "cblas_ssymm",
  45. "Illegal Side setting, %d\n", Side);
  46. CBLAS_CallFromC = 0;
  47. RowMajorStrg = 0;
  48. return;
  49. }
  50. if( Uplo == CblasUpper) UL='U';
  51. else if ( Uplo == CblasLower ) UL='L';
  52. else
  53. {
  54. cblas_xerbla(3, "cblas_ssymm",
  55. "Illegal Uplo setting, %d\n", Uplo);
  56. CBLAS_CallFromC = 0;
  57. RowMajorStrg = 0;
  58. return;
  59. }
  60. #ifdef F77_CHAR
  61. F77_UL = C2F_CHAR(&UL);
  62. F77_SD = C2F_CHAR(&SD);
  63. #endif
  64. F77_ssymm(F77_SD, F77_UL, &F77_M, &F77_N, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
  65. } else if (layout == CblasRowMajor)
  66. {
  67. RowMajorStrg = 1;
  68. if( Side == CblasRight) SD='L';
  69. else if ( Side == CblasLeft ) SD='R';
  70. else
  71. {
  72. cblas_xerbla(2, "cblas_ssymm",
  73. "Illegal Side setting, %d\n", Side);
  74. CBLAS_CallFromC = 0;
  75. RowMajorStrg = 0;
  76. return;
  77. }
  78. if( Uplo == CblasUpper) UL='L';
  79. else if ( Uplo == CblasLower ) UL='U';
  80. else
  81. {
  82. cblas_xerbla(3, "cblas_ssymm",
  83. "Illegal Uplo setting, %d\n", Uplo);
  84. CBLAS_CallFromC = 0;
  85. RowMajorStrg = 0;
  86. return;
  87. }
  88. #ifdef F77_CHAR
  89. F77_UL = C2F_CHAR(&UL);
  90. F77_SD = C2F_CHAR(&SD);
  91. #endif
  92. F77_ssymm(F77_SD, F77_UL, &F77_N, &F77_M, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
  93. } else cblas_xerbla(1, "cblas_ssymm",
  94. "Illegal layout setting, %d\n", layout);
  95. CBLAS_CallFromC = 0;
  96. RowMajorStrg = 0;
  97. return;
  98. }