Fixed PF numerics, very fast now. But gives weird curves, need to check