diff options
| author | tradke <tradke@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2001-04-14 17:13:33 +0000 |
|---|---|---|
| committer | tradke <tradke@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2001-04-14 17:13:33 +0000 |
| commit | 2cbf9679f3326c80e65c0a71fdaa764b10af6d4f (patch) | |
| tree | ef4aecd5c8033febc07760c8af6d1df766e81f35 /src/RobinBoundary.c | |
| parent | 4c3b2a577a2935d87c55548c1b1a02d569e16b25 (diff) | |
Fixed wrong indexing for upper boundaries.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@135 6a38eb6e-646e-4a02-a296-d141613ad6c4
Diffstat (limited to 'src/RobinBoundary.c')
| -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); |