convert tabs to (8) spaces. No other changes.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@485 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 62 insertions, 62 deletions

diff --git a/src/GRHydro_Con2PrimM_polytype_pt.c b/src/GRHydro_Con2PrimM_polytype_pt.c
index 97d7648..4931415 100644
--- a/src/GRHydro_Con2PrimM_polytype_pt.c
+++ b/src/GRHydro_Con2PrimM_polytype_pt.c
@@ -60,9 +60,9 @@
 #define NEWT_TOL2       (1.0e-15)  /* TOL of new 1D^*_{v^2} gnr2 method */
 #define MIN_NEWT_TOL2   (1.0e-10)  /* TOL of new 1D^*_{v^2} gnr2 method */
 
-#define W_TOO_BIG	(1.e20)    /* \gamma^2 (\rho_0 + u + p) is assumed
+#define W_TOO_BIG        (1.e20)    /* \gamma^2 (\rho_0 + u + p) is assumed
                                       to always be smaller than this.  This
-				      is used to detect solver failures */
+                                      is used to detect solver failures */
 
 #define FAIL_VAL        (1.e30)    /* Generic value to which we set variables when a problem arises */
 
@@ -80,23 +80,23 @@ struct LocGlob {
 
 // Declarations: 
 static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocGlob *lgp,
-				     void (*funcd) (CCTK_REAL [], CCTK_REAL [], 
-						    CCTK_REAL [], CCTK_REAL [][1], 
-						    CCTK_REAL *, CCTK_REAL *, CCTK_INT, CCTK_REAL, 
-						    struct LocGlob *lgp) );
+                                     void (*funcd) (CCTK_REAL [], CCTK_REAL [], 
+                                                    CCTK_REAL [], CCTK_REAL [][1], 
+                                                    CCTK_REAL *, CCTK_REAL *, CCTK_INT, CCTK_REAL, 
+                                                    struct LocGlob *lgp) );
 
 static  void func_W( CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], CCTK_REAL [][1], 
