Płonący statek

Dzisiaj koniec serii nt. różnych fraktali. Skończymy w płomieniach.

Ponad 2 lata temu na kursie asemblera na PWr mieliśmy za zadanie napisać program w dowolnym języku wysokiego poziomu wykorzystując przy tym zarówno wstawki asemblerowe jak i wywołania procedur napisanych w czystym asemblerze, dołączanych jako pliki .OBJ. Dodatkowo w ramach ćwiczenia trzeba było wykorzystać koprocesor numeryczny. Mój wybór na Pascala z racji młodzieńczych sentymentów. A czy jest coś fajniejszego i prostszego do napisania na koprocesor niż fraktale i tryb 13h? Raczej nie ;)

W początkowej wersji program obsługiwał tylko zbiory Julii, wczoraj i dzisiaj dopisałem zbiór Mandelbrota oraz bardzo ciekawy fraktal, "Płonący statek".

Płonący statek jest niezwykle prostą modyfikacją zbioru Mandelbrota. Fraktal ten to zbiór takich liczb zespolonych c, że ciąg z_n definiowany jako
 z₀ = 0
 z_n+1 = (|Re(z_n)| + i * |Im(z_n)|)² + c
nie dąży do nieskończoności przy n -> ∞.

Jak widać jedyną różnicą między płonącym statkiem o zbiorem Mandelbrota jest wartość bezwzględna z części rzeczywistej i części urojonej z_n. Jedna taka mała modyfikacja (cztery dodatkowe linijki kodu) znacząco zmienia wygląd fraktala.

Dodany tutaj programik działa dość sprawnie - oczywiście nie na sprzęcie z epoki 386, aż tak dobrze nie ma ;) Wykorzystałem tryb 13h udostępniany przez karty VGA. Dzięki temu wyświetlanie grafiki było proste i przyjemne, wadą jest oczywiście niska rozdzielczość (320x200), prostokątne piksele, tylko 256 kolorów (i tak wykorzystuję tylko szarości) i brak możliwości zrobienia print screena ;) Teraz chyba już nikt nie rysuje po ekranie odwołując się bezpośrednio do pamięci karty graficznej ;)

Klawiszologia:
W/S/A/D - przesuwanie fraktala
Q/E - pomniejszenie/powiększenie
+/- - zwiększenie/zmniejszenie dokładności obliczeń
1/2/3 - wybór fraktala: zbiór Julii/Mandelbrota/Płonący statek
U/J/I/K - dla zbioru Julii modyfikacja parametru c
Z - wyjście


Program do pobrania:
http://wojtekm.ovh.org/fractals.exe
Kody źródłowe:
http://pastebin.com/kumyrRum - cały program.
http://pastebin.com/RrJb6cP2 - procedury do obsługi trybu 13h.

Więcej o samych fraktalach:
http://en.wikipedia.org/wiki/Julia_set
http://en.wikipedia.org/wiki/Mandelbrot_set
http://en.wikipedia.org/wiki/Burning_Ship_fractal

I trybie 13h: http://en.wikipedia.org/wiki/Mode_13h