Inplace operators on complex numbers

[dagss]

This is not allowed:

cdef complex z = 1
z += 1j

Trivially done in a transform.

[dagss]

If #441 doesn't make it, this shouldn't :-)

[robertwb]

We need to be sure that

x[f(i)] += 1

has the right semantics, which are different than

x[f(i)] = x[f(i)] + 1

(Same for attribute access.) }}}

[robertwb]

 http://hg.cython.org/cython-devel/rev/c8e87c5b2b58

[dagss]

Mike T on ML complained that it still doesn't work with buffers:

no problem-- it occurred to me that the compiler flags might help. They're pasted below, too

import numpy as np
cimport numpy as cnp
def copy_vec(
    cnp.ndarray[cnp.complex128_t] z_arr
    ):

    cdef cnp.ndarray[cnp.complex128_t] b = np.zeros((len(z_arr),),np.complex128)
    for i in xrange(len(z_arr)):
        b[i] += z_arr[i]


C compiler: /usr/bin/gcc-4.2 -fno-strict-aliasing -fno-common -dynamic -pipe -O2 -DNDEBUG -g -fwrapv -O3 -Wall -Wstrict-prototypes

The error raised in the C compiler is (another source file)

pmri_opt/sensemap.c: In function ‘__pyx_pf_8pmri_opt_8sensemap_2poly_vecc’:
pmri_opt/sensemap.c:2781: error: invalid operands to binary + (have ‘__pyx_t_double_complex’ and ‘__pyx_t_double_complex’)

[mark]

Replying to dagss:

Mike T on ML complained that it still doesn't work with buffers: {{{ no problem-- it occurred to me that the compiler flags might help. They're pasted below, too import numpy as np cimport numpy as cnp def copy_vec( cnp.ndarray[cnp.complex128_t] z_arr ): cdef cnp.ndarray[cnp.complex128_t] b = np.zeros((len(z_arr),),np.complex128) for i in xrange(len(z_arr)): b[i] += z_arr[i] C compiler: /usr/bin/gcc-4.2 -fno-strict-aliasing -fno-common -dynamic -pipe -O2 -DNDEBUG -g -fwrapv -O3 -Wall -Wstrict-prototypes }}} The error raised in the C compiler is (another source file) {{{ pmri_opt/sensemap.c: In function ‘pyx_pf_8pmri_opt_8sensemap_2poly_vecc’: pmri_opt/sensemap.c:2781: error: invalid operands to binary + (have ‘pyx_t_double_complex’ and ‘pyx_t_double_complex’) }}}

Inplace buffer operators were fixed in one of the commits for buffers and memoryviews, this should work now.