BUG: Allos linalg.solve to handle zero equations.

This was previously correct, but the special case was missing
since the gufunc code cannot handle it.
Merge pull request #3854 from charris/backport-1.8-3848

Backport 1.8 3848