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

fix sbgemm bfloat16 conversion errors introduced in PR 4488

This commit is contained in:
Martin Kroeker 2024-02-07 19:57:18 +01:00 committed by GitHub
parent 22b487b622
commit e9f480111e
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 6 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
test/compare_sgemm_sbgemm.c
View File

@ -81,16 +81,6 @@ float16to32 (bfloat16_bits f16)
return f32.v; return f32.v;
} }
float
float32to16 (float32_bits f32)
{
bfloat16_bits f16;
f16.bits.s = f32.bits.s;
f16.bits.e = f32.bits.e;
f16.bits.m = (uint32_t) f32.bits.m >> 16;
return f32.v;
}
int int
main (int argc, char *argv[]) main (int argc, char *argv[])
{ {
@ -110,6 +100,8 @@ main (int argc, char *argv[])
float C[m * n]; float C[m * n];
bfloat16_bits AA[m * k], BB[k * n]; bfloat16_bits AA[m * k], BB[k * n];
float DD[m * n], CC[m * n]; float DD[m * n], CC[m * n];
bfloat16 atmp,btmp;
blasint one=1;
for (j = 0; j < m; j++) for (j = 0; j < m; j++)
{ {
@ -118,8 +110,10 @@ main (int argc, char *argv[])
A[j * k + i] = ((FLOAT) rand () / (FLOAT) RAND_MAX) + 0.5; A[j * k + i] = ((FLOAT) rand () / (FLOAT) RAND_MAX) + 0.5;
B[j * k + i] = ((FLOAT) rand () / (FLOAT) RAND_MAX) + 0.5; B[j * k + i] = ((FLOAT) rand () / (FLOAT) RAND_MAX) + 0.5;
C[j * k + i] = 0; C[j * k + i] = 0;
AA[j * k + i].v = float32to16( A[j * k + i] ); sbstobf16_(&one, &A[j*k+i], &one, &atmp, &one);
BB[j * k + i].v = float32to16( B[j * k + i] ); sbstobf16_(&one, &B[j*k+i], &one, &btmp, &one);
AA[j * k + i].v = atmp;
BB[j * k + i].v = btmp;
CC[j * k + i] = 0; CC[j * k + i] = 0;
DD[j * k + i] = 0; DD[j * k + i] = 0;
} }