ECE 5760: Laboratory 4

Mandelbrot set visualization.

Introduction.

In this assignment you will implement a system to compute and draw the Mandelbrot set, and to zoom in on pieces of the set. A matlab program shows how the computation is done. Note that the numbers are all complex. The image below is from the program output. The x,y axis units are pixels, not natural units. The x axis range is [-2,1] and the y axis [-1,1]. And Shiva Rajagopal found a video for Mandelbrot Set made by Cornell students! https://www.youtube.com/watch?v=ES-yKOYaXq0 Lyrics: http://www.jonathancoulton.com/wiki/Mandelbrot_Set/Lyrics

Procedures:

• There may be one or more serial processors to control computation and serial communication.
• You can use additional NiosII processors, pancake processors or build custom hardware to compute the image. The computation of a 640x480 image is compute intensive, so you are going to want to implement multiple computational units to speed up the process. You may want to consider load balancing since some regions of the complex plane are much faster to compute than others.
• If you use custom hardware, I suggest using 3.33 fixed point notation for interating the quadratic complex number calculation. The top bit will be the usual 2's comp sign bit. This format give the numbers a dynamic range of +/-3. A matlab program to check the effect of limiting the range shows that the computation works. If you use NiosII cpus to do the calcuations use 3.29 format.

Examples from class:

• Skyler Schneider 1, 2, 3, 4
• Venkatesh Santhanagopalan and Rick Wong 1, 2, 3
• Jinda Cui and Jiawei Yang, Full set.;
  Satellite double spiral: x_mid=-0.743 643 900 055; y_mid=0.131 825 890 901; zoom level=3000;
  neural network: range mode: x=(-0.37432239 -0.373655726 ); y=(0.65938433 0.66005100); zoom mode: x_mid=-0.37398906; y_mid=0.65971767; zoom level=3000;
  devil's eye: range mode: x=(-1.87 -1.85); y=(-0.01 0.01); zoom mode: x_mid=-1.86; y_mid=0; zoom level=100;
  alien's egg: range mode: x=(-1.863 -1.85); y=(-0.00153846 0.00153846); zoom mode: x_mid=-1.8615; y_mid=0; zoom level=650;
  lava: range mode: x=(-1.38 -1.36);y=(-0.026 -0.006); zoom mode: x_mid=-1.37; y_mid=-0.016; zoom level=100;
• Sandeep Gangundi and Sion Wang 600x_zoom
  

  image name	x coord	y coord
  flower	[-0.37465401, -0.37332411]	[0.659227668, 0.66020767]
  asymmetric mandelbrot	[0.435396403, 0.451687191]	[0.367981352, 0.380210061]
  pattern	[-0.758,-0.75]	[0.05,0.06]
  spiral	[-0.403, -0.399]	[-0.600, -0.603]

Other links

• A student report on implementing a mandelbrot generator on a DE1 board.
• Mark Bowers FPGA-based Mandelbrot generator
• Mandelbrot set
• MIke Field Mandelbrot on FPGA
• Who needs VRAM anyway? Or, calculating Mandelbrot on an FPGA
• Properties of Mandelbrot set
• Optimizations

• There may be one or more serial processors to control computation and serial communication,
  but you do not need to use a serial processor if you do not want to.
• The corners of the region being plotted should be displayed either on the VGA screen, on a serial console, or on the LCD display.
• Upon pressing KEY3, the system should compute and display the full Mandelbrot set, as shown above.
• The VGA graphics interface should run the VGA at 640x480 resolution with 8-bit color per pixel. Colors may aid in debugging and certainly look cool. You may want to map color to the number of iterations to divergence, or to the log of that number. A reasonable approximation of the log is to just use the position of the largest non-zero bit in the count (in matlab fix(log2(count)). A modified matlab program and image show the effect of log-compressing the numerical range. The zoomed image (see below) also shows different detail.
• There will be a zoom interface using serial communication to a terminal or a graphical interface. It is sufficient to pick a center coordinate and zoom in/out by a factor of two. When you trigger the zoom calculation, the specified corners should be reassigned to the corners of the display and 640x480 new points in those intervals should be iterated again.
• Part of your grade will depend upon how fast you can render the full 640x480 set and two other regions picked on demo day. The elapsed time in seconds and fractions of a second to draw should be shown in decimal on the 7-segment LEDs, or LCD, or console or on the VGA screen.
• You should be able to choose the maximum number of iterations using serial input, or toggle switches, or however you want to do it.

Be prepared to demo your design to your TA in lab.

Your written lab report should include the sections mentioned in the policy page, and :

• A table of times to compute various regions of the Mandelbrot set close to the areas specified in the list.
  • x=[-2 1], y=[-1 1]
  • x=[-0.758,-0.75], y=[0.05,0.06]
  • x=[-1.45, -1.3], y=[-0.07, 0.07]
• A collection of photographs of your favorite regions, and their coordinates. The camera should be good enough to resolve individual pixels.
• A detailed dsecription of your SOPC design.
• A heavily commented listing of your Verilog design and GCC code (if you use NiosII) or Syrup code (if you use Pancake).