coded FormJacobian, compiles and runs without crashing, added psivec to FormFunction to update psi from v (thought about getting rid of psi entirely in favor of v, but that is not a good idea because the psi complex vector is used many many times in the derivative calculations, and also psi must be a field of mode in multimode calculations with only one v vector)