narray 0.5.9.7 → 0.5.9.8

data/src/ChangeLog

@@ -1,9 +1,44 @@
+2010-09-13  David MacMahon <davidm at astro.berkeley...>
+
+	* narray.c (na_where2):
+	  - Convert to NA_BYTE by calling "obj.ne(0)"
+	* ver 0.5.9p8
+
+2010-07-12  Masahiro TANAKA  <masa16.tanaka at gmail.com>
+
+	* narray.c (na_s_to_na):
+	  - improve argument check in NArray.to_na
+
+2010-04-30  Masahiro TANAKA  <masa16.tanaka at gmail.com>
+
+	* narray.c (na_s_to_na): 
+	  - the number of arguments was not checked.
+	    Thanks to S D for report.
+
+2010-01-16  Masahiro TANAKA  <masa16.tanaka at gmail.com>
+
+	* na_index.c (na_index_analysis): 
+	  - access out-of-bounds "struct slice" array.
+	* SPEC.en, SPEC.ja:
+	  - fix wrong explanation in count_false/true.
+	* na_random.c (na_random): 
+	  - modify random method not to change itself.
+	* lib/narray_ext.rb:
+	  - modify randomn! method to replace itself.
+	    Thank Yusuke ENDOH for report and comments.
+
+2009-06-14  Masahiro TANAKA  <masa16.tanaka at gmail.com>
+
+	* narray.h:
+	  include sys/stdint.h rather than sys/types.h
+	  for C99 confomance
+	* ver 0.5.9p7
+
 2009-05-20  Masahiro TANAKA  <masa16.tanaka at gmail.com>
 
 	* lib/narray_ext.rb (mean, stddev,rms,rsmdev,covariance): 
 	  - swith to DFLOAT type if integer array.
 	  - thanks to Philip Silva for report.
-	* ver 0.5.9p7
 
 2009-03-11  Masahiro TANAKA  <masa16.tanaka at gmail.com>
 

data/src/README.ja

@@ -1,53 +1,53 @@
 	Ruby/NArray  ver 0.5.9 (2006-08-09)	by Masahiro TANAKA
 