-		     CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
-		     struct LocGlob * lgp);
+                     CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
+                     struct LocGlob * lgp);
 
 static CCTK_INT gnr2( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocGlob *lgp,
-			    void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
-					   CCTK_REAL [][1],CCTK_REAL *,CCTK_REAL *,CCTK_INT, CCTK_REAL,
-					   struct LocGlob *lgp) );
+                            void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
+                                           CCTK_REAL [][1],CCTK_REAL *,CCTK_REAL *,CCTK_INT, CCTK_REAL,
+                                           struct LocGlob *lgp) );
 
 static void func_rho(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[], 
-		     CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
-		     struct LocGlob *lgp);
+                     CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
+                     struct LocGlob *lgp);
 
 void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt) (
     CCTK_INT  *handle, CCTK_REAL *gamma_eos, 
@@ -171,11 +171,11 @@ RETURN VALUE: of retval = (i*100 + j)  where
                    given tolerances;
              2 ->  Newton-Raphson procedure encountered a numerical divergence 
                    (occurrence of "nan" or "+/-inf" ;
-	     
+             
          j = 0 -> success 
              1 -> failure: some sort of failure in Newton-Raphson; 
              2 -> failure: unphysical vsq = v^2  value at initial guess;
-	     3 -> failure: W<0 or W>W_TOO_BIG
+             3 -> failure: W<0 or W>W_TOO_BIG
              4 -> failure: v^2 > 1 
              ( used to be  5 -> failure: rho,uu <= 0   but now sets epsnegative to non-zero )
 
@@ -226,37 +226,37 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt)  (
 
 #if(DEBUG_CON2PRIMM)
   fprintf(stdout," *dens        = %26.16e \n", *dens_in         );
-  fprintf(stdout," *sx 	        = %26.16e \n", *sx_in 	  );    
-  fprintf(stdout," *sy  	= %26.16e \n", *sy_in  	  );    
-  fprintf(stdout," *sz 	        = %26.16e \n", *sz_in 	  );    
-  fprintf(stdout," *Sc 	        = %26.16e \n", *sc_in 	  );    
-  fprintf(stdout," *Bconsx 	= %26.16e \n", *Bconsx_in );    
-  fprintf(stdout," *Bconsy 	= %26.16e \n", *Bconsy_in );    
-  fprintf(stdout," *Bconsz 	= %26.16e \n", *Bconsz_in );    
-  fprintf(stdout," *rho         = %26.16e \n", *rho 	  );    
-  fprintf(stdout," *velx        = %26.16e \n", *velx 	  );    
-  fprintf(stdout," *vely        = %26.16e \n", *vely	  );    
-  fprintf(stdout," *velz        = %26.16e \n", *velz	  );    
+  fprintf(stdout," *sx                 = %26.16e \n", *sx_in           );    
+  fprintf(stdout," *sy          = %26.16e \n", *sy_in            );    
+  fprintf(stdout," *sz                 = %26.16e \n", *sz_in           );    
+  fprintf(stdout," *Sc                 = %26.16e \n", *sc_in           );    
+  fprintf(stdout," *Bconsx         = %26.16e \n", *Bconsx_in );    
+  fprintf(stdout," *Bconsy         = %26.16e \n", *Bconsy_in );    
+  fprintf(stdout," *Bconsz         = %26.16e \n", *Bconsz_in );    
+  fprintf(stdout," *rho         = %26.16e \n", *rho           );    
+  fprintf(stdout," *velx        = %26.16e \n", *velx           );    
+  fprintf(stdout," *vely        = %26.16e \n", *vely          );    
+  fprintf(stdout," *velz        = %26.16e \n", *velz          );    
   fprintf(stdout," *epsilon     = %26.16e \n", *epsilon     );
   fprintf(stdout," *pressure    = %26.16e \n", *pressure    );
-  fprintf(stdout," *Bx 	        = %26.16e \n", *Bx 	  );    
-  fprintf(stdout," *By 	        = %26.16e \n", *By 	  );    
-  fprintf(stdout," *Bz 	        = %26.16e \n", *Bz 	  );    
-  fprintf(stdout," *bsq	        = %26.16e \n", *bsq	  );    
+  fprintf(stdout," *Bx                 = %26.16e \n", *Bx           );    
+  fprintf(stdout," *By                 = %26.16e \n", *By           );    
+  fprintf(stdout," *Bz                 = %26.16e \n", *Bz           );    
+  fprintf(stdout," *bsq                = %26.16e \n", *bsq          );    
   fprintf(stdout," *w_lorentz   = %26.16e \n", *w_lorentz   );
-  fprintf(stdout," *gxx 	= %26.16e \n", *gxx 	  );    
-  fprintf(stdout," *gxy 	= %26.16e \n", *gxy 	  );    
-  fprintf(stdout," *gxz 	= %26.16e \n", *gxz 	  );    
-  fprintf(stdout," *gyy 	= %26.16e \n", *gyy 	  );    
-  fprintf(stdout," *gyz 	= %26.16e \n", *gyz 	  );    
-  fprintf(stdout," *gzz 	= %26.16e \n", *gzz 	  );    
-  fprintf(stdout," *uxx 	= %26.16e \n", *uxx 	  );    
-  fprintf(stdout," *uxy 	= %26.16e \n", *uxy 	  );    
-  fprintf(stdout," *uxz	        = %26.16e \n", *uxz	  );    
-  fprintf(stdout," *uyy 	= %26.16e \n", *uyy 	  );    
-  fprintf(stdout," *uyz 	= %26.16e \n", *uyz 	  );    
-  fprintf(stdout," *uzz	        = %26.16e \n", *uzz	  );    
-  fprintf(stdout," *det	        = %26.16e \n", *det	  );    
+  fprintf(stdout," *gxx         = %26.16e \n", *gxx           );    
+  fprintf(stdout," *gxy         = %26.16e \n", *gxy           );    
+  fprintf(stdout," *gxz         = %26.16e \n", *gxz           );    
+  fprintf(stdout," *gyy         = %26.16e \n", *gyy           );    
+  fprintf(stdout," *gyz         = %26.16e \n", *gyz           );    
+  fprintf(stdout," *gzz         = %26.16e \n", *gzz           );    
+  fprintf(stdout," *uxx         = %26.16e \n", *uxx           );    
+  fprintf(stdout," *uxy         = %26.16e \n", *uxy           );    
+  fprintf(stdout," *uxz                = %26.16e \n", *uxz          );    
+  fprintf(stdout," *uyy         = %26.16e \n", *uyy           );    
+  fprintf(stdout," *uyz         = %26.16e \n", *uyz           );    
+  fprintf(stdout," *uzz                = %26.16e \n", *uzz          );    
+  fprintf(stdout," *det                = %26.16e \n", *det          );    
   fprintf(stdout," *epsnegative = %10d    \n", *epsnegative );
   fprintf(stdout," *retval      = %26.16e \n", *retval      );
   fflush(stdout);
@@ -326,7 +326,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt)  (
 
   gammasq = 1. / (1. - vsq);
   gamma  = sqrt(gammasq);
-	
+        
   // Always calculate rho from D and gamma so that using D in EOS remains consistent
   //   i.e. you don't get positive values for dP/d(vsq) . 
   rho0 = lg.D / gamma ;
@@ -342,8 +342,8 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt)  (
   // Make sure that W is large enough so that v^2 < 1 : 
   i_increase = 0;
   while( (( W_last*W_last*W_last * ( W_last + 2.*lg.Bsq ) 
-	    - lg.QdotBsq*(2.*W_last + lg.Bsq) ) <= W_last*W_last*(lg.Qtsq-lg.Bsq*lg.Bsq))
-	 && (i_increase < 10) ) {
+            - lg.QdotBsq*(2.*W_last + lg.Bsq) ) <= W_last*W_last*(lg.Qtsq-lg.Bsq*lg.Bsq))
+         && (i_increase < 10) ) {
     W_last *= 10.;
     i_increase++;
   }
