diff options
| author | tradke <tradke@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2001-08-29 14:39:44 +0000 |
|---|---|---|
| committer | tradke <tradke@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2001-08-29 14:39:44 +0000 |
| commit | 7e4e6fd2ee1f92984a7538558a7f94515bbf08b1 (patch) | |
| tree | 9531de7a047e22768fbbd09dde9feae885e0d4db | |
| parent | b989e639d042e8cf2948de347f721b67ef20181f (diff) | |
Fixed bug in Robin BC.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@154 6a38eb6e-646e-4a02-a296-d141613ad6c4
| -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 e4dd291..c80e006 100644 --- a/src/RobinBoundary.c +++ b/src/RobinBoundary.c @@ -454,10 +454,10 @@ void CCTK_FCALL CCTK_FNAME (BndRobinVN) u_src = sqrt (SQR (dx * x[src]) + SQR (dy * y[src]) + SQR (dz * z[src])); \ } \ \ - aux = decay * distance * (u_src / SQR (r[src]) + u_dst / SQR (r[dst])); \ + aux = decay * distance * (u_src + u_dst) / SQR (r[src] + r[dst]); \ \ data = (cctk_type *) GH->data[var][0]; \ - data[dst] = (cctk_type) ((2*aux* data[src] + finf*(1 - aux)) / (1 + aux)); \ + data[dst] = (cctk_type) ((2*aux*finf + data[src]*(1 - aux)) / (1 + aux)); \ } @@ -696,7 +696,7 @@ static int ApplyBndRobin (const cGH *GH, sGHex = (SymmetryGHex *) CCTK_GHExtension (GH, "Symmetry"); /* get the decay rate as a double */ - decay = 0.25 * (double) npow; + decay = (double) npow; /* precompute the distance to all 8 neighbors in a 3D grid */ dist[0] = 0; /* not used */ |