apps/plugins/mandelbrot.c

   1 /***************************************************************************
   2  *             __________               __   ___.
   3  *   Open      \______   \ ____   ____ |  | _\_ |__   _______  ___
   4  *   Source     |       _//  _ \_/ ___\|  |/ /| __ \ /  _ \  \/  /
   5  *   Jukebox    |    |   (  <_> )  \___|    < | \_\ (  <_> > <  <
   6  *   Firmware   |____|_  /\____/ \___  >__|_ \|___  /\____/__/\_ \
   7  *                     \/            \/     \/    \/            \/
   8  * $Id$
   9  *
  10  * Copyright (C) 2004 Matthias Wientapper
  11  * Heavily extended 2005 Jens Arnold
  12  *
  13  *
  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.
  18  *
  19  * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY OF ANY
  20  * KIND, either express or implied.
  21  *
  22  ****************************************************************************/
  23 #include "plugin.h"
  24
  25 #ifdef HAVE_LCD_BITMAP
  26 #include "lib/grey.h"
  27 #include "lib/xlcd.h"
  28
  29 PLUGIN_HEADER
  30
  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
  43
  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
  55
  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)
  69
  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
  82
  83 #define MANDELBROT_RC_QUIT BUTTON_RC_STOP
  84
  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)
  98
  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
 114
 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
 128
 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
 140
 141 #elif CONFIG_KEYPAD == SANSA_FUZE_PAD
 142 #define MANDELBROT_QUIT         BUTTON_POWER
 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)
 152
 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
 164
 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
 176
 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)
 188
 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)
 202
 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
 216
 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
 228
 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
 242
 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
 258
 259 #elif CONFIG_KEYPAD == COWOND2_PAD
 260 #define MANDELBROT_QUIT          BUTTON_POWER
 261
 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
 273
 274 #else
 275 #error No keymap defined!
 276 #endif
 277
 278 #ifdef HAVE_TOUCHSCREEN
 279 #ifndef MANDELBROT_UP
 280 #define MANDELBROT_UP            BUTTON_TOPMIDDLE
 281 #endif
 282 #ifndef MANDELBROT_DOWN
 283 #define MANDELBROT_DOWN          BUTTON_BOTTOMMIDDLE
 284 #endif
 285 #ifndef MANDELBROT_LEFT
 286 #define MANDELBROT_LEFT          BUTTON_MIDLEFT
 287 #endif
 288 #ifndef MANDELBROT_RIGHT
 289 #define MANDELBROT_RIGHT         BUTTON_MIDRIGHT
 290 #endif
 291 #ifndef MANDELBROT_ZOOM_IN_PRE
 292 #define MANDELBROT_ZOOM_IN_PRE   BUTTON_TOPRIGHT
 293 #endif
 294 #ifndef MANDELBROT_ZOOM_IN
 295 #define MANDELBROT_ZOOM_IN      (BUTTON_TOPRIGHT | BUTTON_REL)
 296 #endif
 297 #ifndef MANDELBROT_ZOOM_OUT_PRE
 298 #define MANDELBROT_ZOOM_OUT_PRE  BUTTON_TOPLEFT
 299 #endif
 300 #ifndef MANDELBROT_ZOOM_OUT
 301 #define MANDELBROT_ZOOM_OUT     (BUTTON_TOPLEFT | BUTTON_REL)
 302 #endif
 303 #ifndef MANDELBROT_MAXITER_INC
 304 #define MANDELBROT_MAXITER_INC   BUTTON_BOTTOMRIGHT
 305 #endif
 306 #ifndef MANDELBROT_MAXITER_DEC
 307 #define MANDELBROT_MAXITER_DEC   BUTTON_BOTTOMLEFT
 308 #endif
 309 #ifndef MANDELBROT_RESET
 310 #define MANDELBROT_RESET         BUTTON_CENTER
 311 #endif
 312 #endif
 313
 314 #if LCD_DEPTH < 8
 315 #define USEGSLIB
 316 #define MYLCD(fn) grey_ub_ ## fn
 317 #define MYLCD_UPDATE()
 318 #define MYXLCD(fn) grey_ub_ ## fn
 319 #else
 320 #define UPDATE_FREQ (HZ/50)
 321 #define MYLCD(fn) rb->lcd_ ## fn
 322 #define MYLCD_UPDATE() rb->lcd_update();
 323 #define MYXLCD(fn) xlcd_ ## fn
 324 #endif
 325
 326 static const struct plugin_api* rb;
 327
 328 /* Fixed point format s5.26: sign, 5 bits integer part, 26 bits fractional part */
 329 static long x_min;
 330 static long x_max;
 331 static long x_step;
 332 static long x_delta;
 333 static long y_min;
 334 static long y_max;
 335 static long y_step;
 336 static long y_delta;
 337
 338 static int px_min = 0;
 339 static int px_max = LCD_WIDTH;
 340 static int py_min = 0;
 341 static int py_max = LCD_HEIGHT;
 342
 343 static int step_log2;
 344 static unsigned max_iter;
 345
 346 #ifdef USEGSLIB
 347 GREY_INFO_STRUCT
 348 static unsigned char *gbuf;
 349 static size_t  gbuf_size = 0;
 350 static unsigned char imgbuffer[LCD_HEIGHT];
 351 #else
 352 static fb_data imgbuffer[LCD_HEIGHT];
 353 #endif
 354
 355 /* 8 entries cyclical, last entry is black (convergence) */
 356 #ifdef HAVE_LCD_COLOR
 357 static const fb_data color[9] = {
 358     LCD_RGBPACK(255, 0, 159), LCD_RGBPACK(159, 0, 255), LCD_RGBPACK(0, 0, 255),
 359     LCD_RGBPACK(0, 159, 255), LCD_RGBPACK(0, 255, 128), LCD_RGBPACK(128, 255, 0),
 360     LCD_RGBPACK(255, 191, 0), LCD_RGBPACK(255, 0, 0),   LCD_RGBPACK(0, 0, 0)
 361 };
 362 #else /* greyscale */
 363 static const unsigned char color[9] = {
 364     255, 223, 191, 159, 128, 96, 64, 32, 0
 365 };
 366 #endif
 367
 368 #if CONFIG_CPU == SH7034
 369
 370 #define MULS16_ASR10(a, b) muls16_asr10(a, b)
 371 static inline short muls16_asr10(short a, short b)
 372 {
 373     short r;
 374     asm (
 375         "muls    %[a],%[b]   \n"
 376         "sts     macl,%[r]   \n"
 377         "shlr8   %[r]        \n"
 378         "shlr2   %[r]        \n"
 379         : /* outputs */
 380         [r]"=r"(r)
 381         : /* inputs */
 382         [a]"r"(a),
 383         [b]"r"(b)
 384     );
 385     return r;
 386 }
 387
 388 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
 389 static inline long muls32_asr26(long a, long b)
 390 {
 391     long r, t1, t2, t3;
 392     asm (
 393         /* Signed 32bit * 32bit -> 64bit multiplication.
 394            Notation: xxab * xxcd, where each letter represents 16 bits.
 395            xx is the 64 bit sign extension.  */
 396         "swap.w  %[a],%[t1]  \n" /* t1 = ba */
 397         "mulu    %[t1],%[b]  \n" /* a * d */
 398         "swap.w  %[b],%[t3]  \n" /* t3 = dc */
 399         "sts     macl,%[t2]  \n" /* t2 = a * d */
 400         "mulu    %[t1],%[t3] \n" /* a * c */
 401         "sts     macl,%[r]   \n" /* hi = a * c */
 402         "mulu    %[a],%[t3]  \n" /* b * c */
 403         "clrt                \n"
 404         "sts     macl,%[t3]  \n" /* t3 = b * c */
 405         "addc    %[t2],%[t3] \n" /* t3 += t2, carry -> t2 */
 406         "movt    %[t2]       \n"
 407         "mulu    %[a],%[b]   \n" /* b * d */
 408         "mov     %[t3],%[t1] \n" /* t2t3 <<= 16 */
 409         "xtrct   %[t2],%[t1] \n"
 410         "mov     %[t1],%[t2] \n"
 411         "shll16  %[t3]       \n"
 412         "sts     macl,%[t1]  \n" /* lo = b * d */
 413         "clrt                \n" /* hi.lo += t2t3 */
 414         "addc    %[t3],%[t1] \n"
 415         "addc    %[t2],%[r]  \n"
 416         "cmp/pz  %[a]        \n" /* ab >= 0 ? */
 417         "bt      1f          \n"
 418         "sub     %[b],%[r]   \n" /* no: hi -= cd (sign extension of ab is -1) */
 419     "1:                      \n"
 420         "cmp/pz  %[b]        \n" /* cd >= 0 ? */
 421         "bt      2f          \n"
 422         "sub     %[a],%[r]   \n" /* no: hi -= ab (sign extension of cd is -1) */
 423     "2:                      \n"
 424         /* Shift right by 26 and return low 32 bits */
 425         "shll2   %[r]        \n" /* hi <<= 6 */
 426         "shll2   %[r]        \n"
 427         "shll2   %[r]        \n"
 428         "shlr16  %[t1]       \n" /* (unsigned)lo >>= 26 */
 429         "shlr8   %[t1]       \n"
 430         "shlr2   %[t1]       \n"
 431         "or      %[t1],%[r]  \n" /* combine result */
 432         : /* outputs */
 433         [r] "=&r"(r),
 434         [t1]"=&r"(t1),
 435         [t2]"=&r"(t2),
 436         [t3]"=&r"(t3)
 437         : /* inputs */
 438         [a] "r"  (a),
 439         [b] "r"  (b)
 440     );
 441     return r;
 442 }
 443
 444 #elif defined CPU_COLDFIRE
 445
 446 #define MULS16_ASR10(a, b) muls16_asr10(a, b)
 447 static inline short muls16_asr10(short a, short b)
 448 {
 449     asm (
 450         "muls.w  %[a],%[b]   \n"
 451         "asr.l   #8,%[b]     \n"
 452         "asr.l   #2,%[b]     \n"
 453         : /* outputs */
 454         [b]"+d"(b)
 455         : /* inputs */
 456         [a]"d" (a)
 457     );
 458     return b;
 459 }
 460
 461 /* Needs the EMAC initialised to fractional mode w/o rounding and saturation */
 462 #define MULS32_INIT() coldfire_set_macsr(EMAC_FRACTIONAL)
 463 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
 464 static inline long muls32_asr26(long a, long b)
 465 {
 466     long r, t1;
 467     asm (
 468         "mac.l   %[a], %[b], %%acc0  \n" /* multiply */
 469         "move.l  %%accext01, %[t1]   \n" /* get low part */
 470         "movclr.l %%acc0, %[r]       \n" /* get high part */
 471         "asl.l   #5, %[r]            \n" /* hi <<= 5, plus one free */
 472         "lsr.l   #3, %[t1]           \n" /* lo >>= 3 */
 473         "and.l   #0x1f, %[t1]        \n" /* mask out unrelated bits */
 474         "or.l    %[t1], %[r]         \n" /* combine result */
 475         : /* outputs */
 476         [r] "=d"(r),
 477         [t1]"=d"(t1)
 478         : /* inputs */
 479         [a] "d" (a),
 480         [b] "d" (b)
 481     );
 482     return r;
 483 }
 484
 485 #elif defined CPU_ARM
 486
 487 #define MULS32_ASR26(a, b) muls32_asr26(a, b)
 488 static inline long muls32_asr26(long a, long b)
 489 {
 490     long r, t1;
 491     asm (
 492         "smull   %[r], %[t1], %[a], %[b]     \n"
 493         "mov     %[r], %[r], lsr #26         \n"
 494         "orr     %[r], %[r], %[t1], lsl #6   \n"
 495         : /* outputs */
 496         [r] "=&r,&r,&r"(r),
 497         [t1]"=&r,&r,&r"(t1)
 498         : /* inputs */
 499         [a] "%r,%r,%r" (a),
 500         [b] "r,0,1" (b)
 501     );
 502     return r;
 503 }
 504
 505 #endif /* CPU */
 506
 507 /* default macros */
 508 #ifndef MULS16_ASR10
 509 #define MULS16_ASR10(a, b) ((short)(((long)(a) * (long)(b)) >> 10))
 510 #endif
 511 #ifndef MULS32_ASR26
 512 #define MULS32_ASR26(a, b) ((long)(((long long)(a) * (long long)(b)) >> 26))
 513 #endif
 514 #ifndef MULS32_INIT
 515 #define MULS32_INIT()
 516 #endif
 517
 518 int ilog2_fp(long value) /* calculate integer log2(value_fp_6.26) */
 519 {
 520     int i = 0;
 521
 522     if (value <= 0) {
 523         return -32767;
 524     } else if (value > (1L<<26)) {
 525         while (value >= (2L<<26)) {
 526             value >>= 1;
 527             i++;
 528         }
 529     } else {
 530         while (value < (1L<<26)) {
 531             value <<= 1;
 532             i--;
 533         }
 534     }
 535     return i;
 536 }
 537
 538 void recalc_parameters(void)
 539 {
 540     x_step = (x_max - x_min) / LCD_WIDTH;
 541     x_delta = (x_step * LCD_WIDTH) / 8;
 542     y_step = (y_max - y_min) / LCD_HEIGHT;
 543     y_delta = (y_step * LCD_HEIGHT) / 8;
 544     step_log2 = ilog2_fp(MIN(x_step, y_step));
 545     max_iter = MAX(15, -15 * step_log2 - 45);
 546 }
 547
 548 #if CONFIG_LCD == LCD_SSD1815
 549 /* Recorder, Ondio: pixel_height == 1.25 * pixel_width */
 550 #define MB_HEIGHT (LCD_HEIGHT*5/4)
 551 #else
 552 /* square pixels */
 553 #define MB_HEIGHT LCD_HEIGHT
 554 #endif
 555
 556 #define MB_XOFS (-0x03000000L)       /* -0.75 (s5.26) */
 557 #if 3000*MB_HEIGHT/LCD_WIDTH >= 2400 /* width is limiting factor */
 558 #define MB_XFAC (0x06000000LL)       /* 1.5 (s5.26) */
 559 #define MB_YFAC (MB_XFAC*MB_HEIGHT/LCD_WIDTH)
 560 #else                                /* height is limiting factor */
 561 #define MB_YFAC (0x04cccccdLL)       /* 1.2 (s5.26) */
 562 #define MB_XFAC (MB_YFAC*LCD_WIDTH/MB_HEIGHT)
 563 #endif
 564
 565 void init_mandelbrot_set(void)
 566 {
 567     x_min = MB_XOFS-MB_XFAC;
 568     x_max = MB_XOFS+MB_XFAC;
 569     y_min = -MB_YFAC;
 570     y_max = MB_YFAC;
 571     recalc_parameters();
 572 }
 573
 574 void calc_mandelbrot_low_prec(void)
 575 {
 576     long start_tick, last_yield;
 577 #ifndef USEGSLIB
 578     long next_update = *rb->current_tick;
 579     int last_px = px_min;
 580 #endif
 581     unsigned n_iter;
 582     long a32, b32;
 583     short x, x2, y, y2, a, b;
 584     int p_x, p_y;
 585
 586     start_tick = last_yield = *rb->current_tick;
 587
 588     for (p_x = 0, a32 = x_min; p_x < px_max; p_x++, a32 += x_step) {
 589         if (p_x < px_min)
 590             continue;
 591         a = a32 >> 16;
 592         for (p_y = LCD_HEIGHT-1, b32 = y_min; p_y >= py_min; p_y--, b32 += y_step) {
 593             if (p_y >= py_max)
 594                 continue;
 595             b = b32 >> 16;
 596             x = a;
 597             y = b;
 598             n_iter = 0;
 599
 600             while (++n_iter <= max_iter) {
 601                 x2 = MULS16_ASR10(x, x);
 602                 y2 = MULS16_ASR10(y, y);
 603
 604                 if (x2 + y2 > (4<<10)) break;
 605
 606                 y = 2 * MULS16_ASR10(x, y) + b;
 607                 x = x2 - y2 + a;
 608             }
 609
 610             if  (n_iter > max_iter)
 611                 imgbuffer[p_y] = color[8];
 612             else
 613                 imgbuffer[p_y] = color[n_iter & 7];
 614
 615             /* be nice to other threads:
 616              * if at least one tick has passed, yield */
 617             if  (*rb->current_tick > last_yield) {
 618                 rb->yield();
 619                 last_yield = *rb->current_tick;
 620             }
 621         }
 622 #ifdef USEGSLIB
 623         grey_ub_gray_bitmap_part(imgbuffer, 0, py_min, 1,
 624                                  p_x, py_min, 1, py_max - py_min);
 625 #else
 626         rb->lcd_bitmap_part(imgbuffer, 0, py_min, 1,
 627                             p_x, py_min, 1, py_max - py_min);
 628         if ((p_x == px_max - 1) || TIME_AFTER(*rb->current_tick, next_update))
 629         {
 630             next_update = *rb->current_tick + UPDATE_FREQ;
 631             rb->lcd_update_rect(last_px, py_min, p_x - last_px + 1,
 632                                 py_max - py_min);
 633             last_px = p_x;
 634         }
 635 #endif
 636     }
 637 }
 638
 639 void calc_mandelbrot_high_prec(void)
 640 {
 641     long start_tick, last_yield;
 642 #ifndef USEGSLIB
 643     long next_update = *rb->current_tick;
 644     int last_px = px_min;
 645 #endif
 646     unsigned n_iter;
 647     long x, x2, y, y2, a, b;
 648     int p_x, p_y;
 649
 650     MULS32_INIT();
 651     start_tick = last_yield = *rb->current_tick;
 652
 653     for (p_x = 0, a = x_min; p_x < px_max; p_x++, a += x_step) {
 654         if (p_x < px_min)
 655             continue;
 656         for (p_y = LCD_HEIGHT-1, b = y_min; p_y >= py_min; p_y--, b += y_step) {
 657             if (p_y >= py_max)
 658                 continue;
 659             x = a;
 660             y = b;
 661             n_iter = 0;
 662
 663             while (++n_iter <= max_iter) {
 664                 x2 = MULS32_ASR26(x, x);
 665                 y2 = MULS32_ASR26(y, y);
 666
 667                 if (x2 + y2 > (4L<<26)) break;
 668
 669                 y = 2 * MULS32_ASR26(x, y) + b;
 670                 x = x2 - y2 + a;
 671             }
 672
 673             if  (n_iter > max_iter)
 674                 imgbuffer[p_y] = color[8];
 675             else
 676                 imgbuffer[p_y] = color[n_iter & 7];
 677
 678             /* be nice to other threads:
 679              * if at least one tick has passed, yield */
 680             if  (*rb->current_tick > last_yield) {
 681                 rb->yield();
 682                 last_yield = *rb->current_tick;
 683             }
 684         }
 685 #ifdef USEGSLIB
 686         grey_ub_gray_bitmap_part(imgbuffer, 0, py_min, 1,
 687                                  p_x, py_min, 1, py_max - py_min);
 688 #else
 689         rb->lcd_bitmap_part(imgbuffer, 0, py_min, 1,
 690                             p_x, py_min, 1, py_max-py_min);
 691         if ((p_x == px_max - 1) || TIME_AFTER(*rb->current_tick, next_update))
 692         {
 693             next_update = *rb->current_tick + UPDATE_FREQ;
 694             rb->lcd_update_rect(last_px, py_min, p_x - last_px + 1,
 695                                 py_max - py_min);
 696             last_px = p_x;
 697         }
 698 #endif
 699     }
 700 }
 701
 702 void cleanup(void *parameter)
 703 {
 704     (void)parameter;
 705 #ifdef USEGSLIB
 706     grey_release();
 707 #endif
 708 }
 709
 710 #define REDRAW_NONE    0
 711 #define REDRAW_PARTIAL 1
 712 #define REDRAW_FULL    2
 713
 714 enum plugin_status plugin_start(const struct plugin_api* api, const void* parameter)
 715 {
 716     int button;
 717     int lastbutton = BUTTON_NONE;
 718     int redraw = REDRAW_FULL;
 719
 720     rb = api;
 721     (void)parameter;
 722
 723 #ifdef USEGSLIB
 724     /* get the remainder of the plugin buffer */
 725     gbuf = (unsigned char *) rb->plugin_get_buffer(&gbuf_size);
 726
 727     /* initialize the greyscale buffer.*/
 728     if (!grey_init(rb, gbuf, gbuf_size, GREY_ON_COP,
 729                    LCD_WIDTH, LCD_HEIGHT, NULL))
 730     {
 731         rb->splash(HZ, "Couldn't init greyscale display");
 732         return 0;
 733     }
 734     grey_show(true); /* switch on greyscale overlay */
 735 #else
 736     xlcd_init(rb);
 737 #endif
 738
 739 #if LCD_DEPTH > 1
 740     rb->lcd_set_backdrop(NULL);
 741 #endif
 742
 743     init_mandelbrot_set();
 744
 745     /* main loop */
 746     while (true) {
 747         if (redraw > REDRAW_NONE) {
 748 #ifdef HAVE_ADJUSTABLE_CPU_FREQ
 749             rb->cpu_boost(true);
 750 #endif
 751             if (redraw == REDRAW_FULL) {
 752                 MYLCD(clear_display)();
 753                 MYLCD_UPDATE();
 754             }
 755
 756             if (step_log2 <= -10) /* select precision */
 757                 calc_mandelbrot_high_prec();
 758             else
 759                 calc_mandelbrot_low_prec();
 760
 761 #ifdef HAVE_ADJUSTABLE_CPU_FREQ
 762             rb->cpu_boost(false);
 763 #endif
 764             px_min = 0;
 765             px_max = LCD_WIDTH;
 766             py_min = 0;
 767             py_max = LCD_HEIGHT;
 768             redraw = REDRAW_NONE;
 769         }
 770
 771         button = rb->button_get(true);
 772         switch (button) {
 773 #ifdef MANDELBROT_RC_QUIT
 774         case MANDELBROT_RC_QUIT:
 775 #endif
 776         case MANDELBROT_QUIT:
 777 #ifdef USEGSLIB
 778             grey_release();
 779 #endif
 780             return PLUGIN_OK;
 781
 782         case MANDELBROT_ZOOM_OUT:
 783 #ifdef MANDELBROT_ZOOM_OUT_PRE
 784             if (lastbutton != MANDELBROT_ZOOM_OUT_PRE)
 785                 break;
 786 #endif
 787             x_min -= x_delta;
 788             x_max += x_delta;
 789             y_min -= y_delta;
 790             y_max += y_delta;
 791             recalc_parameters();
 792             redraw = REDRAW_FULL;
 793             break;
 794
 795
 796         case MANDELBROT_ZOOM_IN:
 797 #ifdef MANDELBROT_ZOOM_IN_PRE
 798             if (lastbutton != MANDELBROT_ZOOM_IN_PRE)
 799                 break;
 800 #endif
 801 #ifdef MANDELBROT_ZOOM_IN2
 802         case MANDELBROT_ZOOM_IN2:
 803 #endif
 804             x_min += x_delta;
 805             x_max -= x_delta;
 806             y_min += y_delta;
 807             y_max -= y_delta;
 808             recalc_parameters();
 809             redraw = REDRAW_FULL;
 810             break;
 811
 812         case MANDELBROT_UP:
 813             y_min += y_delta;
 814             y_max += y_delta;
 815             MYXLCD(scroll_down)(LCD_HEIGHT/8);
 816             MYLCD_UPDATE();
 817             py_max = (LCD_HEIGHT/8);
 818             redraw = REDRAW_PARTIAL;
 819             break;
 820
 821         case MANDELBROT_DOWN:
 822             y_min -= y_delta;
 823             y_max -= y_delta;
 824             MYXLCD(scroll_up)(LCD_HEIGHT/8);
 825             MYLCD_UPDATE();
 826             py_min = (LCD_HEIGHT-LCD_HEIGHT/8);
 827             redraw = REDRAW_PARTIAL;
 828             break;
 829
 830         case MANDELBROT_LEFT:
 831             x_min -= x_delta;
 832             x_max -= x_delta;
 833             MYXLCD(scroll_right)(LCD_WIDTH/8);
 834             MYLCD_UPDATE();
 835             px_max = (LCD_WIDTH/8);
 836             redraw = REDRAW_PARTIAL;
 837             break;
 838
 839         case MANDELBROT_RIGHT:
 840             x_min += x_delta;
 841             x_max += x_delta;
 842             MYXLCD(scroll_left)(LCD_WIDTH/8);
 843             MYLCD_UPDATE();
 844             px_min = (LCD_WIDTH-LCD_WIDTH/8);
 845             redraw = REDRAW_PARTIAL;
 846             break;
 847
 848         case MANDELBROT_MAXITER_DEC:
 849 #ifdef MANDELBROT_MAXITER_DEC_PRE
 850             if (lastbutton != MANDELBROT_MAXITER_DEC_PRE)
 851                 break;
 852 #endif
 853             if (max_iter >= 15) {
 854                 max_iter -= max_iter / 3;
 855                 redraw = REDRAW_FULL;
 856             }
 857             break;
 858
 859         case MANDELBROT_MAXITER_INC:
 860 #ifdef MANDELBROT_MAXITER_INC_PRE
 861             if (lastbutton != MANDELBROT_MAXITER_INC_PRE)
 862                 break;
 863 #endif
 864             max_iter += max_iter / 2;
 865             redraw = REDRAW_FULL;
 866             break;
 867
 868         case MANDELBROT_RESET:
 869             init_mandelbrot_set();
 870             redraw = REDRAW_FULL;
 871             break;
 872
 873         default:
 874             if (rb->default_event_handler_ex(button, cleanup, NULL)
 875                 == SYS_USB_CONNECTED)
 876                 return PLUGIN_USB_CONNECTED;
 877             break;
 878         }
 879         if (button != BUTTON_NONE)
 880             lastbutton = button;
 881     }
 882 #ifdef USEGSLIB
 883     grey_release();
 884 #endif
 885     return PLUGIN_OK;
 886 }
 887 #endif