Ppm

Preconditioner for F, size(P)=[Nj x Nk]. Fu = uhat, (PF)u = Puhat, so ubar==Puhat