-* Ruby/NArray$B$NFCD'!#(B
-
-  + Ruby$B$G!"9bB.$J?tCM7W;;$,2DG=!#(B
-  + $BMWAG$K$O!"(B8,16,32 bit $B@0?t!"C1@:EY(B/$BG\@:EY$N<B?t(B/$BJ#AG?t!"(B
-    $B$*$h$S(B Ruby$B%*%V%8%'%/%H$r%5%]!<%H!#(B
-  + $BItJ,G[Ns$N<h=P$7!"BeF~$bMF0W!#(B
-    $BMWAG0LCV$N;XDj$K$O!"?tCM!"HO0O!"%$%s%G%C%/%9$NG[Ns$,;HMQ2DG=!#(B
-  + +,-,*,/,%,** $B$d(B $B;;=Q4X?t$N1i;;$O!"MWAG(B-$BBP(B-$BMWAG$G$*$3$J$&!#(B
-  + $BG[NsF1;N$N1i;;!&BeF~$O!"3F<!85$N%5%$%:$,F1$8$G$"$k$3$H$,I,MW!#(B
-    $B$?$@$7!"%5%$%:$,(B1$B$N<!85$O!"B>J}$NG[Ns$N%5%$%:$K9g$o$;$F(B
-    $B!V7+$jJV$7!WF1$8MWAG$,E,MQ$5$l$k!#(B
-  + FFTW ($B9bB.%U!<%j%(JQ49(B) version 3 $B$O<!$N%b%8%e!<%k$G%5%]!<%H!#(B
-  + NImage (X11$B%$%a!<%8I=<((B) $B%/%i%9ImB0!#(B(nimage/ $B%G%#%l%/%H%j;2>H(B)
-  + Ruby/PGPLOT ($B%0%i%U%#%C%/%9%i%$%V%i%j!"JLW}(B) $B$K$F(B
-     XY$B%0%i%U!"%R%9%H%0%i%`!"Ey9b@~!"%$%a!<%8I=<(2DG=!#(B
-  + $B?tCM7W;;!&2hA|=hM}!&%G!<%?2r@O$J$II}9-$$1~MQ$,2DG=(B($B$H;W$&(B)$B!#(B
-
-* $BN`;wIJ(B
+* Ruby/NArrayの特徴。
+
+  + Rubyで、高速な数値計算が可能。
+  + 要素には、8,16,32 bit 整数、単精度/倍精度の実数/複素数、
+    および Rubyオブジェクトをサポート。
+  + 部分配列の取出し、代入も容易。
+    要素位置の指定には、数値、範囲、インデックスの配列が使用可能。
+  + +,-,*,/,%,** や 算術関数の演算は、要素-対-要素でおこなう。
+  + 配列同士の演算・代入は、各次元のサイズが同じであることが必要。
+    ただし、サイズが1の次元は、他方の配列のサイズに合わせて
+    「繰り返し」同じ要素が適用される。
+  + FFTW (高速フーリエ変換) version 3 は次のモジュールでサポート。
+  + NImage (X11イメージ表示) クラス附属。(nimage/ ディレクトリ参照)
+  + Ruby/PGPLOT (グラフィックスライブラリ、別悃) にて
+     XYグラフ、ヒストグラム、等高線、イメージ表示可能。
+  + 数値計算・画像処理・データ解析など幅広い応用が可能(と思う)。
+
+* 類似品
 
   + Python/NumPy, Perl/PDL, Yorick, IDL
 
-* $BIT==J,$JE@(B
+* 不十分な点
 
-  + $B%a%=%C%I$,ITB-!#(B
-  + $B%P%0=P$7$,IT==J,!#(B
-  + $B%I%-%e%a%s%H$,$J$$!#(B
+  + メソッドが不足。
+  + バグ出しが不十分。
+  + ドキュメントがない。
 
-* $B%$%s%9%H!<%kJ}K!(B
+* インストール方法
 
-  Ruby$B$NI8=`E*$J3HD%%i%$%V%i%j$HF1$8$G$9!#%=!<%9$rE83+$7$?%G%#%l%/%H%j$G!"(B
+  Rubyの標準的な拡張ライブラリと同じです。ソースを展開したディレクトリで、
 
     ruby extconf.rb
     make
     make site-install
-    ($B$^$?$O(B make install)
+    (または make install)
 
-  $B$H<B9T$7$^$9!#(B
+  と実行します。
 
-* $BF0:n3NG'(B
+* 動作確認
 
   ruby 1.8.4 (2005-12-24) [i686-linux]
   gcc version 4.1.1 20060525 (Red Hat 4.1.1-1)
 
-* $BG[I[>r7o(B
+* 配布条件
 
-  Ruby$B$N%i%$%;%s%9$K=`5r$7$^$9!#(B
+  Rubyのライセンスに準拠します。
 
-* $BCx<T(B
+* 著者
 
-  $BEDCf>;9((B
-  $B$40U8+!$%P%0%l%]!<%H$=$NB>$O(B masa16.tanaka@gmail.com $B$^$G$*4j$$$7$^$9!#(B
+  田中昌宏
+  ご意見，バグレポートその他は masa16.tanaka@gmail.com までお願いします。

data/src/SPEC.en

@@ -1,5 +1,5 @@
 
-	Ruby/NArray  ver 0.5.9p5 (2008-05-28)	by Masahiro TANAKA
+	Ruby/NArray  ver 0.5.9p8 (2010-12-06)	by Masahiro TANAKA
 
 
 Class method:
@@ -218,9 +218,9 @@ Byte swap
 
 Boolean / mask related
     self.count_false	count # of elements whose value == 0 (only for
-			byte, sint, and int)
+			byte type)
     self.count_true	count # of elements whose value != 0 (only for
-			byte, sint, and int)
+			byte type)
     self.mask( mask )	same as self[ mask ], but exclusively for masking.
 			Unlike [], a int or sint mask is accepted.
 

data/src/SPEC.ja

@@ -1,9 +1,9 @@
 
-	Ruby/NArray  ver 0.5.9p5 (2008-05-28)	by Masahiro TANAKA
+	Ruby/NArray  ver 0.5.9p8 (2010-12-06)	by Masahiro TANAKA
 
 
-$B%/%i%9%a%=%C%I(B:
-    NArray.new(typecode, size, ...)	$BG[Ns$r@8@.$9$k!#MWAG$O(B0$B$G=i4|2=!#(B
+クラスメソッド:
+    NArray.new(typecode, size, ...)	配列を生成する。要素は0で初期化。
 
     NArray.byte(size,...)		1 byte unsigned integer
     NArray.sint(size,...)		2 byte signed integer
@@ -13,88 +13,88 @@
     NArray.scomplex(size,...)		single precision complex
     NArray.complex(size,...)		double precision complex
     NArray.object(size,...)		Ruby object
-					$B0J>eMWAG$O(B0$B$^$?$O(Bnil$B$G=i4|2=!#(B
+					以上要素は0またはnilで初期化。
 
-    NArray.to_na(array)			NArray$B$KJQ49(B
+    NArray.to_na(array)			NArrayに変換
     NArray.to_na(string,type[,size,..])
     NArray[...]
       NArray[1,5,10.0]  #=>  NArray.float(3):[1.0, 5.0, 10.0]
       NArray[1..10]     #=>  NArray.int(10):[1,2,3,4,5,6,7,8,9,10]
 
 
-$B%/%i%9JQ?t(B:
-    CLASS_DIMENSION	$B%G!<%?$H$7$F07$o$l$k<!85!#(B
-			NArray$B$O(B0$B!#(BNVector$B$O(B1$B!#(BNMatrix$B$O(B2$B!#(B
+クラス変数:
+    CLASS_DIMENSION	データとして扱われる次元。
+			NArrayは0。NVectorは1。NMatrixは2。
 
-$BG[Ns>pJs;2>H(B
-    self.dim	    $B<!85(B($B%$%s%G%C%/%9$N?t(B)$B$rJV$9!#(B
-    self.rank	    $B<!85(B($B%$%s%G%C%/%9$N?t(B)$B$rJV$9!#(B
-    self.shape	    $B<!85$4$H$N%5%$%:$rJV$9!#(B
-    self.total	    $BA4MWAG?t$rJV$9!#(B
+配列情報参照
+    self.dim	    次元(インデックスの数)を返す。
+    self.rank	    次元(インデックスの数)を返す。
+    self.shape	    次元ごとのサイズを返す。
+    self.total	    全要素数を返す。
 
-$B%$%s%G%C%/%9;2>H(B
+インデックス参照
     self[ dim0, dim1, ... ]
 
-  -- $B%$%s%G%C%/%90z?t$K;XDj$G$-$k$b$N(B:  $B?tCM!"HO0O!"G[Ns!"(Btrue, false
-  -- $B%$%s%G%C%/%9$N=g=x(B:		FORTRAN $B7?(B
-  -- $BE:;z0z?t$,(B1$B$D$N>l9g!"B?<!85G[Ns$O(Bflatten$B$5$l$?(B1$B<!85G[Ns$H$_$J$5$l$k!#(B
-     $BNc(B: a $B$,(B 3x3 $BG[Ns$N$H$-!"(Ba[3] $B$O(B a[0,1] $B$NMWAG$r;X$9!#(B
-
-    a[ 1, 2, -1 ]	$BMWAG$N<h$j=P$7!#Ii?t$O:G8e$+$i?t$($k(B(-1$B$,:G8e(B)
-			$BMWAG;XDj$N<!85$O=LLs$5$l$k!#(B
-    a[ 0..3, 4..1 ]	$BHO0O<h$j=P$7!#HO0O$N:G8e$,:G=i$h$jA0$J$i$P5U=g$K$J$k!#(B
-    a[ [1,4,2] ]	$B%$%s%G%C%/%9G[Ns!#MWAG$,(B[a[1],a[4],a[2]]$B$N(BNArray$B$,JV$k!#(B
-    a[]			a.dup $B$HF1$8!#(B
-    a[ 0, true ]	a[0, 0..-1] $B$HF1$8!#(B
-    a[ 0, false ]	a$B$,(B3$B<!85$N$H$-!"(Ba[0,true,true] $B$HF1$8!#(B
-			$B>JN,$5$l$?<!85$9$Y$F$K(Btrue$B$r;XDj$7$?$N$HF1$8!#(B
-    a[ mask ]		$B%^%9%-%s%0(B. mask $B$OD9$5$,(B a $B$HEy$7$$(B byte $B7?(B 
-			NArray. mask $B$N3FMWAG$NCM$K1~$8$F!"(Ba $B$N$=$l$>(B
-			$B$l$OMn$5$l$k(B(0$B$N>l9g(B)$B$+!"J];}$5$l$k(B(0$B0J30$N>l9g(B)$B!#(B
-			$BNc(B: 
+  -- インデックス引数に指定できるもの:  数値、範囲、配列、true, false
+  -- インデックスの順序:		FORTRAN 型
+  -- 添字引数が1つの場合、多次元配列はflattenされた1次元配列とみなされる。
+     例: a が 3x3 配列のとき、a[3] は a[0,1] の要素を指す。
+
+    a[ 1, 2, -1 ]	要素の取り出し。負数は最後から数える(-1が最後)
+			要素指定の次元は縮約される。
+    a[ 0..3, 4..1 ]	範囲取り出し。範囲の最後が最初より前ならば逆順になる。
+    a[ [1,4,2] ]	インデックス配列。要素が[a[1],a[4],a[2]]のNArrayが返る。
+    a[]			a.dup と同じ。
+    a[ 0, true ]	a[0, 0..-1] と同じ。
+    a[ 0, false ]	aが3次元のとき、a[0,true,true] と同じ。
+			省略された次元すべてにtrueを指定したのと同じ。
+    a[ mask ]		マスキング. mask は長さが a と等しい byte 型 
+			NArray. mask の各要素の値に応じて、a のそれぞ
+			れは落される(0の場合)か、保持される(0以外の場合)。
+			例: 
 				a=NArray.float(2,2).indgen!
 				p a[ a.lt 3 ]
 				--> [ 0.0, 1.0, 2.0 ]
-				(a.lt 3 $B$O(B byte $B7?(B NArray $B$rJV$9(B)
-				($BF1$8$3$H$O(B a[ (a.lt 3).where ] $B$G$b=PMh$k(B)
+				(a.lt 3 は byte 型 NArray を返す)
+				(同じことは a[ (a.lt 3).where ] でも出来る)
 
-  -- self.slice(...)	self[...] $B$HF1$8$@$,!"D9$5$,#1$K$J$C$?<!85$rMn(B
-			$B$5$:(B(self[]$B$OMn$9(B)$B!"$b$H$N%i%s%/$rJ]$D!#C"$7!"(B
-			1$B<!85%$%s%G%C%/%9IU$1$H%^%9%-%s%0$ONc30(B([]$B$HF1(B
-			$B$8(B)$B!#(B
+  -- self.slice(...)	self[...] と同じだが、長さが１になった次元を落
+			さず(self[]は落す)、もとのランクを保つ。但し、
+			1次元インデックス付けとマスキングは例外([]と同
+			じ)。
 
-$B%$%s%G%C%/%9BeF~!#(B-- $B<h=P$7$H$[$\F1$8%k!<%k!#(B
+インデックス代入。-- 取出しとほぼ同じルール。
 
     a[ 1, 2, 3 ] = 1
     a[ 0..3, 1..4, 2..5 ] = 2
     a[ [1,3,2,4], true ] = 3
-    a[] = 4			a.fill!(4) $B$HF1$8!#(B
+    a[] = 4			a.fill!(4) と同じ。
   
-    a[0..2]  = b[1..5]		--> $BMWAG?t$,0[$J$k$N$G%(%i!<!#(B
-    a[1,2]   = b[0..2,1..3]	[1,2]$B$r;OE@$H$7$FBeF~!#(B
-    a[0..2,0..3]  = b[0..2,1]	$B7+$jJV$7BeF~!#(B
+    a[0..2]  = b[1..5]		--> 要素数が異なるのでエラー。
+    a[1,2]   = b[0..2,1..3]	[1,2]を始点として代入。
+    a[0..2,0..3]  = b[0..2,1]	繰り返し代入。
 				( a[0,0]=b[0,1],..,a[0,3]=b[0,1] )
 
-$BCM$N%;%C%H!#(B
-    self.indgen!([start[,step]]) start$B$+$i(Bstep$B$E$DA}2C$7$?CM$r%;%C%H!#(B
-    self.fill!(value)		 $B$9$Y$F$NMWAG$K(Bvalue$B$r%;%C%H!#(B
-    self.random!(max)		 0<=x<max $B$N0lMM$J%i%s%@%`CM$r@8@.!#(B
+値のセット。
+    self.indgen!([start[,step]]) startからstepづつ増加した値をセット。
+    self.fill!(value)		 すべての要素にvalueをセット。
+    self.random!(max)		 0<=x<max の一様なランダム値を生成。
 				 using MT19337
-    self.randomn		 $BJ?6Q(B0$B!"J,;6(B1$B$N@55,J,I[$N%i%s%@%`CM$r@8@.!#(B
+    self.randomn		 平均0、分散1の正規分布のランダム値を生成。
 				 (Box-Muller)
-    NArray.srand([seed])	 $BMp?t$N%7!<%I$r@_Dj!#(B
-				 $B>JN,;~$O;~9o$+$i<+F0@8@.!#(B
+    NArray.srand([seed])	 乱数のシードを設定。
+				 省略時は時刻から自動生成。
 
-$B1i;;(B: $BMWAG$4$H$K$*$3$J$&!#(B
+演算: 要素ごとにおこなう。
     a = NArray.float(3,3).indgen
     b = NArray.float(3,3).fill(10)
     c = a*b	# --> NArray.float(3,3)
 
     a = NArray.float(3,1).indgen
     b = NArray.float(1,3).fill(10)
-    c = a*b	# --> NArray.float(3,3) -- size=1$B$N<!85$O3HD%$9$k!#(B
+    c = a*b	# --> NArray.float(3,3) -- size=1の次元は拡張する。
 
-$B;;=Q1i;;;R(B
+算術演算子
     -self
     self + other
     self - other
@@ -110,19 +110,19 @@
     self.div! other
     self.mod! other
 
-    self.mul_add(other,dim,...)   (self * other).sum(dim,...)$B$H$[$\F1$8!#(B
-				  $B$?$@$7ESCf$GG[Ns$r:n$i$J$$!#(B
+    self.mul_add(other,dim,...)   (self * other).sum(dim,...)とほぼ同じ。
+				  ただし途中で配列を作らない。
 
-$B%S%C%H1i;;;R(B($B@0?t$N$_2DG=(B)
+ビット演算子(整数のみ可能)
     ~self
     self & other
     self | other
     self ^ other
 
-$BHf3S(B
-  -- $BMWAG$4$H$KCM$rHf3S$7!"7k2L$r(BBYTE$B7?(B NArray$B$rJV$9!#(B
-     true/false$B$G$J$$$3$H$KCm0U!#(B
-    self.eq other  ( ver 0.5.4 $B0J9_!"(B== $B1i;;;R$OGQ;_$7$^$7$?!#(B)
+比較
+  -- 要素ごとに値を比較し、結果をBYTE型 NArrayを返す。
+     true/falseでないことに注意。
+    self.eq other  ( ver 0.5.4 以降、== 演算子は廃止しました。)
     self.ne other
     self.gt other
     self >  other
@@ -133,89 +133,89 @@
     self.le other
     self <= other
 
-    self.and other  $BMWAG$4$H$N>r7oHf3S!#(B
+    self.and other  要素ごとの条件比較。
     self.or  other
     self.xor other
     self.not other
 
-    self.all?     $BMWAG$,$9$Y$F??$J$i$P??!#(B
-    self.any?     $BMWAG$N$I$l$+$,??$J$i$P??!#(B
-    self.none?    $BMWAG$N$I$l$+$,??$J$i$P??!#(B
-    self.where	  $BMWAG$,??$N%$%s%G%C%/%9G[Ns$rJV$9!#(B
-    self.where2	  $BMWAG$,??$H56$N%$%s%G%C%/%9G[Ns$r4^$`(B(Ruby)$BG[Ns$rJV$9!#(B
-
-	$BNc(B: idx_t,idx_f = (a>12).where2
-
-$BE}7W(B
-    self.sum(dim,..)	    $B;XDj$7$?<!85$NOB!#(B
-    self.mean(dim,..)       $B;XDj$7$?<!85$NJ?6Q!#(B
-    self.stddev(dim,..)     $B;XDj$7$?<!85$NI8=`JP:9(B($BI8K\I8=`JP:9(B)$B!#(B
-    self.rms(dim,..)        $B;XDj$7$?<!85$N(Broot mean square$B!#(B
-    self.rmsdev(dim,..)     $B;XDj$7$?<!85$N(Broot mean square deviation$B!#(B
-    self.min(dim,..)        $B;XDj$7$?<!85$N:G>.!#(B
-    self.max(dim,..)        $B;XDj$7$?<!85$N:GBg!#(B
-			    ($B>JN,;~$OA4$F$N<!85!#(BRange$B;XDj2D!#(B)
-    self.median(dim)	    0..dim$B$N<!85$NCf1{CM!#>JN,;~$O$9$Y$F$N<!85!#(B
-
-$B%=!<%H(B
-    self.sort(dim)	     0..dim$B$N<!85$G%=!<%H!#>JN,;~$O$9$Y$F$N<!85!#(B
-    self.sort_index(dim)     $B%=!<%H$7$?%$%s%G%C%/%9$rJV$9!#(B
-			     self[self.sort_index] $B$O(B self.sort $B$HF1Ey!#(B
-
-$BE>CV(B
+    self.all?     要素がすべて真ならば真。
+    self.any?     要素のどれかが真ならば真。
+    self.none?    要素のどれかが真ならば真。
+    self.where	  要素が真のインデックス配列を返す。
+    self.where2	  要素が真と偽のインデックス配列を含む(Ruby)配列を返す。
+
+	例: idx_t,idx_f = (a>12).where2
+
+統計
+    self.sum(dim,..)	    指定した次元の和。
+    self.mean(dim,..)       指定した次元の平均。
+    self.stddev(dim,..)     指定した次元の標準偏差(標本標準偏差)。
+    self.rms(dim,..)        指定した次元のroot mean square。
+    self.rmsdev(dim,..)     指定した次元のroot mean square deviation。
+    self.min(dim,..)        指定した次元の最小。
+    self.max(dim,..)        指定した次元の最大。
+			    (省略時は全ての次元。Range指定可。)
+    self.median(dim)	    0..dimの次元の中央値。省略時はすべての次元。
+
+ソート
+    self.sort(dim)	     0..dimの次元でソート。省略時はすべての次元。
+    self.sort_index(dim)     ソートしたインデックスを返す。
+			     self[self.sort_index] は self.sort と同等。
+
+転置
     self.transpose( dim0, dim1, .. )
-	$BG[Ns$NE>CV!#(Bself$B$NBh(B(dim0)$B<!85$,?7$7$$G[Ns$NBh(B0$B<!85$K$J$k!#(B
-	$BIi?t$O8e$+$i$N=gHV!#(B
-	transpose(-1,1..-2,0) $B$G:G=i$H:G8e$rF~$l49$(!#(B
+	配列の転置。selfの第(dim0)次元が新しい配列の第0次元になる。
+	負数は後からの順番。
+	transpose(-1,1..-2,0) で最初と最後を入れ換え。
 
-$B%$%s%G%C%/%9$NJQ99(B ($BMWAG?t$OITJQ(B)
+インデックスの変更 (要素数は不変)
     self.reshape!(size,...)
     self.shape= size,...
-    self.newdim!(dim,...)	$B;XDj0LCV$K%5%$%:(B1$B$N<!85$rA^F~$9$k!#(B
-
-$B%G!<%?$N;2>H(B
-    self.refer			self$B$N%G!<%?$r;2>H$9$kJL$N%*%V%8%'%/%H$r:n@.!#(B
-    self.reshape(size,...)	self.refer.reshape! $B$HF1MM!#(B
-    self.newdim(dim,...)	self.refer.newdim! $B$HF1MM!#(B
-
-$B7?JQ49(B
-    self.floor self$B$h$j>.$5$$:GBg$N@0?t$rJV$9!#(B
-    self.ceil  self$B$h$jBg$-$$:G>.$N@0?t$rJV$9!#(B
-    self.round self$B$K$b$C$H$b6a$$@0?t$rJV$9!#(B
-    self.to_f  $BCM$rIbF0>.?tE@?t$KJQ49$9$k!#(B
-    self.to_i  $BCM$r@0?t$KJQ49$9$k!#(B 
-    self.to_a  $BCM$r(BRuby$B$NG[Ns$KJQ49$9$k!#(B
-    self.to_s  $B%P%$%J%j%G!<%?$r$=$N$^$^(BRuby$B$NJ8;zNs%G!<%?$KJQ49$9$k!#(B
-    self.to_string  $B3FMWAG$rJ8;zNs$KJQ49$9$k!#(B
-
-$B%$%F%l!<%?(B
+    self.newdim!(dim,...)	指定位置にサイズ1の次元を挿入する。
+
+データの参照
+    self.refer			selfのデータを参照する別のオブジェクトを作成。
+    self.reshape(size,...)	self.refer.reshape! と同様。
+    self.newdim(dim,...)	self.refer.newdim! と同様。
+
+型変換
+    self.floor selfより小さい最大の整数を返す。
+    self.ceil  selfより大きい最小の整数を返す。
+    self.round selfにもっとも近い整数を返す。
+    self.to_f  値を浮動小数点数に変換する。
+    self.to_i  値を整数に変換する。 
+    self.to_a  値をRubyの配列に変換する。
+    self.to_s  バイナリデータをそのままRubyの文字列データに変換する。
+    self.to_string  各要素を文字列に変換する。
+
+イテレータ
     self.each {|i| ...}
     self.collect {|i| ...}
     self.collect! {|i| ...}
 
-$B%P%$%H%9%o%C%W(B
-    self.swap_byte	$B%P%$%H%9%o%C%W(B
-    self.hton		$B%M%C%H%o!<%/%P%$%H%*!<%@!<$KJQ49(B
+バイトスワップ
+    self.swap_byte	バイトスワップ
+    self.hton		ネットワークバイトオーダーに変換
     self.ntoh		
-    self.htov		VAX$B%P%$%H%*!<%@!<$KJQ49(B
+    self.htov		VAXバイトオーダーに変換
     self.vtoh		
 
-Boolean / $B%^%9%/4X78(B
-    self.count_false	$BCM(B == 0 $B$NMWAG?t(B (byte, sint, int $B$N$_(B)
-    self.count_true	$BCM(B == 1 $B$NMWAG?t(B (byte, sint, int $B$N$_(B)
-    self.mask( mask )	self[ mask ] $B$HF1$8$@$+%^%9%-%s%0@lMQ(B.
-			[] $B$H0c$$(B int, sint $B$N%^%9%/$b;H$($k(B.
+Boolean / マスク関係
+    self.count_false	値 == 0 の要素数 (byte型のみ)
+    self.count_true	値 == 1 の要素数 (byte型のみ)
+    self.mask( mask )	self[ mask ] と同じだかマスキング専用.
+			[] と違い int, sint のマスクも使える.
 
-$BJ#AG?t(B
+複素数
     self.real
     self.imag
     self.conj
     self.conj!
     self.angle		atan2(self.imag, self.real)
-    self.imag= other	$B5u?tItJ,$K(Bother$B$r%;%C%H!#(B
-    self.im		$B5u?tG\!#(B
+    self.imag= other	虚数部分にotherをセット。
+    self.im		虚数倍。
 
-NMath $B%b%8%e!<%k(B
+NMath モジュール
     sqrt(x)
     exp(x)
     log(x)
@@ -233,55 +233,55 @@ NMath $B%b%8%e!<%k(B
     covariance
 
 
-FFTW $B%b%8%e!<%k(B (fftw-2.1.3$B$r(Bshared lib$B$G%3%s%Q%$%k$7$?$b$N$G3NG'(B)
- ($BJL%b%8%e!<%k(B)
+FFTW モジュール (fftw-2.1.3をshared libでコンパイルしたもので確認)
+ (別モジュール)
     fftw(x,[1|-1])
-    convol(a,b)		FFTW$B$rMQ$$$?>v$_9~$_!#(B
+    convol(a,b)		FFTWを用いた畳み込み。
 
 
 NMatrix
 
-  NArray$B$N%5%V%/%i%9!#:G=i$N(B2$B<!85$r(BMatrix$B$H$7$FMQ$$$k!#(B
-  $B;D$j$N<!85$OB?<!85G[Ns$H$7$F07$o$l$k!#(B
-  $B<!85$N=g=x$O!"?t3X$G$NI=5-$H$O5U(B:  a_ij => a[j,i]
+  NArrayのサブクラス。最初の2次元をMatrixとして用いる。
+  残りの次元は多次元配列として扱われる。
+  次元の順序は、数学での表記とは逆:  a_ij => a[j,i]
 
-  $B%a%=%C%I(B:
-    +,-	 $BAj<j$,(B NMatrix $B$N$H$-$K1i;;2D!#(B
-    *	 $BAj<j$,(B NMatrix $B$^$?$O(B NVector $B$N$H$-$O(B Matrix$B@Q!#(B
-	 $BAj<j$,(B Numeric $B$^$?$O(B NArray $B$N$H$-$O(B Scalar$B@Q!#(B
-	 $BNc(B: NMatrix[[1,2],[3,4]] * [1,10]
+  メソッド:
+    +,-	 相手が NMatrix のときに演算可。
+    *	 相手が NMatrix または NVector のときは Matrix積。
+	 相手が Numeric または NArray のときは Scalar積。
+	 例: NMatrix[[1,2],[3,4]] * [1,10]
 	       == NMatrix[ [[1,2],[3,4]], [[10,20],[30,40]] ]
-    /	 $BAj<j$,(B Numeric $B$^$?$O(B NArray $B$N$H$-$O(BScalar$B=|;;!#(B
-	 $BAj<j$,(B square NMatrix $B$N$H$-$O(BLU$B$K$h$j@~7AJ}Dx<0$r2r$/!#(B
+    /	 相手が Numeric または NArray のときはScalar除算。
+	 相手が square NMatrix のときはLUにより線形方程式を解く。
 	 a/b == b.lu.solve(a)
 
-    transpose		$B0z?t$r>JN,$7$?>l9g$O!":G=i$N(BMatrix$B<!85$r8r49!#(B
+    transpose		引数を省略した場合は、最初のMatrix次元を交換。
     diagonal(other)
-    diagonal!(other)	$BBP3QMWAG$KCM$r%;%C%H!#0z?t>JN,;~$O(B1$B$r%;%C%H!#(B
-    I			$BBP3QMWAG$KCM$K(B1$B$r%;%C%H!#(B
-    inverse		$B5U9TNs(B
-    lu			LU$BJ,2r$r7W;;!#(BNMatrixLU $B%/%i%9$N%$%s%9%?%s%9$rJV$9!#(B
+    diagonal!(other)	対角要素に値をセット。引数省略時は1をセット。
+    I			対角要素に値に1をセット。
+    inverse		逆行列
+    lu			LU分解を計算。NMatrixLU クラスのインスタンスを返す。
 
 
 NVector
 
-  NArray$B$N%5%V%/%i%9!#:G=i$N(B1$B<!85$r(BVector$B$H$7$FMQ$$$k!#(B
-  $B;D$j$N<!85$OB?<!85G[Ns$H$7$F07$o$l$k!#(B
+  NArrayのサブクラス。最初の1次元をVectorとして用いる。
+  残りの次元は多次元配列として扱われる。
 
-  $B%a%=%C%I(B:
-    +,-	 $BAj<j$,(B NVector $B$N$H$-$K1i;;2D!#(B
-    *	 $BAj<j$,(B NMatrix $B$N$H$-$O(B Matrix$B@Q!#(B
-    	 $BAj<j$,(B NVector $B$N$H$-$O(B $BFb@Q!#(B
-	 $BAj<j$,(B Numeric $B$^$?$O(B NArray $B$N$H$-$O(B Scalar$B@Q!#(B
-    /	 $BAj<j$,(B Numeric $B$^$?$O(B NArray $B$N$H$-$O(B Scalar$B=|;;!#(B
-	 $BAj<j$,(B square NMatrix $B$N$H$-$O(BLU$B$K$h$j@~7AJ}Dx<0$r2r$/!#(B
+  メソッド:
+    +,-	 相手が NVector のときに演算可。
+    *	 相手が NMatrix のときは Matrix積。
+    	 相手が NVector のときは 内積。
+	 相手が Numeric または NArray のときは Scalar積。
+    /	 相手が Numeric または NArray のときは Scalar除算。
+	 相手が square NMatrix のときはLUにより線形方程式を解く。
 	 v/m == m.lu.solve(v)
 
 NMatrixLU
 
-  NMatrix#lu $B%a%=%C%I$K$h$j:n$i$l$k!#(B
-  LU (NMatrix) $B$H(B pivot (NVector) $B$r4^$`!#(B
+  NMatrix#lu メソッドにより作られる。
+  LU (NMatrix) と pivot (NVector) を含む。
 
-  $B%a%=%C%I(B:
-    solve(other)	LU$BJ,2r$N7k2L$r;H$C$F(B other $B$r2r$/!#(B
-			other $B$O(B NMatrix $B$^$?$O(B NVector $B$N%$%s%9%?%s%9!#(B
+  メソッド:
+    solve(other)	LU分解の結果を使って other を解く。
+			other は NMatrix または NVector のインスタンス。

data/src/extconf.rb

@@ -83,16 +83,20 @@ na_math
 na_linalg
 )
 
-if have_header("sys/types.h")
+header = "stdint.h"
+unless have_header(header)
   header = "sys/types.h"
-else
-  header = nil
+  unless have_header(header)
+    header = nil
+  end
 end
 
 have_type("u_int8_t", header)
+have_type("uint8_t", header)
 have_type("int16_t", header)
 have_type("int32_t", header)
 have_type("u_int32_t", header)
+have_type("uint32_t", header)
 #have_library("m")
 #have_func("sincos")
 #have_func("asinh")

data/src/lib/narray_ext.rb

@@ -160,7 +160,12 @@ class NArray
       self.class.refer(rr)
     end
   end
-  alias randomn! randomn
+  #alias randomn! randomn
+
+  def randomn!
+    self[]= random
+    self
+  end
 
   #SFloatOne = NArray.sfloat(1).fill!(1)
 end

data/src/na_index.c

@@ -230,9 +230,11 @@ static int
       }
       multi_ellip = 1;
     } else {
-      size = na_index_test( idx[i], ary->shape[j], &sl[j] );
-      if (size != 1)
-	total *= size;
+      if (j < ary->rank) {
+         size = na_index_test( idx[i], ary->shape[j], &sl[j] );
+         if (size != 1)
+             total *= size;
+      }
       ++j;
     }
   }

data/src/na_random.c

@@ -370,8 +370,9 @@ static void RndC(int n, char *p1, int i1, double rmax)
 na_func_t RndFuncs =
   { TpErr, RndB, RndI, RndL, RndF, RndD, RndX, RndC, TpErr };
 
+
 static VALUE
- na_random(int argc, VALUE *argv, VALUE self)
+ na_random_bang(int argc, VALUE *argv, VALUE self)
 {
   VALUE  vmax;
   struct NARRAY *ary;
@@ -397,10 +398,16 @@ static VALUE
   return self;
 }
 
+static VALUE
+ na_random(int argc, VALUE *argv, VALUE self)
+{
+  return na_random_bang(argc, argv, na_clone(self));
+}
+
 void
 Init_na_random()
 {
     rb_define_singleton_method(cNArray,"srand",na_s_srand,-1);
-    rb_define_method(cNArray, "random!", na_random,-1);
-    rb_define_alias(cNArray,  "random","random!");
+    rb_define_method(cNArray, "random!", na_random_bang,-1);
+    rb_define_method(cNArray, "random",  na_random,-1);
 }

data/src/narray.c

@@ -18,7 +18,7 @@ VALUE cNArray, cNArrayScalar, cComplex;
 
 ID na_id_beg, na_id_end, na_id_exclude_end;
 ID na_id_minus, na_id_abs, na_id_power;
-ID na_id_compare, na_id_and, na_id_or;
+ID na_id_compare, na_id_ne, na_id_and, na_id_or;
 ID na_id_class_dim;
 ID na_id_add, na_id_sbt, na_id_mul, na_id_div, na_id_mod;
 ID na_id_real, na_id_imag;
@@ -722,7 +722,7 @@ static VALUE
 
 /* singleton method: 
    NArray.to_na( string, type, size1,size2,...,sizeN )
-   NArray.to_na( object )
+   NArray.to_na( array )
 */
 static VALUE
  na_s_to_na(int argc, VALUE *argv, VALUE klass)
@@ -731,19 +731,23 @@ static VALUE
   VALUE v;
   struct NARRAY *ary;
 
-  switch(TYPE(argv[0])) {
-  case T_STRING:
+  if (argc < 1) {
+    rb_raise(rb_eArgError, "Argument is required");
+  }
+  if (TYPE(argv[0]) == T_STRING) {
     return na_str_to_na(argc-1,argv+1,argv[0]);
-  case T_ARRAY:
-    if (argc>1)
-      rb_raise(rb_eArgError,"extra arguments");
+  }
+  if (argc > 1) {
+    rb_raise(rb_eArgError, "Only one array argument must be provided");
+  }
+  if (TYPE(argv[0]) == T_ARRAY) {
     return na_ary_to_nary( argv[0], klass );
-  default:
-    v = na_make_object(na_object_type(argv[0]),1,&shape,klass);
-    GetNArray(v,ary);
-    SetFuncs[ary->type][NA_ROBJ](1, ary->ptr,0, argv,0);
-    return v;
   }
+  if (NA_IsNArray(argv[0])) {
+    return argv[0];
+  }
+  rb_raise(rb_eTypeError, "Argument must be Array or String (or NArray)");
+  return Qnil;
 }
 
 
@@ -1024,8 +1028,12 @@ static VALUE
   int32_t *idx1, *idx0;
   struct NARRAY *ary, *a1, *a0; /* a1=true, a0=false */
 
-  obj = na_cast_object(obj,NA_BYTE);
   GetNArray(obj,ary);
+  /* Convert to NA_BYTE by calling "obj.ne(0)", if needed */
+  if(ary->type != NA_BYTE) {
+    obj = rb_funcall(obj, na_id_ne, 1, INT2FIX(0));
+    GetNArray(obj,ary);
+  }
   n = ary->total;
 
   /* Count true */
@@ -1278,6 +1286,7 @@ void
     na_id_usec 	= rb_intern("usec");
     na_id_now 	= rb_intern("now");
     na_id_compare = rb_intern("<=>");
+    na_id_ne    = rb_intern("ne");
     na_id_and   = rb_intern("&&");
     na_id_or    = rb_intern("||");
     na_id_minus = rb_intern("-@");

data/src/narray.h

@@ -15,19 +15,30 @@
 
 #include "narray_config.h"
 
+#ifdef HAVE_STDINT_H
+# include <stdint.h>
+#endif
+
 #ifdef HAVE_SYS_TYPES_H
 # include <sys/types.h>
 #endif
 
-#define NARRAY_VERSION "0.5.9p6"
-#define NARRAY_VERSION_CODE 596
+#define NARRAY_VERSION "0.5.9p8"
+#define NARRAY_VERSION_CODE 598
 
 /*
   Data types used in NArray :
   Please modify these types if your system has any different type.
 */
+
+
+/* NA_BYTE : unsigned 8-bit integer */
 #ifndef HAVE_U_INT8_T
-typedef unsigned char          u_int8_t; /* NA_BYTE */
+# ifdef HAVE_UINT8_T
+typedef uint8_t			u_int8_t;
+# else
+typedef unsigned char		u_int8_t;
+# endif
 #endif
 
 #ifndef HAVE_INT16_T
@@ -50,14 +61,19 @@ typedef int                    int32_t;  /* NA_LINT */
 # endif
 #endif /* HAVE_INT32_T */
 
+/* unsigned 32-bit integer */
 #ifndef HAVE_U_INT32_T
-# if SIZEOF_LONG == 4
-typedef unsigned long                   u_int32_t;  /* NA_LINT */
+# ifdef HAVE_UINT32_T
+typedef uint32_t			u_int32_t;
 # else
-#  if SIZEOF_INT == 4
-typedef unsigned int                    u_int32_t;  /* NA_LINT */
+#  if SIZEOF_LONG == 4
+typedef unsigned long                   u_int32_t;  
 #  else
+#   if SIZEOF_INT == 4
+typedef unsigned int                    u_int32_t;
+#   else
 ---->> Please define u_int32_t manually because sizeof(long) != 4. <<----
+#   endif
 #  endif
 # endif
 #endif /* HAVE_U_INT32_T */

metadata

@@ -1,33 +1,26 @@
 --- !ruby/object:Gem::Specification 
-rubygems_version: 0.9.4
-specification_version: 1
 name: narray
 version: !ruby/object:Gem::Version 
-  version: 0.5.9.7
-date: 2009-06-14 00:00:00 +09:00
-summary: N-dimensional Numerical Array class for Ruby
-require_paths: 
-- .
-email: masa16.tanaka@gmail.com
-homepage: http://narray.rubyforge.org/
-rubyforge_project: narray
-description: 
-autorequire: 
-default_executable: 
-bindir: bin
-has_rdoc: false
-required_ruby_version: !ruby/object:Gem::Version::Requirement 
-  requirements: 
-  - - ">"
-    - !ruby/object:Gem::Version 
-      version: 0.0.0
-  version: 
+  version: 0.5.9.8
 platform: ruby
-signing_key: 
-cert_chain: 
-post_install_message: 
 authors: 
 - Masahiro Tanaka
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2010-12-06 00:00:00 +09:00
+default_executable: 
+dependencies: []
+
+description: 
+email: masa16.tanaka@gmail.com
+executables: []
+
+extensions: 
+- src/extconf.rb
+extra_rdoc_files: []
+
 files: 
 - src/ChangeLog
 - src/MANIFEST
@@ -82,17 +75,33 @@ files:
 - src/bench/dummy.m
 - src/bench/dummy.py
 - src/bench/dummy.rb
-test_files: []
+has_rdoc: true
+homepage: http://narray.rubyforge.org/
+licenses: []
 
+post_install_message: 
 rdoc_options: []
 
-extra_rdoc_files: []
-
-executables: []
-
-extensions: 
-- src/extconf.rb
+require_paths: 
+- .
+required_ruby_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
 requirements: []
 
-dependencies: []
+rubyforge_project: narray
+rubygems_version: 1.3.5
+signing_key: 
+specification_version: 3
+summary: N-dimensional Numerical Array class for Ruby
+test_files: []