ruby-dev

<199802190520.OAA31625@picachu.netlab.co.jp> の、
   "[ruby-dev:1495] Re: tkmandel" において、
   "matz@netlab.co.jp (Yukihiro Matsumoto)"さんは書きました：

> |ズーム機能のバグ修正と手直し（ゆっくり深まるように）です：
> 
> このパッチ↓でdevided by zeroがでるようになってしまいました．
> どうしてなんでしょう？
> 
> |-  size = (d - a).abs if (size < (d - a).abs)
> |+  size = (d - b).abs if (size < (d - b).abs)


sizeが0になるような場合に、1に上げてやる必要がありますね。


> 
> さて，数回ズームすると落っこちる件ですが，毎回イメージを生成
> していた部分を流用するようにしたら，とりあえず落ちないように
> なりました．恐らくTkの内部でリークしているのだと思います．
> 
>                                 まつもと ゆきひろ /:|)

す、すごい。すごいきれいに治ってる。

すばらしいです。あと、一応、倍率が少しずつ深まるように
もう一段sqrtをかませたのと、初期値を少し上げてみました。

それから、

#   z = c = Complex(re, im)
#   for i in 0 .. $max_depth
#     z = (z * z) + c
#     break if z.abs2 > 5  ← ココ。
#   end
#   return i

の部分、以前は z.abs2 > 4 相当だったのですが、拡張モジュールになりますが
少し上げて 5 にしてやると、より深まるようです。

ありがとうございました。

M.Kanemitsu

Attachments (1)

tkmandel.rb (4.16 KB, text/x-ruby)

# require "complex"
require "mandel"
require "tkclass"
require "GD"

DefaultMaxDepth = 50
DefaultSX = -2.25
DefaultSY = 1.75
DefaultEX = 1.25
DefaultEY = -1.75

def reset
  $max_depth = DefaultMaxDepth
  $s_re = DefaultSX
  $s_im = DefaultSY
  $e_re = DefaultEX
  $e_im = DefaultEY
  $dx = ($e_re - $s_re).abs / Width
  $dy = ($e_im - $s_im).abs / Height
end


Width = 400
Height = 400

$c = Canvas.new {
       width Width
       height Height
     }
$c.pack

$c_rect = Rectangle.new($c, 0, 0, Width+1, Height+1)
$c_rect.fill "white"



$colors = []

def gd_setup
  $im = GD::Image.new(Width, Height)

  $white = $im.colorAllocate(255, 255, 255)
  $black = $im.colorAllocate(0, 0, 0)

  $im.transparent($white)
  $im.interlace = TRUE

  $color_max = 0
  for i in 0 .. 125
    $colors[$color_max] = $im.colorAllocate(250 - (i*2), i*2, 0)
    $color_max += 1
  end
  for i in 0 .. 125
    $colors[$color_max] = $im.colorAllocate(0, 250 - (i*2), i*2)
    $color_max += 1
  end
  $color_max -= 1
end


def zoom(a, b, c, d)
  center_x = (a + c) / 2
  center_y = (b + d) / 2
  size = (c - a).abs
  size = (d - b).abs if (size < (d - b).abs)
  size = 1 if (size < 1)
  zoom_rate = ((Width + Height) / 2).to_f / size
  $max_depth = ($max_depth.to_f * Math.sqrt(Math.sqrt(Math.sqrt(zoom_rate)))).to_i

  move_x_rate = (center_x - (Width / 2)).to_f / (Width / 2)
  move_y_rate = (center_y - (Height / 2)).to_f / (Height / 2)

  center_re = ($s_re + $e_re) / 2
  center_im = ($s_im + $e_im) / 2
  c_size_re = ($e_re - $s_re).abs
  c_size_im = ($e_im - $s_im).abs

  center_re = center_re + (move_x_rate * (c_size_re / 2))
  center_im = center_im - (move_y_rate * (c_size_im / 2))

  $s_re = center_re - ((c_size_re / 2) / zoom_rate)
  $s_im = center_im + ((c_size_im / 2) / zoom_rate)
  $e_re = center_re + ((c_size_re / 2) / zoom_rate)
  $e_im = center_im - ((c_size_im / 2) / zoom_rate)

  $dx = ($e_re - $s_re).abs / Width
  $dy = ($e_im - $s_im).abs / Height
end


def mandel(x, y)
  re = $s_re + ($dx * x)
  im = $s_im - ($dy * y)
#   z = c = Complex(re, im)
#   for i in 0 .. $max_depth
#     z = (z * z) + c
#     break if z.abs2 > 5
#   end
#   return i
  return Mandel.loop(re, im, $max_depth)
end

