Complex support in dolfin

1

Hello!

I’m trying to solve a generalized eigenvalue problem by specifying complex values of shift for the shift-invert algorithm. My matrices A and B are real.

Does any version of dolfin (not dolfinx) currently support complex shifts to the SLEPcEigensolver object? I see that complex support is possible according to https://fenicsproject.org/fenics-project-roadmap-2019/
, is this a reference to dolfinx?

Given that A and B are real themselves, if dolfin can be built from source against complex PETSc this should solve the problem right?

Thanks!

2

There is a bit of discussion on this question at Solving complex-valued PDE

1 Like
3

Thanks, I did see that but that talks more about splitting the problem into real and imaginary parts. My question is more about applying a complex shift to real matrices A and B in the eigenvalue problem.

The SLEPcEigensolver object is just a wrapper around the PETSc EPS object which can support complex shifts if built against complex PETSc right?

Is it at all possible to build dolfin against complex PETSc? I dont want to solve the complex GEVP just apply a complex shift.

4

As discussed in that thread, you’ll need dolfinx if you want complex number support in the software.

1 Like