How use scalar type and complex type at the same time?

[Randall Mackie]

I have a code where I'm solving Maxwell's equations in 3D, so this is complex,
but then later I'm also solving a system of equations on the model space, and
this is real. I just simply use complex for both - it's just easier to use complex
for solving the real systems and not worry about the additional overhead, at least
to me.

Randy


Zhifeng Sheng wrote:
> Dear all
> 
> I am working on a FEM program, which can be used to solve time domain 
> and frequency domain EM problems. I finished the time domain solver part 
> and now I need to implement the frequency domain solver.... So I built 
> Petsc with --scalar-type = complex, and then my time domain solver can 
> not be compiled, it says " can not convert double to PetscScalar".
> 
> Is there any way that I can handle both real matrix and complex matrix 
> at the same time? (without changing my old code too much?)
> 
> Thanks
> Best regards
> Zhifeng Sheng
>