[dennou-ruby:003305] Re: ruby-lapack の ruby 1.9.1 対応.

[Seiya Nishizawa <seiya@xxxxxxxxxxxxxx>]

佐々木様

CLAPACK 同梱の F2C を見てみると、
typedef long int integer;
となっていて、8バイトですが、

debian パッケージの libf2c-dev に入っている f2c.h をみると
typedef int integer;
となっていて 4バイトです。

ruby-lapack は integer が 8バイトでであるとして作っているので、
配列を壊してしまっているのだと思います。


西澤誠也

2011/3/24 Seiya Nishizawa <seiya@xxxxxxxxxxxxxx>:
> 佐々木様
>
> debian 6.0 amd64
> libatlas-base-3.8.3-27 (package)
>
> で以下のようになります。
>
> 同じ 6.0 で ruby, atlas のバージョンが違うのはなぜでしょう
>
>
>
> kurtosis% ruby1.9.1 -r narray samples/dsyevr.rb
> 3
> NArray.float(3):
> [ 2.66454e-15, 1.0, 3.0 ]
> NArray.float(3,3):
> [ [ 0.57735, -0.57735, 0.57735 ],
>  [ -0.707107, -9.42055e-16, 0.707107 ],
>  [ 0.408248, 0.816497, 0.408248 ] ]
> NArray.int(12):
> [ 1, 3, 1, 3, 1, 3, 0, 0, 0, 0, 0, 0 ]
> NArray.float(78):
> [ 78.0, 0.0, 0.0, 1.0, 2.0, 1.0, 6.0, 1.0, 0.0, -2.0, -0.5, -4.15668e-13, ... ]
> NArray.int(60):
> [ 30, 0, 0, 1, 1, 1, 1, 2, 3, 3, 3, 2, 1, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, ... ]
> 0
> NArray.float(3,3):
> [ [ 1.0, 1.0, 0.0 ],
>  [ 1.0, 2.0, 1.0 ],
>  [ 0.0, 1.0, 1.0 ] ]
> kurtosis% ruby1.9.1 -v
> ruby 1.9.2p0 (2010-08-18 revision 29036) [x86_64-linux]
>
>
> 西澤誠也
>
>
>>> At Thu, 24 Mar 2011 16:14:34 +0900,
>>> Seiya Nishizawa <seiya@xxxxxxxxxxxxxx> wrote:
>>>> ruby 最新版で試してみましたが、
>>>> 結果は同じでした。
>>>
>>> うげ.やな塩梅ですね...
>>> # どこ見たら良いんだろう...
>>>
>>> 手元の環境ですが,
>>>
>>>  OS: Debian 6.0(squeeze) amd64
>>>  Ruby1.9: ruby 1.9.2p180 (2011-02-18 revision 30909) [x86_64-linux]
>>>  NArray: libnarray-ruby1.9.1 0.5.9p7-1 (package)
>>>  (C)LAPACK: libatlas-base-dev 3.8.3-30 (package)
>>>
>>> です.  これは (C)LAPACK の問題かなぁ...
>
>
>
>
>
> --
> Seiya Nishizawa
> Center for Planetary Science, Kobe University
> TEL : 078-803-6472
> FAX : 078-803-5731
>



-- 
Seiya Nishizawa
Center for Planetary Science, Kobe University
TEL : 078-803-6472
FAX : 078-803-5731