[Numpy-discussion] ANN: PyCPX 0.01, a numpy/cython wrapper for CPlex.

2011-09-18 Thread Hoyt Koepke
. Thanks, --Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ ___ NumPy-Discussion

Re: [Numpy-discussion] Linear assignment problem: Jonker-Volgenant algorithm

2011-05-17 Thread Hoyt Koepke
, I don't know if I'm an expert in this area, but I rub shoulders with a few and I'm glad to be of help. -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com

Re: [Numpy-discussion] Linear assignment problem: Jonker-Volgenant algorithm

2011-05-16 Thread Hoyt Koepke
) has solvers for this problem and has python bindings. It's under the boost license (is that ok?). It might be a bit heavyweight for this, though, but it's great software. -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics

Re: [Numpy-discussion] Linear assignment problem: Jonker-Volgenant algorithm

2011-05-16 Thread Hoyt Koepke
of implementations out there, as it's often a standard coding project in undergrad algorithms courses. http://en.wikipedia.org/wiki/Hungarian_algorithm has links to a few. -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http

Re: [Numpy-discussion] Linear assignment problem: Jonker-Volgenant algorithm

2011-05-16 Thread Hoyt Koepke
be coded up efficiently in a few hundred lines of code. Hope that answers your questions :-). --Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com

Re: [Numpy-discussion] Compiling numpy using icc gets missing library error

2011-03-27 Thread Hoyt Koepke
a patch together for the distutils. However, I'd want someone to review it since I'm not that confident in my knowledge of the distutils code. I can also try to turn this into a more complete description for the wiki. -- Hoyt + Hoyt Koepke

Re: [Numpy-discussion] Compiling numpy using icc gets missing library error

2011-03-25 Thread Hoyt Koepke
wrong there? Thanks! -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ setupenv.sh Description: Bourne shell script

[Numpy-discussion] ANN: TreeDict 0.12, a lightweight library for hierarchical parameter/data management

2011-03-25 Thread Hoyt Koepke
the existing python representations of the matlab objects. Discussion/comments on this point are welcome. + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com

Re: [Numpy-discussion] Compiling numpy using icc gets missing library error

2011-03-25 Thread Hoyt Koepke
to my real work. -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ ___ NumPy

Re: [Numpy-discussion] Compiling numpy using icc gets missing library error

2011-03-25 Thread Hoyt Koepke
, 2011 at 2:30 PM, Pauli Virtanen p...@iki.fi wrote: On Fri, 25 Mar 2011 13:46:41 -0700, Hoyt Koepke wrote: [clip] Also, I could still not get the CloughTocher2DInterpolator to not segfault. Backtrace would be useful here. It's probably best to recompile with -O0 and some debug flags enabled

[Numpy-discussion] Compiling numpy using icc gets missing library error

2011-03-22 Thread Hoyt Koepke
something could have gone wrong? Thanks! -- Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ build.log.gz Description: GNU

[Numpy-discussion] ndarray special methods should raise NotImplemented instead of TypeError

2011-03-13 Thread Hoyt Koepke
+ Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http

Re: [Numpy-discussion] ndarray special methods should raise NotImplemented instead of TypeError

