Preliminary 1 Simulations

The code appears to be running so it is now possible to obtain some prelimiary results for the base set of paramters to investigate, $E_{\text{APPLIED}}$=-0.16, -0.18, -0.20, -0.22, -0.24, -0.26, -0.28, -0.30; $c_{\theta}^{\infty}$=0.006, 0.012


In [1]:
from distributed import LocalCluster
from distributed import Client
from extremefill2D.fextreme import init_sim, restart_sim, iterate_sim, multi_init_sim
from extremefill2D.fextreme.plot import vega_plot_treants, vega_plot_treant
import vega
from extremefill2D.fextreme.tools import get_by_uuid, outer_dict, pmap
from toolz.curried import map, pipe, curry
import itertools

%reload_ext yamlmagic

In [2]:
cluster = LocalCluster(nanny=True, n_workers=8, threads_per_worker=1)
client = Client(cluster)
client


Out[2]:
<Client: scheduler="127.0.0.1:8786" processes=8 cores=8>

In [6]:
client.shutdown()

In [4]:
treants = multi_init_sim('../../scripts/params1.json',
                         '../../data',
                         pmap(client),
                         dict(appliedPotential=(-0.16, -0.18, -0.20, -0.22, -0.24, -0.26, -0.28, -0.30),
                              bulkSuppressor=(0.006, 0.012)),
                         tags=['prelim1'])

In [5]:
print(treants)


[<Treant: '550e1938-5843-4e72-8f2c-5209e68bc53c'>, <Treant: '766e328f-cecd-4fc7-aaa6-9c59efc528bb'>, <Treant: '63124df8-95ea-4cf5-9f64-ae2040f6dc5a'>, <Treant: '41727a23-eb07-485d-9587-3f97f94a3f70'>, <Treant: 'd3255296-2d7a-4857-a49c-2e21259b0959'>, <Treant: '068a2809-eb5a-4a00-bc11-8d5820dce979'>, <Treant: '35a13610-ac65-4c5b-b24e-cd16b912fed4'>, <Treant: '5950f67d-fb09-4292-8633-fe2e46c9b738'>, <Treant: '1eb186bc-ea19-485a-bc8c-31933dbf0980'>, <Treant: '2431195c-f15f-423b-a5b3-5b112affa977'>, <Treant: '506cee39-53ec-4ce7-b57e-56c9c528859c'>, <Treant: 'f4e978e0-ef90-4856-bc3f-f28ad898138a'>, <Treant: 'd7edd8f4-71e6-4d93-bcfd-644d715b1ac2'>, <Treant: '95499a89-a7fd-4aae-af34-dc91cf0234ce'>, <Treant: 'dbd3ed44-449f-4b48-be0a-730c3718f2eb'>, <Treant: '17d5767e-c57d-4c90-b015-e767ad66e495'>]

In [6]:
treant_and_errors = pmap(client)(iterate_sim(iterations=20, steps=100), treants)

In [ ]: