Wrzlprmft
commented
2 years ago I tried to reproduce your problem by running the following code:
import symengine
import numpy as np
from jitcdde import jitcdde, y, t
symbol_names = 'τ, βH, βC, βD, βN, η, α, σ, λm, λp, κC, κD, κE, fc'
param_symbols = τ, βH, βC, βD, βN, η, α, σ, λm, λp, κC, κD, κE, fc = symengine.symbols(symbol_names)
f1 = 1
f2 = 1
def model():
h1, c1, d1, e1, n1, s1 = [y(i) for i in range(6)]
h2, c2, d2, e2, n2, s2 = [y(i+6) for i in range(6)]
h1d = y(0, t-τ)
s1d = y(5, t-τ)
h2d = y(6, t-τ)
s2d = y(11, t-τ)
h1c = y(0, t-τ*fc)
h2c = y(6, t-τ*fc)
yield -h1 + βH * (1 + α * (σ * s1d)**η ) / ( (1 + (σ * s1d)**η) * (1 + h1d**η) )
yield -c1 + βC / (1 + h1c**η) + λm * e1 - λp * c1 * d1 - κC * c1 * f2
yield -d1 + βD + λm * e1 - λp * c1 * d1 - κD * d1 * f2
yield -e1 - λm * e1 + λp * c1 * d1 - κE * e1 * f2
yield -f1 + βN - f1 * (κC * c2 + κD * d2 + κE * e2)
yield -s1 + f1 * (κC * c2 + κD * d2 + κE * e2)
yield -h2 + βH * (1 + α * (σ * s2d)**η ) / ( (1 + (σ * s2d)**η) * (1 + h2d**η) )
yield -c2 + βC / (1 + h2c**η) + λm * e2 - λp * c2 * d2 - κC * c2 * f1
yield -d2 + βD + λm * e2 - λp * c2 * d2 - κD * d2 * f1
yield -e2 - λm * e2 + λp * c2 * d2 - κE * e2 * f1
yield -f2 + βN - f2 * (κC * c1 + κD * d1 + κE * e1)
yield -s2 + f2 * (κC * c1 + κD * d1 + κE * e1)
# initialize integrator
DDE = jitcdde(model, verbose=True, max_delay=3, n=12, control_pars=param_symbols)
model_location = DDE.save_compiled('model.so', overwrite=True)… which works fine for me.
Note that I had to pick some values for f1 and f2 here, so my wild guess would be that those have a problematic type, structure, value, or similar in your case. How you wrap things in functions and similar should not affect the kind of error you are observing.
mwappner
I'm getting a long and cryptic message when trying to compile a module to later save it to a file to make multiple short integrations in parallel. The error in question reads
with a really long traceback, which I can try to provide if needed. The system in question is 12-D composed of two 6-D coupled systems. Here's the code, which I fear I need to provide in full, since a more minimal example does work (see below). Assume
get_maximum_delayis a function defined somewhere that retrieves the maximum valueτandτ*fcwill have. You can setmaxdelay=3for testing purposes.After running this, which raises the error, the program is supposed to launch a bunch of processes to load the model and integrate it in parallel for a bunch of different parameter values.
Just as a sanity check, here's a very minimal example of a 2D system that does work when doing this. Ignore of course the fact that we are not multiprocessing here, since the error arises before calling the integration anyway.