fix wrong score recording.
[kugel-rb.git] / apps / plugins / mandelbrot.c
blob61469d759781d96cb8f6f10f2d21b62a5a495f6a
1 /***************************************************************************
2 * __________ __ ___.
3 * Open \______ \ ____ ____ | | _\_ |__ _______ ___
4 * Source | _// _ \_/ ___\| |/ /| __ \ / _ \ \/ /
5 * Jukebox | | ( <_> ) \___| < | \_\ ( <_> > < <
6 * Firmware |____|_ /\____/ \___ >__|_ \|___ /\____/__/\_ \
7 * \/ \/ \/ \/ \/
8 * $Id$
10 * Copyright (C) 2004 Matthias Wientapper
11 * Heavily extended 2005 Jens Arnold
14 * This program is free software; you can redistribute it and/or
15 * modify it under the terms of the GNU General Public License
16 * as published by the Free Software Foundation; either version 2
17 * of the License, or (at your option) any later version.
19 * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY OF ANY
20 * KIND, either express or implied.
22 ****************************************************************************/
23 #include "plugin.h"
25 #ifdef HAVE_LCD_BITMAP
26 #include "lib/grey.h"
27 #include "lib/xlcd.h"
29 PLUGIN_HEADER
31 /* variable button definitions */
32 #if CONFIG_KEYPAD == RECORDER_PAD
33 #define MANDELBROT_QUIT BUTTON_OFF
34 #define MANDELBROT_UP BUTTON_UP
35 #define MANDELBROT_DOWN BUTTON_DOWN
36 #define MANDELBROT_LEFT BUTTON_LEFT
37 #define MANDELBROT_RIGHT BUTTON_RIGHT
38 #define MANDELBROT_ZOOM_IN BUTTON_PLAY
39 #define MANDELBROT_ZOOM_OUT BUTTON_ON
40 #define MANDELBROT_MAXITER_INC BUTTON_F2
41 #define MANDELBROT_MAXITER_DEC BUTTON_F1
42 #define MANDELBROT_RESET BUTTON_F3
44 #elif CONFIG_KEYPAD == ARCHOS_AV300_PAD
45 #define MANDELBROT_QUIT BUTTON_OFF
46 #define MANDELBROT_UP BUTTON_UP
47 #define MANDELBROT_DOWN BUTTON_DOWN
48 #define MANDELBROT_LEFT BUTTON_LEFT
49 #define MANDELBROT_RIGHT BUTTON_RIGHT
50 #define MANDELBROT_ZOOM_IN BUTTON_SELECT
51 #define MANDELBROT_ZOOM_OUT BUTTON_ON
52 #define MANDELBROT_MAXITER_INC BUTTON_F2
53 #define MANDELBROT_MAXITER_DEC BUTTON_F1
54 #define MANDELBROT_RESET BUTTON_F3
56 #elif CONFIG_KEYPAD == ONDIO_PAD
57 #define MANDELBROT_QUIT BUTTON_OFF
58 #define MANDELBROT_UP BUTTON_UP
59 #define MANDELBROT_DOWN BUTTON_DOWN
60 #define MANDELBROT_LEFT BUTTON_LEFT
61 #define MANDELBROT_RIGHT BUTTON_RIGHT
62 #define MANDELBROT_ZOOM_IN_PRE BUTTON_MENU
63 #define MANDELBROT_ZOOM_IN (BUTTON_MENU | BUTTON_REL)
64 #define MANDELBROT_ZOOM_IN2 (BUTTON_MENU | BUTTON_UP)
65 #define MANDELBROT_ZOOM_OUT (BUTTON_MENU | BUTTON_DOWN)
66 #define MANDELBROT_MAXITER_INC (BUTTON_MENU | BUTTON_RIGHT)
67 #define MANDELBROT_MAXITER_DEC (BUTTON_MENU | BUTTON_LEFT)
68 #define MANDELBROT_RESET (BUTTON_MENU | BUTTON_OFF)
70 #elif (CONFIG_KEYPAD == IRIVER_H100_PAD) || \
71 (CONFIG_KEYPAD == IRIVER_H300_PAD)
72 #define MANDELBROT_QUIT BUTTON_OFF
73 #define MANDELBROT_UP BUTTON_UP
74 #define MANDELBROT_DOWN BUTTON_DOWN
75 #define MANDELBROT_LEFT BUTTON_LEFT
76 #define MANDELBROT_RIGHT BUTTON_RIGHT
77 #define MANDELBROT_ZOOM_IN BUTTON_SELECT
78 #define MANDELBROT_ZOOM_OUT BUTTON_MODE
79 #define MANDELBROT_MAXITER_INC (BUTTON_ON | BUTTON_RIGHT)
80 #define MANDELBROT_MAXITER_DEC (BUTTON_ON | BUTTON_LEFT)
81 #define MANDELBROT_RESET BUTTON_REC
83 #define MANDELBROT_RC_QUIT BUTTON_RC_STOP
85 #elif (CONFIG_KEYPAD == IPOD_4G_PAD) || \
86 (CONFIG_KEYPAD == IPOD_3G_PAD) || \
87 (CONFIG_KEYPAD == IPOD_1G2G_PAD)
88 #define MANDELBROT_QUIT (BUTTON_SELECT | BUTTON_MENU)
89 #define MANDELBROT_UP BUTTON_MENU
90 #define MANDELBROT_DOWN BUTTON_PLAY
91 #define MANDELBROT_LEFT BUTTON_LEFT
92 #define MANDELBROT_RIGHT BUTTON_RIGHT
93 #define MANDELBROT_ZOOM_IN BUTTON_SCROLL_FWD
94 #define MANDELBROT_ZOOM_OUT BUTTON_SCROLL_BACK
95 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
96 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
97 #define MANDELBROT_RESET (BUTTON_SELECT | BUTTON_PLAY)
99 #elif CONFIG_KEYPAD == IAUDIO_X5M5_PAD
100 #define MANDELBROT_QUIT BUTTON_POWER
101 #define MANDELBROT_UP BUTTON_UP
102 #define MANDELBROT_DOWN BUTTON_DOWN
103 #define MANDELBROT_LEFT BUTTON_LEFT
104 #define MANDELBROT_RIGHT BUTTON_RIGHT
105 #define MANDELBROT_ZOOM_IN_PRE BUTTON_SELECT
106 #define MANDELBROT_ZOOM_IN (BUTTON_SELECT | BUTTON_REL)
107 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_SELECT
108 #define MANDELBROT_ZOOM_OUT (BUTTON_SELECT | BUTTON_REPEAT)
109 #define MANDELBROT_MAXITER_INC_PRE BUTTON_PLAY
110 #define MANDELBROT_MAXITER_INC (BUTTON_PLAY | BUTTON_REL)
111 #define MANDELBROT_MAXITER_DEC_PRE BUTTON_PLAY
112 #define MANDELBROT_MAXITER_DEC (BUTTON_PLAY | BUTTON_REPEAT)
113 #define MANDELBROT_RESET BUTTON_REC
115 #elif CONFIG_KEYPAD == GIGABEAT_PAD
116 #define MANDELBROT_QUIT BUTTON_POWER
117 #define MANDELBROT_UP BUTTON_UP
118 #define MANDELBROT_DOWN BUTTON_DOWN
119 #define MANDELBROT_LEFT BUTTON_LEFT
120 #define MANDELBROT_RIGHT BUTTON_RIGHT
121 #define MANDELBROT_ZOOM_IN_PRE BUTTON_SELECT
122 #define MANDELBROT_ZOOM_IN (BUTTON_SELECT | BUTTON_REL)
123 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_SELECT
124 #define MANDELBROT_ZOOM_OUT (BUTTON_SELECT | BUTTON_REPEAT)
125 #define MANDELBROT_MAXITER_INC BUTTON_VOL_UP
126 #define MANDELBROT_MAXITER_DEC BUTTON_VOL_DOWN
127 #define MANDELBROT_RESET BUTTON_A
129 #elif CONFIG_KEYPAD == SANSA_E200_PAD
130 #define MANDELBROT_QUIT BUTTON_POWER
131 #define MANDELBROT_UP BUTTON_UP
132 #define MANDELBROT_DOWN BUTTON_DOWN
133 #define MANDELBROT_LEFT BUTTON_LEFT
134 #define MANDELBROT_RIGHT BUTTON_RIGHT
135 #define MANDELBROT_ZOOM_IN BUTTON_SCROLL_FWD
136 #define MANDELBROT_ZOOM_OUT BUTTON_SCROLL_BACK
137 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
138 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
139 #define MANDELBROT_RESET BUTTON_REC
141 #elif CONFIG_KEYPAD == SANSA_FUZE_PAD
142 #define MANDELBROT_QUIT (BUTTON_HOME|BUTTON_REPEAT)
143 #define MANDELBROT_UP BUTTON_UP
144 #define MANDELBROT_DOWN BUTTON_DOWN
145 #define MANDELBROT_LEFT BUTTON_LEFT
146 #define MANDELBROT_RIGHT BUTTON_RIGHT
147 #define MANDELBROT_ZOOM_IN BUTTON_SCROLL_FWD
148 #define MANDELBROT_ZOOM_OUT BUTTON_SCROLL_BACK
149 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
150 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
151 #define MANDELBROT_RESET (BUTTON_SELECT | BUTTON_REPEAT)
153 #elif CONFIG_KEYPAD == SANSA_C200_PAD
154 #define MANDELBROT_QUIT BUTTON_POWER
155 #define MANDELBROT_UP BUTTON_UP
156 #define MANDELBROT_DOWN BUTTON_DOWN
157 #define MANDELBROT_LEFT BUTTON_LEFT
158 #define MANDELBROT_RIGHT BUTTON_RIGHT
159 #define MANDELBROT_ZOOM_IN BUTTON_VOL_UP
160 #define MANDELBROT_ZOOM_OUT BUTTON_VOL_DOWN
161 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
162 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
163 #define MANDELBROT_RESET BUTTON_REC
165 #elif CONFIG_KEYPAD == SANSA_CLIP_PAD
166 #define MANDELBROT_QUIT BUTTON_POWER
167 #define MANDELBROT_UP BUTTON_UP
168 #define MANDELBROT_DOWN BUTTON_DOWN
169 #define MANDELBROT_LEFT BUTTON_LEFT
170 #define MANDELBROT_RIGHT BUTTON_RIGHT
171 #define MANDELBROT_ZOOM_IN BUTTON_VOL_UP
172 #define MANDELBROT_ZOOM_OUT BUTTON_VOL_DOWN
173 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
174 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
175 #define MANDELBROT_RESET BUTTON_HOME
177 #elif CONFIG_KEYPAD == SANSA_M200_PAD
178 #define MANDELBROT_QUIT BUTTON_POWER
179 #define MANDELBROT_UP BUTTON_UP
180 #define MANDELBROT_DOWN BUTTON_DOWN
181 #define MANDELBROT_LEFT BUTTON_LEFT
182 #define MANDELBROT_RIGHT BUTTON_RIGHT
183 #define MANDELBROT_ZOOM_IN BUTTON_VOL_UP
184 #define MANDELBROT_ZOOM_OUT BUTTON_VOL_DOWN
185 #define MANDELBROT_MAXITER_INC (BUTTON_SELECT | BUTTON_RIGHT)
186 #define MANDELBROT_MAXITER_DEC (BUTTON_SELECT | BUTTON_LEFT)
187 #define MANDELBROT_RESET (BUTTON_SELECT | BUTTON_UP)
189 #elif CONFIG_KEYPAD == IRIVER_H10_PAD
190 #define MANDELBROT_QUIT BUTTON_POWER
191 #define MANDELBROT_UP BUTTON_SCROLL_UP
192 #define MANDELBROT_DOWN BUTTON_SCROLL_DOWN
193 #define MANDELBROT_LEFT BUTTON_LEFT
194 #define MANDELBROT_RIGHT BUTTON_RIGHT
195 #define MANDELBROT_ZOOM_IN_PRE BUTTON_PLAY
196 #define MANDELBROT_ZOOM_IN (BUTTON_PLAY | BUTTON_REL)
197 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_PLAY
198 #define MANDELBROT_ZOOM_OUT (BUTTON_PLAY | BUTTON_REPEAT)
199 #define MANDELBROT_MAXITER_INC BUTTON_FF
200 #define MANDELBROT_MAXITER_DEC BUTTON_REW
201 #define MANDELBROT_RESET (BUTTON_PLAY | BUTTON_REW)
203 #elif CONFIG_KEYPAD == IRIVER_IFP7XX_PAD
204 #define MANDELBROT_QUIT BUTTON_EQ
205 #define MANDELBROT_UP BUTTON_UP
206 #define MANDELBROT_DOWN BUTTON_DOWN
207 #define MANDELBROT_LEFT BUTTON_LEFT
208 #define MANDELBROT_RIGHT BUTTON_RIGHT
209 #define MANDELBROT_ZOOM_IN_PRE BUTTON_SELECT
210 #define MANDELBROT_ZOOM_IN (BUTTON_SELECT | BUTTON_REL)
211 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_SELECT
212 #define MANDELBROT_ZOOM_OUT (BUTTON_SELECT | BUTTON_REPEAT)
213 #define MANDELBROT_MAXITER_INC (BUTTON_PLAY | BUTTON_RIGHT)
214 #define MANDELBROT_MAXITER_DEC (BUTTON_PLAY | BUTTON_LEFT)
215 #define MANDELBROT_RESET BUTTON_MODE
217 #elif CONFIG_KEYPAD == GIGABEAT_S_PAD
218 #define MANDELBROT_QUIT BUTTON_BACK
219 #define MANDELBROT_UP BUTTON_UP
220 #define MANDELBROT_DOWN BUTTON_DOWN
221 #define MANDELBROT_LEFT BUTTON_LEFT
222 #define MANDELBROT_RIGHT BUTTON_RIGHT
223 #define MANDELBROT_ZOOM_IN BUTTON_VOL_UP
224 #define MANDELBROT_ZOOM_OUT BUTTON_VOL_DOWN
225 #define MANDELBROT_MAXITER_INC BUTTON_PREV
226 #define MANDELBROT_MAXITER_DEC BUTTON_NEXT
227 #define MANDELBROT_RESET BUTTON_MENU
229 #elif CONFIG_KEYPAD == MROBE100_PAD
230 #define MANDELBROT_QUIT BUTTON_POWER
231 #define MANDELBROT_UP BUTTON_UP
232 #define MANDELBROT_DOWN BUTTON_DOWN
233 #define MANDELBROT_LEFT BUTTON_LEFT
234 #define MANDELBROT_RIGHT BUTTON_RIGHT
235 #define MANDELBROT_ZOOM_IN_PRE BUTTON_SELECT
236 #define MANDELBROT_ZOOM_IN (BUTTON_SELECT | BUTTON_REL)
237 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_SELECT
238 #define MANDELBROT_ZOOM_OUT (BUTTON_SELECT | BUTTON_REPEAT)
239 #define MANDELBROT_MAXITER_INC BUTTON_MENU
240 #define MANDELBROT_MAXITER_DEC BUTTON_PLAY
241 #define MANDELBROT_RESET BUTTON_DISPLAY
243 #elif CONFIG_KEYPAD == IAUDIO_M3_PAD
244 #define MANDELBROT_QUIT BUTTON_RC_REC
245 #define MANDELBROT_UP BUTTON_RC_VOL_UP
246 #define MANDELBROT_DOWN BUTTON_RC_VOL_DOWN
247 #define MANDELBROT_LEFT BUTTON_RC_REW
248 #define MANDELBROT_RIGHT BUTTON_RC_FF
249 #define MANDELBROT_ZOOM_IN_PRE BUTTON_RC_PLAY
250 #define MANDELBROT_ZOOM_IN (BUTTON_RC_PLAY | BUTTON_REL)
251 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_RC_PLAY
252 #define MANDELBROT_ZOOM_OUT (BUTTON_RC_PLAY | BUTTON_REPEAT)
253 #define MANDELBROT_MAXITER_INC_PRE BUTTON_RC_MODE
254 #define MANDELBROT_MAXITER_INC (BUTTON_RC_MODE|BUTTON_REL)
255 #define MANDELBROT_MAXITER_DEC_PRE BUTTON_RC_MODE
256 #define MANDELBROT_MAXITER_DEC (BUTTON_RC_MODE|BUTTON_REPEAT)
257 #define MANDELBROT_RESET BUTTON_RC_MENU
259 #elif CONFIG_KEYPAD == COWOND2_PAD
260 #define MANDELBROT_QUIT BUTTON_POWER
262 #elif CONFIG_KEYPAD == CREATIVEZVM_PAD
263 #define MANDELBROT_QUIT BUTTON_BACK
264 #define MANDELBROT_UP BUTTON_UP
265 #define MANDELBROT_DOWN BUTTON_DOWN
266 #define MANDELBROT_LEFT BUTTON_LEFT
267 #define MANDELBROT_RIGHT BUTTON_RIGHT
268 #define MANDELBROT_ZOOM_IN BUTTON_PLAY
269 #define MANDELBROT_ZOOM_OUT BUTTON_MENU
270 #define MANDELBROT_MAXITER_INC (BUTTON_UP | BUTTON_CUSTOM)
271 #define MANDELBROT_MAXITER_DEC (BUTTON_DOWN | BUTTON_CUSTOM)
272 #define MANDELBROT_RESET BUTTON_SELECT
274 #elif CONFIG_KEYPAD == PHILIPS_HDD1630_PAD
275 #define MANDELBROT_QUIT BUTTON_POWER
276 #define MANDELBROT_UP BUTTON_UP
277 #define MANDELBROT_DOWN BUTTON_DOWN
278 #define MANDELBROT_LEFT BUTTON_LEFT
279 #define MANDELBROT_RIGHT BUTTON_RIGHT
280 #define MANDELBROT_ZOOM_IN BUTTON_VIEW
281 #define MANDELBROT_ZOOM_OUT BUTTON_PLAYLIST
282 #define MANDELBROT_MAXITER_INC BUTTON_VOL_UP
283 #define MANDELBROT_MAXITER_DEC BUTTON_VOL_DOWN
284 #define MANDELBROT_RESET BUTTON_MENU
286 #elif CONFIG_KEYPAD == ONDAVX747_PAD || CONFIG_KEYPAD == MROBE500_PAD
287 #define MANDELBROT_QUIT BUTTON_POWER
289 #else
290 #error No keymap defined!
291 #endif
293 #ifdef HAVE_TOUCHSCREEN
294 #ifndef MANDELBROT_UP
295 #define MANDELBROT_UP BUTTON_TOPMIDDLE
296 #endif
297 #ifndef MANDELBROT_DOWN
298 #define MANDELBROT_DOWN BUTTON_BOTTOMMIDDLE
299 #endif
300 #ifndef MANDELBROT_LEFT
301 #define MANDELBROT_LEFT BUTTON_MIDLEFT
302 #endif
303 #ifndef MANDELBROT_RIGHT
304 #define MANDELBROT_RIGHT BUTTON_MIDRIGHT
305 #endif
306 #ifndef MANDELBROT_ZOOM_IN_PRE
307 #define MANDELBROT_ZOOM_IN_PRE BUTTON_TOPRIGHT
308 #endif
309 #ifndef MANDELBROT_ZOOM_IN
310 #define MANDELBROT_ZOOM_IN (BUTTON_TOPRIGHT | BUTTON_REL)
311 #endif
312 #ifndef MANDELBROT_ZOOM_OUT_PRE
313 #define MANDELBROT_ZOOM_OUT_PRE BUTTON_TOPLEFT
314 #endif
315 #ifndef MANDELBROT_ZOOM_OUT
316 #define MANDELBROT_ZOOM_OUT (BUTTON_TOPLEFT | BUTTON_REL)
317 #endif
318 #ifndef MANDELBROT_MAXITER_INC
319 #define MANDELBROT_MAXITER_INC BUTTON_BOTTOMRIGHT
320 #endif
321 #ifndef MANDELBROT_MAXITER_DEC
322 #define MANDELBROT_MAXITER_DEC BUTTON_BOTTOMLEFT
323 #endif
324 #ifndef MANDELBROT_RESET
325 #define MANDELBROT_RESET BUTTON_CENTER
326 #endif
327 #endif
329 #if LCD_DEPTH < 8
330 #define USEGSLIB
331 #define MYLCD(fn) grey_ub_ ## fn
332 #define MYLCD_UPDATE()
333 #define MYXLCD(fn) grey_ub_ ## fn
334 #else
335 #define UPDATE_FREQ (HZ/50)
336 #define MYLCD(fn) rb->lcd_ ## fn
337 #define MYLCD_UPDATE() rb->lcd_update();
338 #define MYXLCD(fn) xlcd_ ## fn
339 #endif
341 /* Fixed point format s5.26: sign, 5 bits integer part, 26 bits fractional part */
342 static long x_min;
343 static long x_max;
344 static long x_step;
345 static long x_delta;
346 static long y_min;
347 static long y_max;
348 static long y_step;
349 static long y_delta;
351 static int px_min = 0;
352 static int px_max = LCD_WIDTH;
353 static int py_min = 0;
354 static int py_max = LCD_HEIGHT;
356 static int step_log2;
357 static unsigned max_iter;
359 #ifdef USEGSLIB
360 GREY_INFO_STRUCT
361 static unsigned char *gbuf;
362 static size_t gbuf_size = 0;
363 static unsigned char imgbuffer[LCD_HEIGHT];
364 #else
365 static fb_data imgbuffer[LCD_HEIGHT];
366 #endif
368 /* 8 entries cyclical, last entry is black (convergence) */
369 #ifdef HAVE_LCD_COLOR
370 static const fb_data color[9] = {
371 LCD_RGBPACK(255, 0, 159), LCD_RGBPACK(159, 0, 255), LCD_RGBPACK(0, 0, 255),
372 LCD_RGBPACK(0, 159, 255), LCD_RGBPACK(0, 255, 128), LCD_RGBPACK(128, 255, 0),
373 LCD_RGBPACK(255, 191, 0), LCD_RGBPACK(255, 0, 0), LCD_RGBPACK(0, 0, 0)
375 #else /* greyscale */
376 static const unsigned char color[9] = {
377 255, 223, 191, 159, 128, 96, 64, 32, 0
379 #endif
381 #if CONFIG_CPU == SH7034
383 #define MULS16_ASR10(a, b) muls16_asr10(a, b)
384 static inline short muls16_asr10(short a, short b)
386 short r;
387 asm (
388 "muls %[a],%[b] \n"
389 "sts macl,%[r] \n"
390 "shlr8 %[r] \n"
391 "shlr2 %[r] \n"
392 : /* outputs */
393 [r]"=r"(r)
394 : /* inputs */
395 [a]"r"(a),
396 [b]"r"(b)
398 return r;
401 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
402 static inline long muls32_asr26(long a, long b)
404 long r, t1, t2, t3;
405 asm (
406 /* Signed 32bit * 32bit -> 64bit multiplication.
407 Notation: xxab * xxcd, where each letter represents 16 bits.
408 xx is the 64 bit sign extension. */
409 "swap.w %[a],%[t1] \n" /* t1 = ba */
410 "mulu %[t1],%[b] \n" /* a * d */
411 "swap.w %[b],%[t3] \n" /* t3 = dc */
412 "sts macl,%[t2] \n" /* t2 = a * d */
413 "mulu %[t1],%[t3] \n" /* a * c */
414 "sts macl,%[r] \n" /* hi = a * c */
415 "mulu %[a],%[t3] \n" /* b * c */
416 "clrt \n"
417 "sts macl,%[t3] \n" /* t3 = b * c */
418 "addc %[t2],%[t3] \n" /* t3 += t2, carry -> t2 */
419 "movt %[t2] \n"
420 "mulu %[a],%[b] \n" /* b * d */
421 "mov %[t3],%[t1] \n" /* t1t3 = t2t3 << 16 */
422 "xtrct %[t2],%[t1] \n"
423 "shll16 %[t3] \n"
424 "sts macl,%[t2] \n" /* lo = b * d */
425 "clrt \n" /* hi.lo += t1t3 */
426 "addc %[t3],%[t2] \n"
427 "addc %[t1],%[r] \n"
428 "cmp/pz %[a] \n" /* ab >= 0 ? */
429 "bt 1f \n"
430 "sub %[b],%[r] \n" /* no: hi -= cd (sign extension of ab is -1) */
431 "1: \n"
432 "cmp/pz %[b] \n" /* cd >= 0 ? */
433 "bt 2f \n"
434 "sub %[a],%[r] \n" /* no: hi -= ab (sign extension of cd is -1) */
435 "2: \n"
436 /* Shift right by 26 and return low 32 bits */
437 "shll2 %[r] \n" /* hi <<= 6 */
438 "shll2 %[r] \n"
439 "shll2 %[r] \n"
440 "shlr16 %[t2] \n" /* (unsigned)lo >>= 26 */
441 "shlr8 %[t2] \n"
442 "shlr2 %[t2] \n"
443 "or %[t2],%[r] \n" /* combine result */
444 : /* outputs */
445 [r] "=&r"(r),
446 [t1]"=&r"(t1),
447 [t2]"=&r"(t2),
448 [t3]"=&r"(t3)
449 : /* inputs */
450 [a] "r" (a),
451 [b] "r" (b)
453 return r;
456 #elif defined CPU_COLDFIRE
458 #define MULS16_ASR10(a, b) muls16_asr10(a, b)
459 static inline short muls16_asr10(short a, short b)
461 asm (
462 "muls.w %[a],%[b] \n"
463 "asr.l #8,%[b] \n"
464 "asr.l #2,%[b] \n"
465 : /* outputs */
466 [b]"+d"(b)
467 : /* inputs */
468 [a]"d" (a)
470 return b;
473 /* Needs the EMAC initialised to fractional mode w/o rounding and saturation */
474 #define MULS32_INIT() coldfire_set_macsr(EMAC_FRACTIONAL)
475 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
476 static inline long muls32_asr26(long a, long b)
478 long r, t1;
479 asm (
480 "mac.l %[a], %[b], %%acc0 \n" /* multiply */
481 "move.l %%accext01, %[t1] \n" /* get low part */
482 "movclr.l %%acc0, %[r] \n" /* get high part */
483 "asl.l #5, %[r] \n" /* hi <<= 5, plus one free */
484 "lsr.l #3, %[t1] \n" /* lo >>= 3 */
485 "and.l #0x1f, %[t1] \n" /* mask out unrelated bits */
486 "or.l %[t1], %[r] \n" /* combine result */
487 : /* outputs */
488 [r] "=d"(r),
489 [t1]"=d"(t1)
490 : /* inputs */
491 [a] "d" (a),
492 [b] "d" (b)
494 return r;
497 #elif defined CPU_ARM
499 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
500 static inline long muls32_asr26(long a, long b)
502 long r, t1;
503 asm (
504 "smull %[r], %[t1], %[a], %[b] \n"
505 "mov %[r], %[r], lsr #26 \n"
506 "orr %[r], %[r], %[t1], lsl #6 \n"
507 : /* outputs */
508 [r] "=&r,&r,&r"(r),
509 [t1]"=&r,&r,&r"(t1)
510 : /* inputs */
511 [a] "%r,%r,%r" (a),
512 [b] "r,0,1" (b)
514 return r;
517 #endif /* CPU */
519 /* default macros */
520 #ifndef MULS16_ASR10
521 #define MULS16_ASR10(a, b) ((short)(((long)(a) * (long)(b)) >> 10))
522 #endif
523 #ifndef MULS32_ASR26
524 #define MULS32_ASR26(a, b) ((long)(((long long)(a) * (long long)(b)) >> 26))
525 #endif
526 #ifndef MULS32_INIT
527 #define MULS32_INIT()
528 #endif
530 int ilog2_fp(long value) /* calculate integer log2(value_fp_6.26) */
532 int i = 0;
534 if (value <= 0) {
535 return -32767;
536 } else if (value > (1L<<26)) {
537 while (value >= (2L<<26)) {
538 value >>= 1;
539 i++;
541 } else {
542 while (value < (1L<<26)) {
543 value <<= 1;
544 i--;
547 return i;
550 void recalc_parameters(void)
552 x_step = (x_max - x_min) / LCD_WIDTH;
553 x_delta = (x_step * LCD_WIDTH) / 8;
554 y_step = (y_max - y_min) / LCD_HEIGHT;
555 y_delta = (y_step * LCD_HEIGHT) / 8;
556 step_log2 = ilog2_fp(MIN(x_step, y_step));
557 max_iter = MAX(15, -15 * step_log2 - 45);
560 #if CONFIG_LCD == LCD_SSD1815
561 /* Recorder, Ondio: pixel_height == 1.25 * pixel_width */
562 #define MB_HEIGHT (LCD_HEIGHT*5/4)
563 #else
564 /* square pixels */
565 #define MB_HEIGHT LCD_HEIGHT
566 #endif
568 #define MB_XOFS (-0x03000000L) /* -0.75 (s5.26) */
569 #if 3000*MB_HEIGHT/LCD_WIDTH >= 2400 /* width is limiting factor */
570 #define MB_XFAC (0x06000000LL) /* 1.5 (s5.26) */
571 #define MB_YFAC (MB_XFAC*MB_HEIGHT/LCD_WIDTH)
572 #else /* height is limiting factor */
573 #define MB_YFAC (0x04cccccdLL) /* 1.2 (s5.26) */
574 #define MB_XFAC (MB_YFAC*LCD_WIDTH/MB_HEIGHT)
575 #endif
577 void init_mandelbrot_set(void)
579 x_min = MB_XOFS-MB_XFAC;
580 x_max = MB_XOFS+MB_XFAC;
581 y_min = -MB_YFAC;
582 y_max = MB_YFAC;
583 recalc_parameters();
586 void calc_mandelbrot_low_prec(void)
588 long start_tick, last_yield;
589 #ifndef USEGSLIB
590 long next_update = *rb->current_tick;
591 int last_px = px_min;
592 #endif
593 unsigned n_iter;
594 long a32, b32;
595 short x, x2, y, y2, a, b;
596 int p_x, p_y;
598 start_tick = last_yield = *rb->current_tick;
600 for (p_x = 0, a32 = x_min; p_x < px_max; p_x++, a32 += x_step) {
601 if (p_x < px_min)
602 continue;
603 a = a32 >> 16;
604 for (p_y = LCD_HEIGHT-1, b32 = y_min; p_y >= py_min; p_y--, b32 += y_step) {
605 if (p_y >= py_max)
606 continue;
607 b = b32 >> 16;
608 x = a;
609 y = b;
610 n_iter = 0;
612 while (++n_iter <= max_iter) {
613 x2 = MULS16_ASR10(x, x);
614 y2 = MULS16_ASR10(y, y);
616 if (x2 + y2 > (4<<10)) break;
618 y = 2 * MULS16_ASR10(x, y) + b;
619 x = x2 - y2 + a;
622 if (n_iter > max_iter)
623 imgbuffer[p_y] = color[8];
624 else
625 imgbuffer[p_y] = color[n_iter & 7];
627 /* be nice to other threads:
628 * if at least one tick has passed, yield */
629 if (*rb->current_tick > last_yield) {
630 rb->yield();
631 last_yield = *rb->current_tick;
634 #ifdef USEGSLIB
635 grey_ub_gray_bitmap_part(imgbuffer, 0, py_min, 1,
636 p_x, py_min, 1, py_max - py_min);
637 #else
638 rb->lcd_bitmap_part(imgbuffer, 0, py_min, 1,
639 p_x, py_min, 1, py_max - py_min);
640 if ((p_x == px_max - 1) || TIME_AFTER(*rb->current_tick, next_update))
642 next_update = *rb->current_tick + UPDATE_FREQ;
643 rb->lcd_update_rect(last_px, py_min, p_x - last_px + 1,
644 py_max - py_min);
645 last_px = p_x;
647 #endif
651 void calc_mandelbrot_high_prec(void)
653 long start_tick, last_yield;
654 #ifndef USEGSLIB
655 long next_update = *rb->current_tick;
656 int last_px = px_min;
657 #endif
658 unsigned n_iter;
659 long x, x2, y, y2, a, b;
660 int p_x, p_y;
662 MULS32_INIT();
663 start_tick = last_yield = *rb->current_tick;
665 for (p_x = 0, a = x_min; p_x < px_max; p_x++, a += x_step) {
666 if (p_x < px_min)
667 continue;
668 for (p_y = LCD_HEIGHT-1, b = y_min; p_y >= py_min; p_y--, b += y_step) {
669 if (p_y >= py_max)
670 continue;
671 x = a;
672 y = b;
673 n_iter = 0;
675 while (++n_iter <= max_iter) {
676 x2 = MULS32_ASR26(x, x);
677 y2 = MULS32_ASR26(y, y);
679 if (x2 + y2 > (4L<<26)) break;
681 y = 2 * MULS32_ASR26(x, y) + b;
682 x = x2 - y2 + a;
685 if (n_iter > max_iter)
686 imgbuffer[p_y] = color[8];
687 else
688 imgbuffer[p_y] = color[n_iter & 7];
690 /* be nice to other threads:
691 * if at least one tick has passed, yield */
692 if (*rb->current_tick > last_yield) {
693 rb->yield();
694 last_yield = *rb->current_tick;
697 #ifdef USEGSLIB
698 grey_ub_gray_bitmap_part(imgbuffer, 0, py_min, 1,
699 p_x, py_min, 1, py_max - py_min);
700 #else
701 rb->lcd_bitmap_part(imgbuffer, 0, py_min, 1,
702 p_x, py_min, 1, py_max-py_min);
703 if ((p_x == px_max - 1) || TIME_AFTER(*rb->current_tick, next_update))
705 next_update = *rb->current_tick + UPDATE_FREQ;
706 rb->lcd_update_rect(last_px, py_min, p_x - last_px + 1,
707 py_max - py_min);
708 last_px = p_x;
710 #endif
714 void cleanup(void *parameter)
716 (void)parameter;
717 #ifdef USEGSLIB
718 grey_release();
719 #endif
722 #define REDRAW_NONE 0
723 #define REDRAW_PARTIAL 1
724 #define REDRAW_FULL 2
726 enum plugin_status plugin_start(const void* parameter)
728 int button;
729 int lastbutton = BUTTON_NONE;
730 int redraw = REDRAW_FULL;
732 (void)parameter;
734 #ifdef USEGSLIB
735 /* get the remainder of the plugin buffer */
736 gbuf = (unsigned char *) rb->plugin_get_buffer(&gbuf_size);
738 /* initialize the greyscale buffer.*/
739 if (!grey_init(gbuf, gbuf_size, GREY_ON_COP,
740 LCD_WIDTH, LCD_HEIGHT, NULL))
742 rb->splash(HZ, "Couldn't init greyscale display");
743 return 0;
745 grey_show(true); /* switch on greyscale overlay */
746 #endif
748 #if LCD_DEPTH > 1
749 rb->lcd_set_backdrop(NULL);
750 #endif
752 init_mandelbrot_set();
754 /* main loop */
755 while (true) {
756 if (redraw > REDRAW_NONE) {
757 #ifdef HAVE_ADJUSTABLE_CPU_FREQ
758 rb->cpu_boost(true);
759 #endif
760 if (redraw == REDRAW_FULL) {
761 MYLCD(clear_display)();
762 MYLCD_UPDATE();
765 if (step_log2 <= -10) /* select precision */
766 calc_mandelbrot_high_prec();
767 else
768 calc_mandelbrot_low_prec();
770 #ifdef HAVE_ADJUSTABLE_CPU_FREQ
771 rb->cpu_boost(false);
772 #endif
773 px_min = 0;
774 px_max = LCD_WIDTH;
775 py_min = 0;
776 py_max = LCD_HEIGHT;
777 redraw = REDRAW_NONE;
780 button = rb->button_get(true);
781 switch (button) {
782 #ifdef MANDELBROT_RC_QUIT
783 case MANDELBROT_RC_QUIT:
784 #endif
785 case MANDELBROT_QUIT:
786 #ifdef USEGSLIB
787 grey_release();
788 #endif
789 return PLUGIN_OK;
791 case MANDELBROT_ZOOM_OUT:
792 #ifdef MANDELBROT_ZOOM_OUT_PRE
793 if (lastbutton != MANDELBROT_ZOOM_OUT_PRE)
794 break;
795 #endif
796 x_min -= x_delta;
797 x_max += x_delta;
798 y_min -= y_delta;
799 y_max += y_delta;
800 recalc_parameters();
801 redraw = REDRAW_FULL;
802 break;
805 case MANDELBROT_ZOOM_IN:
806 #ifdef MANDELBROT_ZOOM_IN_PRE
807 if (lastbutton != MANDELBROT_ZOOM_IN_PRE)
808 break;
809 #endif
810 #ifdef MANDELBROT_ZOOM_IN2
811 case MANDELBROT_ZOOM_IN2:
812 #endif
813 x_min += x_delta;
814 x_max -= x_delta;
815 y_min += y_delta;
816 y_max -= y_delta;
817 recalc_parameters();
818 redraw = REDRAW_FULL;
819 break;
821 case MANDELBROT_UP:
822 y_min += y_delta;
823 y_max += y_delta;
824 MYXLCD(scroll_down)(LCD_HEIGHT/8);
825 MYLCD_UPDATE();
826 py_max = (LCD_HEIGHT/8);
827 redraw = REDRAW_PARTIAL;
828 break;
830 case MANDELBROT_DOWN:
831 y_min -= y_delta;
832 y_max -= y_delta;
833 MYXLCD(scroll_up)(LCD_HEIGHT/8);
834 MYLCD_UPDATE();
835 py_min = (LCD_HEIGHT-LCD_HEIGHT/8);
836 redraw = REDRAW_PARTIAL;
837 break;
839 case MANDELBROT_LEFT:
840 x_min -= x_delta;
841 x_max -= x_delta;
842 MYXLCD(scroll_right)(LCD_WIDTH/8);
843 MYLCD_UPDATE();
844 px_max = (LCD_WIDTH/8);
845 redraw = REDRAW_PARTIAL;
846 break;
848 case MANDELBROT_RIGHT:
849 x_min += x_delta;
850 x_max += x_delta;
851 MYXLCD(scroll_left)(LCD_WIDTH/8);
852 MYLCD_UPDATE();
853 px_min = (LCD_WIDTH-LCD_WIDTH/8);
854 redraw = REDRAW_PARTIAL;
855 break;
857 case MANDELBROT_MAXITER_DEC:
858 #ifdef MANDELBROT_MAXITER_DEC_PRE
859 if (lastbutton != MANDELBROT_MAXITER_DEC_PRE)
860 break;
861 #endif
862 if (max_iter >= 15) {
863 max_iter -= max_iter / 3;
864 redraw = REDRAW_FULL;
866 break;
868 case MANDELBROT_MAXITER_INC:
869 #ifdef MANDELBROT_MAXITER_INC_PRE
870 if (lastbutton != MANDELBROT_MAXITER_INC_PRE)
871 break;
872 #endif
873 max_iter += max_iter / 2;
874 redraw = REDRAW_FULL;
875 break;
877 case MANDELBROT_RESET:
878 init_mandelbrot_set();
879 redraw = REDRAW_FULL;
880 break;
882 default:
883 if (rb->default_event_handler_ex(button, cleanup, NULL)
884 == SYS_USB_CONNECTED)
885 return PLUGIN_USB_CONNECTED;
886 break;
888 if (button != BUTTON_NONE)
889 lastbutton = button;
891 #ifdef USEGSLIB
892 grey_release();
893 #endif
894 return PLUGIN_OK;
896 #endif