[cig-commits] commit: Change from a rectangle to a square
Mercurial
hg at geodynamics.org
Fri Feb 17 09:47:21 PST 2012
changeset: 51:d8119d27f1e7
user: Walter Landry <wlandry at caltech.edu>
date: Fri Feb 17 09:45:39 2012 -0800
files: constants.hxx
description:
Change from a rectangle to a square
diff -r 5b1d92fd417b -r d8119d27f1e7 constants.hxx
--- a/constants.hxx Fri Feb 17 09:44:44 2012 -0800
+++ b/constants.hxx Fri Feb 17 09:45:39 2012 -0800
@@ -1,15 +1,16 @@ const int NN(512);
-const int NN(512);
+const int NN(128);
const int Nx(NN);
-const int Ny(2*NN);
+const int Ny(NN);
const double min_eta=1;
-const double max_eta=1e4;
+const double max_eta=1e3;
const double eta_jump=max_eta-min_eta;
#include <cmath>
const double log_max_eta=std::log(max_eta);
const double log_min_eta=std::log(min_eta);
-const double middle=(25 + 0.00000001)/64;
+// const double middle=(25 + 0.00000001)/64;
// const double middle=(25 + 32/Nx + 0.00000001)/64;
-// const double middle=0.4;
+// const double middle=(25 - 32/2048.0)/64;
+const double middle=0.4;
const double h(1.0/Nx);
const double pi(atan(1.0)*4);
More information about the CIG-COMMITS
mailing list