def calc
  return if $current_rect

  for x in 0 .. Width - 1
    depth = mandel(x, $calc_y)
    if depth >= $max_depth
      $im.setPixel(x, $calc_y, $black)
    else
      $im.setPixel(x, $calc_y, $colors[ $color_max * depth / $max_depth ])
    end
  end

  $calc_y += 1
  if (($calc_y % 20) == 0)
    print "#{($calc_y * 100 / Height)}% done. -- depth #{$max_depth}\n"
    open("#tmpmandel#.gif", "w") do |f|
      $im.gif(f)
    end
    unless $gif
      $gif = TkPhotoImage.new { file "#tmpmandel#.gif" }
    else
      $gif.file "#tmpmandel#.gif"
    end
    unless $mandel
      $mandel = TkcImage.new($c, Width/2, Height/2) { image $gif }
    else
      $mandel.image($gif)
    end
  end

  if ($calc_y > Height - 1)
    $calc_y = StartCalcY
    $calc_on = false
    open("#tmpmandel#.gif", "w") do |f|
      $im.gif(f)
    end
    $gif.file "#tmpmandel#.gif"
    $mandel.image($gif)
#    $gif.destroy
#    $gin = nil
  end

  if $calc_on
    Tk.after(1) { calc() }
  end
end

reset()
$calc_y = StartCalcY = 0
$calc_on = true
gd_setup()
calc()

def clear
#  $mandel.destroy if $mandel
  $calc_y = StartCalcY
end

$start_x = $start_y = 0
$current_rect = nil

def do_press(x, y)
  $start_x = x
  $start_y = y
  $current_rect = Rectangle.new($c, x, y, x, y) { outline "white" }
end

def do_motion(x, y)
  if $current_rect
    $current_rect.coords $start_x, $start_y, x, y
  end
end

def do_release(x, y)
  if $current_rect
    $current_rect.coords $start_x, $start_y, x, y
    $current_rect.destroy
    $current_rect = nil
    clear()
    $calc_on = true
    gd_setup()
    zoom($start_x, $start_y, x, y)
    calc()
  end
end

$c.bind("1", proc{|e| do_press e.x, e.y})
$c.bind("B1-Motion", proc{|x, y| do_motion x, y}, "%x %y")
$c.bind("ButtonRelease-1", proc{|x, y| do_release x, y}, "%x %y")

begin
  Tk.mainloop
ensure
  File.delete("#tmpmandel#.gif")
end

Thread

Prev Next

In This Thread

Prev Next

[#1371] Re: [ruby-list:6104] Re: パースエラー等の嵐を Dos窓で受けるには ...? — Masao Kanemitsu <masao-k@...>

[#1377] Email your AD to 57 MILLION People ONLY $99 — 64687116@...

[#1381] bound method — shugo@... (Shugo Maeda)

[#1386] downcase — Shin-ichiro HARA <sinara@...>

[#1397] Tk.mainloopと Thread — Masao Kanemitsu <masao-k@...>

[#1406] Inclusion test — Kazuhisa YANAGAWA <katze@...>

[#1408] GD — Masao Kanemitsu <masao-k@...>

[#1410] GD — Masao Kanemitsu <masao-k@...>

[#1412] GD — Masao Kanemitsu <masao-k@...>

[#1424] Queue#waiters — Kazuhisa YANAGAWA <katze@...>

[#1449] perse error — Masao Kanemitsu <masao-k@...>

[#1455] parse.y — 助田 雅紀 <masaki.suketa@...>

[#1457] RE: ruby.def(mswin32) — 助田 雅紀 <masaki.suketa@...>

[#1463] iterator block passing (Re: rgrep) — shugo@... (Shugo Maeda)

[#1468] tcltk library — Masao Kanemitsu <masao-k@...>

[#1473] argumant type — Masao Kanemitsu <masao-k@...>

[#1481] tkmandel — Masao Kanemitsu <masao-k@...>

[#1483] 拡張モジュール on Windows — Masao Kanemitsu <masao-k@...>

[#1501] Re: 拡張モジュールon Windows — Masao Kanemitsu <masao-k@...>

[#1505] final.rb:35: [BUG] bug in variable assignment — shugo@... (Shugo Maeda)

[#1506] [BUG] instance_of? — keiju@... (Keiju ISHITSUKA)

[#1507] syscall — shugo@... (Shugo Maeda)

[#1508] flock -> lock — shugo@... (Shugo Maeda)

[#1512] Re: ruby 1.1b8 released — matz@... (Yukihiro Matsumoto)

[#1520] SizedQueue and some question — keiju@... (Keiju ISHITSUKA)

[#1528] Mutex#unlock — shugo@... (Shugo Maeda)

[#1536] return_value patch (Re: Mutex/ConditionVariable/Queue) — matz@... (Yukihiro Matsumoto)

[#1540] OS/2 port — KENN@...

[#1541] Re: OS/2 port — MORIOKA Toru <vette@...>

[ruby-dev:1500] Re: tkmandel

Attachments (1)

Thread

In This Thread

[#1455] parse.y — 助田　雅紀 <masaki.suketa@...>

[#1457] RE: ruby.def(mswin32) — 助田　雅紀 <masaki.suketa@...>