Restricting trial function to converge the solution

variational formulation
6

Nowhere in the post does it state that you can do pk_space(0), pk_space(1).
These should be fe.FunctionAssigner(pk_space, [phiP_0.function_space(), ksi_0.function_space()])