not recalculating true anomaly for post-compiled phase curves

[ben-cassese]

Although polynomial_limb_darkened_transit calculates f internally, the phase curve functions all take 2D/3D coefficients as inputs, not times. Need to be sure to re-compute f in the actual OblateSystem.lightcurve() method before computing those coefficients.