numpy-discussion — Discussion list for all users of Numerical Python

[Numpy-discussion] A Job offer

[Casandra <Ca...@ta...>]

Dear Friend,

Taken into account the expansion of a field of
activity of our firm, 
which is located in Europe and engaged in modeling and
fashion
industry in the world markets. 
We require Project Manager. The job is of a partial
employment, which means that
you will not spend a lot of energy and time on it. The
salary
depends on the price of projects executed by you 
(about 30000-50000 dollars a year). If our offer has
interested you and you are a serious and hardworking
person, who is not afraid to open a new field of
activity for yourself, you are probably the person
that we are looking for. Also, we will be glad to see you
in our company. To suit this job you have to be at least
21 years of age or older and be a US or Australian
citizen.
Tahnk you.

gal...@ga...

[Numpy-discussion] 1.0b5 Installation on OS X

[David A. <irb...@gm...>]

Hi,

I'm unable to install the 1.0b5 release from the .mpkg on OS X - when
it comes to the two checkboxes to select the platlib and the scripts,
both are greyed out.  The installer then claims 'Cannot continue:
Nothing to install'.

Any suggestions?

Cheers,

Dave

[Numpy-discussion] PHxoeARMA

[Sib R. <phi...@al...>]

Hi
=20
QUIT OVE k RPA y YIN r G FOR Y c OUR PH l AR g MAC n Y
=20
S w AV i E up t h o 50 z % wit p h http://www.criaboedri.info

Re: [Numpy-discussion] (no subject)

[<yuk...@ya...>]

:―― INFORMATION  ―――――――――――――――――――――――――:
 不正・悪質なサイトを一切排除しておりますので
 安心してご利用ください。   http://love-match.bz/pc/?02
:――――――――――――――――――――――――――――――――――:

*・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・*


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』
　-----------------------------------------------------------------
　　登録料・利用料　・・・・・・・・・【無料】
　　メールの送受信　・・・・・・・・・【無料】
　　ユーザーの検索　・・・・・・・・・【無料】
　　掲示板の閲覧・書込み　・・・・・・【無料】
　　画像交換・アップロード　・・・・・【無料】
　　アドレス交換・電話番号交換　・・・【無料】
　-----------------------------------------------------------------
　　どれだけ使っても全て無料！　http://love-match.bz/pc/?02

━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ いつでも女性ユーザーがいっぱい。その理由は？

　-----------------------------------------------------------------
　ＰＣ＆モバイルに対応！いつでもどこでも気軽に楽しめる！
　-----------------------------------------------------------------
　仕事中は携帯電話から、プライベートは自宅からのんびりと。
　気になる相手といつでも繋がっているから、新密度も急速にＵＰ。
　http://love-match.bz/pc/?02

　-----------------------------------------------------------------
　ＰＣから簡単プロフィール作成。ネット初心者でもラクラク参加ＯＫ
　-----------------------------------------------------------------
　面倒な登録は一切不要。パソコンから簡単なプロフィールを作成して
　初心者の方や女性でもすぐに参加できます。
　http://love-match.bz/pc/?02
　-----------------------------------------------------------------
　自由恋愛対応！直電・直メ交換支援ツール
　-----------------------------------------------------------------
　基本的にメールアドレスや電話番号は非公開ですが
　仲良くなった人だけにメールアドレスや電話番号を教える事ができます。
　http://love-match.bz/pc/?02

　-----------------------------------------------------------------
　写真アップロードに対応！好みの相手を素早くＣＨＥＣＫ！
　-----------------------------------------------------------------
　待ち合わせ場所にイメージとまったく違う人が来たら…。
　ピュアックスなら会う前に写真交換ができるから、そんな不安も解消。
　http://love-match.jp/pc/?06

　-----------------------------------------------------------------
　スレッド掲示板で秘密のパートナー検索も効率ＵＰ！
　-----------------------------------------------------------------
　メインの掲示板のほかにスレッド型の掲示板を設置。
　メル友から秘密のパートナーまで目的別のユーザーが集う掲示板です。
　http://love-match.jp/pc/?06


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 毎日500人近くのユーザーが続々参加中！！

□-----------------------------------------------------------------

