From 51dead090e5df59aa97d7ce42dde94f1c1bf7407 Mon Sep 17 00:00:00 2001
From: CGH0S7 <776459475@qq.com>
Date: Sat, 28 Feb 2026 13:49:45 +0800
Subject: [PATCH] =?UTF-8?q?bssn=5Frhs=5Fc:=20=E8=9E=8D=E5=90=88=E6=9C=80?=
 =?UTF-8?q?=E7=BB=88RHS=E4=B8=A4=E5=BE=AA=E7=8E=AF=E4=B8=BA=E4=B8=80?=
 =?UTF-8?q?=E5=BE=AA=E7=8E=AF=EF=BC=8C=E7=94=A8=E5=B1=80=E9=83=A8=E5=8F=98?=
 =?UTF-8?q?=E9=87=8F=E4=BC=A0=E9=80=92fij=E4=B8=AD=E9=97=B4=E5=80=BC=20(Mo?=
 =?UTF-8?q?dify=206)?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
---
 AMSS_NCKU_source/bssn_rhs_c.C | 27 ++++++++++++---------------
 1 file changed, 12 insertions(+), 15 deletions(-)

diff --git a/AMSS_NCKU_source/bssn_rhs_c.C b/AMSS_NCKU_source/bssn_rhs_c.C
index ea2351d..ea88692 100644
--- a/AMSS_NCKU_source/bssn_rhs_c.C
+++ b/AMSS_NCKU_source/bssn_rhs_c.C
@@ -812,23 +812,20 @@ int f_compute_rhs_bssn(int *ex, double &T,
                 + (alpn1[i] / chin1[i]) * f[i]
             );
 
-            fxx[i] = alpn1[i] * (Rxx[i] - EIGHT * PI * Sxx[i]) - fxx[i];
-            fxy[i] = alpn1[i] * (Rxy[i] - EIGHT * PI * Sxy[i]) - fxy[i];
-            fxz[i] = alpn1[i] * (Rxz[i] - EIGHT * PI * Sxz[i]) - fxz[i];
-            fyy[i] = alpn1[i] * (Ryy[i] - EIGHT * PI * Syy[i]) - fyy[i];
-            fyz[i] = alpn1[i] * (Ryz[i] - EIGHT * PI * Syz[i]) - fyz[i];
-            fzz[i] = alpn1[i] * (Rzz[i] - EIGHT * PI * Szz[i]) - fzz[i];
-        }
-        // 8ms //
-        for (int i=0;i<all;i+=1) {
+            double l_fxx = alpn1[i] * (Rxx[i] - EIGHT * PI * Sxx[i]) - fxx[i];
+            double l_fxy = alpn1[i] * (Rxy[i] - EIGHT * PI * Sxy[i]) - fxy[i];
+            double l_fxz = alpn1[i] * (Rxz[i] - EIGHT * PI * Sxz[i]) - fxz[i];
+            double l_fyy = alpn1[i] * (Ryy[i] - EIGHT * PI * Syy[i]) - fyy[i];
+            double l_fyz = alpn1[i] * (Ryz[i] - EIGHT * PI * Syz[i]) - fyz[i];
+            double l_fzz = alpn1[i] * (Rzz[i] - EIGHT * PI * Szz[i]) - fzz[i];
 
             /* Aij_rhs = fij - gij * f */
-            Axx_rhs[i] = fxx[i] - gxx[i] * f[i];
-            Ayy_rhs[i] = fyy[i] - gyy[i] * f[i];
-            Azz_rhs[i] = fzz[i] - gzz[i] * f[i];
-            Axy_rhs[i] = fxy[i] - gxy[i] * f[i];
-            Axz_rhs[i] = fxz[i] - gxz[i] * f[i];
-            Ayz_rhs[i] = fyz[i] - gyz[i] * f[i];
+            Axx_rhs[i] = l_fxx - gxx[i] * f[i];
+            Ayy_rhs[i] = l_fyy - gyy[i] * f[i];
+            Azz_rhs[i] = l_fzz - gzz[i] * f[i];
+            Axy_rhs[i] = l_fxy - gxy[i] * f[i];
+            Axz_rhs[i] = l_fxz - gxz[i] * f[i];
+            Ayz_rhs[i] = l_fyz - gyz[i] * f[i];
 
             /* Now: store A_il A^l_j into fij: */
             fxx[i] =