[SciPy-user] returning floats from a function
Perez Kawumi
pepe_kawumi at yahoo.co.uk
Sat Sep 8 07:15:42 EDT 2007
Hi,
im having problems returning floating point numbers in the whitney_function matrices at the end of my program. The answers are being rounded off to whole numbers. I think it might be the way im defining my whitney_functions initially(the first highlighted block). Please help.
Thanks Perez
def whitney(elem_num,x,y):
global ELEMENTS, NODE_COORD
whitney_functions = []
whitney_functions = zeros((3,2))
perez = simplex2D(elem_num,x,y)
trinodes = ELEMENTS[[elem_num],:]-1
x1 = NODE_COORD[trinodes[0,0],0]
y1 = NODE_COORD[trinodes[0,0],1]
x2 = NODE_COORD[trinodes[0,1],0]
y2 = NODE_COORD[trinodes[0,1],1]
x3 = NODE_COORD[trinodes[0,2],0]
y3 = NODE_COORD[trinodes[0,2],1]
#area = 0.5*abs(det([1 x1 y1; 1 x2 y2; 1 x3 y3]));
temp = linalg.inv(([x1,x2,x3],[y1,y2,y3],[1,1,1]))
b = temp[:,[0]]
c = temp[:,[1]]
nabla_lambda = temp[:,0:2]
whitney_functions[[0],:] = perez[0]*nabla_lambda[[1],:] - perez[1]*nabla_lambda[[0],:]
whitney_functions[[1],:] = perez[0]*nabla_lambda[[2],:] - perez[2]*nabla_lambda[[0],:]
whitney_functions[[2],:] = perez[1]*nabla_lambda[[2],:] - perez[2]*nabla_lambda[[1],:]
return whitney_functions
___________________________________________________________
Want ideas for reducing your carbon footprint? Visit Yahoo! For Good http://uk.promotions.yahoo.com/forgood/environment.html
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20070908/7b0214b1/attachment.html>
More information about the SciPy-User
mailing list