<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.0.6487.1">
<TITLE>Problem with trace function?</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P><FONT SIZE=2 FACE="Arial">Hi all,</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">I have some trouble with the new numarray library.</FONT>
<BR><FONT SIZE=2 FACE="Arial">Most of the time the numarray code prints out some negative value whereas the</FONT>
<BR><FONT SIZE=2 FACE="Arial">Numeric code yields positive values.</FONT>
<BR><FONT SIZE=2 FACE="Arial">This behaviour justs shows up when using large arrays.</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">Can anybody point out what I need to change in the numarray version to get the Numeric behaviour?</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">System: WinXP, python2.3.4, Numeric23.5, numarray1.1.1</FONT>
<BR><FONT SIZE=2 FACE="Arial">-----</FONT>
<BR><FONT SIZE=2 FACE="Arial">NUMARRAY-VERSION:</FONT>
<BR><FONT SIZE=2 FACE="Arial">from numarray import greater, reshape, trace</FONT>
<BR><FONT SIZE=2 FACE="Arial">from numarray.random_array import standard_normal</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">n=1000</FONT>
<BR><FONT SIZE=2 FACE="Arial">A = greater(standard_normal(n*n), 0.9)</FONT>
<BR><FONT SIZE=2 FACE="Arial">A = reshape(A, (n, n))</FONT>
<BR><FONT SIZE=2 FACE="Arial">print trace(A)</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">-----</FONT>
<BR><FONT SIZE=2 FACE="Arial">NUMERIC-Version:</FONT>
<BR><FONT SIZE=2 FACE="Arial">from Numeric import greater, reshape, trace</FONT>
<BR><FONT SIZE=2 FACE="Arial">from RandomArray import standard_normal</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">n=1000</FONT>
<BR><FONT SIZE=2 FACE="Arial">A = greater(standard_normal(n*n), 0.9)</FONT>
<BR><FONT SIZE=2 FACE="Arial">A = reshape(A, (n, n))</FONT>
<BR><FONT SIZE=2 FACE="Arial">print trace(A)</FONT>
<BR><FONT SIZE=2 FACE="Arial">-----</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">Thank you for your help.</FONT>
<BR><FONT SIZE=2 FACE="Arial">Gerald</FONT>
</P>
<BR>
</BODY>
</HTML>