Files
lapack/CBLAS/testing/c_zblas1.c
Sébastien Villemot 776a0223b0 Fix CBLAS tests with 64-bit indexing
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.
2022-11-27 13:18:46 +01:00

75 lines
1.8 KiB
C

/*
* c_zblas1.c
*
* The program is a C wrapper for zcblat1.
*
* Written by Keita Teranishi. 2/11/1998
*
*/
#include "cblas_test.h"
#include "cblas.h"
void F77_zaxpy(const CBLAS_INT *N, const void *alpha, void *X,
const CBLAS_INT *incX, void *Y, const CBLAS_INT *incY)
{
cblas_zaxpy(*N, alpha, X, *incX, Y, *incY);
return;
}
void F77_zcopy(const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY)
{
cblas_zcopy(*N, X, *incX, Y, *incY);
return;
}
void F77_zdotc(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX,
const void *Y, const CBLAS_INT *incY,void *dotc)
{
cblas_zdotc_sub(*N, X, *incX, Y, *incY, dotc);
return;
}
void F77_zdotu(const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY,void *dotu)
{
cblas_zdotu_sub(*N, X, *incX, Y, *incY, dotu);
return;
}
void F77_zdscal(const CBLAS_INT *N, const double *alpha, void *X,
const CBLAS_INT *incX)
{
cblas_zdscal(*N, *alpha, X, *incX);
return;
}
void F77_zscal(const CBLAS_INT *N, const void * *alpha, void *X,
const CBLAS_INT *incX)
{
cblas_zscal(*N, alpha, X, *incX);
return;
}
void F77_zswap( const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY)
{
cblas_zswap(*N,X,*incX,Y,*incY);
return;
}
CBLAS_INT F77_izamax(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
{
if (*N < 1 || *incX < 1) return(0);
return(cblas_izamax(*N, X, *incX)+1);
}
double F77_dznrm2(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
{
return cblas_dznrm2(*N, X, *incX);
}
double F77_dzasum(const CBLAS_INT *N, void *X, const CBLAS_INT *incX)
{
return cblas_dzasum(*N, X, *incX);
}