complex support

[bruno.zerbo at gmail.com]

Now I understand, thank you
Bruno


Il Monday 29 September 2008 12:23:48 Jed Brown ha scritto:
> On Mon 2008-09-29 11:56, bruno.zerbo at gmail.com wrote:
> > Hi,
> > I use PetscScalar to specify a complex and PetscReal for a double,
> > 1) is it wrong?
>
> Matrices and vectors have entries of type PetscScalar (which is complex
> if compiled --with-scalar-type=complex).  When solving a complex
> problem, there will be certain values (e.g. coordinates and parameters)
> which are always real.  This is the main purpose of PetscReal.  You
> can't make vectors of type PetscReal.
>
> > 2) is there a constant equal to the imaginary unit?
>
> PETSC_i
>
> Jed