floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Merge pull request #4817 from martin-frbg/fix4807 

Fix SCAL on x86 and RISCV_GENERIC
This commit is contained in:
Martin Kroeker 2024-07-26 14:56:44 +02:00 committed by GitHub
parent fb7c53c5e5 a090011fbf
commit 25e148ec58
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
3 changed files with 28 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
kernel/riscv64/scal.c
View File

@ -43,7 +43,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x, BLAS
if ( (n <= 0) || (inc_x <= 0))
return(0);
if (dummy2 == 0) {
if (dummy2 == 1) {
while(j < n)
{

2
kernel/x86/scal.S
View File

@ -65,7 +65,7 @@
#else
movl 32(%esp),%edi
movl 36(%esp),%esi
movl 54(%esp),%ecx
movl 52(%esp),%ecx
#endif
ftst

48
utest/test_gemv.c
View File

@ -12,47 +12,49 @@
CTEST(sgemv, 0_nan_inf)
{
int i;
blasint N = 17;
blasint incX = 1;
blasint incY = 1;
float alpha = 0.0;
float beta = 0.0;
char trans = 'N';
float A[N * N];
float X[N];
float Y[N];
float A[17 * 17];
float X[17];
float Y[17];
memset(A, 0, sizeof(A));
memset(X, 0, sizeof(X));
for (int i = 0; i < (N - 1); i += 2)
for (i = 0; i < (N - 1); i += 2)
{
Y[i] = NAN;
Y[i + 1] = INFINITY;
}
Y[N - 1] = NAN;
BLASFUNC(sgemv)(&trans, &N, &N, &alpha, A, &N, X, &incX, &beta, Y, &incY);
for (int i = 0; i < N; i ++)
for (i = 0; i < N; i ++)
ASSERT_TRUE(Y[i] == 0.0);
}
CTEST(sgemv, 0_nan_inf_incy_2)
{
blasint N = 17;
int i;
blasint N = 17;
blasint Ny = 33;
blasint incX = 1;
blasint incY = 2;
float alpha = 0.0;
float beta = 0.0;
char trans = 'N';
float A[N * N];
float X[N];
float Y[Ny];
float A[17 * 17];
float X[17];
float Y[33];
float *ay = Y;
memset(A, 0, sizeof(A));
memset(X, 0, sizeof(X));
memset(Y, 0, sizeof(Y));
for (int i = 0; i < (N - 1); i += 2)
for (i = 0; i < (N - 1); i += 2)
{
ay[0] = NAN;
ay += 2;
@ -61,7 +63,7 @@ CTEST(sgemv, 0_nan_inf_incy_2)
}
Y[Ny - 1] = NAN;
BLASFUNC(sgemv)(&trans, &N, &N, &alpha, A, &N, X, &incX, &beta, Y, &incY);
for (int i = 0; i < Ny; i ++)
for (i = 0; i < Ny; i ++)
ASSERT_TRUE(Y[i] == 0.0);
}
@ -70,47 +72,49 @@ CTEST(sgemv, 0_nan_inf_incy_2)
#ifdef BUILD_DOUBLE
CTEST(dgemv, 0_nan_inf)
{
int i;
blasint N = 17;
blasint incX = 1;
blasint incY = 1;
double alpha = 0.0;
double beta = 0.0;
char trans = 'N';
double A[N * N];
double X[N];
double Y[N];
double A[17 * 17];
double X[17];
double Y[17];
memset(A, 0, sizeof(A));
memset(X, 0, sizeof(X));
for (int i = 0; i < (N - 1); i += 2)
for (i = 0; i < (N - 1); i += 2)
{
Y[i] = NAN;
Y[i + 1] = INFINITY;
}
Y[N - 1] = NAN;
BLASFUNC(dgemv)(&trans, &N, &N, &alpha, A, &N, X, &incX, &beta, Y, &incY);
for (int i = 0; i < N; i ++)
for (i = 0; i < N; i ++)
ASSERT_TRUE(Y[i] == 0.0);
}
CTEST(dgemv, 0_nan_inf_incy_2)
{
blasint N = 17;
int i;
blasint N = 17;
blasint Ny = 33;
blasint incX = 1;
blasint incY = 2;
double alpha = 0.0;
double beta = 0.0;
char trans = 'N';
double A[N * N];
double X[N];
double Y[Ny];
double A[17 * 17];
double X[17];
double Y[33];
double *ay = Y;
memset(A, 0, sizeof(A));
memset(X, 0, sizeof(X));
memset(Y, 0, sizeof(Y));
for (int i = 0; i < (N - 1); i += 2)
for (i = 0; i < (N - 1); i += 2)
{
ay[0] = NAN;
ay += 2;
@ -119,7 +123,7 @@ CTEST(dgemv, 0_nan_inf_incy_2)
}
Y[Ny - 1] = NAN;
BLASFUNC(dgemv)(&trans, &N, &N, &alpha, A, &N, X, &incX, &beta, Y, &incY);
for (int i = 0; i < Ny; i ++)
for (i = 0; i < Ny; i ++)
ASSERT_TRUE(Y[i] == 0.0);
}