Hi Timothy,
On Fri, Feb 29, 2008 at 12:16 AM, Timothy Hochberg wrote:
> On Thu, Feb 28, 2008 at 1:47 PM, Andrea Gavana <[EMAIL PROTECTED]>
> wrote:
> > Hi All,
> >
> >I have some problems in figuring out a solution for an issue I am
> > trying to solve. I have a 3D grid of dimension Nx, Ny, Nz
On Thu, Feb 28, 2008 at 06:55:11PM -0600, Robert Kern wrote:
> On Thu, Feb 28, 2008 at 6:43 PM, Geoffrey Irving <[EMAIL PROTECTED]> wrote:
> > > The magic is in In[27]. We reshape the array of vectors to be
> > > compatible with the shape of the array of matrices. When we multiply
> > > the two
On Thu, Feb 28, 2008 at 6:43 PM, Geoffrey Irving <[EMAIL PROTECTED]> wrote:
> > The magic is in In[27]. We reshape the array of vectors to be
> > compatible with the shape of the array of matrices. When we multiply
> > the two together, it is as if we multiplied two (n,3,3) matrices, the
> > la
On Thu, Feb 28, 2008 at 05:57:29PM -0600, Robert Kern wrote:
> On Thu, Feb 28, 2008 at 4:34 PM, Geoffrey Irving <[EMAIL PROTECTED]> wrote:
> > Hello,
> >
> > I have a large number of points (shape (n,3)), and a matching
> > number of 3x3 matrices (shape (n,3,3)), and I want to compute
> > the pr
On Thu, Feb 28, 2008 at 4:34 PM, Geoffrey Irving <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I have a large number of points (shape (n,3)), and a matching
> number of 3x3 matrices (shape (n,3,3)), and I want to compute
> the product of each matrix times the corresponding point.
>
> I can't see a wa
On Thu, Feb 28, 2008 at 1:47 PM, Andrea Gavana <[EMAIL PROTECTED]>
wrote:
> Hi All,
>
>I have some problems in figuring out a solution for an issue I am
> trying to solve. I have a 3D grid of dimension Nx, Ny, Nz; for every
> cell of this grid, I calculate the cell centroids (with the cell
> c
On Thu, Feb 28, 2008 at 3:37 PM, Christian Heimes <[EMAIL PROTECTED]> wrote:
> Stuart Brorson wrote:
> > Also, what do these specs say about 0^?
>
> See for yourself
> http://www.open-std.org/jtc1/sc22/wg14/www/docs/n1124.pdf The
> interesting information are in Annex F.9 and Annex G.6.
Thank
Hello,
I have a large number of points (shape (n,3)), and a matching
number of 3x3 matrices (shape (n,3,3)), and I want to compute
the product of each matrix times the corresponding point.
I can't see a way to do this operation with dot or tensordot,
since these routines either sum across an inde
Stuart Brorson wrote:
> math.pow(0, -1)
>> Traceback (most recent call last):
>> File "", line 1, in
>> ValueError: math domain error
>
> Why isn't this one inf?
The standard says return inf and raise a divide-by-zero floating point
exception. Since we can't do both in Python we sticked to
On Thu, Feb 28, 2008 at 2:11 PM, Robert Kern <[EMAIL PROTECTED]> wrote:
> On Thu, Feb 28, 2008 at 2:00 PM, Sameer DCosta <[EMAIL PROTECTED]> wrote:
> > was going to create a ticket for this, but I don't think I have
> > permissions to do this. Can someone who has the permissions, please
> >
Hi All,
I have some problems in figuring out a solution for an issue I am
trying to solve. I have a 3D grid of dimension Nx, Ny, Nz; for every
cell of this grid, I calculate the cell centroids (with the cell
coordinates x, y, and z) and then I try to find which cell centroid is
the closest to
> Please checkout Mark Dickinson's and my trunk-math branch of Python 2.6.
> We have put lots of effort into fixing edge cases of floats, math and
> cmath functions. The return values are either based on the latest
> revision of IEEE 754 or the last public draft of the C99 standard (1124,
> Annex F
Stuart Brorson wrote:
> I have been poking at the limits of NumPy's handling of powers of
> zero. I find some results which are disturbing, at least to me.
> Here they are:
[SNIPP]
Please checkout Mark Dickinson's and my trunk-math branch of Python 2.6.
We have put lots of effort into fixing ed
On Thu, Feb 28, 2008 at 2:00 PM, Sameer DCosta <[EMAIL PROTECTED]> wrote:
> I think having a record array field renaming api is a good idea.. I
> was going to create a ticket for this, but I don't think I have
> permissions to do this. Can someone who has the permissions, please
> create it? Tha
I think having a record array field renaming api is a good idea.. I
was going to create a ticket for this, but I don't think I have
permissions to do this. Can someone who has the permissions, please
create it? Thanks.
Sameer
On Thu, Feb 28, 2008 at 10:48 AM, Travis E. Oliphant
<[EMAIL PROTECTED]
That's very strange, I've made a local installation of numpy 1.0.4
and scipy 0.6 from official source, and I use the same kind of
installation script as in
http://www.scipy.org/Installing_SciPy/Linux/BuildingFromSource/GCC_1,
so I always erase all the local directory to prevent this kind of
confl
Hello,
I'm using svn numpy and get the following error upon executing
/usr/local/bin/python2.5 setup.py config --noisy
--cc=/opt/intel/cce/10.0.025/bin/icc --compiler=intel --fcompiler=intel
build_clib build_ext
I see:
conv_template:> build/src.linux-x86_64-2.5/numpy/core/src/scalartypes.in
Lisandro Dalcin wrote:
> On 2/27/08, Travis E. Oliphant <[EMAIL PROTECTED]> wrote:
>> Did this discussion resolve with a fix that can go in before 1.0.5 is
>> released?
>
> I believe the answer is yes, but we have to choose:
>
> 1- Use the regepx based solution of David.
A good idea, but a feat
Timothy Hochberg wrote:
> I suppose one could use NaN + j NaN, as some have proposed, but is seems
> unnecessary;
Except don't we want the result to be a complex number (in terms of data
storage) so it can be part of a complex array?
-Chris
--
Christopher Barker, Ph.D.
Oceanographer
Emergen
Loic,
I was not able to reproduce this. Have you tried doing a clean
install by removing your numpy and scipy directories in the
site-packages and reinstalling. I've had old files in my install
directory cause problems in the past.
You are correct, the tests should not effect one another.
Chri
On Wed, Feb 27, 2008 at 4:10 PM, Stuart Brorson <[EMAIL PROTECTED]> wrote:
> I have been poking at the limits of NumPy's handling of powers of
> zero. I find some results which are disturbing, at least to me.
> Here they are:
[SNIP]
>
> ** 0^(x+y*i): This one is tricky; please bear with me
Sameer DCosta wrote:
> Hi,
>
> I'm having trouble renaming record array fields if they contain object
> arrays in them. I followed the solutions posted by Robert Kern and
> Stefan van der Walt (Thanks again) but it doesn't look like this
> method works in all cases. For reference:
> http://project
Hi,
I'm having trouble renaming record array fields if they contain object
arrays in them. I followed the solutions posted by Robert Kern and
Stefan van der Walt (Thanks again) but it doesn't look like this
method works in all cases. For reference:
http://projects.scipy.org/pipermail/numpy-discus
On Thu, Feb 28, 2008 at 7:35 AM, Neal Becker <[EMAIL PROTECTED]> wrote:
> I tried numpyx.pyx with cython-0.9.6.12.
These were written for and still work with Pyrex. If it doesn't work
with Cython then that is either a bug in Cython or an intentional
incompatibility of Cython.
--
Robert Kern
"I
I'd drive to Austin and wash and wax their cars if there was an OSX
distribution.
OK, mild exaggeration but I'd buy them a beer or two.
On Feb 28, 2008, at 9:18 AM, Steve Lianoglou wrote:
>> On Wed, 27 Feb 2008, Travis Vaught apparently wrote:
>>> http://www.enthought.com/epd
>>
>> Looks good.
> On Wed, 27 Feb 2008, Travis Vaught apparently wrote:
>> http://www.enthought.com/epd
>
> Looks good.
> An increasing number of my students are buying Macs,
> so the OSX support will be very welcome.
Yeah ... agreed, this is great!
-steve
___
Numpy-disc
On Thu, Feb 28, 2008 at 3:41 PM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> > Arnar wrote
>
> > from scipy import linalg
> > facearray-=facearray.mean(0) #mean centering
> > u, s, vt = linalg.svd(facearray, 0)
> > scores = u*s
> > facespace = vt.T
>
> hi Arnar
> when i do this i get these
> Arnar wrote
> from scipy import linalg
> facearray-=facearray.mean(0) #mean centering
> u, s, vt = linalg.svd(facearray, 0)
> scores = u*s
> facespace = vt.T
hi Arnar
when i do this i get these
u =< 'numpy.core.defmatrix.matrix'> (4, 4)
that matches the eigenvectors matrix in my previous data
s=
On 2/27/08, Travis E. Oliphant <[EMAIL PROTECTED]> wrote:
> Did this discussion resolve with a fix that can go in before 1.0.5 is
> released?
I believe the answer is yes, but we have to choose:
1- Use the regepx based solution of David.
2- Move to use 'index' instead of 'find' as proposed by Al
On Thu, Feb 28, 2008 at 8:17 AM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> i all
> I am learning PCA method by reading up Turk&Petland papers etc
> while trying out PCA on a set of greyscale images using python, and
> numpy I tried to create eigenvectors and facespace.
>
> i have
> faces
I tried numpyx.pyx with cython-0.9.6.12. Here's what I got:
In [2]: import numpyx
In [3]: numpyx.test()
-=-=-=-=-=-=-=-=-=-=
printing array info for ndarray at 0x0
print number of dimensions: 0
address of strides: 0x0
strides:
memory dump:
-1e-30
-=-=-=-=-=-=-=-=-=-=
-=-=-=-=-=-=-=-=-=-=
Travis E. Oliphant wrote:
> Neal Becker wrote:
>> Travis E. Oliphant wrote:
>>
>>
>>
>>
>> The code for this is a bit hard to understand. It does appear that it
>> only
>> searches for a conversion on the 2nd argument. I don't think that's
>> desirable behavior.
>>
...
>
> Thus, what is miss
>> ** 0^0: This is problematic.
>
> Accessible discussion:
> http://en.wikipedia.org/wiki/Exponentiation#Zero_to_the_zero_power>
Thanks. That was quite informative. Indeed, I communicated with a
math professor at MIT who also more or less convinced me that 0^0 = 1.
Stuart
OK, what you are getting are not the eigenvectors of you data, but the
eigenvectors of the transposition of your data (I suppose).
You have two options :
- either you make an eigen analysis of your data and get 12 eigenvectors
- either you make an eigen analysis of the transposition of your data an
On Feb 28, 1:27 pm, "Matthieu Brucher" wrote
> If your images are 4x3, your eigenvector must be 12 long.
hi
thanx for reply
i am using 4 images each of size 4X3
the covariance matrix obtained from adjfaces*faces_trans is 4X4 in
size and that produces the evalues and eigenvectors given here
eva
Hi,
If your images are 4x3, your eigenvector must be 12 long.
Matthieu
2008/2/28, [EMAIL PROTECTED] <[EMAIL PROTECTED]>:
>
> i all
> I am learning PCA method by reading up Turk&Petland papers etc
> while trying out PCA on a set of greyscale images using python, and
> numpy I tried to create eige
36 matches
Mail list logo
