getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Getfem-users] Getfem python interface questions


From: Julien Pommier
Subject: Re: [Getfem-users] Getfem python interface questions
Date: Sun, 25 Nov 2007 22:51:41 +0100
User-agent: Internet Messaging Program (IMP) 3.2.6

Hi ShuntTim,

I think the problem is elsewhere as your expression is working just fine. An
example of python script using eval is the demo_plate.py one.

About the equivalence of types, yes you are correct. Do you run into specific
problems on AMD64 ?

Quoting LUK ShunTim <address@hidden>:

> Hello,
>
> I'm learning to use the python interface by trying to translate the
> M-files in the tests/matlab directory into their python counterparts.
>
> I start with tutorial1.m and things seem to go well until this line
>
> R=gf_mesh_fem_get(mf, 'eval', {'(x-.5).^2 + (y-.5).^2 + x/5 - y/3'});
>
> Here's the python equivalent I used:
> R=mf.eval('(x[0] - 0.5)**2 + (x[1] - 0.5)**2 + x[0]/5.0 - x[1]/3.0')
>
> I got this error traceback
>
> </error>
> ([], array([], type=Int32))
> Traceback (most recent call last):
>    File "tutorial1.py", line 50, in ?
>      R=mf.eval('(x[0] - 0.5)**2 + (x[1] - 0.5)**2 + x[0]/5.0 - x[1]/3.0')
>    File "/usr/lib/python2.4/site-packages/getfem/getfem.py", line 865,
> in eval
>      x=P[:,0]; r=numarray.array(eval(expression))
> IndexError: Index out of range
> </error>
>
> I think I don't know how to use the eval() function correctly.
>
> Another question: is my understanding of the equivalence of types in the
> matlab-interface document correct?
>
> Matlab                Python
>
> ivec          numarray rank-1 integer array (of default precision)
> vec           numarray rank-1 float array (of default precision)
> imat          numarray rank-2 integer array (of default precision)
> mat           numarray rank-2 float array (of default precision)
>
> On my debian/sid amd64, the default precision for array elements are
> Int64 and Float64, respectively.
>
> If anyone is interested, after I get them working, I'm happy to share
> them. But don't hold your breadth, :-) I just started my learning and I
> can only do it intermittently. I'm sure I have to seek your help again.
>
> Thanks in advance,
> ST
> --
>
>
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
>






reply via email to

[Prev in Thread] Current Thread [Next in Thread]