(*) Added BOLD_ode_precompute, which uses a neuronal signal stored in the global BOLDParams. Speed-ups of an order of magnitude. Has to be ode1 to get correct sampling of external signal (using resample/decimate)
(*) Added support for the ode_precompute code using ode1, which correctly samples the neuronal signal for simulating the BOLD response
(*) Initially added a rounding function so that the half-steps of ode4 match the pre-computed values, but actually this distorts the input, so I stuck with ode1. Half steps with external drives can be difficult. Don't know the theory of it, but it doesn't seem to be working properly...
(*) Added an option to just sample a subset of time in plot_V (long runs otherwise obscure the local dynamics, which is the point of this function)
(*) Better axis labels on plot_structure_function
(*) ReSampleBOLD added -- resamples the BOLD signal to match an extermental sampling rate. Still playing around with this script.
(*) decimatecolumnwise is a short function to decimate each column of a matrix using the same parameters. Probably a bsxfun way to do this or something, but oh well here this works.