Re: [Numpy-discussion] Large symmetrical matrix

2008-06-11 Thread Peter Skomoroch
You could do something like this, where MyDIstanceFunc is written with weave and some intelligent caching. This only computes the upper diagonal elements in a sparse matrix: self.NumDatapoints = len(self.DataPoints) self.Distances = sparse.lil_matrix((self.NumDatapoints,self.NumDatapoints)) for

Re: [Numpy-discussion] eigenvector and eigenface

2008-03-11 Thread Peter Skomoroch
see this page I found in my del.icio.us links, sorry I forgot to mention it at the time of the thread: http://www.owlnet.rice.edu/~elec301/Projects99/faces/code.html All the best On Mon, Mar 10, 2008 at 1:37 AM, royG <[EMAIL PROTECTED]> wrote: > friends > I am learning eigenfaces using numpy .

Re: [Numpy-discussion] confusion about eigenvector

2008-03-11 Thread Peter Skomoroch
I found this in my del.icio.us links, sorry I forgot to mention it at the time: http://www.owlnet.rice.edu/~elec301/Projects99/faces/code.html All the best On Thu, Mar 6, 2008 at 10:39 AM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > ok..I coded everything again from scratch..looks like i was

Re: [Numpy-discussion] eigenvector and eigenface

2008-03-09 Thread Peter Skomoroch
See this thread: http://www.mail-archive.com/numpy-discussion@scipy.org/msg06877.html On Mon, Mar 10, 2008 at 1:37 AM, royG <[EMAIL PROTECTED]> wrote: > friends > I am learning eigenfaces using numpy . i use data from N images and > create eigenvectors to get a 'sorted eigenvectors' array of si

Re: [Numpy-discussion] PCA on set of face images

2008-02-29 Thread Peter Skomoroch
b^H * vl[:,i] > > where a^H denotes transpose(conjugate(a)). > On Sat, Mar 1, 2008 at 12:41 AM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > > On Mar 1, 12:57 am, "Peter Skomoroch" wrote: > I think > > > matlab example should be easy to transla

Re: [Numpy-discussion] PCA on set of face images

2008-02-29 Thread Peter Skomoroch
Forgot the url: http://www.cis.hut.fi/Opinnot/T-61.2010/harjoitustyo_en07.shtml On Fri, Feb 29, 2008 at 2:56 PM, Peter Skomoroch <[EMAIL PROTECTED]> wrote: > Here is the page I referenced for the octave version ... it includes > examples very similar to what you want. I will be po

Re: [Numpy-discussion] PCA on set of face images

2008-02-29 Thread Peter Skomoroch
nfaces clear Z; for i=1:length(eigensorted) Z(:,:,1,i)=reshape(eigensorted(:,i)-1.5*min(min(min(eigensorted))), 19,19); end montage(Z) %show top 16 eigenfaces clear Z; for i=1:16 Z(:,:,1,i)=reshape(eigensorted(:,i)-min(min(min(eigensorted))), 19,19); end montage(Z) On Fri, Feb 29, 2008 at 2:

Re: [Numpy-discussion] PCA on set of face images

2008-02-29 Thread Peter Skomoroch
RoyG, The timing of your question couldn't be better, I just did an blog post on this (I also plugged scipy and the EPD): http://www.datawrangling.com/python-montage-code-for-displaying-arrays.html The code basically replicates the matlab montage() function and approach to handling grayscale ima

Re: [Numpy-discussion] Best way to run python parallel

2007-03-29 Thread Peter Skomoroch
If you want to use PyMPI or PyPar, I'm writing a series of tutorials on how to get them running on Amazon EC2, http://www.datawrangling.com/on-demand-mpi-cluster-with-python-and-ec2-part-1-of-3.html I'm using PyMPI on a 20 node EC2 cluster and everything seems groovy, but I'm relatively new to