Boudaries clean, ready for 1,2d

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@88 6a38eb6e-646e-4a02-a296-d141613ad6c4

1 files changed, 31 insertions, 29 deletions

diff --git a/src/RadiationBoundary.c b/src/RadiationBoundary.c
index e48df15..bab18a7 100644
--- a/src/RadiationBoundary.c
+++ b/src/RadiationBoundary.c
@@ -74,20 +74,21 @@
 
 #define SQR(a) ((a)*(a))
 
-void RadiativeBCi(cGH *GH, 
-		  int *sw,  
-		  int *doBC, 
-		  int *lsh,
-		  CCTK_REAL *dxyz, 
-		  CCTK_REAL dt, 
-		  CCTK_REAL *var_n, 
-		  CCTK_REAL *var_p,
-		  CCTK_REAL *x, 
-		  CCTK_REAL *y, 
-		  CCTK_REAL *z, 
-		  CCTK_REAL *r,
-		  CCTK_REAL var0, 
-		  CCTK_REAL v0) 
+int BndApplyRadiative3Di(cGH *GH,
+			  int  gdim,
+			  int *sw,  
+			  int *doBC, 
+			  int *lssh,
+			  CCTK_REAL *dxyz, 
+			  CCTK_REAL dt, 
+			  CCTK_REAL *var_n, 
+			  CCTK_REAL *var_p,
+			  CCTK_REAL *x, 
+			  CCTK_REAL *y, 
+			  CCTK_REAL *z, 
+			  CCTK_REAL *r,
+			  CCTK_REAL var0, 
+			  CCTK_REAL v0) 
 {
   int i,j,k;
 
@@ -129,8 +130,8 @@ void RadiativeBCi(cGH *GH,
   /* Lower x-bound: x(2,:,:) --> x[xgp2] */
 
   if (doBC[0]==1) {
-    for (k=0;k<lsh[2];k++) {
-      for (j=0;j<lsh[1];j++) {
+    for (k=0;k<lssh[2];k++) {
+      for (j=0;j<lssh[1];j++) {
 	for (i=sw[0]-1;i>=0;i--) {
 	  
 	  xgp0 = CCTK_GFINDEX3D(GH,i  ,j,k);
@@ -151,9 +152,9 @@ void RadiativeBCi(cGH *GH,
   /* Upper x-bound: x(nx,:,:) --> xgp[xgp0] */
 
   if (doBC[1]==1){
-    for (k=0;k<lsh[2];k++) {
-      for (j=0;j<lsh[1];j++) {
-	for (i=lsh[0]-sw[0];i<lsh[0];i++) {
+    for (k=0;k<lssh[2];k++) {
+      for (j=0;j<lssh[1];j++) {
+	for (i=lssh[0]-sw[0];i<lssh[0];i++) {
 	  
 	  xgp0 = CCTK_GFINDEX3D(GH,i  ,j,k);
 	  xgp1 = CCTK_GFINDEX3D(GH,i-1,j,k);
@@ -173,8 +174,8 @@ void RadiativeBCi(cGH *GH,
   /* Lower y-bound */
 
   if (doBC[2] == 1) {
-    for (k=0;k<lsh[2];k++) {
-      for (i=0;i<lsh[0];i++) {
+    for (k=0;k<lssh[2];k++) {
+      for (i=0;i<lssh[0];i++) {
 	for (j=sw[1]-1;j>=0;j--) {
 	  
 	  ygp0 = CCTK_GFINDEX3D(GH,i,j  ,k);
@@ -195,9 +196,9 @@ void RadiativeBCi(cGH *GH,
   /* Upper y bound */
 
   if (doBC[3] == 1) {
-    for (k=0;k<lsh[2];k++) {
-      for (i=0;i<lsh[0];i++) {
-	for (j=lsh[1]-sw[1];j<lsh[1];j++) {
+    for (k=0;k<lssh[2];k++) {
+      for (i=0;i<lssh[0];i++) {
+	for (j=lssh[1]-sw[1];j<lssh[1];j++) {
 	  
 	  ygp0 = CCTK_GFINDEX3D(GH,i,j  ,k);
 	  ygp1 = CCTK_GFINDEX3D(GH,i,j-1,k);
@@ -217,8 +218,8 @@ void RadiativeBCi(cGH *GH,
   /* Lower z-bound */
 
   if (doBC[4]==1) {
-    for (j=0;j<lsh[1];j++) {
-      for (i=0;i<lsh[0];i++) {
+    for (j=0;j<lssh[1];j++) {
+      for (i=0;i<lssh[0];i++) {
 	for (k=sw[2]-1;k>=0;k--) {
 	  
 	  zgp0 = CCTK_GFINDEX3D(GH,i,j,k  );
@@ -239,9 +240,9 @@ void RadiativeBCi(cGH *GH,
   /* Upper z-bound */
 
   if (doBC[5] == 1) {
-    for (j=0;j<lsh[1];j++) {
-      for (i=0;i<lsh[0];i++) {
-	for (k=lsh[2]-sw[2];k<lsh[2];k++) {
+    for (j=0;j<lssh[1];j++) {
+      for (i=0;i<lssh[0];i++) {
+	for (k=lssh[2]-sw[2];k<lssh[2];k++) {
 	  
 	  zgp0 = CCTK_GFINDEX3D(GH,i,j,k  );
 	  zgp1 = CCTK_GFINDEX3D(GH,i,j,k-1);
@@ -257,4 +258,5 @@ void RadiativeBCi(cGH *GH,
       }
     }
   }
+  return(0);
 }