diff --git a/SRC/cgelsy.f b/SRC/cgelsy.f
index 557742a9c..e8dfe3bbc 100644
--- a/SRC/cgelsy.f
+++ b/SRC/cgelsy.f
@@ -266,9 +266,9 @@
 *
       INFO = 0
       NB1 = ILAENV( 1, 'CGEQRF', ' ', M, N, -1, -1 )
-      NB2 = ILAENV( 1, 'CGERQF', ' ', M, N, -1, -1 )
+      NB2 = ILAENV( 1, 'CTZRZF', ' ', M, N, -1, -1 )
       NB3 = ILAENV( 1, 'CUNMQR', ' ', M, N, NRHS, -1 )
-      NB4 = ILAENV( 1, 'CUNMRQ', ' ', M, N, NRHS, -1 )
+      NB4 = ILAENV( 1, 'CUNMRZ', ' ', M, N, NRHS, -1 )
       NB = MAX( NB1, NB2, NB3, NB4 )
       LWKOPT = MAX( 1, MN+2*N+NB*(N+1), 2*MN+NB*NRHS )
       WORK( 1 ) = CMPLX( LWKOPT )
diff --git a/SRC/dgelsy.f b/SRC/dgelsy.f
index 54bc881c2..7f51af280 100644
--- a/SRC/dgelsy.f
+++ b/SRC/dgelsy.f
@@ -275,9 +275,9 @@
             LWKOPT = 1
          ELSE
             NB1 = ILAENV( 1, 'DGEQRF', ' ', M, N, -1, -1 )
-            NB2 = ILAENV( 1, 'DGERQF', ' ', M, N, -1, -1 )
+            NB2 = ILAENV( 1, 'DTZRZF', ' ', M, N, -1, -1 )
             NB3 = ILAENV( 1, 'DORMQR', ' ', M, N, NRHS, -1 )
-            NB4 = ILAENV( 1, 'DORMRQ', ' ', M, N, NRHS, -1 )
+            NB4 = ILAENV( 1, 'DORMRZ', ' ', M, N, NRHS, -1 )
             NB = MAX( NB1, NB2, NB3, NB4 )
             LWKMIN = MN + MAX( 2*MN, N + 1, MN + NRHS )
             LWKOPT = MAX( LWKMIN,
diff --git a/SRC/sgelsy.f b/SRC/sgelsy.f
index bf8368414..22f65eac1 100644
--- a/SRC/sgelsy.f
+++ b/SRC/sgelsy.f
@@ -276,9 +276,9 @@
             LWKOPT = 1
          ELSE
             NB1 = ILAENV( 1, 'SGEQRF', ' ', M, N, -1, -1 )
-            NB2 = ILAENV( 1, 'SGERQF', ' ', M, N, -1, -1 )
+            NB2 = ILAENV( 1, 'STZRZF', ' ', M, N, -1, -1 )
             NB3 = ILAENV( 1, 'SORMQR', ' ', M, N, NRHS, -1 )
-            NB4 = ILAENV( 1, 'SORMRQ', ' ', M, N, NRHS, -1 )
+            NB4 = ILAENV( 1, 'SORMRZ', ' ', M, N, NRHS, -1 )
             NB = MAX( NB1, NB2, NB3, NB4 )
             LWKMIN = MN + MAX( 2*MN, N + 1, MN + NRHS )
             LWKOPT = MAX( LWKMIN,
diff --git a/SRC/zgelsy.f b/SRC/zgelsy.f
index 6042cd55b..4a0cdb935 100644
--- a/SRC/zgelsy.f
+++ b/SRC/zgelsy.f
@@ -266,9 +266,9 @@
 *
       INFO = 0
       NB1 = ILAENV( 1, 'ZGEQRF', ' ', M, N, -1, -1 )
-      NB2 = ILAENV( 1, 'ZGERQF', ' ', M, N, -1, -1 )
+      NB2 = ILAENV( 1, 'ZTZRZF', ' ', M, N, -1, -1 )
       NB3 = ILAENV( 1, 'ZUNMQR', ' ', M, N, NRHS, -1 )
-      NB4 = ILAENV( 1, 'ZUNMRQ', ' ', M, N, NRHS, -1 )
+      NB4 = ILAENV( 1, 'ZUNMRZ', ' ', M, N, NRHS, -1 )
       NB = MAX( NB1, NB2, NB3, NB4 )
       LWKOPT = MAX( 1, MN+2*N+NB*( N+1 ), 2*MN+NB*NRHS )
       WORK( 1 ) = DCMPLX( LWKOPT )