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 > <[email protected]> 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 > [email protected] > http://mail.scipy.org/mailman/listinfo/numpy-discussion _______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
