Faster log_fac calculation
Tests for log_probability in multinomial