src.scripts.sim module

src.scripts.sim.axes_sim(fig, freqx, freqy, waitx, waity)

Plot histogram of frequency and wait time

Example:

figure(pjoin([dir_out,'images']),func=axes_sim,funcdata=[freq_bins/240,freq_hist,wait_bins,wait_hist])

src.scripts.sim.freq_wait(a_lst)

reshape (T,X^2) to (T,X,X)

calculate excitation pixels, new excitation pixels, frequency and wait time of new excitations

src.scripts.sim.noise0d(X=16384, T=240, DT=3, amp=0.1, frac=0.03)

for each x position, choose a fraction of uniformly distributed time points, apply perturbations of duration DT

src.scripts.sim.solve_params(a_lst)

solve FitzhughNagumo reaction with stochastic perturbation on 1D grid

save the time evolution as tif images with shape (T,X^2)