ENH: special: add xlogy = x*log(y) with xlogy=0 for x=0

This is a somewhat common special case in scipy.stats, and it's nicer to have a special function for it, rather than having to deal with special-casing separately in every instance.

See e.g. gh-467