Merge pull request #3675 from martin-frbg/issue3654
workaround ThunderX2 DNRM2 fault with ssq=inf,scale=0
This commit is contained in:
commit
bf8998a9f4
|
@ -404,6 +404,7 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x)
|
||||||
#else
|
#else
|
||||||
nrm2_compute(n, x, inc_x, &ssq, &scale);
|
nrm2_compute(n, x, inc_x, &ssq, &scale);
|
||||||
#endif
|
#endif
|
||||||
|
if (fabs(scale) <1.e-300) return 0.;
|
||||||
ssq = sqrt(ssq) * scale;
|
ssq = sqrt(ssq) * scale;
|
||||||
|
|
||||||
return ssq;
|
return ssq;
|
||||||
|
|
Loading…
Reference in New Issue