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

Fix division by zero in the non-x86 codepath

This commit is contained in:
Martin Kroeker 2021-04-29 09:47:18 +02:00 committed by GitHub
parent 352efdd13a
commit c5fb91f1bc
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 12 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
interface/zrotg.c
View File

@ -79,8 +79,12 @@ void NAME(FLOAT *DA, FLOAT *DB, FLOAT *C, FLOAT *S){
aa_i = fabs(da_r); aa_i = fabs(da_r);
} }
scale = (aa_i / aa_r); if (aa_r == ZERO) {
ada = aa_r * sqrt(ONE + scale * scale); ada = 0.;
} else {
scale = (aa_i / aa_r);
ada = aa_r * sqrt(ONE + scale * scale);
}
bb_r = fabs(db_r); bb_r = fabs(db_r);
bb_i = fabs(db_i); bb_i = fabs(db_i);
@ -90,9 +94,12 @@ void NAME(FLOAT *DA, FLOAT *DB, FLOAT *C, FLOAT *S){
bb_i = fabs(bb_r); bb_i = fabs(bb_r);
} }
scale = (bb_i / bb_r); if (bb_r == ZERO) {
adb = bb_r * sqrt(ONE + scale * scale); adb = 0.;
} else {
scale = (bb_i / bb_r);
adb = bb_r * sqrt(ONE + scale * scale);
}
scale = ada + adb; scale = ada + adb;
aa_r = da_r / scale; aa_r = da_r / scale;