Dear FreeFEM Developer,
I try to solve the following eigenvalue problem using SLPEc parallelly.
The FE space for
[u, v, p] is
fespace Xh(Th, [P2, P2, P1]);
Note that the variables U and V are independent of spatial variables, so they can not be defined through
I solve the eigenvalue using
EPSSolve(A, B, vectors = vec, values = val, ....);
My question is that U and V are not defined in any finite element space, how can we compose
vectors = vec?
It seems impossible. I checked all the SLEPc examples provided by FreeFEM (link), but I didn’t find similar cases.
Can we store the eigenvector in an array, just like this
real[int, int] EigVec(Xh.ndof+2, num);
num denotes the number of eigenvalues.
Then we get eigenvectors through
Xh [eu, ev, ep]; eu = EigVec[0:Xh.ndof-1];
eu denotes the leading eigenvector.