예제 #1
0
# === MODEL

from spux.models.randomwalk import Randomwalk
dt = 0.1
model = Randomwalk(dt=dt)
framework.assign(model)

# === INIT FRAMEWORK

model.attach(workers=None)

from spux.utils.seed import Seed
from spux.utils.sandbox import Sandbox

sandbox = Sandbox('sandbox-synthesize')
seed = Seed(4)

# setup SPUX framework
framework.setup(sandbox=sandbox, seed=seed, verbosity=1, outputdir='datasets')

# init SPUX framework
framework.init()

# === SYNTHESIS

from error import error
from inputset import inputset
from timeset import timeset
from spux.io.loader import load_dict
import numpy
connector = utils.select ('auto')

# === SAMPLER with attached executors

from script_executors import sampler

# === SAMPLING

# SANDBOX
# use fast tmpfs
from spux.utils.sandbox import Sandbox
sandbox = Sandbox (path = '/dev/shm/spux-sandbox')

# SEED
from spux.utils.seed import Seed
seed = Seed (8)

# init executor
sampler.executor.init (connector)

# setup sampler
sampler.setup (sandbox = sandbox, verbosity = 2, index = 0, seed = seed, lock = 50, informative = True)

# init sampler (use prior for generating starting values)
sampler.init ()

# generate samples from posterior distribution
sampler (12)

# exit executor
sampler.executor.exit ()
예제 #3
0
from hydro import Hydro
from precipitation import precipitation_coefficients
from initial import initial

model = Hydro (precipitation_coefficients, initial, dt=10, integrator='RKSSP', numbify=0)
framework.assign (model)

# === INIT FRAMEWORK

model.attach (workers = None)

from spux.utils.seed import Seed
from spux.utils.sandbox import Sandbox

sandbox = Sandbox ('sandbox-synthesize')
seed = Seed (1)

# setup SPUX framework
framework.setup (sandbox = sandbox, seed = seed, verbosity = 1, outputdir = 'datasets_synthetic')

# init SPUX framework
framework.init ()

# === SYNTHESIS

from exact import exact
# from error import error
error = None

parameters = exact ['parameters']
steps = 10000