　リエ（２１歳/会社員）
　いつも１人でエッチなことを考えてます。
　メールだといろいろ話せるんだけど、実際に会うとあまりしゃべれなく
　なっちゃうので、盛り上げてくれるような楽しい男の人いないかな？
　引っ込み思案のせいか、男性経験はあまり無いんです。
　優しく＆楽しくリードしてくれる男性からのメール待ってます。
　[写真有り] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　真菜（２４歳/フリーター）
　彼氏が浮気してて超アタマきたっ！まなだって遊びたい盛りだし、ずっと
　ガマンしてたのにさ！かっこいい人見つけて思いっきりふってやるつもりで
　登録してみた（笑）
　[写真有り] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　みさ（３４歳/専業主婦）
　殆ど家に帰ってこない仕事人間のだんなさまと二人きりの毎日で、ほんと
　寂しい思いをしています。年下の男の子がいれば仲良くなりたいな。
　年下の人とは付き合ったことがないので興味津々です(^^)
　[写真無し] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　恭子（２８歳/会社員）
　彼氏とはいつも同じようなセックスばかりでかなり冷め気味です。
　誰かあたしと熱いセックスを楽しみませんか？めんどくさい事は
　言いません。ただ、いつもと違うドキドキするような事がしたい
　だけなんです。
　[写真無し] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　ななこ（２８歳/主婦）
　半年前にだんなと別れて今は×1です。
　夜のお仕事なので、昼間まったりと過ごしませんか？
　心身ともに疲れ気味で、今、激しく癒されたいです。
　[写真有り] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　祥子（３１歳/クリエイター）
　平日は18時くらいまでは大体仕事してるので、その後に食事したり
　楽しく飲んだりできるパートナー希望です。年上でも年下でも
　かまいませんので気軽にメールを送って頂けると嬉しいです。
　[写真有り] http://love-match.bz/pc/?02

□-----------------------------------------------------------------

