跳至內容

分形/軟體

Add links
來自華夏公益教科書

分形軟體

[編輯 | 編輯原始碼]

Imagina

[編輯 | 編輯原始碼]

精度

  • 帶有單精度浮點數的 floatexp
  • 帶有雙精度浮點數的 FloatExp

控制

[編輯 | 編輯原始碼]
  • 滑鼠左鍵
    • 拖動 - 平移檢視
    • 雙擊 - 查詢附近的特徵並開啟特徵查詢對話方塊
  • 滑鼠滾輪 - 縮放
  • A - 增加顏色密度
  • S - 減少顏色密度
  • D - 增加迭代限制
  • F - 減少迭代限制
  • E - 顏色向前迴圈
  • R - 顏色向後迴圈

Nanoscope 由 pauldelbrot

[編輯 | 編輯原始碼]
Extra range is needed on a few specific iterations, if zooming past e300 and doing flybys past minibrots below e300.

Nanoscope stores a similar array to KF, 3x double with re, im, and a mag value used for glitch detection. It also stores a pointer that is either null or points to a wide-exponent copy of those three values. If the double precision values underflow (denorm or zero) these get set for that iteration during reference orbit computation. Otherwise the pointer is null. During iteration, if the pointer is not null the next iteration is done using 52-bit mantissa wide-exponent calculations (significantly slower, but for only one iteration).

The circumstances that trigger this are instances like this. Say there's a period-73174 mini at about e380. If the zoom goes very close to that mini, then for images near and deeper than that mini, every 73174th reference orbit iteration is within 1e-380 of zero or so, so it underflows in a double. So wide exponent calculations must be done on those iterations, and those reference orbit entries must be stored with a wide exponent.

For every other iteration below e300, Nanoscope does the calculations with rescaled values that don't need a wide exponent. The rescaling is changed every 500 or so iterations -- this works because until a point is on the verge of escape, its orbit dynamics are dominated by the effects of repelling points in the Julia set nearby, and that typically means its magnitude doubles each iteration. The exponent width of doubles is 11 bits, so from about -1000 to about 1000, representing powers of 2 (not 10), so easily accommodates 500 doublings with plenty of margin for error. Re-rescaling is also done after every iteration that needed a wide-exponent reference point, because the orbit has jumped much closer to 0 again on such iterations.

On "final approach" an escape-bound critical orbit moves faster, with the magnitude squaring each iteration, but by the time this happens the unscaled magnitude is above the e-300 threshold and Nanoscope has switched to bog-standard perturbation calculations without any rescaling or other sneaky tricks. And escape is usually within the next 500 iterations anyway.[1]

NanoMB

[編輯 | 編輯原始碼]
time ./nanomb64_i7.exe --kfr ./test_x1.kfr --maxiters 250000 --period 2489  --width 500 --height 500 --orderM 4 --orderN 4 --output_ppm s.ppm --output_kfb s.kfb --force_type 'floatexp'

週期為 2489

nanomb.exe --kfr test_x1.kfr --period 2449 --width 500 --height 500 --orderM 4 --orderN 4 --output_ppm s.ppm --force_type floatexp


SuperMB

[編輯 | 編輯原始碼]
  • SuperMB
  • superMB 實際上並不是超級的。它只是對 Claude 在 superfractalthing 主題中釋出的程式碼進行的修改。我用它來進行實驗。它現在有一個基本的 GUI。我已經在下面附加了原始碼 + 依賴項


Knighty's SMB which I think is still a bit faster than KF though being more of a testbed than a usable renderer puts the glitches in distinct sets (G1,.., Gn) with same iteration number where glitch was detected. Next references will then be 1 random pixel from each of the G1,..,Gn and is used only to recalculate the pixels in each set G. Secondary glitches simply generate another set G and you just put then in some queue or stack and keep going at it till no more G sets left. When dealing with glitched pixels you can use the same series expansion as a starting point. ( Gerrit)[2]


在 GitHub 上

[編輯 | 編輯原始碼]

中心

[編輯 | 編輯原始碼]
  1. ↑ fractalforums.org：memory-bandwidth-trade-offs-for-perturbation-rendering
  2. ↑ fractalforums.org : how-to-get-second-reference-when-using-perturbation-theory
檢索自 "https://wikibook.tw/w/index.php?title=Fractals/software&oldid=4327085"
華夏公益教科書