complex
: to tell Gfem that complex number will be used.
When it is used it must be located at the beginning of the program
before any function declarations, otherwise the results will be incorrect.
It can appear more than once in the program but only the first occurrence
counts.
There is purposely no conjug function but barz=Re(z)-I*Im(z) will do.
By default all graphics display the real part. To display the imaginary part do plot(Im(f)).
The functions implemented for complex numbers are:
The linear systems for the PDE are solved by a Gauss complex LU factorization.
WARNING: failure to declare complex
in the program implies all
computation will be done in real, even if I
is used.