Return to Answer

UPDATE I thought it would be neat to try and animate the thing, so I let the $a$ parameter run between $-\pi$ and $\pi$. I generated 600 images and put them together using ffmpeg. Check it out on youtube.

UPDATE I thought it would be neat to try and animate the thing, so I let the $a$ parameter run between $-\pi$ and $\pi$. I generated 600 images and put them together using ffmpeg. Check it out on youtube.

It might not be in the spirit of Mathematica Stack Exchange, but allow me an objection - stuff that is slow in Mathematica should be kept out of it. To wit consider how this little C++ nugget does the grunt work:

UPDATE I thought it would be neat to try and animate the thing, so I let the $a$ parameter run between $-\pi$ and $\pi$. I generated 600 images and put them together using ffmpeg. Check it out on youtube.

It might not be in the spirit of Mathematica Stack Exchange, but allow me an objection - stuff that is slow in Mathematica should be kept out of it. To wit consider how this little C++ nugget does the grunt work:

And this Mathematica code that makes an image and colorizes it from the produced data:

buffer = BinaryReadList["image.bin", "Real32"]; dim = Sqrt[Length@buffer]; shape = Partition[buffer, dim]; bigimg = Image[shape, "Real32", ImageSize -> {dim, dim}]; pic = Rasterize[bigimg, ImageSize -> dim/4]; cf = ColorData["SunsetColors"]; cpic = Colorize[pic, ColorFunction -> cf] 

And this Mathematica code makes an image and colorizes it from the produced data:

buffer = BinaryReadList["image.bin", "Real32"]; dim = Sqrt[Length@buffer]; bigimg = Image[Partition[buffer, dim], "Real32"]; Colorize[Rasterize[bigimg, ImageSize -> dim/4],   ColorFunction -> ColorData["SunsetColors"]]