From 37e3c5ff9dff48d647bd278ea28b64c5cede1c9f Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Wed, 31 Jul 2024 14:26:00 +0200 Subject: [PATCH] Update zscal.c --- kernel/zarch/zscal.c | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/kernel/zarch/zscal.c b/kernel/zarch/zscal.c index 9fcf47cc4..63f36e972 100644 --- a/kernel/zarch/zscal.c +++ b/kernel/zarch/zscal.c @@ -381,7 +381,8 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, x[i + 1] = da_i * x[i]; else x[i + 1] = NAN; - x[i] = temp0; + if (isfinite(x[i])) + x[i] = temp0; i += 2; j++; @@ -396,6 +397,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, while (j < n) { temp0 = da_r * x[i]; + if (isinf(x[i]) x[i + 1] = da_r * x[i + 1]; x[i] = temp0; i += 2;