On 10 December 2014 at 20:36, Chris Barker wrote:
> On Tue, Dec 9, 2014 at 11:03 PM, Sturla Molden
> wrote:
>
>> Nathaniel Smith wrote:
>>
>> > @contextmanager
>> > def tmp_zeros(*args, **kwargs):
>> > arr = np.zeros(*args, **kwargs)
>> > try:
>> > yield arr
>> > finally:
>>
On 25 November 2014 at 19:33, David Cournapeau wrote:
>
>
> On Tue, Nov 25, 2014 at 6:10 PM, Sturla Molden
> wrote:
>
>> David Cournapeau wrote:
>> > Shall we consider > > href="https://github.com/scipy/scipy/issues/4168";>
>> https://github.com/scipy/scipy/issues/4168
>> > to be a
>> > blocker
Hi All,
I have the following (very ugly) line of code:
all_results = np.asarray([transm_hist[date_idx, :, idx_main_set[date_idx]
]*main_flow[date_idx, 0:n_fluids] for date_idx in xrange(n_dates)])
where transm_hist.shape = (n_dates, n_fluids, n_nodes), main_flow.shape =
(n_dates, n_fluids) a
On 15 August 2014 20:32, Fernando Perez wrote:
> Hi folks,
>
> [x-posting to numba and numpy lists, discussion should be on the
> python-ideas list where the python core folks actually reside]
>
> just a quick note that Guido has proposed to adopt the mypy model for type
> annotations in the langu
On 28 March 2014 19:56, Sturla Molden wrote:
> Matthew Brett wrote:
>
> > Does anyone know how their performance compares to MKL or the
> > reference implementations?
>
> http://eigen.tuxfamily.org/index.php?title=Benchmark
Very, very funny and twisted approach to legend-ordering-in-a-plot
ap
Hi,
On 22 March 2014 19:13, Nathaniel Smith wrote:
> Hi all,
>
> After 88 emails we don't have a conclusion in the other thread (see
> [1] for background). But we have to come to some conclusion or another
> if we want @ to exist :-). So I'll summarize where the discussion
> stands and let's see
On 6 February 2013 01:55, Chris Barker - NOAA Federal wrote:
> On Tue, Feb 5, 2013 at 4:32 PM, Matthew Brett wrote:
>>> 4) Numpy-MKL requires the Intel runtime DLLs (MKL is linked statically
>>> btw). I ship those with the installers and append the directory
>>> containing the DLLs to os.environ['
Hi All,
once again, my apologies for a (possibly) very ignorant question,
my google-fu is failing me... also because I am not sure of what
exactly I should look for.
My problem is relatively simple. Let's assume I have two Python
objects, A and B, and one of their attributes can assume a valu
On 2 July 2012 22:11, klo uo wrote:
> On Mon, Jul 2, 2012 at 9:26 PM, Fernando Perez wrote:
>>> ANNOUNCING
>>>
>>> eGenix PyRun - One file Python Runtime
>>>
>>> Version 1.0.0
>>>
>>>
>>> An easy-to-use single file relocatable Python run-time
On 26 June 2012 22:39, John Hunter wrote:
> On Tue, Jun 26, 2012 at 3:27 PM, Thouis (Ray) Jones wrote:
>> +1 !
>>
>> Speaking as someone trying to get started in contributing to numpy, I
>> find this discussion extremely off-putting. It's childish,
>> meaningless, and spiteful, and I think it's d
-- Forwarded message --
From: "Andrea Gavana"
Date: Feb 13, 2012 11:31 PM
Subject: Re: [Numpy-discussion] Creating parallel curves
To: "Jonathan Hilmer"
Thank you Jonathan for this, it's exactly what I was looking for. I' ll try
it tomorrow on th
ed or not, but
again every suggestion is most welcome.
Thank you.
Andrea.
>
>
> Jonathan
>
>
> On Sun, Feb 12, 2012 at 1:37 PM, Charles R Harris
> wrote:
>>
>>
>> On Sun, Feb 12, 2012 at 1:26 PM, Andrea Gavana
>> wrote:
>>>
>&
Charles,
On 12 February 2012 21:00, Charles R Harris wrote:
>
>
> On Fri, Feb 10, 2012 at 9:38 AM, Andrea Gavana
> wrote:
>>
>> Hi All,
>>
>> my apologies for my deep ignorance about math stuff; I guess I
>> should be able to find this
Jonathan,
On 12 February 2012 20:53, Jonathan Hilmer wrote:
> Andrea,
>
> Here is how to do it with splines. I would be more standard to return
> an array of normals, rather than two arrays of x and y components, but
> it actually requires less housekeeping this way. As an aside, I would
> prefe
HI Chris and All,
On 10 February 2012 17:53, Chris Barker wrote:
> Andrea,
>
>> Basically I have a set of x, y data (around 1,000 elements each) and I
>> want to create 2 parallel "curves" (offset curves) to the original
>> one; "parallel" means curves which are displaced from the base curve
>> by
Hi All,
my apologies for my deep ignorance about math stuff; I guess I
should be able to find this out but I keep getting impossible results.
Basically I have a set of x, y data (around 1,000 elements each) and I
want to create 2 parallel "curves" (offset curves) to the original
one; "paralle
Hi All,
I was fiddling here and there with some code doing dynamic import of
stuff, and I noticed that this code:
import os
import sys
init_name = r"C:\Python27\Lib\site-packages\numpy\__init__.py"
directory, module_name = os.path.split(init_name)
main = os.path.splitext(module_name)[0]
sy
Hi Dmitrey,
2011/8/15 Dmitrey :
> Hi all,
> I'm glad to inform you that general constraints handling for interalg (free
> solver with guaranteed user-defined precision) now is available. Despite it
> is very premature and requires lots of improvements, it is already capable
> of outperforming comm
Hi Chris and All,
On 12 August 2011 16:53, Christopher Jordan-Squire wrote:
> Hi Andrea--An easy way to get something like this would be
>
> import numpy as np
> import scipy.stats as stats
>
> sigma = #some reasonable standard deviation for your application
> x = stats.norm.rvs(size=1000, loc=125
Hi Chris & Brennan,
On 15 August 2011 00:59, Brennan Williams wrote:
> You can use scipy.stats.truncnorm, can't you? Unless I misread, you want to
> sample a normal distribution but with generated values only being within a
> specified range? However you also say you want to do this with triangula
Hi All,
I am working on something that appeared to be a no-brainer issue (at the
beginning), by my complete ignorance in statistics is overwhelming and I got
stuck.
What I am trying to do can be summarized as follows
Let's assume that I have to generate a sample of a 1,000 values for a
varia
(Resending as numpy-discussion has a 40 Kb message limit)
On 30 March 2010 22:44, Friedrich Romstedt wrote:
> 2010/3/30 Andrea Gavana :
>> However, from the first 100 or so interpolated simulations, I could
>> gather these findings:
>>
>> 1) Interpolations on *cumulat
Hi Friedrich & All,
On 30 March 2010 21:48, Friedrich Romstedt wrote:
> 2010/3/30 Andrea Gavana :
>> On 29 March 2010 23:44, Friedrich Romstedt wrote:
>>> When you have nice results using 40 Rbfs for each time instant, this
>>> procedure means that the values f
HI Friedrich & All,
On 29 March 2010 23:44, Friedrich Romstedt wrote:
> 2010/3/29 Andrea Gavana :
>> If anyone is interested in a follow up, I have tried a time-based
>> interpolation of my oil profile (and gas and gas injection profiles)
>> using those 40 interpolators
On 29 March 2010 23:13, Brennan Williams wrote:
> Andrea Gavana wrote:
>> Hi Chris and All,
>>
>> On 29 March 2010 22:35, Christopher Barker wrote:
>>
>>> Andrea Gavana wrote:
>>>
>>>>>> Scaling each axis by its standard deviation i
Hi Chris and All,
On 29 March 2010 22:35, Christopher Barker wrote:
> Andrea Gavana wrote:
>>>> Scaling each axis by its standard deviation is a typical first start.
>>>> Shifting and scaling the values such that they each go from 0 to 1 is
>>>> another
Hi Brennan & All,
On 29 March 2010 00:46, Brennan Williams wrote:
> Andrea Gavana wrote:
>> As for your question, the parameter are not spread completely
>> randomly, as this is a collection of simulations done over the years,
>> trying manually different scenarios,
Hi Kevin,
On 29 March 2010 01:38, Kevin Dunn wrote:
>> Message: 5
>> Date: Sun, 28 Mar 2010 00:24:01 +0000
>> From: Andrea Gavana
>> Subject: [Numpy-discussion] Interpolation question
>> To: Discussion of Numerical Python
>> Message-ID:
>>
>&
Hi All,
On 29 March 2010 00:59, Andrea Gavana wrote:
> On 29 March 2010 00:34, Robert Kern wrote:
>> Scaling each axis by its standard deviation is a typical first start.
>> Shifting and scaling the values such that they each go from 0 to 1 is
>> another useful thing to t
On 29 March 2010 00:34, Robert Kern wrote:
> On Sun, Mar 28, 2010 at 18:30, Andrea Gavana wrote:
>> Hi Friedrich & All,
>>
>> On 28 March 2010 23:51, Friedrich Romstedt wrote:
>>> 2010/3/28 Andrea Gavana :
>>>> Example 1
>>>>
>>
Hi Friedrich & All,
On 28 March 2010 23:51, Friedrich Romstedt wrote:
> 2010/3/28 Andrea Gavana :
>> Example 1
>>
>> # o2 and o3 are the number of production wells, split into 2
>> # different categories
>> # inj is the number of injection wells
>>
Hi Brennan & All,
On 28 March 2010 23:36, Brennan Williams wrote:
> Andrea Gavana wrote:
>> Let's see a couple of practical examples (I can share the data if
>> someone is interested).
>>
>>
> Definitely interested in helping solve this one so feel free t
HI Brennan,
On 28 March 2010 22:50, Brennan Williams wrote:
> Andrea Gavana wrote:
>> HI All,
>>
>> On 28 March 2010 19:22, Robert Kern wrote:
>>
>>> On Sun, Mar 28, 2010 at 03:26, Anne Archibald
>>> wrote:
>>>
>>>&g
Hi All,
On 28 March 2010 22:14, Pierre GM wrote:
> On Mar 28, 2010, at 4:47 PM, Andrea Gavana wrote:
>> HI All,
>>
>> On 28 March 2010 19:22, Robert Kern wrote:
>>> On Sun, Mar 28, 2010 at 03:26, Anne Archibald
>>> wrote:
>>>> On 27 M
Hi All,
On 28 March 2010 22:14, Pierre GM wrote:
> On Mar 28, 2010, at 4:47 PM, Andrea Gavana wrote:
>> HI All,
>>
>> On 28 March 2010 19:22, Robert Kern wrote:
>>> On Sun, Mar 28, 2010 at 03:26, Anne Archibald
>>> wrote:
>>>> On 27 M
HI All,
On 28 March 2010 19:22, Robert Kern wrote:
> On Sun, Mar 28, 2010 at 03:26, Anne Archibald
> wrote:
>> On 27 March 2010 20:24, Andrea Gavana wrote:
>>> Hi All,
>>>
>>> I have an interpolation problem and I am having some difficulties
>>
On 28 March 2010 08:26, Anne Archibald wrote:
> On 27 March 2010 20:24, Andrea Gavana wrote:
>> Hi All,
>>
>> I have an interpolation problem and I am having some difficulties
>> in tackling it. I hope I can explain myself clearly enough.
>>
>> Basically
Hi All,
I have an interpolation problem and I am having some difficulties
in tackling it. I hope I can explain myself clearly enough.
Basically, I have a whole bunch of 3D fluid flow simulations (close to
1000), and they are a result of different combinations of parameters.
I was planning to
Hi All,
I have tried the solutions proposed in the previous thread and it
looks like Chris' one is the fastest for my purposes. Now, I have a
question which is probably more conceptual than
implementation-related.
I started this little thread as my task is to read medium to
(relatively) big u
Hi All,
On Fri, May 22, 2009 at 11:28 PM, David Warde-Farley wrote:
> On 22-May-09, at 6:13 PM, Christopher Barker wrote:
>
>> that's why I put a sys.maxint at the end of the series...
>
> Oops! I foolishly assumed the sequence was unaltered. That makes a lot
> more sense.
Thank you guys for your
Hi All,
this should be a very easy question but I am trying to make a
script run as fast as possible, so please bear with me if the solution
is easy and I just overlooked it.
I have a list of integers, like this one:
indices = [1,2,3,4,5,6,7,8,9,255,256,257,258,10001,10002,10003,10004]
>Fro
Hi All,
I am tryin to "vectorize" 3 nested for loops but I am not having
much success. Here is the code I use:
import numpy
import numpy.ma as masked
grid = numpy.zeros((nx, ny), dtype=numpy.float32)
xOut = numpy.zeros((nx, ny), dtype=numpy.float32)
yOut = numpy.zeros((nx, ny), dtype=numpy.f
Hi All,
I am tryin to "vectorize" 3 nested for loops but I am not having
much success. Here is the code I use:
import numpy
import numpy.ma as masked
grid = numpy.zeros((nx, ny), dtype=numpy.float32)
xOut = numpy.zeros((nx, ny), dtype=numpy.float32)
yOut = numpy.zeros((nx, ny), dtype=numpy.f
Hi All,
On Mon, Oct 20, 2008 at 12:48 PM, Alan G Isaac wrote:
> On 10/20/2008 5:20 AM Andrea Gavana apparently wrote:
>> this is probably a very silly question, but combinatorial math is
>> not exactly my strength and I am not even sure on how to formulate the
>> questi
Hi All,
this is probably a very silly question, but combinatorial math is
not exactly my strength and I am not even sure on how to formulate the
question. I apologize if it is a very elementary problem.
Let's suppose that I have 60 oil wells and 3 surface facilities. Every
well must be tied-in
Hi Rob & All,
On Sat, Sep 27, 2008 at 4:05 PM, Rob Clewley wrote:
> Hi Andrea,
>
>>I was wondering if someone had any suggestions/references/snippets
>> of code on how to find the minimum distance between 2 paths in 3D.
>> Basically, for every path, I have I series of points (x, y, z) and I
>>
Hi All,
I was wondering if someone had any suggestions/references/snippets
of code on how to find the minimum distance between 2 paths in 3D.
Basically, for every path, I have I series of points (x, y, z) and I
would like to know if there is a better way, other than comparing
point by point th
Hi Stefan & All,
On Sun, May 25, 2008 at 8:59 PM, Stéfan van der Walt wrote:
> Hi Andrea
>
> 2008/5/25 Andrea Gavana <[EMAIL PROTECTED]>:
>>> When you bench the Cython code, you'll have to take out the Python
>>> calls (for checking dtype etc.), otherwi
Hi Stefan & All,
On Sat, May 24, 2008 at 8:11 PM, Stéfan van der Walt wrote:
> Hi Andrea
>
> 2008/5/24 Andrea Gavana <[EMAIL PROTECTED]>:
>> Number Of Cells: 5
>> ---
Hi All,
On Fri, May 23, 2008 at 4:50 PM, Andrea Gavana wrote:
> Hi Peter & All,
>
> On Fri, May 23, 2008 at 4:16 PM, Peter Creasey wrote:
>> Hi Andrea,
>>
>> 2008/5/23 "Andrea Gavana" <[EMAIL PROTECTED]>:
>>> And so on. The probelm with t
Hi Peter & All,
On Fri, May 23, 2008 at 4:16 PM, Peter Creasey wrote:
> Hi Andrea,
>
> 2008/5/23 "Andrea Gavana" <[EMAIL PROTECTED]>:
>> And so on. The probelm with this approach is that I lose the original
>> indices for which I want all the inequality t
Hi Stefan & All,
On Fri, May 23, 2008 at 1:02 AM, Stéfan van der Walt wrote:
> Hi Andrea
>
> 2008/5/23 Andrea Gavana <[EMAIL PROTECTED]>:
>> Thank you very much for this! I am going to try it and time it,
>> comparing it with the other implementations. I think I n
Hi Stefan,
On Thu, May 22, 2008 at 10:23 PM, Stéfan van der Walt wrote:
> Hi Andrea
>
> 2008/5/22 Andrea Gavana <[EMAIL PROTECTED]>:
>> By the way, about the solution Francesc posted:
>>
>> xyzReq = (xCent >= xMin) & (xCent <= xMax) & \
Hi Chris and All,
On Thu, May 22, 2008 at 8:40 PM, Christopher Barker wrote:
> Andrea Gavana wrote:
>> By the way, about the solution Francesc posted:
>>
>> xyzReq = (xCent >= xMin) & (xCent <= xMax) & \
>> (yCent >= yMin) & (y
Hi All,
On Thu, May 22, 2008 at 7:46 PM, Robert Kern wrote:
> On Thu, May 22, 2008 at 12:26 PM, Stéfan van der Walt <[EMAIL PROTECTED]>
> wrote:
>> Just to clarify things in my mind: is VTK *that* slow? I find that
>> surprising, since it is written in C or C++.
>
> Performance can depend more
Hi Francesc & All,
On Thu, May 22, 2008 at 1:04 PM, Francesc Alted wrote:
> I don't know if this is what you want, but you can get the boolean
> arrays separately, do the intersection and finally get the interesting
> values (by using fancy indexing) or coordinates (by using .nonzero()).
> Here it
Hi Stefan & All,
On Thu, May 22, 2008 at 12:29 PM, Stéfan van der Walt wrote:
> Hi Andrea
>
> 2008/5/22 Andrea Gavana <[EMAIL PROTECTED]>:
>>I am building some 3D grids for visualization starting from a much
>> bigger grid. I build these grids by satisfying c
Hi All,
I am building some 3D grids for visualization starting from a much
bigger grid. I build these grids by satisfying certain conditions on
x, y, z coordinates of their cells: up to now I was using VTK to
perform this operation, but VTK is slow as a turtle, so I thought to
use numpy to get
Hi All,
as I didn't get any answer, I thought I might repost with some
more analysis the trouble I am having with interp2d.
The attached script produces 3 different results depending on the
value of the parameter numColumns, i.e.:
1) numColumns = 20
Traceback (most recent call last):
File
Hi All,
On Mon, Apr 28, 2008 at 12:41 PM, Andrea Gavana wrote:
> Hi All,
>
> I have 2 matrices coming from 2 different simulations: the first
> column of the matrices is a date (time) at which all the other results
> in the matrix have been reported (simulation step). In the
Hi All,
I have 2 matrices coming from 2 different simulations: the first
column of the matrices is a date (time) at which all the other results
in the matrix have been reported (simulation step). In these 2
matrices, very often the simulation steps do not coincide, so I just
want to interpolate
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.
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
ocatable arrays as outputs
(whether from a subroutine or from a module).
> On 2/1/08, Andrea Gavana <[EMAIL PROTECTED]> wrote:
> > Hi All,
> >
> > I sent a couple of messages to f2py mailing list, but it seems
> > like my problem has no simple solution so I tho
Hi All,
I sent a couple of messages to f2py mailing list, but it seems
like my problem has no simple solution so I thought to ask for some
suggestions here.
Basically, I read some huge unformatted binary files which contain
time-step data from a reservoir simulation. I don't know the
dimensio
65 matches
Mail list logo