The following code produces two qualitatively different sets of results on about 50% of runs with SpatialKappa v2.1.0. On 50% of runs NMDA transitions from the c0 to the c1 state, as expected. On the other 50% of runs no transitions appear to take place.
The problem seems to be the use of the implicit definition of Glu(). If this is replaced with the explicit definition of the default state Glu(x), the expected behaviour occurs.
%agent: NMDA(b0, state~c0~c1)
%agent: Glu(x)
'C0-C1' NMDA(b0, state~c0), Glu(x) -> NMDA(b0!1, state~c1), Glu(x!1) @ 0.008
'Create Glu' -> Glu(x) @ 31207.554688
%init: 15 NMDA()
# %init: 15 NMDA(b0, state~c0)
%init: 1 Glu() # BREAKS sims about 50% of time
# %init: 1 Glu(x) # DOESN'T BREAK sims at all
%var: 'Glu' Glu(x)
%var: 'NMDAC0' NMDA(state~c0)
%var: 'NMDAC1' NMDA(state~c1)
%plot: 'NMDAC0'
%plot: 'NMDAC1'
The following code produces two qualitatively different sets of results on about 50% of runs with SpatialKappa v2.1.0. On 50% of runs
NMDAtransitions from thec0to thec1state, as expected. On the other 50% of runs no transitions appear to take place.The problem seems to be the use of the implicit definition of
Glu(). If this is replaced with the explicit definition of the default stateGlu(x), the expected behaviour occurs.