[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