BUG: special: multigammaln returned wrong size array (closes ticket #1849)

The only change to the function is to remove a few lines that set axis=-1 if the size of the input array is 1.  I'm guessing there was once a good reason for this, but later changes made this incorrect.