2011-03-13 Thread Hoyt Koepke
what I was hoping for. I assume in my case, I just need __array_priority__ 10 to gain priority over a matrix (from reading http://docs.scipy.org/doc/numpy/reference/arrays.classes.html). Is that correct? Thanks! --Hoyt + Hoyt Koepke

Re: [Numpy-discussion] Slicing slower than matrix multiplication?

2009-12-11 Thread Hoyt Koepke
+ Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com ++ ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org

Re: [Numpy-discussion] Faster way to generate a rotation matrix?

2009-03-03 Thread Hoyt Koepke
at optimizing that kind of thing too. --Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com

Re: [Numpy-discussion] how to get the corresponding eigenvector for a specific eigen value?

2009-02-08 Thread Hoyt Koepke
afterwards. Quick and dirty, but should work. --Hoyt -- + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + hoy...@gmail.com

Re: [Numpy-discussion] glibc error

2009-01-24 Thread Hoyt Koepke
___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- + Hoyt Koepke + University of Washington Department of Statistics + http

Re: [Numpy-discussion] glibc memory corruption when running numpy.test()

2008-11-20 Thread Hoyt Koepke
that it doesn't hang. I still get a number of test failures, however (build + test logs attached). Thanks a lot for the help! --Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + [EMAIL

Re: [Numpy-discussion] glibc memory corruption when running numpy.test()

2008-11-20 Thread Hoyt Koepke
-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- + Hoyt Koepke + University of Washington Department of Statistics + http://www.stat.washington.edu/~hoytak/ + [EMAIL PROTECTED

[Numpy-discussion] glibc memory corruption when running numpy.test()

2008-11-19 Thread Hoyt Koepke
, but they don't seem to be at first glance. I've attached the full log. If there is anything more you want me to do with this, I'd be happy to. Thanks! --Hoyt + Hoyt Koepke + University of Washington Department of Statistics + http

Re: [Numpy-discussion] glibc memory corruption when running numpy.test()

2008-11-19 Thread Hoyt Koepke
Sorry; I also added '-fPIC' to the compile flags, though it may work without having it there. I just got errors related to not having it and so threw it at everything... --Hoyt On Wed, Nov 19, 2008 at 11:41 PM, Hoyt Koepke [EMAIL PROTECTED] wrote: Attached. I had a bunch of issues getting

Re: [Numpy-discussion] numpy, swig and TNT-Arrays

2008-11-11 Thread Hoyt Koepke
___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED

Re: [Numpy-discussion] problem installing numpy using scons

2008-10-25 Thread Hoyt Koepke
Excellent; thank you -- I will test it out soon. -- Hoyt On Sat, Oct 25, 2008 at 2:57 AM, David Cournapeau [EMAIL PROTECTED] wrote: On Sun, Oct 19, 2008 at 7:57 PM, David Cournapeau [EMAIL PROTECTED] wrote: The trunk requires some changes which are not released yet. It will be in the next

[Numpy-discussion] problem installing numpy using scons

2008-10-18 Thread Hoyt Koepke
Hello, I'm trying to install the latest numpy using setupscons.py, but it gives me an error that I can't track down. It's in the latest numpy release (5946). with the latest numpyscons from the 0.9.3 branch. Here's the last part of output: is bootstrapping ? True Executing scons command (pkg

Re: [Numpy-discussion] Slicing without a priori knowledge of dimension

2008-08-06 Thread Hoyt Koepke
a pointer as to what it is? Thanks in advance, Matt ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department

Re: [Numpy-discussion] Large symmetrical matrix

2008-06-11 Thread Hoyt Koepke
-- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED] +++ ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy

Re: [Numpy-discussion] segmentation fault

2008-05-28 Thread Hoyt Koepke
Virtanen ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak

Re: [Numpy-discussion] Quick Question about Optimization

2008-05-20 Thread Hoyt Koepke
,:] and theta[1-curidx] with theta[1-curidx, :]. (I may have missed some others.) weave.blitz is currently limited to just array operations... it doesn't really support the numpy functions. Hope that helps a little -- Hoyt -- +++ Hoyt Koepke UBC Department of Computer

Re: [Numpy-discussion] Quick Question about Optimization

2008-05-19 Thread Hoyt Koepke
, :]) wv.blitz(theta[curidx, :] = expfac_theta * theta[2-curidx] - (1-expfac_theta)) idx_spk = np.where(v = theta) S[n,idx_spk] = 1 theta[curidx, idx_spk] += b # Flop to handle previous stuff curidx = 2 - curidx -- +++ Hoyt Koepke UBC Department

Re: [Numpy-discussion] Quick Question about Optimization

2008-05-19 Thread Hoyt Koepke
then there is no copying. What am I missing? Actually, I think you are correct. My bad. I was mainly thinking in terms of weave.blitz, where it would make a difference, then translating back... --Hoyt +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca

Re: [Numpy-discussion] Faster

2008-05-04 Thread Hoyt Koepke
and then update the others only if the row you're updating has that minimum value in it. Then, when scanning for the min dist, you only need to scan O(n) rows. Sorry, let me clarify -- Update the entries corresponding to entries in the row you're updating if they are the same as the

Re: [Numpy-discussion] Faster

2008-05-03 Thread Hoyt Koepke
is the distance matrix. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http

Re: [Numpy-discussion] Faster

