BUG: specify output types for gmst functions

I noticed that the vectorized greenwich_mean_sidereal_time and greenwich_sidereal_time returned arrays with dtype=object because I forgot to mark the output type of the function. This doesn't impact the scalar version so existing behaviour should be unchanged.

I also found (via testing) that I had missed the long double type for the fused type.