Re: [Numpy-discussion] Numpy/Fortran puzzle (?)

Hi Andrea You should add a dimension argument in your Fortran code, also you should write a f2py header in the same Fortran code. Remember, numpy memory is C order wise. You can specify in numpy the ordering of the matrices you pass when you create them. F2py automatically deals with matrices ,

Re: [Numpy-discussion] what do I get if I build with MKL?

One major advantage you can have using mkl is installing numexpr compiling it with MLK. That's a strong suggestion to easily use mkl and go faster on common operations. Xavier On 20/04/2013 1:16 AM, Matthieu Brucher matthieu.bruc...@gmail.com wrote: The graph is a comparison of the dot calls,

Re: [Numpy-discussion] idea of optimisation?

Actually this can be a good idea. i didn't thought using he sorting. i'll try thanks for yours ideas Xavier 2011/12/7 Tony Yu tsy...@gmail.com On Tue, Dec 6, 2011 at 2:51 AM, Xavier Barthelemy xab...@gmail.comwrote: ok let me be more precise I have an Z array which is the elevation from

Re: [Numpy-discussion] idea of optimisation?

Yes I understood what you said. I know these tools, and I am using them. I was just wandering if someone has a more one-liner-pythonic way to do it. I don't think it's worth importing a new fortran module. Thanks anyway :) Xavier 2011/12/6 David Froger david.fro...@gmail.com Excerpts from

Re: [Numpy-discussion] ignore NAN in numpy.true_divide()

Hi, I don't know if it is the best choice, but this is what I do in my code: for each slice: indexnonNaN=np.isfinite(SliceOf Toto) SliceOf TotoWithoutNan= SliceOf Toto [indexnonNaN] and then perform all operation I want o on the last array. i hope it does answer your question Xavier

Re: [Numpy-discussion] ignore NAN in numpy.true_divide()

the TSFC_avg is:, TSFC_avg On Tue, Dec 6, 2011 at 9:50 AM, Xavier Barthelemy xab...@gmail.comwrote: Hi, I don't know if it is the best choice, but this is what I do in my code: for each slice: indexnonNaN=np.isfinite(SliceOf Toto) SliceOf TotoWithoutNan= SliceOf Toto [indexnonNaN

[Numpy-discussion] idea of optimisation?

Hi everyone I was wondering if there is a more optimal way to write what follows: I am studying waves, so I have an array of wave crests positions, Xcrest and the positions of the ZeroCrossings, Xzeros. The goal is to find between which Xzeros my xcrest are. XXX1=XCrest

Re: [Numpy-discussion] idea of optimisation?

ok let me be more precise I have an Z array which is the elevation from this I extract a discrete array of Zero Crossing, and another discrete array of Crests. len(crest) is different than len(Xzeros). I have a threshold method to detect my valid crests, and sometimes there are 2 crests between

Re: [Numpy-discussion] compiling and linking MKL with numpy/scipy

architecture, in case it might help. I always had trouble getting programs (other than numpy though) to link and execute properly with MKL. You might also try to play with LD_PRELOAD. Good luck, -=- Olivier 2011/11/23 Xavier Barthelemy xab...@gmail.com Dear everyone I have been a bit better

[Numpy-discussion] compiling and linking MKL with numpy/scipy

dear everyone, I am hard trying to compile numpy / and scipy with mkl. unfortunately it does not work. i have tried a lot of solution, and the closest for me to work is: *I have * Linux manning 2.6.32-27-server #49-Ubuntu SMP Thu Dec 2 02:05:21 UTC 2010 x86_64 GNU/Linux python -c 'import

Re: [Numpy-discussion] compiling and linking MKL with numpy/scipy

Dear everyone I have been a bit better patching the files I finally get something in .basrc source /opt/intel/composer_xe_2011_sp1.7.256/bin/compilervars.sh intel64 source /opt/intel/composer_xe_2011_sp1.7.256/mkl/bin/mklvars.sh intel64 lp64 export OMP_NUM_THREADS=4 export