Hello, A while ago someone pointed me towards test_slepc.py for an example on how to solve an eigenvalue problem. I have been using this as a model to compute Quasi-Geostrophic basin modes in a closed basin. We have a version working in FENicS and I am very keen to get it working in Firedrake as well. I have a good start, which can be found on github, https://github.com/francispoulin/firedrakeQG/blob/master/basin_modes_qg.py Unfortunately, I don't know how to do a lot of things. For example specify number of eigenvalues to determine, specify a target eigenvalue, what method to use, or even just look at the eigenvector. Are there any other examples that someone might be able to share that might to able to do this? Cheers, Francis ------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637