On 20-Feb-09, at 6:41 AM, Olivier Grisel wrote:
> Alright, thanks for the reply.
>
> Is there a canonical way /sample code to gain low level access to
> blas / lapack
> atlas routines using ctypes from numpy / scipy code?
>
> I don't mind fixing the dimensions and the ndtype of my array if it
On 20-Feb-09, at 10:39 AM, Robert Kern wrote:
>> Fiddling around I've found other blas/lapack
>> functions where the same arg is offered but the choice actually does
>> something.
>
> Examples?
An even better example is scipy.linalg.fblas.dgemv, the matrix-vector
equivalent of dgemm. overwrite_
On 20-Feb-09, at 10:39 AM, Robert Kern wrote:
>> Fiddling around I've found other blas/lapack
>> functions where the same arg is offered but the choice actually does
>> something.
>
> Examples?
scipy.lib.lapack.flapack.dpotri, for example. I'm not sure of the
proper usage, but when I pass it an
:00 AM, numpy-discussion-requ...@scipy.org wrote:
Message: 2
Date: Fri, 20 Feb 2009 09:52:03 -0600
From: Robert Kern
Subject: Re: [Numpy-discussion] inplace dot products
To: Discussion of Numerical Python
Message-ID:
<3d375d730902200752k6b71efc9gb1e0ac1d260d...@mail.gmail.com>
Conte
On Fri, Feb 20, 2009 at 05:18, David Warde-Farley wrote:
> Hi Olivier,
>
> There was this idea posted on the Scipy-user list a while back:
>
>http://projects.scipy.org/pipermail/scipy-user/2008-August/017954.html
>
> but it doesn't look like he got anywhere with it, or even got a
> respons
On Fri, Feb 20, 2009 at 06:25, David Warde-Farley wrote:
>
> On 20-Feb-09, at 6:39 AM, David Cournapeau wrote:
>
>> You can just use ctypes to access ATLAS, as you would do for any
>> library. Or do you mean something else ?
>
> Say, David... :)
>
> Do you have any idea why the pyf wrapper for fbl
On 20-Feb-09, at 6:39 AM, David Cournapeau wrote:
> You can just use ctypes to access ATLAS, as you would do for any
> library. Or do you mean something else ?
Say, David... :)
Do you have any idea why the pyf wrapper for fblas3 completely ignores
the overwrite_c argument? Fiddling around I'v
Olivier Grisel wrote:
> 2009/2/20 David Warde-Farley :
>
>> Hi Olivier,
>>
>> There was this idea posted on the Scipy-user list a while back:
>>
>>http://projects.scipy.org/pipermail/scipy-user/2008-August/017954.html
>>
>> but it doesn't look like he got anywhere with it, or even got a
2009/2/20 David Warde-Farley :
> Hi Olivier,
>
> There was this idea posted on the Scipy-user list a while back:
>
> http://projects.scipy.org/pipermail/scipy-user/2008-August/017954.html
>
> but it doesn't look like he got anywhere with it, or even got a
> response.
>
> I just tried it and
Hi Olivier,
There was this idea posted on the Scipy-user list a while back:
http://projects.scipy.org/pipermail/scipy-user/2008-August/017954.html
but it doesn't look like he got anywhere with it, or even got a
response.
I just tried it and I observe the same behaviour. A quick look a
Hi numpist people,
I discovered the ufunc and there ability to compute the results on
preallocated arrays:
>>> a = arange(10, dtype=float32)
>>> b = arange(10, dtype=float32) + 1
>>> c = add(a, b, a)
>>> c is a
True
>>> a
array([ 1., 3., 5., 7., 9., 11., 13., 15., 17.,
19.],
11 matches
Mail list logo
