64-BitBrainstorm_2026/AMSS-NCKU
8
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bssn_rhs_c: 融合最终RHS两循环为一循环，用局部变量传递fij中间值 (Modify 6)

Browse Source 
Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
This commit is contained in:
CGH0S7
2026-02-28 13:49:45 +08:00
parent 34d6922a66
commit 51dead090e
1 changed files with 12 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
AMSS_NCKU_source/bssn_rhs_c.C
View File

@@ -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] =