Sébastien Villemot
776a0223b0
This patch replaces all instances of “int” by “CBLAS_INT” in CBLAS/testing/*.c (except for RowMajorStrg which has must remain an “int”, since it is declared as such in CBLAS main code). Interestingly, the bug would not manifest on low endian architectures, because interpreting an int64_t pointer as an int32_t pointer works there (as long as the integers fits in 31 bits, which is the case here). But on big endian architectures, this of course fails badly.
84 lines
2.0 KiB
C
84 lines
2.0 KiB
C
/*
|
|
* c_dblas1.c
|
|
*
|
|
* The program is a C wrapper for dcblat1.
|
|
*
|
|
* Written by Keita Teranishi. 2/11/1998
|
|
*
|
|
*/
|
|
#include "cblas_test.h"
|
|
#include "cblas.h"
|
|
double F77_dasum(const CBLAS_INT *N, double *X, const CBLAS_INT *incX)
|
|
{
|
|
return cblas_dasum(*N, X, *incX);
|
|
}
|
|
|
|
void F77_daxpy(const CBLAS_INT *N, const double *alpha, const double *X,
|
|
const CBLAS_INT *incX, double *Y, const CBLAS_INT *incY)
|
|
{
|
|
cblas_daxpy(*N, *alpha, X, *incX, Y, *incY);
|
|
return;
|
|
}
|
|
|
|
void F77_dcopy(const CBLAS_INT *N, double *X, const CBLAS_INT *incX,
|
|
double *Y, const CBLAS_INT *incY)
|
|
{
|
|
cblas_dcopy(*N, X, *incX, Y, *incY);
|
|
return;
|
|
}
|
|
|
|
double F77_ddot(const CBLAS_INT *N, const double *X, const CBLAS_INT *incX,
|
|
const double *Y, const CBLAS_INT *incY)
|
|
{
|
|
return cblas_ddot(*N, X, *incX, Y, *incY);
|
|
}
|
|
|
|
double F77_dnrm2(const CBLAS_INT *N, const double *X, const CBLAS_INT *incX)
|
|
{
|
|
return cblas_dnrm2(*N, X, *incX);
|
|
}
|
|
|
|
void F77_drotg( double *a, double *b, double *c, double *s)
|
|
{
|
|
cblas_drotg(a,b,c,s);
|
|
return;
|
|
}
|
|
|
|
void F77_drot( const CBLAS_INT *N, double *X, const CBLAS_INT *incX, double *Y,
|
|
const CBLAS_INT *incY, const double *c, const double *s)
|
|
{
|
|
|
|
cblas_drot(*N,X,*incX,Y,*incY,*c,*s);
|
|
return;
|
|
}
|
|
|
|
void F77_dscal(const CBLAS_INT *N, const double *alpha, double *X,
|
|
const CBLAS_INT *incX)
|
|
{
|
|
cblas_dscal(*N, *alpha, X, *incX);
|
|
return;
|
|
}
|
|
|
|
void F77_dswap( const CBLAS_INT *N, double *X, const CBLAS_INT *incX,
|
|
double *Y, const CBLAS_INT *incY)
|
|
{
|
|
cblas_dswap(*N,X,*incX,Y,*incY);
|
|
return;
|
|
}
|
|
|
|
double F77_dzasum(const CBLAS_INT *N, void *X, const CBLAS_INT *incX)
|
|
{
|
|
return cblas_dzasum(*N, X, *incX);
|
|
}
|
|
|
|
double F77_dznrm2(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
|
|
{
|
|
return cblas_dznrm2(*N, X, *incX);
|
|
}
|
|
|
|
CBLAS_INT F77_idamax(const CBLAS_INT *N, const double *X, const CBLAS_INT *incX)
|
|
{
|
|
if (*N < 1 || *incX < 1) return(0);
|
|
return (cblas_idamax(*N, X, *incX)+1);
|
|
}
|