@@ -357,7 +357,7 @@ void CCTK_FCALL CCTK_FNAME(GRHydro_Con2PrimM_Polytype_pt)  (
   *retval = 1.0*oned_newton_raphson( x_1d, 1, gammaeos, &lg, func_W ) ;  
 
   W = x_1d[0];
-	
+        
   /* Problem with solver, so return denoting error before doing anything further */
   if( ((*retval) != 0.) || (W == FAIL_VAL) ) {
     *retval = *retval*100.+1.;
@@ -487,9 +487,9 @@ static void validate_x(CCTK_REAL x[2], CCTK_REAL x0[2] )
 
 *****************************************************************/
 static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocGlob *lgp,
-				   void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
-						  CCTK_REAL [][1], CCTK_REAL *, 
-						  CCTK_REAL *, CCTK_INT, CCTK_REAL, struct LocGlob *) )
+                                   void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
+                                                  CCTK_REAL [][1], CCTK_REAL *, 
+                                                  CCTK_REAL *, CCTK_INT, CCTK_REAL, struct LocGlob *) )
 {
   CCTK_REAL f, df, dx[1], x_old[1], resid[1], 
     jac[1][1];
@@ -528,8 +528,8 @@ static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammae
 //    //METHOD specific:
 //    i_increase = 0;
 //    while( (( x[0]*x[0]*x[0] * ( x[0] + 2.*lgp->Bsq ) - 
-//	      lgp->QdotBsq*(2.*x[0] + lgp->Bsq) ) <= x[0]*x[0]*(lgp->Qtsq-lgp->Bsq*lgp->Bsq))
-//	   && (i_increase < 10) ) {
+//              lgp->QdotBsq*(2.*x[0] + lgp->Bsq) ) <= x[0]*x[0]*(lgp->Qtsq-lgp->Bsq*lgp->Bsq))
+//           && (i_increase < 10) ) {
 //      x[0] -= (1.*i_increase) * dx[0] / 10. ;
 //      i_increase++;
 //    }
@@ -561,7 +561,7 @@ static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammae
     if( doing_extra == 1 ) i_extra++ ;
 
     if( ((fabs(errx) <= NEWT_TOL)&&(doing_extra == 0)) || 
-	(i_extra > EXTRA_NEWT_ITER) || (n_iter >= (MAX_NEWT_ITER-1)) ) {
+        (i_extra > EXTRA_NEWT_ITER) || (n_iter >= (MAX_NEWT_ITER-1)) ) {
       keep_iterating = 0;
     }
 
@@ -574,7 +574,7 @@ static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammae
   if( (!finite(f)) || (!finite(df)) || (!finite(x[0]))  ) {
 #if(DEBUG_CON2PRIMM)
     fprintf(stderr,"\ngnr not finite, f,df,x_o,x,W_o,W,rho_o,rho = %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e \n",
-    	    f,df,x[0],x_old[0],lgp->W_for_gnr2_old,lgp->W_for_gnr2,lgp->rho_for_gnr2_old,lgp->rho_for_gnr2); fflush(stderr); 
+                f,df,x[0],x_old[0],lgp->W_for_gnr2_old,lgp->W_for_gnr2,lgp->rho_for_gnr2_old,lgp->rho_for_gnr2); fflush(stderr); 
 #endif
     return(2);
   }
