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

1 files changed, 34 insertions, 11 deletions

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 <math.h>
-
 #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);
 }