From 779e0b79c1af4f3a1e1e66164dc3d539fd1869bc Mon Sep 17 00:00:00 2001 From: lanfer Date: Tue, 18 Jul 2000 16:56:01 +0000 Subject: Boundary routines: reworked Robin to have conformal calling interfaces fixed two bugs: the more interresting one: coordinates do not have timelevels. Replaced most CCTK_WARNS with CCTK_VWarns to be more informative. gerd git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@115 6a38eb6e-646e-4a02-a296-d141613ad6c4 --- src/RobinBoundary.c | 45 ++++++++++++++++++++++++++++++++++----------- 1 file changed, 34 insertions(+), 11 deletions(-) (limited to 'src/RobinBoundary.c') diff --git a/src/RobinBoundary.c b/src/RobinBoundary.c index b6f3267..bcecc4d 100644 --- a/src/RobinBoundary.c +++ b/src/RobinBoundary.c @@ -1,18 +1,40 @@ + /*@@ + @file RobinBoundary.c + @date Tue Jul 18 18:07:24 2000 + @author Gerd Lanfermann + @desc + internal routines to apply 1d,2d,3d Robin + boundary conditions + @enddesc + @@*/ #include - #include "cctk.h" #include "cctk_Parameters.h" -void BndApplyRobin3Di(cGH *GH, - int gdim, - int *sw, - int *doBC, - int *lsh, - CCTK_REAL *var, - CCTK_REAL *r, - CCTK_REAL finf, - CCTK_INT npow) +/*@@ + @routine BndApplyRobin3Di + @date Tue Jul 18 18:06:53 2000 + @author Gerd Lanfermann + @desc + Apply 3d Robin Boundary condition - internal routine + @enddesc + @calls + @calledby + @history + + @endhistory + +@@*/ + +int BndApplyRobin3Di(cGH *GH, + int *sw, + int *doBC, + int *lsh, + CCTK_REAL *var, + CCTK_REAL *r, + CCTK_REAL finf, + CCTK_INT npow) { int xgp1,xgp2,xgp3,xgp4; @@ -53,7 +75,7 @@ void BndApplyRobin3Di(cGH *GH, xgp4 = CCTK_GFINDEX3D(GH,lsh[0]-4,j,k); var[xgp1] = finf + ( - pow(r[xgp3],npow)*( var[xgp4] - finf ) + pow(r[xgp3],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); @@ -134,4 +156,5 @@ void BndApplyRobin3Di(cGH *GH, } } } + return(0); } -- cgit v1.2.3