HSPでシダを描画する

少し流行に乗り遅れた感がありますが，HSPでシダを描画したので載せます．

元ネタ

「プログラムでシダを描画する」一覧

ソースコード

tytyさんによる「再帰モジュール」が別途必要です．
ダウンロードし，m_rec.hspファイルを同じディレクトリに配置してください．

Nの値を調節すると，いい具合の時間で終わります．

sida.hsp

#include "hsp3dish.as"
#include "m_rec.hsp"

rec_init 100000

#module

#defcfunc W1x double px, double py
return 0.836 * px + 0.044 * py

#defcfunc W1y double px, double py
return -0.044 * px + 0.836 * py + 0.169

#defcfunc W2x double px, double py
return -0.141 * px + 0.302 * py

#defcfunc W2y double px, double py
return 0.302 * px + 0.141 * py + 0.127

#defcfunc W3x double px, double py
return 0.141 * px - 0.302 * py

#defcfunc W3y double px, double py
return 0.302 * px + 0.141 * py + 0.169

#defcfunc W4x double px, double py
return 0.0

#defcfunc W4y double px, double py
return 0.175337 * py

#global

#const N 20

    screen 0, 640, 480

*main

    // 画面の書き換えを止める
    redraw 0

    color 255, 255, 255
    boxf

    pos 10, 10
    color 0, 0, 0
    mes "Drawing..."

    k = N
    x = 0.0
    y = 0.0

    rec_begin k, x, y

        if (0 < k) {
            rec_call k-1, W1x(x, y), W1y(x, y)

            if (rnd(3) == 0) :rec_call k-1, W2x(x, y), W2y(x, y)
            if (rnd(3) == 0) :rec_call k-1, W3x(x, y), W3y(x, y)
            if (rnd(3) == 0) :rec_call k-1, W4x(x, y), W4y(x, y)
        } else {
            color 0, 255, 0
            dx = int(x*ginfo_winx) + ginfo_winx/2
            dy = ginfo_winy - int(y*ginfo_winy)
            boxf dx, dy, dx+1, dy+1
        }
    rec_end

    color 255, 255, 255
    boxf 0, 0, 100, 100
    pos 10, 10
    color 0, 0, 0: mes "Finished"

    // 画面の書き換えを反映
    redraw 1