Skip to content
Snippets Groups Projects
Commit 9d7d9b2a authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Additional changes to the two MC studies

parent 08f59dbb
No related branches found
No related tags found
No related merge requests found
...@@ -74,7 +74,6 @@ for depth in depths: ...@@ -74,7 +74,6 @@ for depth in depths:
'lower': Delta_center-DeltaDelta, 'lower': Delta_center-DeltaDelta,
'upper': Delta_center+DeltaDelta}, 'upper': Delta_center+DeltaDelta},
} }
theta_prior = {'F0': {'upper': 30.000006381121477, 'lower': 29.999993618878523, 'type': 'unif'}, 'F1': {'upper': 1.0143020701400378e-10, 'lower': 9.8569792985996225e-11, 'type': 'unif'}, 'F2': 0, 'Delta': {'upper': -0.20155527961896461, 'lower': -0.24155527961896459, 'type': 'unif'}, 'Alpha': {'upper': 2.8924321897264367, 'lower': 2.8524321897264366, 'type': 'unif'}}
ntemps = 2 ntemps = 2
log10temperature_min = -1 log10temperature_min = -1
......
Source diff could not be displayed: it is too large. Options to address this: view the blob.
...@@ -4,7 +4,6 @@ import os ...@@ -4,7 +4,6 @@ import os
import sys import sys
import time import time
ID = sys.argv[1] ID = sys.argv[1]
outdir = sys.argv[2] outdir = sys.argv[2]
...@@ -31,9 +30,9 @@ VF0 = VF1 = 100 ...@@ -31,9 +30,9 @@ VF0 = VF1 = 100
DeltaF0 = VF0 * np.sqrt(3)/(np.pi*Tspan) DeltaF0 = VF0 * np.sqrt(3)/(np.pi*Tspan)
DeltaF1 = VF1 * np.sqrt(45/4.)/(np.pi*Tspan**2) DeltaF1 = VF1 * np.sqrt(45/4.)/(np.pi*Tspan**2)
depths = np.linspace(100, 400, 9) depths = np.linspace(100, 400, 7)
nsteps = 20 nsteps = 50
run_setup = [((nsteps, 0), 20, False), run_setup = [((nsteps, 0), 20, False),
((nsteps, 0), 7, False), ((nsteps, 0), 7, False),
((nsteps, 0), 2, False), ((nsteps, 0), 2, False),
...@@ -58,17 +57,17 @@ for depth in depths: ...@@ -58,17 +57,17 @@ for depth in depths:
startTime = time.time() startTime = time.time()
theta_prior = {'F0': {'type': 'unif', theta_prior = {'F0': {'type': 'unif',
'lower': F0-DeltaF0/2., 'lower': F0_center-DeltaF0,
'upper': F0+DeltaF0/2.}, 'upper': F0_center+DeltaF0},
'F1': {'type': 'unif', 'F1': {'type': 'unif',
'lower': F1-DeltaF1/2., 'lower': F1_center-DeltaF1,
'upper': F1+DeltaF1/2.}, 'upper': F1_center+DeltaF1},
'F2': F2, 'F2': F2,
'Alpha': Alpha, 'Alpha': Alpha,
'Delta': Delta 'Delta': Delta
} }
ntemps = 1 ntemps = 2
log10temperature_min = -1 log10temperature_min = -1
nwalkers = 100 nwalkers = 100
......
...@@ -9,4 +9,4 @@ Log=CollectedOutput/log.$(Process) ...@@ -9,4 +9,4 @@ Log=CollectedOutput/log.$(Process)
request_cpus = 1 request_cpus = 1
request_memory = 8 GB request_memory = 8 GB
Queue 400 Queue 100
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment