On 10 December 2014 at 20:36, Chris Barker chris.bar...@noaa.gov wrote:
On Tue, Dec 9, 2014 at 11:03 PM, Sturla Molden sturla.mol...@gmail.com
wrote:
Nathaniel Smith n...@pobox.com wrote:
@contextmanager
def tmp_zeros(*args, **kwargs):
arr = np.zeros(*args, **kwargs)
try:
On 25 November 2014 at 19:33, David Cournapeau courn...@gmail.com wrote:
On Tue, Nov 25, 2014 at 6:10 PM, Sturla Molden sturla.mol...@gmail.com
wrote:
David Cournapeau courn...@gmail.com wrote:
Shall we consider a
href=https://github.com/scipy/scipy/issues/4168;
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)
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 language:
On 28 March 2014 19:56, Sturla Molden wrote:
Matthew Brett matthew.br...@gmail.com 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
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 if
On 6 February 2013 01:55, Chris Barker - NOAA Federal wrote:
On Tue, Feb 5, 2013 at 4:32 PM, Matthew Brett matthew.br...@gmail.com 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
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
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 -
available
On 26 June 2012 22:39, John Hunter wrote:
On Tue, Jun 26, 2012 at 3:27 PM, Thouis (Ray) Jones tho...@gmail.com 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
-- Forwarded message --
From: Andrea Gavana andrea.gav...@gmail.com
Date: Feb 13, 2012 11:31 PM
Subject: Re: [Numpy-discussion] Creating parallel curves
To: Jonathan Hilmer jkhil...@gmail.com
Thank you Jonathan for this, it's exactly what I was looking for. I' ll try
it tomorrow
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 a constant
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
prefer to
Charles,
On 12 February 2012 21:00, Charles R Harris wrote:
On Fri, Feb 10, 2012 at 9:38 AM, Andrea Gavana andrea.gav...@gmail.com
wrote:
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
Harris
charlesr.har...@gmail.com wrote:
On Sun, Feb 12, 2012 at 1:26 PM, Andrea Gavana andrea.gav...@gmail.com
wrote:
Charles,
On 12 February 2012 21:00, Charles R Harris wrote:
On Fri, Feb 10, 2012 at 9:38 AM, Andrea Gavana andrea.gav...@gmail.com
wrote:
Hi All,
my
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; parallel
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 = rC:\Python27\Lib\site-packages\numpy\__init__.py
directory, module_name = os.path.split(init_name)
main = os.path.splitext(module_name)[0]
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 triangular
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 Dmitrey,
2011/8/15 Dmitrey tm...@ukr.net:
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
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
(Resending as numpy-discussion has a 40 Kb message limit)
On 30 March 2010 22:44, Friedrich Romstedt friedrichromst...@gmail.com wrote:
2010/3/30 Andrea Gavana andrea.gav...@gmail.com:
However, from the first 100 or so interpolated simulations, I could
gather these findings:
1
Hi Friedrich All,
On 30 March 2010 21:48, Friedrich Romstedt wrote:
2010/3/30 Andrea Gavana andrea.gav...@gmail.com:
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 for one time instant
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 try.
Ah, magnifico! Thank you
Hi Kevin,
On 29 March 2010 01:38, Kevin Dunn wrote:
Message: 5
Date: Sun, 28 Mar 2010 00:24:01 +
From: Andrea Gavana andrea.gav...@gmail.com
Subject: [Numpy-discussion] Interpolation question
To: Discussion of Numerical Python numpy-discussion@scipy.org
Message-ID
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, without having in mind a proper
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 useful thing to try.
Ah, magnifico! Thank you Robert
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 is a typical first start.
Shifting and scaling the values such that they each go from 0
HI Friedrich All,
On 29 March 2010 23:44, Friedrich Romstedt wrote:
2010/3/29 Andrea Gavana andrea.gav...@gmail.com:
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 (and even
On 28 March 2010 08:26, Anne Archibald wrote:
On 27 March 2010 20:24, Andrea Gavana andrea.gav...@gmail.com 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, I have a whole bunch of 3D
HI All,
On 28 March 2010 19:22, Robert Kern wrote:
On Sun, Mar 28, 2010 at 03:26, Anne Archibald peridot.face...@gmail.com
wrote:
On 27 March 2010 20:24, Andrea Gavana andrea.gav...@gmail.com wrote:
Hi All,
I have an interpolation problem and I am having some difficulties
in tackling
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 peridot.face...@gmail.com
wrote:
On 27 March 2010 20:24, Andrea Gavana andrea.gav
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 peridot.face...@gmail.com
wrote:
On 27 March 2010 20:24, Andrea Gavana andrea.gav
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 peridot.face...@gmail.com
wrote:
On 27 March 2010 20:24, Andrea Gavana andrea.gav...@gmail.com wrote:
Hi
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 to email
the data (obviously not 1,000 Eclipse smspec
Hi Friedrich All,
On 28 March 2010 23:51, Friedrich Romstedt wrote:
2010/3/28 Andrea Gavana andrea.gav...@gmail.com:
Example 1
# o2 and o3 are the number of production wells, split into 2
# different categories
# inj is the number of injection wells
# fomts is the final oil recovery
rbf
On 29 March 2010 00:34, Robert Kern wrote:
On Sun, Mar 28, 2010 at 18:30, Andrea Gavana andrea.gav...@gmail.com wrote:
Hi Friedrich All,
On 28 March 2010 23:51, Friedrich Romstedt wrote:
2010/3/28 Andrea Gavana andrea.gav...@gmail.com:
Example 1
# o2 and o3 are the number of production
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,
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 help,
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]
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),
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),
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
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
question. I apologize
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
would
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
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
-
| Rank | Method Name | Execution Time
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.), otherwise you're comparing apples and
oranges. After I
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 this approach is that I lose the original
indices for which I want all
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 need to study a
bit your code as I know
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 tests to succeed:
To have the original indices you
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
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 certain conditions on
x, y, z coordinates
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 is
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 on the
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) (yCent = yMax) \
(zCent = zMin) (zCent = zMax)
xyzReq
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) \
(yCent = yMin) (yCent = yMax) \
(zCent = zMin
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):
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
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
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
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 thought to ask for some
suggestions here.
Basically, I read some huge unformatted
62 matches
Mail list logo
