trying simple fixed point iteration now, x <- mix(x, F(x-1`mod`n), g). g = 1/2 in my current test, convergence metric fluctuating around 0.005, want it to be around 1000x smaller to consider it ready....

strobing gif Show more

but these are solutions to the fixed point iteration, not necessarily solutions to the reaction-diffusion time crystal problem...

found a typo in my earlier jacobian code, no wonder it didn't work properly

links to full screen stroboscopic HTML Show more

I fixed all the Jacobian issues I hope (there was some i,j vs j,i confusion as well as the U vs V typos...)

Now I'm redetecting the optimum loop period at every iteration in case it changes, a little slower but gives smoother loops it seems...

link to stroboscopic HTML Show more


The O(S^4) asymptotics are harsh.

I can run 8 copies at 64x64 without running out of memory, 4 copies at 80x80, and conjecture only 1 copy at 112x112.

If I used single precision floating point instead of double precision, or paged out half of the data (accessed infrequently) to disk, I would be able to run 1 copy at 134x134. But I don't know if there is a version of GSL that supports single precision.

If I did both, I could run 1 copy at 160x160, which is getting into the realm of interesting sizes...

Sign in to participate in the conversation

Welcome to, an instance for discussions around cultural freedom, experimental, new media art, net and computational culture, and things like that. This is part of a family of services that include mailing lists, group chat, and XMPP.