Work around old numpy version on server

np.full doesn’t exist there.
Avoid overflows in var-bound cased by log(gamma())