workaround fault with ssq=inf,scale=0

diff --git a/kernel/arm64/dznrm2_thunderx2t99.c b/kernel/arm64/dznrm2_thunderx2t99.c
index fba2fe8..e342b0b 100644 (file)
--- a/kernel/arm64/dznrm2_thunderx2t99.c
+++ b/kernel/arm64/dznrm2_thunderx2t99.c
@@ -404,6 +404,7 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x)
 #else
        nrm2_compute(n, x, inc_x, &ssq, &scale);
 #endif
+       if (fabs(scale) <1.e-300) return 0.;
        ssq = sqrt(ssq) * scale;
 
        return ssq;