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

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106
  1. /*
  2. *
  3. * cblas_csymm.c
  4. * This program is a C interface to csymm.
  5. * Written by Keita Teranishi
  6. * 4/8/1998
  7. *
  8. */
  9. #include "cblas.h"
  10. #include "cblas_f77.h"
  11. void cblas_csymm(const CBLAS_LAYOUT layout, const CBLAS_SIDE Side,
  12. const CBLAS_UPLO Uplo, const int M, const int N,
  13. const void *alpha, const void *A, const int lda,
  14. const void *B, const int ldb, const void *beta,
  15. void *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_csymm", "Illegal Side setting, %d\n", Side);
  45. CBLAS_CallFromC = 0;
  46. RowMajorStrg = 0;
  47. return;
  48. }
  49. if( Uplo == CblasUpper) UL='U';
  50. else if ( Uplo == CblasLower ) UL='L';
  51. else
  52. {
  53. cblas_xerbla(3, "cblas_csymm", "Illegal Uplo setting, %d\n", Uplo);
  54. CBLAS_CallFromC = 0;
  55. RowMajorStrg = 0;
  56. return;
  57. }
  58. #ifdef F77_CHAR
  59. F77_UL = C2F_CHAR(&UL);
  60. F77_SD = C2F_CHAR(&SD);
  61. #endif
  62. F77_csymm(F77_SD, F77_UL, &F77_M, &F77_N, alpha, A, &F77_lda,
  63. B, &F77_ldb, beta, C, &F77_ldc);
  64. } else if (layout == CblasRowMajor)
  65. {
  66. RowMajorStrg = 1;
  67. if( Side == CblasRight) SD='L';
  68. else if ( Side == CblasLeft ) SD='R';
  69. else
  70. {
  71. cblas_xerbla(2, "cblas_csymm", "Illegal Side setting, %d\n", Side);
  72. CBLAS_CallFromC = 0;
  73. RowMajorStrg = 0;
  74. return;
  75. }
  76. if( Uplo == CblasUpper) UL='L';
  77. else if ( Uplo == CblasLower ) UL='U';
  78. else
  79. {
  80. cblas_xerbla(3, "cblas_csymm", "Illegal Uplo setting, %d\n", Uplo);
  81. CBLAS_CallFromC = 0;
  82. RowMajorStrg = 0;
  83. return;
  84. }
  85. #ifdef F77_CHAR
  86. F77_UL = C2F_CHAR(&UL);
  87. F77_SD = C2F_CHAR(&SD);
  88. #endif
  89. F77_csymm(F77_SD, F77_UL, &F77_N, &F77_M, alpha, A, &F77_lda,
  90. B, &F77_ldb, beta, C, &F77_ldc);
  91. }
  92. else cblas_xerbla(1, "cblas_csymm", "Illegal layout setting, %d\n", layout);
  93. CBLAS_CallFromC = 0;
  94. RowMajorStrg = 0;
  95. return;
  96. }