I followed R, c, rx, cx = eigensolver.get_eigenpair(i) to be able to plot my eigenfunctions. However, in my case rx and cx are vector valued since my trial functions are also vector valued. In particular
(u1, u2) = TrialFunctions(Vh)
Do you know how can I extract the first and second component of rx and cx?
Thanks for the help.