Previous Up Next
Previous: 3.9.1 Iter
Up: 3.9 Other features
Next: 3.9.3 Scal()

3.9.2 Complex numbers

Gfem can handle complex coefficients with 4 dedicated keywords:

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.


Christophe Prud'homme

Previous Up Next