2008-05-03 Thread Hoyt Koepke
-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED

Re: [Numpy-discussion] types as functions convert 1 elm arrays to scalars

2008-05-01 Thread Hoyt Koepke
report? --Hoyt On Mon, Apr 28, 2008 at 11:51 PM, Charles R Harris [EMAIL PROTECTED] wrote: On Tue, Apr 29, 2008 at 12:28 AM, Hoyt Koepke [EMAIL PROTECTED] wrote: Hello, I have a quick question that I'm hoping will improve my numpy understanding. I noticed some behavior when using

Re: [Numpy-discussion] Distribution Functions Change Behavior

2008-04-28 Thread Hoyt Koepke
/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED] +++ ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http

Re: [Numpy-discussion] Distribution Functions Change Behavior

2008-04-28 Thread Hoyt Koepke
. --Hoyt On Mon, Apr 28, 2008 at 8:29 AM, Hoyt Koepke [EMAIL PROTECTED] wrote: I may not understand what you are asking, Rich, but I'm not sure I agree with Alan. A Gaussian fit to data x should fit exactly as well as data fit to ax, a 0, just with a variance a^2 times the original. The only

Re: [Numpy-discussion] Generating Bell Curves (was: Using normal() )

2008-04-25 Thread Hoyt Koepke
-- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED] +++ ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman

Re: [Numpy-discussion] access ndarray in C++

2008-04-23 Thread Hoyt Koepke
/index.html [2] http://www.riverbankcomputing.co.uk/sip/index.php ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC

Re: [Numpy-discussion] #734: interactive docstring search (lookfor)

2008-04-10 Thread Hoyt Koepke
Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion -- +++ Hoyt Koepke UBC Department of Computer Science http://www.cs.ubc.ca/~hoytak/ [EMAIL PROTECTED

Re: [Numpy-discussion] #734: interactive docstring search (lookfor)

2008-04-10 Thread Hoyt Koepke
In [1]: import numpy In [2]: numpy.info('eigvals') *** Found in numpy.linalg *** eigvals(a) Fair enough Don't know why I missed that, prob relied too much on Google search :-) Having it as part of iPython does make sense. --Hoyt ___

Re: [Numpy-discussion] Quikest way to create a symetric (diagonal???) matrix ?

2008-03-26 Thread Hoyt Koepke
If the rest of the matrix is already zeros and memory wasn't a problem, you could just use A_sym = A + A.T - diag(diag(A)) If memory was an issue, I'd suggest weave.inline (if that's a viable option) or pyrex to do the loop, which would be about as fast as you could get. --Hoyt On Wed, Mar

Re: [Numpy-discussion] dunno what array operation I'm looking for...

2008-03-21 Thread Hoyt Koepke
Try result = A[1:] - A[:-1] --Hoyt On Fri, Mar 21, 2008 at 7:43 PM, Chris Withers [EMAIL PROTECTED] wrote: Hi All, Say I have an array like: measurements = array([100,109,115,117]) What do I do to it to get: array([9, 6, 2]) Is the following really the best way? result =

[Numpy-discussion] numpy.random.RandomState threadsafe?

2008-03-11 Thread Hoyt Koepke
This should be a really quick question. Is a RandomState object thread safe? I'm wanting to use a common RandomState object in a multithreaded program, and I need to know if it's necessary to protect it with a lock (which wouldn't be difficult). Thanks! --Hoyt

Re: [Numpy-discussion] numpy.random.RandomState threadsafe?

2008-03-11 Thread Hoyt Koepke
Okay, thanks! I won't be using the multivariate_normal function in my code, so this should work fine. --Hoyt On Tue, Mar 11, 2008 at 3:00 PM, Robert Kern [EMAIL PROTECTED] wrote: On Tue, Mar 11, 2008 at 4:18 PM, Hoyt Koepke [EMAIL PROTECTED] wrote: This should be a really quick question

Re: [Numpy-discussion] Optimize speed of for loop using numpy

2008-02-25 Thread Hoyt Koepke
I would definitely suggest using scipy's weave.inline for this. It seems like this particular function can be translated into C code really easily, which would give you a HUGE speed up. Look at some of the examples in scipy/weave/examples to see how to do this. The numpy book also has a section