scipy/scipy — pull request #455 by WarrenWeckesser (2013-02-28 18:26:22+00:00)
Title: 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.