implement and add this function to math_base.h
It is useful for various machine learning models, notably Latent Dirichlet Allocation.
lindahua
