cblas_cher2k.c File Reference

#include "cblas.h"
#include "cblas_f77.h"

Go to the source code of this file.

Macros
#define	F77_TR   &TR
#define	F77_UL   &UL
#define	F77_N   N
#define	F77_K   K
#define	F77_lda   lda
#define	F77_ldb   ldb
#define	F77_ldc   ldc

Functions
void	cblas_cher2k (const CBLAS_LAYOUT layout, const CBLAS_UPLO Uplo, const CBLAS_TRANSPOSE Trans, const CBLAS_INT N, const CBLAS_INT K, const void *alpha, const void *A, const CBLAS_INT lda, const void *B, const CBLAS_INT ldb, const float beta, void *C, const CBLAS_INT ldc)

Macro Definition Documentation

F77_K

#define F77_K   K

F77_lda

#define F77_lda   lda

F77_ldb

#define F77_ldb   ldb

F77_ldc

#define F77_ldc   ldc

F77_N

#define F77_N   N

F77_TR

#define F77_TR   &TR

F77_UL

#define F77_UL   &UL

Function Documentation

cblas_cher2k()

void cblas_cher2k	(	const CBLAS_LAYOUT	layout,
		const CBLAS_UPLO	Uplo,
		const CBLAS_TRANSPOSE	Trans,
		const CBLAS_INT	N,
		const CBLAS_INT	K,
		const void *	alpha,
		const void *	A,
		const CBLAS_INT	lda,
		const void *	B,
		const CBLAS_INT	ldb,
		const float	beta,
		void *	C,
		const CBLAS_INT	ldc )

Definition at line 12 of file cblas_cher2k.c.

{
   char UL, TR;
#ifdef F77_CHAR
   F77_CHAR F77_TR, F77_UL;
#else
   #define F77_TR &TR
   #define F77_UL &UL
#endif
 
#ifdef F77_INT
   F77_INT F77_N=N, F77_K=K, F77_lda=lda, F77_ldb=ldb;
   F77_INT F77_ldc=ldc;
#else
   #define F77_N N
   #define F77_K K
   #define F77_lda lda
   #define F77_ldb ldb
   #define F77_ldc ldc
#endif
 
   extern int CBLAS_CallFromC;
   extern int RowMajorStrg;
   float ALPHA[2];
   const float *alp=(float *)alpha;
 
   CBLAS_CallFromC = 1;
   RowMajorStrg = 0;
 
   if( layout == CblasColMajor )
   {
 
      if( Uplo == CblasUpper) UL='U';
      else if ( Uplo == CblasLower ) UL='L';
      else
      {
         cblas_xerbla(2, "cblas_cher2k", "Illegal Uplo setting, %d\n", Uplo);
         CBLAS_CallFromC = 0;
         RowMajorStrg = 0;
         return;
      }
 
      if( Trans == CblasTrans) TR ='T';
      else if ( Trans == CblasConjTrans ) TR='C';
      else if ( Trans == CblasNoTrans )   TR='N';
      else
      {
         cblas_xerbla(3, "cblas_cher2k", "Illegal Trans setting, %d\n", Trans);
         CBLAS_CallFromC = 0;
         RowMajorStrg = 0;
         return;
      }
 
      #ifdef F77_CHAR
         F77_UL = C2F_CHAR(&UL);
         F77_TR = C2F_CHAR(&TR);
      #endif
 
      F77_cher2k(F77_UL, F77_TR, &F77_N, &F77_K, alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
   } else if (layout == CblasRowMajor)
   {
      RowMajorStrg = 1;
 
      if( Uplo == CblasUpper) UL='L';
      else if ( Uplo == CblasLower ) UL='U';
      else
      {
         cblas_xerbla(2, "cblas_cher2k", "Illegal Uplo setting, %d\n", Uplo);
         CBLAS_CallFromC = 0;
         RowMajorStrg = 0;
         return;
      }
      if( Trans == CblasTrans) TR ='N';
      else if ( Trans == CblasConjTrans ) TR='N';
      else if ( Trans == CblasNoTrans )   TR='C';
      else
      {
         cblas_xerbla(3, "cblas_cher2k", "Illegal Trans setting, %d\n", Trans);
         CBLAS_CallFromC = 0;
         RowMajorStrg = 0;
         return;
      }
      #ifdef F77_CHAR
         F77_UL = C2F_CHAR(&UL);
         F77_TR = C2F_CHAR(&TR);
      #endif
 
      ALPHA[0]= *alp;
      ALPHA[1]= -alp[1];
      F77_cher2k(F77_UL,F77_TR, &F77_N, &F77_K, ALPHA, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
   }
   else  cblas_xerbla(1, "cblas_cher2k", "Illegal layout setting, %d\n", layout);
   CBLAS_CallFromC = 0;
   RowMajorStrg = 0;
   return;
}