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 >
