lib/numru/gphys/interpolate.rb in gphys-1.4.3.2 vs lib/numru/gphys/interpolate.rb in gphys-1.5.0
- old
+ new
@@ -1,5 +1,6 @@
+# -*- coding: utf-8 -*-
require "numru/gphys/gphys"
require "numru/dcl" # math1/gt2dlib is used for pure-2D interpolation.
# Also for dcl_fig_cut
require "numru/dcl_mouse" # for mouse_cut, mouse_cut_repeat
require "narray_miss"
@@ -474,11 +475,11 @@
ids1 = NArray.int(len)
f1 = NArray.float(len)
ids2 = NArray.int(len)
f2 = NArray.float(len)
for j in 0...len
- ## [開発メモ] (高速化) このループはCにしたほうがいい（1Dなのでまあいいけど，下の2Dのをするなら一緒に）．その際，g2ictr は探索を高速化してないので，interpo_find_loc_1D の2D版の探索をした上で g2ibl2 を直接呼ぶ方がいい．
+ ## [髢狗匱繝｡繝｢] (鬮倬溷喧) 縺薙ｮ繝ｫ繝ｼ繝励ｯC縺ｫ縺励◆縺ｻ縺縺後＞縺ｼ1D縺ｪ縺ｮ縺ｧ縺ｾ縺ゅ＞縺縺代←ｼ御ｸ九ｮ2D縺ｮ繧偵☆繧九↑繧我ｸ邱偵↓ｼ会ｼ弱◎縺ｮ髫幢ｼ携2ictr 縺ｯ謗｢邏｢繧帝ｫ倬溷喧縺励※縺ｪ縺縺ｮ縺ｧｼ景nterpo_find_loc_1D 縺ｮ2D迚医ｮ謗｢邏｢繧偵＠縺滉ｸ翫〒 g2ibl2 繧堤峩謗･蜻ｼ縺ｶ譁ｹ縺後＞縺ｼ
begin
uxi, uyi = DCL.g2ictr(txi[j], tyi[j])
ids1[j] = [ [uxi.floor,0].max, uxg.length-2 ].min
ids2[j] = [ [uyi.floor,0].max, uyg.length-2 ].min
f1[j] = uxi - ids1[j]
@@ -501,10 +502,10 @@
f1 = NArray.float(lenx,leny)
ids2 = NArray.int(lenx,leny)
f2 = NArray.float(lenx,leny)
for k in 0...leny
for j in 0...lenx
- ## [開発メモ] (高速化) このループはCにしたほうがいい(2Dだし特に)．その際，g2ictr は探索を高速化してないので，interpo_find_loc_1D の2D版の探索をした上で g2ibl2 を直接呼ぶ方がいい．
+ ## [髢狗匱繝｡繝｢] (鬮倬溷喧) 縺薙ｮ繝ｫ繝ｼ繝励ｯC縺ｫ縺励◆縺ｻ縺縺後＞縺(2D縺縺礼音縺ｫ)ｼ弱◎縺ｮ髫幢ｼ携2ictr 縺ｯ謗｢邏｢繧帝ｫ倬溷喧縺励※縺ｪ縺縺ｮ縺ｧｼ景nterpo_find_loc_1D 縺ｮ2D迚医ｮ謗｢邏｢繧偵＠縺滉ｸ翫〒 g2ibl2 繧堤峩謗･蜻ｼ縺ｶ譁ｹ縺後＞縺ｼ
begin
uxi, uyi = DCL.g2ictr(txi[j], tyi[k])
ids1[j,k] = [ [uxi.floor,0].max, uxg.length-2 ].min
ids2[j,k] = [ [uyi.floor,0].max, uyg.length-2 ].min
f1[j,k] = uxi - ids1[j,k]