/ concept-collection / turing-surface
Sign in
concept-collection / turing-surface
23 lines · 546 BCodeBlameHistory
62e6cc4Simplify UI text and built-in .m script commentsJeremy Magland 1% Allen-Cahn on a closed surface: interfaces form, then coarsen.
3% du/dt = eps2*lap_g(u) + u - u^3
4%
62e6cc4Simplify UI text and built-in .m script commentsJeremy Magland 5% Same scheme as models/schnakenberg.m.
7function [U, u] = init(noise)
8 U = analys(noise);
9 u = synth(U);
10end
12function [Un, u] = step(U, lam, gx, gy, gz, eps2, dt, niter)
13 u = synth(U);
15 Bu = U + dt * analys(u - u.^3);
16 Un = Bu ./ (1 + (dt * eps2) * lam);
18 for k = 1:niter
62e6cc4Simplify UI text and built-in .m script commentsJeremy Magland 19 % Placeholder: dlap = lap_g - lap_s (see models/schnakenberg.m).
0af3386Reaction-diffusion on spherical-harmonic surfacesJeremy Magland 20 dLu = 0 * Un;
21 Un = (Bu + (dt * eps2) * dLu) ./ (1 + (dt * eps2) * lam);
22 end
23end
moveopenescclose