Hi Francis, On 06/12/16 04:24, Francis Poulin wrote:
Hello,
Sorry for the bother but I am wondering if anyone had any advice on my problem. The code can be found at the following:
https://github.com/francispoulin/firedrakeQG/blob/master/basin_modes_sw.py
There are two problems:
1) I can't seem to plot the eigenfunctions. Maybe I didn't split them correctly? I'm not sure.
From your submitted demo, it seems that you can now plot these, is that right?
2) The eigenvalues are not at all the ones I'm looking for. I know from the fenics version that there is an eigenvalue with an imaginary part of 3.14. These might be valid but I suspect that they are very noisy and not physically interesting.
I'm not sure about that. If I run a small problem and compute all the eigenvalues using scipy.linalg.eig, the maximum value of the computed real and imaginary parts are 0.108... which agrees, I think with your firedrake and fenics code. Thanks, Lawrence