diff options
-rw-r--r-- | src/RobinBoundary.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/RobinBoundary.c b/src/RobinBoundary.c index de601a2..1ce8cfa 100644 --- a/src/RobinBoundary.c +++ b/src/RobinBoundary.c @@ -75,7 +75,7 @@ int BndApplyRobin3Di(cGH *GH, xgp4 = CCTK_GFINDEX3D(GH,lsh[0]-4,j,k); var[xgp1] = finf + ( - pow(r[xgp3],npow)*( var[xgp4] - finf ) + pow(r[xgp4],npow)*( var[xgp4] - finf ) -3.0*pow(r[xgp3],npow)*( var[xgp3] - finf ) +3.0*pow(r[xgp2],npow)*( var[xgp2] - finf ) )/pow(r[xgp1],npow); @@ -112,7 +112,7 @@ int BndApplyRobin3Di(cGH *GH, ygp4 = CCTK_GFINDEX3D(GH,i,lsh[1]-4,k); var[ygp1] = finf + ( - pow(r[ygp3],npow)*( var[ygp4] - finf ) + pow(r[ygp4],npow)*( var[ygp4] - finf ) -3.0*pow(r[ygp3],npow)*( var[ygp3] - finf ) +3.0*pow(r[ygp2],npow)*( var[ygp2] - finf ) )/pow(r[ygp1],npow); @@ -148,7 +148,7 @@ int BndApplyRobin3Di(cGH *GH, zgp4 = CCTK_GFINDEX3D(GH,i,j,lsh[2]-4); var[zgp1] = finf + ( - pow(r[zgp3],npow)*( var[zgp4] - finf ) + pow(r[zgp4],npow)*( var[zgp4] - finf ) -3.0*pow(r[zgp3],npow)*( var[zgp3] - finf ) +3.0*pow(r[zgp2],npow)*( var[zgp2] - finf ) )/pow(r[zgp1],npow); |