0+ComplexScalar(5+8j, Function()) should be a complex data object

[AndreaCodd]

0+ComplexScalar(5+8j, Function()) should be a complex data object

[AndreaCodd]

when gEx1 is complex

        print("zero Area before",integrate(whereZero(length(g))))
        g +=  whereZero(length(g))*1e-10
        print("zero Area after",integrate(whereZero(length(g))))

       Does not give a zero area for "zero Area after",

It is necessary to write
print("zero Area before",integrate(whereZero(length(g))))
g += (1+1j)*whereZero(length(g))*1e-10
print("zero Area after",integrate(whereZero(length(g))))

It should be possible to add a complex Scalar function to a real Scalar function to get a complex Scalar function.

[LutzGross]

This is an in-place addition. One can argue that an update of a real data object with complex numbers is not allowed (not necessarily saying it should not be possible) but an exception should be thrown if this is not allowed.

[aellery]

I am having some trouble replicating this bug. With the latest commit, the following code prints "x is complex" to the console.

from esys.escript import *
from esys.finley import Rectangle
domain=Rectangle(l0=1, l1=1, n0=10, n1=10)
x=0+ComplexScalar(5+8j, Function(domain))
if x.isComplex() is True:
print("x is complex")

[esys-escript]

Closing the issue as nobody can reproduce the bug.