181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/* 281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_cgemm.c 481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * This program is a C interface to cgemm. 581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * Written by Keita Teranishi 681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 4/8/1998 781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray */ 981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 1081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas.h" 1181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas_f77.h" 1281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murrayvoid cblas_cgemm(const enum CBLAS_ORDER Order, const enum CBLAS_TRANSPOSE TransA, 1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const enum CBLAS_TRANSPOSE TransB, const int M, const int N, 1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const int K, const void *alpha, const void *A, 1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const int lda, const void *B, const int ldb, 1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const void *beta, void *C, const int ldc) 1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{ 1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray char TA, TB; 1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_CHAR 2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_CHAR F77_TA, F77_TB; 2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_TA &TA 2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_TB &TB 2481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 2581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 2681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT 2781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_M=M, F77_N=N, F77_K=K, F77_lda=lda, F77_ldb=ldb; 2881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_ldc=ldc; 2981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 3081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_M M 3181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_N N 3281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_K K 3381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_lda lda 3481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_ldb ldb 3581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_ldc ldc 3681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 3781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 3881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int CBLAS_CallFromC; 3981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int RowMajorStrg; 4081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 4181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 1; 4281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 4381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Order == CblasColMajor ) 4481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(TransA == CblasTrans) TA='T'; 4681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransA == CblasConjTrans ) TA='C'; 4781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransA == CblasNoTrans ) TA='N'; 4881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 4981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 5081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_cgemm", "Illegal TransA setting, %d\n", TransA); 5181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 5281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 5381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 5481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 5681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(TransB == CblasTrans) TB='T'; 5781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransB == CblasConjTrans ) TB='C'; 5881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransB == CblasNoTrans ) TB='N'; 5981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 6081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 6181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(3, "cblas_cgemm", "Illegal TransB setting, %d\n", TransB); 6281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 6381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 6481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 6581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 6681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 6781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 6881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 6981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TB = C2F_CHAR(&TB); 7081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 7181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 7281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_cgemm(F77_TA, F77_TB, &F77_M, &F77_N, &F77_K, alpha, A, 7381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &F77_lda, B, &F77_ldb, beta, C, &F77_ldc); 7481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } else if (Order == CblasRowMajor) 7581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 1; 7781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(TransA == CblasTrans) TB='T'; 7881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransA == CblasConjTrans ) TB='C'; 7981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransA == CblasNoTrans ) TB='N'; 8081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 8181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 8281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_cgemm", "Illegal TransA setting, %d\n", TransA); 8381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 8481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 8581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 8681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 8781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(TransB == CblasTrans) TA='T'; 8881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransB == CblasConjTrans ) TA='C'; 8981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( TransB == CblasNoTrans ) TA='N'; 9081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 9181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 9281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_cgemm", "Illegal TransB setting, %d\n", TransB); 9381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 9481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 9581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 9681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 9781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 9881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 9981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TB = C2F_CHAR(&TB); 10081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 10181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 10281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_cgemm(F77_TA, F77_TB, &F77_N, &F77_M, &F77_K, alpha, B, 10381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &F77_ldb, A, &F77_lda, beta, C, &F77_ldc); 10481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 10581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else cblas_xerbla(1, "cblas_cgemm", "Illegal Order setting, %d\n", Order); 10681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 10781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 10881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 10981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray} 110