From f4f0f44703585c0b8fad7cb506a49a168dc0e2c1 Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Tue, 17 Jun 2025 08:36:50 +0200
Subject: [PATCH] Update zscal.c

---
 kernel/arm/zscal.c | 43 ++++++++++++++++++++++++++++++-------------
 1 file changed, 30 insertions(+), 13 deletions(-)

diff --git a/kernel/arm/zscal.c b/kernel/arm/zscal.c
index d00e748f1..51509ecda 100644
--- a/kernel/arm/zscal.c
+++ b/kernel/arm/zscal.c
@@ -51,22 +51,39 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F
 
     inc_x2 = 2 * inc_x;
     if (dummy2 == 0) {
-        for (i = 0; i < n; i++)
+        for ( i=0; i<n; i++ )
         {
-            if (da_r == 0.0 && da_i == 0.0)
-            {
-                x[ip] = 0.0;
-                x[ip+1] = 0.0;
-            }
-            else
-            {
-                temp    = da_r * x[ip]   - da_i * x[ip+1];
-                x[ip+1] = da_r * x[ip+1] + da_i * x[ip]  ;
-                x[ip] = temp;
-            }
+                if ( da_r == 0.0 )
+                {
+                        if ( da_i == 0.0 )
+                        {
+                                temp = 0.0;
+                                x[ip+1] = 0.0 ;
+                        }
+                        else
+                        {
+                                temp = - da_i * x[ip+1] ;
+                                x[ip+1] = da_i * x[ip]  ;
+                        }
+                }
+                else
+                {
+                        if ( da_i == 0.0 )
+                        {
+                                temp    = da_r * x[ip]  ;
+                                x[ip+1] = da_r * x[ip+1];
+                        }
+                        else
+                        {
+                                temp    = da_r * x[ip]   - da_i * x[ip+1] ;
+                                x[ip+1] = da_r * x[ip+1] + da_i * x[ip]   ;
+                        }
+                }
+                x[ip]   = temp;
 
-            ip += inc_x2;
+                ip += inc_x2;
         }
+
         return(0);
     }
     for (i = 0; i < n; i++)