@@ -611,8 +611,8 @@ static CCTK_INT oned_newton_raphson( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammae
          n    = dimension of x[];
  *********************************************************************************/
 static void func_W(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[], 
-		   CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
-		   struct LocGlob *lgp)
+                   CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
+                   struct LocGlob *lgp)
 {
   CCTK_INT retval, ntries;
   CCTK_REAL  W, x_rho[1], rho, rho_g ;
@@ -689,8 +689,8 @@ static void func_W(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[],
 
 *****************************************************************/
 static CCTK_INT gnr2( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocGlob *lgp,
-			    void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
-					   CCTK_REAL [][1],CCTK_REAL *,CCTK_REAL *,CCTK_INT,CCTK_REAL, struct LocGlob *lgp) )
+                            void (*funcd) (CCTK_REAL [], CCTK_REAL [], CCTK_REAL [], 
+                                           CCTK_REAL [][1],CCTK_REAL *,CCTK_REAL *,CCTK_INT,CCTK_REAL, struct LocGlob *lgp) )
 {
   CCTK_REAL f, df, dx[1], x_old[1], resid[1], 
     jac[1][1];
@@ -746,7 +746,7 @@ static CCTK_INT gnr2( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocG
 
     // See if we've done the extra iterations, or have done too many iterations:
     if( ((fabs(errx) <= NEWT_TOL2)&&(doing_extra == 0)) || 
-	(i_extra > EXTRA_NEWT_ITER) || (n_iter >= (MAX_NEWT_ITER-1)) ) {
+        (i_extra > EXTRA_NEWT_ITER) || (n_iter >= (MAX_NEWT_ITER-1)) ) {
       keep_iterating = 0;
     }
 
@@ -759,7 +759,7 @@ static CCTK_INT gnr2( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocG
   if( (!finite(f)) || (!finite(df)) || (!finite(x[0]))  ) {
 #if(DEBUG_CON2PRIMM)
     fprintf(stderr,"\ngnr2 not finite, f,df,x_o,x,W_o,W,rho_o,rho = %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e %26.20e \n",
-    	    f,df,x[0],x_old[0],lgp->W_for_gnr2_old,lgp->W_for_gnr2,lgp->rho_for_gnr2_old,lgp->rho_for_gnr2); fflush(stderr); 
+                f,df,x[0],x_old[0],lgp->W_for_gnr2_old,lgp->W_for_gnr2,lgp->rho_for_gnr2_old,lgp->rho_for_gnr2); fflush(stderr); 
 #endif
     return(2);
   }
@@ -797,8 +797,8 @@ static CCTK_INT gnr2( CCTK_REAL x[], CCTK_INT n, CCTK_REAL gammaeos, struct LocG
  *********************************************************************************/
 // for the isentropic version:   eq.  (27)
 static void func_rho(CCTK_REAL x[], CCTK_REAL dx[], CCTK_REAL resid[], 
-		     CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
-		     struct LocGlob *lgp)
+                     CCTK_REAL jac[][1], CCTK_REAL *f, CCTK_REAL *df, CCTK_INT n, CCTK_REAL gammaeos, 
+                     struct LocGlob *lgp)
 {
 
   CCTK_REAL rho, W;