diff --git a/AMSS_NCKU_source/bssn_rhs_c.C b/AMSS_NCKU_source/bssn_rhs_c.C index 4354866..3eaf764 100644 --- a/AMSS_NCKU_source/bssn_rhs_c.C +++ b/AMSS_NCKU_source/bssn_rhs_c.C @@ -283,9 +283,6 @@ int f_compute_rhs_bssn(int *ex, double &T, + ( gupxy[i]*gupyz[i] + gupyy[i]*gupxz[i] ) * Axy[i] + ( gupxy[i]*gupzz[i] + gupyz[i]*gupxz[i] ) * Axz[i] + ( gupyy[i]*gupzz[i] + gupyz[i]*gupyz[i] ) * Ayz[i]; - Rxx[i] = axx; Ryy[i] = ayy; Rzz[i] = azz; - Rxy[i] = axy; Rxz[i] = axz; Ryz[i] = ayz; - Gamx_rhs[i] = - TWO * ( Lapx[i]*axx + Lapy[i]*axy + Lapz[i]*axz ) + TWO * alpn1[i] * ( -F3o2/chin1[i] * ( chix[i]*axx + chiy[i]*axy + chiz[i]*axz ) - @@ -332,7 +329,6 @@ int f_compute_rhs_bssn(int *ex, double &T, double lfxx = gxxx[i] + gxyy[i] + gxzz[i]; double lfxy = gxyx[i] + gyyy[i] + gyzz[i]; double lfxz = gxzx[i] + gyzy[i] + gzzz[i]; - fxx[i] = lfxx; fxy[i] = lfxy; fxz[i] = lfxz; double gxa = gupxx[i]*Gamxxx[i] + gupyy[i]*Gamxyy[i] + gupzz[i]*Gamxzz[i] + TWO * ( gupxy[i]*Gamxxy[i] + gupxz[i]*Gamxxz[i] + gupyz[i]*Gamxyz[i] );