MBR model (modified Beeler-Reuter).
2D spiral waves

"the BR model (mod = 1) produces breakup because it has a very long wave length and adapts to very short at high frequencies and leads to breakup via slow recovery fronts (the wave back is slow and the wave front is too fast and crashes on it). By speeding up the calcium dynamics by two (mod = 2) we can decrease the wave length and avoid breakup."   Flavio H. Fenton

dtms dxcm   it delay fps
On 2D plane Vm values are displayed by the red color. 2 x values are displayed by the green color. This script uses simple explicit algorithm and 512512 float32 textures with the CLAMP_TO_EDGE boundary conditions. It makes 2 it time steps per frame.

Now we can trace spiral tip dynamics.

Heart rhythms     updated 22 Jan 2012