Wrap xgelsy and change default least square solver from xgelsd to xgelsy.
Merge pull request #3298 from JuliaLang/anj/gelsy

Wrap xgelsy and change default least square solver from xgelsd to xgelsy