cblas_dsymv.c

/*
 *
 * cblas_dsymv.c
 * This program is a C interface to dsymv.
 * Written by Keita Teranishi
 * 4/6/1998
 *
 */

#include"cblas.h"
#include"cblas_f77.h"
voidcblas_dsymv(constCBLAS_LAYOUTlayout,
constCBLAS_UPLOUplo, constintN,
constdoublealpha, constdouble*A, constintlda,
constdouble*X, constintincX, constdoublebeta,
double*Y, constintincY)
{
charUL;
#ifdefF77_CHAR
F77_CHARF77_UL;
#else
#defineF77_UL &UL
#endif
#ifdefF77_INT
F77_INTF77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY;
#else
#defineF77_N N
#defineF77_lda lda
#defineF77_incX incX
#defineF77_incY incY
#endif
externintCBLAS_CallFromC;
externintRowMajorStrg;
RowMajorStrg=0;

CBLAS_CallFromC=1;
if (layout==CblasColMajor)
 {
if (Uplo==CblasUpper) UL='U';
elseif (Uplo==CblasLower) UL='L';
else
 {
cblas_xerbla(2, "cblas_dsymv","Illegal Uplo setting, %d\n",Uplo );
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
 }
#ifdefF77_CHAR
F77_UL=C2F_CHAR(&UL);
#endif
F77_dsymv(F77_UL, &F77_N, &alpha, A, &F77_lda, X,
&F77_incX, &beta, Y, &F77_incY);
 }
elseif (layout==CblasRowMajor)
 {
RowMajorStrg=1;
if (Uplo==CblasUpper) UL='L';
elseif (Uplo==CblasLower) UL='U';
else
 {
cblas_xerbla(2, "cblas_dsymv","Illegal Uplo setting, %d\n", Uplo);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
 }
#ifdefF77_CHAR
F77_UL=C2F_CHAR(&UL);
#endif
F77_dsymv(F77_UL, &F77_N, &alpha,
A ,&F77_lda, X,&F77_incX, &beta, Y, &F77_incY);
 }
elsecblas_xerbla(1, "cblas_dsymv", "Illegal layout setting, %d\n", layout);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}