Work around issue with truncated Float32 representation of lwork in

sgesdd by using nextfloat. See

http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4587&p=11036&hilit=sgesdd#p11036

and

scipy/scipy#5401

Fixes #15784

Author: andreasnoack

base/linalg/lapack.jl

@@ -1572,7 +1572,7 @@ for (geev, gesvd, gesdd, ggsvd, elty, relty) in
                 end
                 chklapackerror(info[])
                 if i == 1
-                    lwork = BlasInt(real(work[1]))
+                    lwork = round(BlasInt, nextfloat(real(work[1])))
                     work = Array($elty, lwork)
                 end
             end