181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/* 281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_zherk.c 481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * This program is a C interface to zherk. 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_zherk(const enum CBLAS_ORDER Order, const enum CBLAS_UPLO Uplo, 1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const enum CBLAS_TRANSPOSE Trans, const int N, const int K, 1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const double alpha, const void *A, const int lda, 1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const double beta, void *C, const int ldc) 1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{ 1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray char UL, TR; 1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_CHAR 1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_CHAR F77_TR, F77_UL; 2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_TR &TR 2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_UL &UL 2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 2481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 2581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT 2681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_N=N, F77_K=K, F77_lda=lda; 2781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_ldc=ldc; 2881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_N N 3081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_K K 3181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_lda lda 3281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_ldc ldc 3381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 3481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 3581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int CBLAS_CallFromC; 3681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int RowMajorStrg; 3781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 3881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 1; 3981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 4081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Order == CblasColMajor ) 4181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Uplo == CblasUpper) UL='U'; 4381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Uplo == CblasLower ) UL='L'; 4481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 4581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_zherk", "Illegal Uplo setting, %d\n", Uplo); 4781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 4881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 4981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 5081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 5281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Trans == CblasTrans) TR ='T'; 5381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Trans == CblasConjTrans ) TR='C'; 5481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Trans == CblasNoTrans ) TR='N'; 5581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 5681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 5781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(3, "cblas_zherk", "Illegal Trans setting, %d\n", Trans); 5881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 5981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 6081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 6181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 6281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 6381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 6481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_UL = C2F_CHAR(&UL); 6581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TR = C2F_CHAR(&TR); 6681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 6781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 6881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_zherk(F77_UL, F77_TR, &F77_N, &F77_K, &alpha, A, &F77_lda, 6981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &beta, C, &F77_ldc); 7081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } else if (Order == CblasRowMajor) 7181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 1; 7381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Uplo == CblasUpper) UL='L'; 7481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Uplo == CblasLower ) UL='U'; 7581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 7681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(3, "cblas_zherk", "Illegal Uplo setting, %d\n", Uplo); 7881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 7981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 8081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 8181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 8281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( Trans == CblasTrans) TR ='N'; 8381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Trans == CblasConjTrans ) TR='N'; 8481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if ( Trans == CblasNoTrans ) TR='C'; 8581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 8681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 8781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(3, "cblas_zherk", "Illegal Trans setting, %d\n", Trans); 8881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 8981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 9081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 9181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 9281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 9381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 9481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_UL = C2F_CHAR(&UL); 9581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_SD = C2F_CHAR(&SD); 9681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 9781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 9881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_zherk(F77_UL, F77_TR, &F77_N, &F77_K, &alpha, A, &F77_lda, 9981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &beta, C, &F77_ldc); 10081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 10181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else cblas_xerbla(1, "cblas_zherk", "Illegal Order setting, %d\n", Order); 10281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 10381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 10481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 10581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray} 106