[Numpy-discussion] question about in-place operations

Massimo DiPierro massimo.dipierro at gmail.com
Tue May 22 11:06:31 EDT 2012


Thank you this does it.



On May 22, 2012, at 9:59 AM, Robert Kern wrote:

> On Tue, May 22, 2012 at 3:47 PM, Massimo DiPierro
> <massimo.dipierro at gmail.com> wrote:
>> Thank you. I will look into numexpr.
>> 
>> Anyway, I do not need arbitrary expressions. If there were something like
>> 
>> numpy.add_scaled(a,scale,b)
>> 
>> with support for scale in int, float, complex, this would be sufficient for me.
> 
> BLAS has the xAXPY functions, which will do this for float and complex.
> 
> import numpy as np
> from scipy.linalg import fblas
> 
> def add_scaled_inplace(a, scale, b):
>    if np.issubdtype(a.dtype, complex):
>        fblas.zaxpy(b, a, a=scale)
>    else:
>        fblas.daxpy(b, a, a=scale)
> 
> -- 
> Robert Kern
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion




More information about the NumPy-Discussion mailing list