[cig-commits] commit: Remove some testing stuff

Mercurial hg at geodynamics.org
Wed Feb 29 19:20:13 PST 2012


changeset:   64:57ee741e9644
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Feb 29 15:31:48 2012 -0800
files:       compute_Cx.cxx compute_Cy.cxx main.cxx
description:
Remove some testing stuff


diff -r bd123bf91338 -r 57ee741e9644 compute_Cx.cxx
--- a/compute_Cx.cxx	Wed Feb 29 15:30:17 2012 -0800
+++ b/compute_Cx.cxx	Wed Feb 29 15:31:48 2012 -0800
@@ -67,30 +67,7 @@ void compute_Cx(const Model &model, cons
 
               dCx_dvx+= -(eta_jump*dvy_y_dv + dx*dzy_yx_dv)/h
                 + (dp_dv + dx*dp_x_dv)/h;
-
-              if(j==Ny/8)
-                std::cout << "Cx "
-                          << i << " "
-                          << j << " "
-                          << x << " "
-                          << middle << " "
-                          << zx[i][j] << " "
-                          << Cx << " "
-                          << dCx_dvx/eta << " "
-                          << "\n";
-
             }
-          if(j==Ny/8)
-            std::cout << "Cx out "
-                      << i << " "
-                      << j << " "
-                      << x << " "
-                      << middle << " "
-                      << zx[i][j] << " "
-                      << Cx << " "
-                      << dCx_dvx/eta << " "
-                      << "\n";
-
           dCx_dzx=dCx_dvx/eta;
         }
     }
diff -r bd123bf91338 -r 57ee741e9644 compute_Cy.cxx
--- a/compute_Cy.cxx	Wed Feb 29 15:30:17 2012 -0800
+++ b/compute_Cy.cxx	Wed Feb 29 15:31:48 2012 -0800
@@ -62,29 +62,7 @@ void compute_Cy(const Model &model, cons
               Cy-=dzx_yx_correction;
 
               dCy_dvy-=(eta_jump*dvx_y_dv + dx*dzx_yx_dv)/h;
-
-              if(j==Ny/8)
-                std::cout << "Cy "
-                          << i << " "
-                          << j << " "
-                          << x << " "
-                          << middle << " "
-                          << zy[i][j] << " "
-                          << Cy << " "
-                          << dCy_dvy/eta << " "
-                          << "\n";
             }
-
-          if(j==Ny/8)
-            std::cout << "Cy out "
-                      << i << " "
-                      << j << " "
-                      << x << " "
-                      << middle << " "
-                      << zy[i][j] << " "
-                      << Cy << " "
-                      << dCy_dvy/eta << " "
-                      << "\n";
           dCy_dzy=dCy_dvy/eta;
         }
     }
diff -r bd123bf91338 -r 57ee741e9644 main.cxx
--- a/main.cxx	Wed Feb 29 15:30:17 2012 -0800
+++ b/main.cxx	Wed Feb 29 15:31:48 2012 -0800
@@ -58,18 +58,6 @@ int main()
   initial(model,zx,zy,log_etax,log_etay,p,fx,fy);
 
   double Resid_p[Nx][Ny], Resid_x[Nx+1][Ny], Resid_y[Nx][Ny+1];
-
-
-  for(int i=0;i<Nx+1;++i)
-    for(int j=0;j<Ny;++j)
-      zx[i][j]=0;
-
-  for(int i=0;i<Nx;++i)
-    for(int j=0;j<Ny+1;++j)
-      zy[i][j]=0;
-
-  // zx[static_cast<int>(middle/h)][Ny/8]=1;
-  zy[static_cast<int>(middle/h)-1][Ny/8]=1;
 
   for(int iteration=0;iteration<max_iterations;++iteration)
     {



More information about the CIG-COMMITS mailing list