　ゅヵ`（２０歳/学生）
　まずゎ会ってみないとはじまらなぃょね？！
　横浜近辺の人で、いろんな意味でオトナな人は
　プロフ付きでめぇる送って☆
　[写真有り] http://love-match.bz/pc/?02

□-----------------------------------------------------------------


   出会いサイトのサクラに騙されないように↓
　┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓
　　【裏】無料の出会い情報
　　-------------------------------------------------------------
　　お金と時間を持て余している人妻の間で、噂になってるあのサイト
　　[登録・利用料全て無料]　http://love-match.bz/pc/?02
　　-------------------------------------------------------------
　　彼女達が求めているのはこんな男性です。
　　?年上女性にリードしてもらいたい、経験少なめの男性
　　?体力・テクニックに自信が有る男性
　　男性会員が不足しています。我こそは、と思う方は今すぐ参加！
　　[登録・利用料全て無料]　http://love-match.bz/pc/?02
　┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛

お問い合わせ先はこちら↓
広東省茂名市人民大街3-6-4-533
友誼網絡公司
139-3668-7892

Re: [Numpy-discussion] how to get info about internals of an array object ?

[Sebastian H. <ha...@ms...>]

On Friday 15 September 2006 12:49, Travis Oliphant wrote:
> Sebastian Haase wrote:
> >On Friday 15 September 2006 10:00, Travis Oliphant wrote:
> >>Sebastian Haase wrote:
> >>>Hi,
> >>>
> >>>what I'm asking is if numpy has an equivalent to numarray's info()
> >
> >function:
> >>>>>>na.arange(10).info()
> >>
> >>numpy.numarray.info(numpy.arange(10))
> >>
> >>(get recent SVN as there were some bugs just fixed.
> >>
> >>-Travis
> >
> >Thanks,
> >
> >should this maybe also be added somewhere in the numpy module itself !?
> >I guess the question is, what the original intent was for numarray to put
> > it in (only for debugging ?) -- then we could decide if this would also
> > make sense for numpy.
> >
> >I have never used numpy.info(obj) - I don't know what is does (compared to
> >__doc__)
>
> It prints the doc string of an object searching for it better than help
> seems to do and without a "pager".
>
> Compare numpy.info(numpy.sin) versus help(numpy.sin)
In PyShell I would miss the output of numpy.info(numpy.sin) since it goes to 
the (C) stdout and not to the (sometimes redirected ) python sys.stdout.
So it seems it prints N.sin.__doc__

The output from help(...) goes to the (python) sys.stdout (or stderr !?)  and 
is quite long (and mostly just generec info about ufuncs ...)

>
> I don't know what the info method was for other than debugging.
>
> What about having the __doc__ attribute of an ndarray return the info?
>
Maybe,
but honestly I would *not* expect any kind of "live inspection" to be done by 
__doc__
I would expect more like 
"ndarray is an efficient way to handle large data sets in python - it is the 
bases for Numerical Python, see www.scipy.org "

Maybe info() should just be an array attribute - just like it was it numarray.

-Sebastian

Re: [Numpy-discussion] how to get info about internals of an array object ?

[Travis O. <oli...@ee...>]

Sebastian Haase wrote:

>On Friday 15 September 2006 10:00, Travis Oliphant wrote:
>  
>
>>Sebastian Haase wrote:
>>    
>>
>>>Hi,
>>>
>>>what I'm asking is if numpy has an equivalent to numarray's info() 
>>>      
>>>
>function:
>  
>
>>>>>>na.arange(10).info()
>>>>>>            
>>>>>>
>>numpy.numarray.info(numpy.arange(10))
>>
>>(get recent SVN as there were some bugs just fixed.
>>
>>-Travis
>>    
>>
>
>Thanks,
>
>should this maybe also be added somewhere in the numpy module itself !?
>I guess the question is, what the original intent was for numarray to put it 
>in (only for debugging ?) -- then we could decide if this would also make 
>sense for numpy.
>
>I have never used numpy.info(obj) - I don't know what is does (compared to 
>__doc__)   
>
It prints the doc string of an object searching for it better than help 
seems to do and without a "pager".

Compare numpy.info(numpy.sin) versus help(numpy.sin)

I don't know what the info method was for other than debugging.  

What about having the __doc__ attribute of an ndarray return the info?


-Travis

Re: [Numpy-discussion] how to get info about internals of an array object ?

[Sebastian H. <ha...@ms...>]

On Friday 15 September 2006 10:00, Travis Oliphant wrote:
> Sebastian Haase wrote:
> > Hi,
> >
> > what I'm asking is if numpy has an equivalent to numarray's info() 
function:
> >>>> na.arange(10).info()
>
> numpy.numarray.info(numpy.arange(10))
>
> (get recent SVN as there were some bugs just fixed.
>
> -Travis

Thanks,

should this maybe also be added somewhere in the numpy module itself !?
I guess the question is, what the original intent was for numarray to put it 
in (only for debugging ?) -- then we could decide if this would also make 
sense for numpy.

I have never used numpy.info(obj) - I don't know what is does (compared to 
__doc__)   but maybe it could be extended like
if isinstance( obj , N.ndarray):
    print  N.numarray.info( obj )

right now it prints nothings (and returns None)

Thanks,
Sebastian Haase

Re: [Numpy-discussion] Experience with Visit?

[Rob H. <he...@ta...>]

I have looked at it, and I even know somebody who is using python to  
put in a data reader for FVCOM (a finite volume ocean model).  He has  
a prototype working.  I have no experience with the python bindings  
personally.

I looked into compiling it from source (for an intel mac), but it was  
daunting so I gave up.

I also contacted them to try to put in lower level python bindings,  
so that we could create 3D graphics from the command line, but they  
didn't seem too interested in that.

Otherwise, it is a very nice workable app.

-Rob


On Sep 15, 2006, at 11:23 AM, O'Keefe, Michael wrote:

> I haven't tried VisIT before but thanks for the link. I also  
> downloaded and am checking it out.
>
> Along this same line of discussion, has anyone tried OOF2 which is  
> an FEA package that also has some strong python connections?
>
> http://www.ctcms.nist.gov/oof/oof2.html
>
> I'm working on a Windows machine and it the current code-base  
> doesn't seem to support Windows out of the box (if at all). Looks  
> like you can put it together for *nix or Mac OS X, though...
>
> --Michael
>
>> -----Original Message-----
>> From: num...@li...
>> [mailto:num...@li...] On
>> Behalf Of Robert Cimrman
>> Sent: Friday, September 15, 2006 8:29
>> To: Discussion of Numerical Python
>> Subject: Re: [Numpy-discussion] Experience with Visit?
>>
>> Travis Oliphant wrote:
>>> Has anybody had any experience with the 3-D visualization software
>>> VISIT?   It has Python bindings and seems to be pretty
>> sophisticated.
>>> I'm wondering why I haven't heard more about it.
>>>
>>> http://www.llnl.gov/visit/
>>
>> No reaction up to now, so...
>>
>> I have just tried the 'getting started' part and was quite impressed,
>> thanks for posting the link! Up to now I have used ParaView
>> and was very
>> satisfied, but the Python bindings of VisIt are a great lure.
>>
>> r.
>>
>> --------------------------------------------------------------
>> -----------
>> Using Tomcat but need to do more? Need to support web
>> services, security?
>> Get stuff done quickly with pre-integrated technology to make
>> your job easier
>> Download IBM WebSphere Application Server v.1.0.1 based on
>> Apache Geronimo
>> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&
>> dat=121642
>> _______________________________________________
>> Numpy-discussion mailing list
>> Num...@li...
>> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>>
>
> ---------------------------------------------------------------------- 
> ---
> Using Tomcat but need to do more? Need to support web services,  
> security?
> Get stuff done quickly with pre-integrated technology to make your  
> job easier
> Download IBM WebSphere Application Server v.1.0.1 based on Apache  
> Geronimo
> http://sel.as-us.falkag.net/sel? 
> cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion

----
Rob Hetland, Associate Professor
Dept. of Oceanography, Texas A&M University
http://pong.tamu.edu/~rob
phone: 979-458-0096, fax: 979-845-6331

Re: [Numpy-discussion] buggy buggy bugyy: format and casting ?

[Robert K. <rob...@gm...>]

Eric Emsellem wrote:
> Hi,
> 
> I am facing a rather frustrating problem with numpy/scipy: after
> upgrading to svn numpy and scipy, and trying to remove most of the now
> unnecessary casting (floats) the program I wrote does not give me the
> right answer. It seems that the answer is somewhat scaled down (but not
> in a simple way). I had this kind of behaviour some time ago when scipy
> and numpy were not fully compatible and then I had to cast the floats
> explicitely in my program to make it work. But now it SHOULD be fully
> compatible as far as I understand the issue. So my question is:
> is there any remaining flaws in terms of casts/format between scipy and
> numpy? I am specifically using the scipy functions:
> - special.erf, special.erfc, orthogonal.ps_roots, integrate.quad...
> and the numpy functions : sum, sqrt, exp, sin, cos, arctan.
> I am doing all the calculations using numpy.float32.

Could you give us a short piece of code that demonstrates the problem? You don't 
state what flaws you were working around before. Note that some of those 
functions in scipy are unavoidably implemented in double precision at the 
C/FORTRAN level.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Martin W. <mar...@gm...>]

On Friday 15 September 2006 20:27, Travis Oliphant wrote:
> Martin Wiechert wrote:
> > On Friday 15 September 2006 19:14, Travis Oliphant wrote:
> >> Martin Wiechert wrote:
> >>> Thanks Albert! Do you also know the corresponding C-API function? It
> >>> cannot be PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose
> >>> signature has no "align", right?
> >>
> >> The DescrConverter function is meant for "O&"-style conversions.  It
> >> can't accept an align function.  We could possibly add something to the
> >> converter to allow specification of alignment through the object to be
> >> converted.
> >
> > I begin to see the light....
> >
> > For dictionaries one could maybe just add an optional key "align".
> > Also an optional key "elsize" or "itemsize" to force the total size of
> > the record may sometimes be useful. E.g. one may want to faithfully map a
> > given C struct. (That's why I was asking for trailing unused bytes.)
> >
> > But I understand that other things have higher priority.
> >
> >> Or, you can just call the __new__ method of the PyArrayDescr_Type object
> >>
> >> res = PyObject_CallMethod((PyObject *)&PyArrayDescr_Type, "__new__",
> >> "Oi", dict_object, 1))
> >>
> >> or call the tp->new method directly:
> >>
> >> args = Py_BuildValue("Oi", dict_object, 1);
> >> PyArrayDescr_Type->tp_new(&PyArrayDescr_Type, args, NULL);
> >> Py_DECREF(args);
> >
> > Thank you! I'll try this.
> >
> >> (I think passing in NULL for the keywords is O.K., but I haven't checked
> >> it).
> >>
> >> -Travis
> >
> > One final question. To me the repr of a dtype with gaps is a little bit
> >
> > puzzling:
> >>>> dtype ({'names': ['a', 'b', 'c'], 'formats': ['<a4', '<f8', '<f4'],
> >
> > 'offsets': [0, 16, 24]})
> > dtype([('a', '|S4'), ('', '|V12'), ('b', '<f8'), ('', '|V12'), ('c',
> > '<f4')])
> >
> > There should be no gap between "b" and "c" but still the repr has ('',
> > '|V12') between them. Am I missing something?
>
> There was a bug I just fixed in the representation of these structures
> with gaps.  It should be fixed in SVN, now.
>
>
> -Travis
>
>

Thanks again.

> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job
> easier Download IBM WebSphere Application Server v.1.0.1 based on Apache
> Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Travis O. <oli...@ie...>]

Martin Wiechert wrote:
> On Friday 15 September 2006 19:14, Travis Oliphant wrote:
>   
>> Martin Wiechert wrote:
>>     
>>> Thanks Albert! Do you also know the corresponding C-API function? It
>>> cannot be PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose
>>> signature has no "align", right?
>>>       
>> The DescrConverter function is meant for "O&"-style conversions.  It
>> can't accept an align function.  We could possibly add something to the
>> converter to allow specification of alignment through the object to be
>> converted.
>>
>>     
>
> I begin to see the light....
>
> For dictionaries one could maybe just add an optional key "align".
> Also an optional key "elsize" or "itemsize" to force the total size of the 
> record may sometimes be useful. E.g. one may want to faithfully map a given C 
> struct. (That's why I was asking for trailing unused bytes.)
>
> But I understand that other things have higher priority.
>
>   
>> Or, you can just call the __new__ method of the PyArrayDescr_Type object
>>
>> res = PyObject_CallMethod((PyObject *)&PyArrayDescr_Type, "__new__",
>> "Oi", dict_object, 1))
>>
>> or call the tp->new method directly:
>>
>> args = Py_BuildValue("Oi", dict_object, 1);
>> PyArrayDescr_Type->tp_new(&PyArrayDescr_Type, args, NULL);
>> Py_DECREF(args);
>>
>>     
>
> Thank you! I'll try this.
>
>   
>> (I think passing in NULL for the keywords is O.K., but I haven't checked
>> it).
>>
>> -Travis
>>
>>     
>
> One final question. To me the repr of a dtype with gaps is a little bit 
> puzzling:
>
>   
>>>> dtype ({'names': ['a', 'b', 'c'], 'formats': ['<a4', '<f8', '<f4'], 
>>>>         
> 'offsets': [0, 16, 24]})
> dtype([('a', '|S4'), ('', '|V12'), ('b', '<f8'), ('', '|V12'), ('c', '<f4')])
>
> There should be no gap between "b" and "c" but still the repr has ('', '|V12') 
> between them. Am I missing something?
>   
There was a bug I just fixed in the representation of these structures 
with gaps.  It should be fixed in SVN, now.


-Travis

Re: [Numpy-discussion] PyArray_DescrConverter - alignment/ trailingunused bytes

[Albert S. <fu...@gm...>]

Argh

> > <snip>
> 
> > One final question. To me the repr of a dtype with gaps is a little bit
> > puzzling:
> >
> > >>> dtype ({'names': ['a', 'b', 'c'], 'formats': ['<a4', '<f8', '<f4'],
> > 'offsets': [0, 16, 24]})
> > dtype([('a', '|S4'), ('', '|V12'), ('b', '<f8'), ('', '|V12'), ('c',
> > '<f4')])
> >
> > There should be no gap between "b" and "c" but still the repr has ('',
> > '|V12')
> > between them. Am I missing something?

I see you're not specifying align=True here. Ignore my last message.

Cheers,

Albert

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Albert S. <fu...@gm...>]

Hello all

> -----Original Message-----
> From: num...@li... [mailto:numpy-
> dis...@li...] On Behalf Of Martin Wiechert
> Sent: 15 September 2006 19:53
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion]PyArray_DescrConverter - alignment /
> trailingunused bytes
> 
> On Friday 15 September 2006 19:14, Travis Oliphant wrote:
> > Martin Wiechert wrote:
> > > Thanks Albert! Do you also know the corresponding C-API function? It
> > > cannot be PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose
> > > signature has no "align", right?
> >
> > The DescrConverter function is meant for "O&"-style conversions.  It
> > can't accept an align function.  We could possibly add something to the
> > converter to allow specification of alignment through the object to be
> > converted.
> >

> <snip>

> One final question. To me the repr of a dtype with gaps is a little bit
> puzzling:
> 
> >>> dtype ({'names': ['a', 'b', 'c'], 'formats': ['<a4', '<f8', '<f4'],
> 'offsets': [0, 16, 24]})
> dtype([('a', '|S4'), ('', '|V12'), ('b', '<f8'), ('', '|V12'), ('c',
> '<f4')])
> 
> There should be no gap between "b" and "c" but still the repr has ('',
> '|V12')
> between them. Am I missing something?

For performance reasons, compilers will typically align integers (and
probably floats) on 4-byte boundaries and apparently, doubles on 16-byte
boundaries.

Because compilers align like this, so does NumPy. This allows you to:

1. Take any kind of C struct definition
2. Convert it to a dtype
3. Create a NumPy array with this dtype
4. Pass the array's data pointer to C code
5. Cast the data pointer to a pointer to your C struct
6. Operate on the pointer to struct as if it were allocated in C

Cheers,

Albert

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailing unused bytes

[Martin W. <mar...@gm...>]

On Friday 15 September 2006 20:05, Travis Oliphant wrote:
> Martin Wiechert wrote:
> > Hi list,
> >
> > I'm using PyArray_DescrConverter with a dict object to create a
> > "struct-like" dtype from C.
> > As the struct contains different data types I run into "unaligned access"
> > problems.
> > Is there a way to force alignment or to get trailing unused bytes in the
> > dtpye?
>
> Besides calling the tp_new function (or the method __new__),  I added
> PyArray_DescrAlignConverter to the C-API.
>
> It is at the end and so does not require re-compilation of extension
> modules (unless you want to use it...)
>

Thanks Travis. I appreciate your work.

> -Travis
>
>
> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job
> easier Download IBM WebSphere Application Server v.1.0.1 based on Apache
> Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailing unused bytes

[Travis O. <oli...@ie...>]

Martin Wiechert wrote:
> Hi list,
>
> I'm using PyArray_DescrConverter with a dict object to create a "struct-like" 
> dtype from C.
> As the struct contains different data types I run into "unaligned access" 
> problems.
> Is there a way to force alignment or to get trailing unused bytes in the 
> dtpye?
>   
Besides calling the tp_new function (or the method __new__),  I added 
PyArray_DescrAlignConverter to the C-API.

It is at the end and so does not require re-compilation of extension 
modules (unless you want to use it...)

-Travis

Re: [Numpy-discussion] =?iso-8859-1?q?PyArray=5FDescrConverter_-_alig?= =?iso-8859-1?q?nment_/=09trailingunused_bytes?=

[Martin W. <mar...@gm...>]

On Friday 15 September 2006 19:14, Travis Oliphant wrote:
> Martin Wiechert wrote:
> > Thanks Albert! Do you also know the corresponding C-API function? It
> > cannot be PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose
> > signature has no "align", right?
>
> The DescrConverter function is meant for "O&"-style conversions.  It
> can't accept an align function.  We could possibly add something to the
> converter to allow specification of alignment through the object to be
> converted.
>

I begin to see the light....

For dictionaries one could maybe just add an optional key "align".
Also an optional key "elsize" or "itemsize" to force the total size of the 
record may sometimes be useful. E.g. one may want to faithfully map a given C 
struct. (That's why I was asking for trailing unused bytes.)

But I understand that other things have higher priority.

> Or, you can just call the __new__ method of the PyArrayDescr_Type object
>
> res = PyObject_CallMethod((PyObject *)&PyArrayDescr_Type, "__new__",
> "Oi", dict_object, 1))
>
> or call the tp->new method directly:
>
> args = Py_BuildValue("Oi", dict_object, 1);
> PyArrayDescr_Type->tp_new(&PyArrayDescr_Type, args, NULL);
> Py_DECREF(args);
>

Thank you! I'll try this.

> (I think passing in NULL for the keywords is O.K., but I haven't checked
> it).
>
> -Travis
>

One final question. To me the repr of a dtype with gaps is a little bit 
puzzling:

>>> dtype ({'names': ['a', 'b', 'c'], 'formats': ['<a4', '<f8', '<f4'], 
'offsets': [0, 16, 24]})
dtype([('a', '|S4'), ('', '|V12'), ('b', '<f8'), ('', '|V12'), ('c', '<f4')])

There should be no gap between "b" and "c" but still the repr has ('', '|V12') 
between them. Am I missing something?


Anyway, thanks again for your help.

Cheers,
Martin

>
> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job
> easier Download IBM WebSphere Application Server v.1.0.1 based on Apache
> Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Travis O. <oli...@ie...>]

Martin Wiechert wrote:
> Thanks Albert! Do you also know the corresponding C-API function? It cannot be 
> PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose signature has no 
> "align", right?
>   
The DescrConverter function is meant for "O&"-style conversions.  It 
can't accept an align function.  We could possibly add something to the 
converter to allow specification of alignment through the object to be 
converted.

Or, you can just call the __new__ method of the PyArrayDescr_Type object

res = PyObject_CallMethod((PyObject *)&PyArrayDescr_Type, "__new__", 
"Oi", dict_object, 1))

or call the tp->new method directly:

args = Py_BuildValue("Oi", dict_object, 1);
PyArrayDescr_Type->tp_new(&PyArrayDescr_Type, args, NULL);
Py_DECREF(args);

(I think passing in NULL for the keywords is O.K., but I haven't checked 
it).

-Travis

Re: [Numpy-discussion] how to get info about internals of an array object ?

[Travis O. <oli...@ie...>]

Sebastian Haase wrote:
> Hi,
> what I'm asking is if numpy has an equivalent to numarray's info() function:
>   
>>>> na.arange(10).info()
>>>>         
numpy.numarray.info(numpy.arange(10))

(get recent SVN as there were some bugs just fixed.

-Travis

[Numpy-discussion] buggy buggy bugyy: format and casting ?

[Eric E. <ems...@ob...>]

Hi,

I am facing a rather frustrating problem with numpy/scipy: after
upgrading to svn numpy and scipy, and trying to remove most of the now
unnecessary casting (floats) the program I wrote does not give me the
right answer. It seems that the answer is somewhat scaled down (but not
in a simple way). I had this kind of behaviour some time ago when scipy
and numpy were not fully compatible and then I had to cast the floats
explicitely in my program to make it work. But now it SHOULD be fully
compatible as far as I understand the issue. So my question is:
is there any remaining flaws in terms of casts/format between scipy and
numpy? I am specifically using the scipy functions:
- special.erf, special.erfc, orthogonal.ps_roots, integrate.quad...
and the numpy functions : sum, sqrt, exp, sin, cos, arctan.
I am doing all the calculations using numpy.float32.

Did anybody witness weird behaviour recently in this context?
(I don't really know what to do now since I already spent hours trying
to find the problem: the point is that the final answer is wrong for
sure, but the intermediate calculations are not easily checked since I
do not know what to really expect there. The program was working before
I touched the casting + upgraded numpy/scipy)

thanks for any help here, and cheers.

Eric

Re: [Numpy-discussion] how to get info about internals of an arrayobject ?

[Travis O. <oli...@ie...>]

Albert Strasheim wrote:
>> -----Original Message-----
>> From: num...@li... [mailto:numpy-
>> dis...@li...] On Behalf Of Sebastian Haase
>> Sent: 15 September 2006 03:21
>> To: numpy-discussion
>> Subject: [Numpy-discussion] how to get info about internals of an
>> arrayobject ?
>>
>> Hi,
>> what I'm asking is if numpy has an equivalent to numarray's info()
>> function:
>>     
>
>   
>> <snip>
>>     
>
>   
>> This was always helpful to me when debugging C binding code.
>>
>> Especially I'm asking if there is any way to get the memory address of an
>> array - for debugging purposes only - of course ;-)
>>     
>
> numpy.array([]).__array_interface__['data'][0]
>   

Or

numpy.ctypes._as_parameter_

(should work even if you don't have ctypes installed)

-Travis

Re: [Numpy-discussion] Experience with Visit?

[O'Keefe, M. <Mic...@nr...>]

I haven't tried VisIT before but thanks for the link. I also downloaded =
and am checking it out.

Along this same line of discussion, has anyone tried OOF2 which is an =
FEA package that also has some strong python connections?

http://www.ctcms.nist.gov/oof/oof2.html

I'm working on a Windows machine and it the current code-base doesn't =
seem to support Windows out of the box (if at all). Looks like you can =
put it together for *nix or Mac OS X, though...

--Michael=20

> -----Original Message-----
> From: num...@li...=20
> [mailto:num...@li...] On=20
> Behalf Of Robert Cimrman
> Sent: Friday, September 15, 2006 8:29
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] Experience with Visit?
>=20
> Travis Oliphant wrote:
> > Has anybody had any experience with the 3-D visualization software=20
> > VISIT?   It has Python bindings and seems to be pretty=20
> sophisticated. =20
> > I'm wondering why I haven't heard more about it.
> >=20
> > http://www.llnl.gov/visit/
>=20
> No reaction up to now, so...
>=20
> I have just tried the 'getting started' part and was quite impressed,=20
> thanks for posting the link! Up to now I have used ParaView=20
> and was very=20
> satisfied, but the Python bindings of VisIt are a great lure.
>=20
> r.
>=20
> --------------------------------------------------------------
> -----------
> Using Tomcat but need to do more? Need to support web=20
> services, security?
> Get stuff done quickly with pre-integrated technology to make=20
> your job easier
> Download IBM WebSphere Application Server v.1.0.1 based on=20
> Apache Geronimo
> http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&
> dat=3D121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>=20

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Martin W. <mar...@gm...>]

On Friday 15 September 2006 17:53, Albert Strasheim wrote:
> Hello all
>
> In [1]: import numpy as N
>
> In [3]: N.dtype({'names' : ['x', 'y'],
>                  'formats' : [N.intc, N.float64]},
>                 align=True)
> Out[3]: dtype([('x', '<i4'), ('', '|V4'), ('y', '<f8')])
>
> The reason you might not have discovered this:
>
> In [2]: N.dtype?
> Type:           type
> Base Class:     <type 'type'>
> String Form:    <type 'numpy.dtype'>
> Namespace:      Interactive
> File:           c:\python24\lib\site-packages\numpy\__init__.py
> Docstring:
>     <no docstring>
>

Thanks Albert! Do you also know the corresponding C-API function? It cannot be 
PyArray_DescrConverter (PyObject *, PyArray_Descr **), whose signature has no 
"align", right?

> Cheers,
>
> Albert
>
> > -----Original Message-----
> > From: num...@li... [mailto:numpy-
> > dis...@li...] On Behalf Of Martin Wiechert
> > Sent: 15 September 2006 16:14
> > To: numpy-discussion
> > Subject: [Numpy-discussion] PyArray_DescrConverter - alignment /
> > trailingunused bytes
> >
> > Hi list,
> >
> > I'm using PyArray_DescrConverter with a dict object to create a "struct-
> > like"
> > dtype from C.
> > As the struct contains different data types I run into "unaligned access"
> > problems.
> > Is there a way to force alignment or to get trailing unused bytes in the
> > dtpye?
> >
> > Thanks, Martin
>
> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job
> easier Download IBM WebSphere Application Server v.1.0.1 based on Apache
> Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailingunused bytes

[Albert S. <fu...@gm...>]

Hello all

In [1]: import numpy as N

In [3]: N.dtype({'names' : ['x', 'y'],
                 'formats' : [N.intc, N.float64]},
                align=True)
Out[3]: dtype([('x', '<i4'), ('', '|V4'), ('y', '<f8')])

The reason you might not have discovered this:

In [2]: N.dtype?
Type:           type
Base Class:     <type 'type'>
String Form:    <type 'numpy.dtype'>
Namespace:      Interactive
File:           c:\python24\lib\site-packages\numpy\__init__.py
Docstring:
    <no docstring>

Cheers,

Albert

> -----Original Message-----
> From: num...@li... [mailto:numpy-
> dis...@li...] On Behalf Of Martin Wiechert
> Sent: 15 September 2006 16:14
> To: numpy-discussion
> Subject: [Numpy-discussion] PyArray_DescrConverter - alignment /
> trailingunused bytes
> 
> Hi list,
> 
> I'm using PyArray_DescrConverter with a dict object to create a "struct-
> like"
> dtype from C.
> As the struct contains different data types I run into "unaligned access"
> problems.
> Is there a way to force alignment or to get trailing unused bytes in the
> dtpye?
> 
> Thanks, Martin

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailing unused bytes

[Martin W. <mar...@gm...>]

On Friday 15 September 2006 16:40, Francesc Altet wrote:
> A Divendres 15 Setembre 2006 16:13, Martin Wiechert va escriure:
> > Hi list,
> >
> > I'm using PyArray_DescrConverter with a dict object to create a
> > "struct-like" dtype from C.
> > As the struct contains different data types I run into "unaligned access"
> > problems.
> > Is there a way to force alignment or to get trailing unused bytes in the
> > dtpye?
>
> One possible solution is to declare void ('V' charcode) types for filling
> the gaps. For example:
>
>
> In [118]: ra=numpy.rec.array("1"*300,
> dtype=[('sval','<a4'),('unused','|V4'), ('dval','<f8')], shape=2)
>
> In [119]: ra['sval']
> Out[119]:
> recarray([1111, 1111],
>       dtype='|S4')
>
> In [120]: ra['dval']
> Out[120]: array([  9.73041595e-72,   9.73041595e-72])
>

Hm, not exactly elegant but should work. Thanks Francesc!

Note that all but trailing gaps can be controlled using a dictionary 
descriptor with an offsets key:

>>> ra = zeros ((10,), dtype = {'names': ['sval', 'dval'], 'formats': ['<a4', 
'<f8'], 'offsets': [0, 8]})
>>> ra.dtype
dtype([('sval', '|S4'), ('', '|V4'), ('dval', '<f8')])

Still if anybody knows a clean solution I would be very interested.

> You can still access the empty spaces if you want (although it is
> nonsense):
>
> In [121]: ra['unused']
> Out[121]:
> recarray([1111, 1111],
>       dtype='|V4')
>
> Cheers,

Re: [Numpy-discussion] PyArray_DescrConverter - alignment / trailing unused bytes

[Francesc A. <fa...@ca...>]

A Divendres 15 Setembre 2006 16:13, Martin Wiechert va escriure:
> Hi list,
>
> I'm using PyArray_DescrConverter with a dict object to create a
> "struct-like" dtype from C.
> As the struct contains different data types I run into "unaligned access"
> problems.
> Is there a way to force alignment or to get trailing unused bytes in the
> dtpye?

One possible solution is to declare void ('V' charcode) types for filling t=
he=20
gaps. For example:


In [118]: ra=3Dnumpy.rec.array("1"*300, dtype=3D[('sval','<a4'),('unused','=
|V4'),
('dval','<f8')], shape=3D2)

In [119]: ra['sval']
Out[119]:
recarray([1111, 1111],
      dtype=3D'|S4')

In [120]: ra['dval']
Out[120]: array([  9.73041595e-72,   9.73041595e-72])

You can still access the empty spaces if you want (although it is nonsense):

In [121]: ra['unused']
Out[121]:
recarray([1111, 1111],
      dtype=3D'|V4')

Cheers,

=2D-=20
>0,0<   Francesc Altet =A0 =A0 http://www.carabos.com/
V   V   C=E1rabos Coop. V. =A0=A0Enjoy Data
 "-"