mesosciencelab / mesohops Goto Github PK

When specifying an integration_param input dictionary manually for a HOPS trajectory calculation, the default dictionary is completely overwritten, removing any non-specified settings entirely. For example, in the quickstart guide, integration parameters are set in the following dictionary:

integration_param = {"INTEGRATOR": "RUNGE_KUTTA"}

However, if running a hops trajectory:

from mesohops.dynamics.hops_trajectory import HopsTrajectory as HOPS

hops = HOPS(
    sys_param,
    noise_param=noise_param,
    hierarchy_param=hierarchy_param,
    eom_param=eom_param,
    integration_param=integration_param,
    )

in which all other dictionaries are satisfied, this will result in a key error for a missing "EARLY_ADAPTIVE_INTEGRATOR" dictionary entry rather than maintaining the default dictionary

integration_param=dict(
    INTEGRATOR="RUNGE_KUTTA",
    EARLY_ADAPTIVE_INTEGRATOR='INCH_WORM',
    EARLY_INTEGRATOR_STEPS=5,
    INCHWORM_CAP=5,
    STATIC_BASIS=None,
    EFFECTIVE_NOISE_INTEGRATION=False,
)

in any portion.