Bug fix for the lib.frame_order.pseudo_ellipse.tmax_pseudo_ellipse_array() function.
The problem was that when theta_x or theta_y were zero, the floating point value of 0.0 would be
returned. This is the incorrect behaviour as the returned value must be an array matching the
dimensions of the phi angle array argument.
