fixed another huge typo in interpolate, now it works for 2d phc