[cig-commits] [commit] inversion: update inversion file, add script (1098a3b)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/burnman

On branch  : inversion
Link       : https://github.com/geodynamics/burnman/compare/9cf237852c30754d89efcc731fd5d788251e3d96...1098a3bf1b7cf8d21e6a06c1cca7639f1b95764a

>---------------------------------------------------------------

commit 1098a3bf1b7cf8d21e6a06c1cca7639f1b95764a
Author: Timo Heister <timo.heister at gmail.com>
Date:   Sat Dec 13 01:03:10 2014 -0500

    update inversion file, add script


>---------------------------------------------------------------

1098a3bf1b7cf8d21e6a06c1cca7639f1b95764a
 inversion/run_orig.sh                                    | 14 ++++++++++++++
 inversion/setup_isochemical_newmisfit_moduli_realdata.py |  6 +++---
 2 files changed, 17 insertions(+), 3 deletions(-)

diff --git a/inversion/run_orig.sh b/inversion/run_orig.sh
new file mode 100755
index 0000000..8d22241
--- /dev/null
+++ b/inversion/run_orig.sh
@@ -0,0 +1,14 @@
+#!/bin/bash
+
+for i in "1" "2" "3" "4";
+do
+  f=orig_$i
+  if [ -e $f ]
+  then
+    python setup_isochemical_newmisfit_moduli_realdata.py continue $f &
+  else
+    python setup_isochemical_newmisfit_moduli_realdata.py run $f &
+  fi
+done
+
+wait
diff --git a/inversion/setup_isochemical_newmisfit_moduli_realdata.py b/inversion/setup_isochemical_newmisfit_moduli_realdata.py
index ace1357..cb9056f 100644
--- a/inversion/setup_isochemical_newmisfit_moduli_realdata.py
+++ b/inversion/setup_isochemical_newmisfit_moduli_realdata.py
@@ -156,20 +156,20 @@ if whattodo=="run":
     whattodo="continue"
 
 if whattodo=="continue":
-    n_runs = 2 # results in 1000 new points
+    n_runs = 100
     for l in range(0,n_runs):
         db = pymc.database.pickle.load(dbname)
         print "*** run=%d/%d, # samples: %d" % (l, n_runs, db.trace('fraction_pv').stats()['n'] )
         S = pymc.MCMC(model, db=db)
         #S.sample(iter=100, burn=10, thin=1)
-        S.sample(iter=10000, burn=0, thin=10) # Search space for 100000 acceptable steps, forget first 1000 and save every 10.
+        S.sample(iter=1000, burn=0, thin=10) # Search space for 100000 acceptable steps, forget first 1000 and save every 10.
         S.db.close()
 
 if whattodo=="plot":
 	files=sys.argv[2:]
 	print "files:",files
         
-        b=0#10000 # burn number
+        b=100#10000 # burn number
         i=1
 
         for t in things: