diff --git a/lapack-netlib/TESTING/MATGEN/clagge.c b/lapack-netlib/TESTING/MATGEN/clagge.c
index 838df0c02..f05905bd7 100644
--- a/lapack-netlib/TESTING/MATGEN/clagge.c
+++ b/lapack-netlib/TESTING/MATGEN/clagge.c
@@ -627,7 +627,7 @@ static integer c__1 = 1;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clagge_(integer *m, integer *n, integer *kl, integer *ku,
+/* Subroutine */ void clagge_(integer *m, integer *n, integer *kl, integer *ku,
 	 real *d__, complex *a, integer *lda, integer *iseed, complex *work, 
 	integer *info)
 {
@@ -638,16 +638,17 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer i__, j;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *),
 	     cscal_(integer *, complex *, complex *, integer *), cgemv_(char *
 	    , integer *, integer *, complex *, complex *, integer *, complex *
 	    , integer *, complex *, complex *, integer *);
     extern real scnrm2_(integer *, complex *, integer *);
     complex wa, wb;
-    extern /* Subroutine */ int clacgv_(integer *, complex *, integer *);
+    extern /* Subroutine */ void clacgv_(integer *, complex *, integer *);
     real wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), clarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void clarnv_(
 	    integer *, integer *, integer *, complex *);
     complex tau;
 
@@ -686,8 +687,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("CLAGGE", &i__1);
-	return 0;
+	xerbla_("CLAGGE", &i__1, 6);
+	return;
     }
 
 /*     initialize A to diagonal matrix */
@@ -713,7 +714,7 @@ static integer c__1 = 1;
 /*     Quick exit if the user wants a diagonal matrix */
 
     if (*kl == 0 && *ku == 0) {
-	return 0;
+	return;
     }
 
 /*     pre- and post-multiply A by random unitary matrices */
@@ -1024,7 +1025,7 @@ static integer c__1 = 1;
 	}
 /* L70: */
     }
-    return 0;
+    return;
 
 /*     End of CLAGGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/claghe.c b/lapack-netlib/TESTING/MATGEN/claghe.c
index 084a3264e..77ed94912 100644
--- a/lapack-netlib/TESTING/MATGEN/claghe.c
+++ b/lapack-netlib/TESTING/MATGEN/claghe.c
@@ -615,7 +615,7 @@ static integer c__1 = 1;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int claghe_(integer *n, integer *k, real *d__, complex *a, 
+/* Subroutine */ void claghe_(integer *n, integer *k, real *d__, complex *a, 
 	integer *lda, integer *iseed, complex *work, integer *info)
 {
     /* System generated locals */
@@ -624,17 +624,17 @@ static integer c__1 = 1;
     complex q__1, q__2, q__3, q__4;
 
     /* Local variables */
-    extern /* Subroutine */ int cher2_(char *, integer *, complex *, complex *
+    extern /* Subroutine */ void cher2_(char *, integer *, complex *, complex *
 	    , integer *, complex *, integer *, complex *, integer *);
     integer i__, j;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *);
     complex alpha;
-    extern /* Subroutine */ int cscal_(integer *, complex *, complex *, 
+    extern /* Subroutine */ void cscal_(integer *, complex *, complex *, 
 	    integer *);
     extern /* Complex */ VOID cdotc_(complex *, integer *, complex *, integer 
 	    *, complex *, integer *);
-    extern /* Subroutine */ int cgemv_(char *, integer *, integer *, complex *
+    extern /* Subroutine */ void cgemv_(char *, integer *, integer *, complex *
 	    , complex *, integer *, complex *, integer *, complex *, complex *
 	    , integer *), chemv_(char *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, complex *, 
@@ -643,7 +643,8 @@ static integer c__1 = 1;
     extern real scnrm2_(integer *, complex *, integer *);
     complex wa, wb;
     real wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), clarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void clarnv_(
 	    integer *, integer *, integer *, complex *);
     complex tau;
 
@@ -678,8 +679,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("CLAGHE", &i__1);
-	return 0;
+	xerbla_("CLAGHE", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -857,7 +858,7 @@ static integer c__1 = 1;
 	}
 /* L80: */
     }
-    return 0;
+    return;
 
 /*     End of CLAGHE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clagsy.c b/lapack-netlib/TESTING/MATGEN/clagsy.c
index 908bdc4a8..de215bc63 100644
--- a/lapack-netlib/TESTING/MATGEN/clagsy.c
+++ b/lapack-netlib/TESTING/MATGEN/clagsy.c
@@ -615,7 +615,7 @@ static integer c__1 = 1;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clagsy_(integer *n, integer *k, real *d__, complex *a, 
+/* Subroutine */ void clagsy_(integer *n, integer *k, real *d__, complex *a, 
 	integer *lda, integer *iseed, complex *work, integer *info)
 {
     /* System generated locals */
@@ -626,14 +626,14 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer i__, j;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *);
     complex alpha;
-    extern /* Subroutine */ int cscal_(integer *, complex *, complex *, 
+    extern /* Subroutine */ void cscal_(integer *, complex *, complex *, 
 	    integer *);
     extern /* Complex */ VOID cdotc_(complex *, integer *, complex *, integer 
 	    *, complex *, integer *);
-    extern /* Subroutine */ int cgemv_(char *, integer *, integer *, complex *
+    extern /* Subroutine */ void cgemv_(char *, integer *, integer *, complex *
 	    , complex *, integer *, complex *, integer *, complex *, complex *
 	    , integer *), caxpy_(integer *, complex *, complex *, 
 	    integer *, complex *, integer *), csymv_(char *, integer *, 
@@ -642,9 +642,10 @@ static integer c__1 = 1;
     extern real scnrm2_(integer *, complex *, integer *);
     integer ii, jj;
     complex wa, wb;
-    extern /* Subroutine */ int clacgv_(integer *, complex *, integer *);
+    extern /* Subroutine */ void clacgv_(integer *, complex *, integer *);
     real wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), clarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void clarnv_(
 	    integer *, integer *, integer *, complex *);
     complex tau;
 
@@ -679,8 +680,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("CLAGSY", &i__1);
-	return 0;
+	xerbla_("CLAGSY", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -910,7 +911,7 @@ static integer c__1 = 1;
 	}
 /* L120: */
     }
-    return 0;
+    return;
 
 /*     End of CLAGSY */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clahilb.c b/lapack-netlib/TESTING/MATGEN/clahilb.c
index 4c4fe1bb9..95b00fb3e 100644
--- a/lapack-netlib/TESTING/MATGEN/clahilb.c
+++ b/lapack-netlib/TESTING/MATGEN/clahilb.c
@@ -645,7 +645,7 @@ static complex c_b6 = {0.f,0.f};
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clahilb_(integer *n, integer *nrhs, complex *a, integer *
+/* Subroutine */ void clahilb_(integer *n, integer *nrhs, complex *a, integer *
 	lda, complex *x, integer *ldx, complex *b, integer *ldb, real *work, 
 	integer *info, char *path)
 {
@@ -670,9 +670,9 @@ static complex c_b6 = {0.f,0.f};
     integer i__, j, m, r__;
     char c2[2];
     integer ti, tm;
-    extern /* Subroutine */ int claset_(char *, integer *, integer *, complex 
-	    *, complex *, complex *, integer *), xerbla_(char *, 
-	    integer *);
+    extern /* Subroutine */ void claset_(char *, integer *, integer *, complex 
+	    *, complex *, complex *, integer *);
+    extern int xerbla_(char *, integer *, ftnlen);
     extern logical lsamen_(integer *, char *, char *);
     complex tmp;
 
@@ -722,8 +722,8 @@ static complex c_b6 = {0.f,0.f};
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("CLAHILB", &i__1);
-	return 0;
+	xerbla_("CLAHILB", &i__1, 7);
+	return;
     }
     if (*n > 6) {
 	*info = 1;
@@ -830,6 +830,6 @@ static complex c_b6 = {0.f,0.f};
 	    }
 	}
     }
-    return 0;
+    return;
 } /* clahilb_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clakf2.c b/lapack-netlib/TESTING/MATGEN/clakf2.c
index e24f47f2d..e15bc58ab 100644
--- a/lapack-netlib/TESTING/MATGEN/clakf2.c
+++ b/lapack-netlib/TESTING/MATGEN/clakf2.c
@@ -615,7 +615,7 @@ static complex c_b1 = {0.f,0.f};
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clakf2_(integer *m, integer *n, complex *a, integer *lda,
+/* Subroutine */ void clakf2_(integer *m, integer *n, complex *a, integer *lda,
 	 complex *b, complex *d__, complex *e, complex *z__, integer *ldz)
 {
     /* System generated locals */
@@ -625,7 +625,7 @@ static complex c_b1 = {0.f,0.f};
 
     /* Local variables */
     integer i__, j, l, ik, jk, mn;
-    extern /* Subroutine */ int claset_(char *, integer *, integer *, complex 
+    extern /* Subroutine */ void claset_(char *, integer *, integer *, complex 
 	    *, complex *, complex *, integer *);
     integer mn2;
 
@@ -737,7 +737,7 @@ static complex c_b1 = {0.f,0.f};
 /* L90: */
     }
 
-    return 0;
+    return;
 
 /*     End of CLAKF2 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clarge.c b/lapack-netlib/TESTING/MATGEN/clarge.c
index 299ec9308..0215436b7 100644
--- a/lapack-netlib/TESTING/MATGEN/clarge.c
+++ b/lapack-netlib/TESTING/MATGEN/clarge.c
@@ -600,7 +600,7 @@ static integer c__1 = 1;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clarge_(integer *n, complex *a, integer *lda, integer *
+/* Subroutine */ void clarge_(integer *n, complex *a, integer *lda, integer *
 	iseed, complex *work, integer *info)
 {
     /* System generated locals */
@@ -610,7 +610,7 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer i__;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *),
 	     cscal_(integer *, complex *, complex *, integer *), cgemv_(char *
 	    , integer *, integer *, complex *, complex *, integer *, complex *
@@ -618,7 +618,8 @@ static integer c__1 = 1;
     extern real scnrm2_(integer *, complex *, integer *);
     complex wa, wb;
     real wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), clarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void clarnv_(
 	    integer *, integer *, integer *, complex *);
     complex tau;
 
@@ -650,8 +651,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("CLARGE", &i__1);
-	return 0;
+	xerbla_("CLARGE", &i__1, 6);
+	return;
     }
 
 /*     pre- and post-multiply A by random unitary matrix */
@@ -702,7 +703,7 @@ static integer c__1 = 1;
 		* a_dim1 + 1], lda);
 /* L10: */
     }
-    return 0;
+    return;
 
 /*     End of CLARGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/claror.c b/lapack-netlib/TESTING/MATGEN/claror.c
index 1fb287af3..cd0d15300 100644
--- a/lapack-netlib/TESTING/MATGEN/claror.c
+++ b/lapack-netlib/TESTING/MATGEN/claror.c
@@ -672,7 +672,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int claror_(char *side, char *init, integer *m, integer *n, 
+/* Subroutine */ void claror_(char *side, char *init, integer *m, integer *n, 
 	complex *a, integer *lda, integer *iseed, complex *x, integer *info)
 {
     /* System generated locals */
@@ -683,23 +683,23 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     integer kbeg, jcol;
     real xabs;
     integer irow, j;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *),
 	     cscal_(integer *, complex *, complex *, integer *);
     extern logical lsame_(char *, char *);
-    extern /* Subroutine */ int cgemv_(char *, integer *, integer *, complex *
+    extern /* Subroutine */ void cgemv_(char *, integer *, integer *, complex *
 	    , complex *, integer *, complex *, integer *, complex *, complex *
 	    , integer *);
     complex csign;
     integer ixfrm, itype, nxfrm;
     real xnorm;
     extern real scnrm2_(integer *, complex *, integer *);
-    extern /* Subroutine */ int clacgv_(integer *, complex *, integer *);
+    extern /* Subroutine */ void clacgv_(integer *, complex *, integer *);
     //extern /* Complex */ VOID clarnd_(complex *, integer *, integer *);
     extern complex clarnd_(integer *, integer *);
-    extern /* Subroutine */ int claset_(char *, integer *, integer *, complex 
-	    *, complex *, complex *, integer *), xerbla_(char *, 
-	    integer *);
+    extern /* Subroutine */ void claset_(char *, integer *, integer *, complex 
+	    *, complex *, complex *, integer *);
+    extern int xerbla_(char *, integer *, ftnlen);
     real factor;
     complex xnorms;
 
@@ -723,7 +723,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     /* Function Body */
     *info = 0;
     if (*n == 0 || *m == 0) {
-	return 0;
+	return;
     }
 
     itype = 0;
@@ -750,8 +750,8 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     }
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLAROR", &i__1);
-	return 0;
+	xerbla_("CLAROR", &i__1, 6);
+	return;
     }
 
     if (itype == 1) {
@@ -815,8 +815,8 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 	if (abs(factor) < 1e-20f) {
 	    *info = 1;
 	    i__2 = -(*info);
-	    xerbla_("CLAROR", &i__2);
-	    return 0;
+	    xerbla_("CLAROR", &i__2, 6);
+	    return;
 	} else {
 	    factor = 1.f / factor;
 	}
@@ -899,7 +899,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 /* L90: */
 	}
     }
-    return 0;
+    return;
 
 /*     End of CLAROR */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clarot.c b/lapack-netlib/TESTING/MATGEN/clarot.c
index 9bc295062..f1e73ed5b 100644
--- a/lapack-netlib/TESTING/MATGEN/clarot.c
+++ b/lapack-netlib/TESTING/MATGEN/clarot.c
@@ -739,7 +739,7 @@ static integer c__8 = 8;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clarot_(logical *lrows, logical *lleft, logical *lright, 
+/* Subroutine */ void clarot_(logical *lrows, logical *lleft, logical *lright, 
 	integer *nl, complex *c__, complex *s, complex *a, integer *lda, 
 	complex *xleft, complex *xright)
 {
@@ -752,7 +752,7 @@ static integer c__8 = 8;
     complex tempx;
     integer ix, iy, nt;
     complex xt[2], yt[2];
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     integer iyt;
 
 
@@ -804,12 +804,12 @@ static integer c__8 = 8;
 /*     Check for errors */
 
     if (*nl < nt) {
-	xerbla_("CLAROT", &c__4);
-	return 0;
+	xerbla_("CLAROT", &c__4, 6);
+	return;
     }
     if (*lda <= 0 || ! (*lrows) && *lda < *nl - nt) {
-	xerbla_("CLAROT", &c__8);
-	return 0;
+	xerbla_("CLAROT", &c__8, 6);
+	return;
     }
 
 /*     Rotate */
@@ -887,7 +887,7 @@ static integer c__8 = 8;
 	a[i__1].r = yt[i__2].r, a[i__1].i = yt[i__2].i;
     }
 
-    return 0;
+    return;
 
 /*     End of CLAROT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatm1.c b/lapack-netlib/TESTING/MATGEN/clatm1.c
index a2a62a5f4..665b16057 100644
--- a/lapack-netlib/TESTING/MATGEN/clatm1.c
+++ b/lapack-netlib/TESTING/MATGEN/clatm1.c
@@ -647,7 +647,7 @@ static integer c__3 = 3;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatm1_(integer *mode, real *cond, integer *irsign, 
+/* Subroutine */ void clatm1_(integer *mode, real *cond, integer *irsign, 
 	integer *idist, integer *iseed, complex *d__, integer *n, integer *
 	info)
 {
@@ -664,9 +664,9 @@ static integer c__3 = 3;
     complex ctemp;
     //extern /* Complex */ VOID clarnd_(complex *, integer *, integer *);
     extern complex clarnd_(integer *, integer *);
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern real slaran_(integer *);
-    extern /* Subroutine */ int clarnv_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void clarnv_(integer *, integer *, integer *, 
 	    complex *);
 
 
@@ -691,7 +691,7 @@ static integer c__3 = 3;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -711,8 +711,8 @@ static integer c__3 = 3;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLATM1", &i__1);
-	return 0;
+	xerbla_("CLATM1", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -848,7 +848,7 @@ L120:
 
     }
 
-    return 0;
+    return;
 
 /*     End of CLATM1 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatm5.c b/lapack-netlib/TESTING/MATGEN/clatm5.c
index dbd1ee642..c2b81ccf3 100644
--- a/lapack-netlib/TESTING/MATGEN/clatm5.c
+++ b/lapack-netlib/TESTING/MATGEN/clatm5.c
@@ -778,7 +778,7 @@ static complex c_b5 = {20.f,0.f};
 /* > \endverbatim */
 /* > */
 /*  ===================================================================== */
-/* Subroutine */ int clatm5_(integer *prtype, integer *m, integer *n, complex 
+/* Subroutine */ void clatm5_(integer *prtype, integer *m, integer *n, complex 
 	*a, integer *lda, complex *b, integer *ldb, complex *c__, integer *
 	ldc, complex *d__, integer *ldd, complex *e, integer *lde, complex *f,
 	 integer *ldf, complex *r__, integer *ldr, complex *l, integer *ldl, 
@@ -793,7 +793,7 @@ static complex c_b5 = {20.f,0.f};
 
     /* Local variables */
     integer i__, j, k;
-    extern /* Subroutine */ int cgemm_(char *, char *, integer *, integer *, 
+    extern /* Subroutine */ void cgemm_(char *, char *, integer *, integer *, 
 	    integer *, complex *, complex *, integer *, complex *, integer *, 
 	    complex *, complex *, integer *);
     complex imeps, reeps;
@@ -1277,6 +1277,6 @@ static complex c_b5 = {20.f,0.f};
 
 /*     End of CLATM5 */
 
-    return 0;
+    return;
 } /* clatm5_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatm6.c b/lapack-netlib/TESTING/MATGEN/clatm6.c
index 54c0af9da..96a3bd290 100644
--- a/lapack-netlib/TESTING/MATGEN/clatm6.c
+++ b/lapack-netlib/TESTING/MATGEN/clatm6.c
@@ -686,7 +686,7 @@ static integer c__24 = 24;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatm6_(integer *type__, integer *n, complex *a, integer 
+/* Subroutine */ void clatm6_(integer *type__, integer *n, complex *a, integer 
 	*lda, complex *b, complex *x, integer *ldx, complex *y, integer *ldy, 
 	complex *alpha, complex *beta, complex *wx, complex *wy, real *s, 
 	real *dif)
@@ -702,10 +702,10 @@ static integer c__24 = 24;
     complex work[26];
     integer i__, j;
     complex z__[64]	/* was [8][8] */;
-    extern /* Subroutine */ int clakf2_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void clakf2_(integer *, integer *, complex *, 
 	    integer *, complex *, complex *, complex *, complex *, integer *);
     real rwork[50];
-    extern /* Subroutine */ int cgesvd_(char *, char *, integer *, integer *, 
+    extern /* Subroutine */ void cgesvd_(char *, char *, integer *, integer *, 
 	    complex *, integer *, real *, complex *, integer *, complex *, 
 	    integer *, complex *, integer *, real *, integer *), clacpy_(char *, integer *, integer *, complex *, integer 
 	    *, complex *, integer *);
@@ -931,7 +931,7 @@ static integer c__24 = 24;
 	    &c__1, &work[2], &c__24, &rwork[8], &info);
     dif[5] = rwork[7];
 
-    return 0;
+    return;
 
 /*     End of CLATM6 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatme.c b/lapack-netlib/TESTING/MATGEN/clatme.c
index 72effbf9a..a905f5608 100644
--- a/lapack-netlib/TESTING/MATGEN/clatme.c
+++ b/lapack-netlib/TESTING/MATGEN/clatme.c
@@ -811,7 +811,7 @@ static integer c__5 = 5;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatme_(integer *n, char *dist, integer *iseed, complex *
+/* Subroutine */ void clatme_(integer *n, char *dist, integer *iseed, complex *
 	d__, integer *mode, real *cond, complex *dmax__, char *rsign, char *
 	upper, char *sim, real *ds, integer *modes, real *conds, integer *kl, 
 	integer *ku, real *anorm, complex *a, integer *lda, complex *work, 
@@ -827,22 +827,22 @@ static integer c__5 = 5;
     integer isim;
     real temp;
     integer i__, j;
-    extern /* Subroutine */ int cgerc_(integer *, integer *, complex *, 
+    extern /* Subroutine */ void cgerc_(integer *, integer *, complex *, 
 	    complex *, integer *, complex *, integer *, complex *, integer *);
     complex alpha;
-    extern /* Subroutine */ int cscal_(integer *, complex *, complex *, 
+    extern /* Subroutine */ void cscal_(integer *, complex *, complex *, 
 	    integer *);
     extern logical lsame_(char *, char *);
-    extern /* Subroutine */ int cgemv_(char *, integer *, integer *, complex *
+    extern /* Subroutine */ void cgemv_(char *, integer *, integer *, complex *
 	    , complex *, integer *, complex *, integer *, complex *, complex *
 	    , integer *);
     integer iinfo;
     real tempa[1];
     integer icols, idist;
-    extern /* Subroutine */ int ccopy_(integer *, complex *, integer *, 
+    extern /* Subroutine */ void ccopy_(integer *, complex *, integer *, 
 	    complex *, integer *);
     integer irows;
-    extern /* Subroutine */ int clatm1_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void clatm1_(integer *, real *, integer *, integer 
 	    *, integer *, complex *, integer *, integer *), slatm1_(integer *,
 	     real *, integer *, integer *, integer *, real *, integer *, 
 	    integer *);
@@ -850,17 +850,18 @@ static integer c__5 = 5;
     extern real clange_(char *, integer *, integer *, complex *, integer *, 
 	    real *);
     integer ir;
-    extern /* Subroutine */ int clarge_(integer *, complex *, integer *, 
+    extern /* Subroutine */ void clarge_(integer *, complex *, integer *, 
 	    integer *, complex *, integer *), clarfg_(integer *, complex *, 
 	    complex *, integer *, complex *), clacgv_(integer *, complex *, 
 	    integer *);
     //extern /* Complex */ VOID clarnd_(complex *, integer *, integer *);
     extern complex clarnd_(integer *, integer *);
     real ralpha;
-    extern /* Subroutine */ int csscal_(integer *, real *, complex *, integer 
+    extern /* Subroutine */ void csscal_(integer *, real *, complex *, integer 
 	    *), claset_(char *, integer *, integer *, complex *, complex *, 
-	    complex *, integer *), xerbla_(char *, integer *),
-	     clarnv_(integer *, integer *, integer *, complex *);
+	    complex *, integer *);
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void clarnv_(integer *, integer *, integer *, complex *);
     integer irsign, iupper;
     complex xnorms;
     integer jcr;
@@ -894,7 +895,7 @@ static integer c__5 = 5;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -986,8 +987,8 @@ static integer c__5 = 5;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLATME", &i__1);
-	return 0;
+	xerbla_("CLATME", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1008,7 +1009,7 @@ static integer c__5 = 5;
     clatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], n, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && abs(*mode) != 6) {
 
@@ -1028,7 +1029,7 @@ static integer c__5 = 5;
 	    alpha.r = q__1.r, alpha.i = q__1.i;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	cscal_(n, &alpha, &d__[1], &c__1);
@@ -1065,7 +1066,7 @@ static integer c__5 = 5;
 	slatm1_(modes, conds, &c__0, &c__0, &iseed[1], &ds[1], n, &iinfo);
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by V and V' */
@@ -1073,7 +1074,7 @@ static integer c__5 = 5;
 	clarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by S and (1/S) */
@@ -1086,7 +1087,7 @@ static integer c__5 = 5;
 		csscal_(n, &r__1, &a[j * a_dim1 + 1], &c__1);
 	    } else {
 		*info = 5;
-		return 0;
+		return;
 	    }
 /* L50: */
 	}
@@ -1096,7 +1097,7 @@ static integer c__5 = 5;
 	clarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1210,7 +1211,7 @@ static integer c__5 = 5;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of CLATME */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatmr.c b/lapack-netlib/TESTING/MATGEN/clatmr.c
index 1265052ab..4a117a9d4 100644
--- a/lapack-netlib/TESTING/MATGEN/clatmr.c
+++ b/lapack-netlib/TESTING/MATGEN/clatmr.c
@@ -998,7 +998,7 @@ static integer c__1 = 1;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatmr_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void clatmr_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, complex *d__, integer *mode, real *cond, complex *
 	dmax__, char *rsign, char *grade, complex *dl, integer *model, real *
 	condl, complex *dr, integer *moder, real *condr, char *pivtng, 
@@ -1022,7 +1022,7 @@ static integer c__1 = 1;
     integer mnsub;
     real onorm;
     integer mxsub, npvts;
-    extern /* Subroutine */ int clatm1_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void clatm1_(integer *, real *, integer *, integer 
 	    *, integer *, complex *, integer *, integer *);
     extern /* Complex */ VOID clatm2_(complex *, integer *, integer *, 
 	    integer *, integer *, integer *, integer *, integer *, integer *, 
@@ -1038,10 +1038,10 @@ static integer c__1 = 1;
     integer igrade;
     extern real clansb_(char *, char *, integer *, integer *, complex *, 
 	    integer *, real *);
-    extern /* Subroutine */ int csscal_(integer *, real *, complex *, integer 
+    extern /* Subroutine */ void csscal_(integer *, real *, complex *, integer 
 	    *);
     logical fulbnd;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     logical badpvt;
     extern real clansp_(char *, char *, integer *, complex *, real *), clansy_(char *, char *, integer *, complex *, integer *, 
 	    real *);
@@ -1077,7 +1077,7 @@ static integer c__1 = 1;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1278,8 +1278,8 @@ static integer c__1 = 1;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLATMR", &i__1);
-	return 0;
+	xerbla_("CLATMR", &i__1, 6);
+	return;
     }
 
 /*     Decide if we can pivot consistently */
@@ -1305,7 +1305,7 @@ static integer c__1 = 1;
     clatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, info);
     if (*info != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && *mode != -6 && *mode != 6) {
 
@@ -1321,7 +1321,7 @@ static integer c__1 = 1;
 	}
 	if (temp == 0.f && (dmax__->r != 0.f || dmax__->i != 0.f)) {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 	if (temp != 0.f) {
 	    q__1.r = dmax__->r / temp, q__1.i = dmax__->i / temp;
@@ -1361,7 +1361,7 @@ static integer c__1 = 1;
 	clatm1_(model, condl, &c__0, &idist, &iseed[1], &dl[1], m, info);
 	if (*info != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1371,7 +1371,7 @@ static integer c__1 = 1;
 	clatm1_(moder, condr, &c__0, &idist, &iseed[1], &dr[1], n, info);
 	if (*info != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2026,7 +2026,7 @@ static integer c__1 = 1;
 /*           Desired scaling impossible */
 
 	    *info = 5;
-	    return 0;
+	    return;
 
 	} else if (*anorm > 1.f && onorm < 1.f || *anorm < 1.f && onorm > 1.f)
 		 {
@@ -2099,6 +2099,6 @@ static integer c__1 = 1;
 
 /*     End of CLATMR */
 
-    return 0;
+    return;
 } /* clatmr_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatms.c b/lapack-netlib/TESTING/MATGEN/clatms.c
index 6f878e400..3f7011166 100644
--- a/lapack-netlib/TESTING/MATGEN/clatms.c
+++ b/lapack-netlib/TESTING/MATGEN/clatms.c
@@ -845,7 +845,7 @@ static logical c_false = FALSE_;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatms_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void clatms_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, real *d__, integer *mode, real *cond, real *dmax__, 
 	integer *kl, integer *ku, char *pack, complex *a, integer *lda, 
 	complex *work, integer *info)
@@ -870,14 +870,14 @@ static logical c_false = FALSE_;
     integer ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     complex ctemp;
     integer idist, mnmin, iskew;
     complex extra, dummy;
-    extern /* Subroutine */ int slatm1_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void slatm1_(integer *, real *, integer *, integer 
 	    *, integer *, real *, integer *, integer *);
     integer ic, jc, nc;
-    extern /* Subroutine */ int clagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void clagge_(integer *, integer *, integer *, 
 	    integer *, real *, complex *, integer *, integer *, complex *, 
 	    integer *), claghe_(integer *, integer *, real *, complex *, 
 	    integer *, integer *, complex *, integer *);
@@ -888,13 +888,14 @@ static logical c_false = FALSE_;
     extern complex clarnd_(integer *, integer *);
     integer minlda;
     complex st;
-    extern /* Subroutine */ int claset_(char *, integer *, integer *, complex 
+    extern /* Subroutine */ void claset_(char *, integer *, integer *, complex 
 	    *, complex *, complex *, integer *), clartg_(complex *, 
-	    complex *, real *, complex *, complex *), xerbla_(char *, integer 
-	    *), clagsy_(integer *, integer *, real *, complex *, 
+	    complex *, real *, complex *, complex *);
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void clagsy_(integer *, integer *, real *, complex *, 
 	    integer *, integer *, complex *, integer *);
     extern real slarnd_(integer *, integer *);
-    extern /* Subroutine */ int clarot_(logical *, logical *, logical *, 
+    extern /* Subroutine */ void clarot_(logical *, logical *, logical *, 
 	    integer *, complex *, complex *, complex *, integer *, complex *, 
 	    complex *);
     logical iltemp, givens;
@@ -929,7 +930,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1068,8 +1069,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLATMS", &i__1);
-	return 0;
+	xerbla_("CLATMS", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1090,7 +1091,7 @@ static logical c_false = FALSE_;
     slatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1119,7 +1120,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	sscal_(&mnmin, &alpha, &d__[1], &c__1);
@@ -2033,7 +2034,7 @@ static logical c_false = FALSE_;
 
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2208,7 +2209,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of CLATMS */
 
diff --git a/lapack-netlib/TESTING/MATGEN/clatmt.c b/lapack-netlib/TESTING/MATGEN/clatmt.c
index 04dc34697..3bb7ca4d9 100644
--- a/lapack-netlib/TESTING/MATGEN/clatmt.c
+++ b/lapack-netlib/TESTING/MATGEN/clatmt.c
@@ -853,7 +853,7 @@ static logical c_false = FALSE_;
 /* > \ingroup complex_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int clatmt_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void clatmt_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, real *d__, integer *mode, real *cond, real *dmax__, 
 	integer *rank, integer *kl, integer *ku, char *pack, complex *a, 
 	integer *lda, complex *work, integer *info)
@@ -876,16 +876,16 @@ static logical c_false = FALSE_;
     integer ipack, ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     complex ctemp;
     integer idist, mnmin;
     complex extra;
     integer iskew;
     complex dummy;
-    extern /* Subroutine */ int slatm7_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void slatm7_(integer *, real *, integer *, integer 
 	    *, integer *, real *, integer *, integer *, integer *);
     integer ic, jc, nc;
-    extern /* Subroutine */ int clagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void clagge_(integer *, integer *, integer *, 
 	    integer *, real *, complex *, integer *, integer *, complex *, 
 	    integer *), claghe_(integer *, integer *, real *, complex *, 
 	    integer *, integer *, complex *, integer *);
@@ -896,13 +896,14 @@ static logical c_false = FALSE_;
     extern complex clarnd_(integer *, integer *);
     integer minlda;
     complex st;
-    extern /* Subroutine */ int claset_(char *, integer *, integer *, complex 
+    extern /* Subroutine */ void claset_(char *, integer *, integer *, complex 
 	    *, complex *, complex *, integer *), clartg_(complex *, 
-	    complex *, real *, complex *, complex *), xerbla_(char *, integer 
-	    *), clagsy_(integer *, integer *, real *, complex *, 
+	    complex *, real *, complex *, complex *);
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void clagsy_(integer *, integer *, real *, complex *, 
 	    integer *, integer *, complex *, integer *);
     extern real slarnd_(integer *, integer *);
-    extern /* Subroutine */ int clarot_(logical *, logical *, logical *, 
+    extern /* Subroutine */ void clarot_(logical *, logical *, logical *, 
 	    integer *, complex *, complex *, complex *, integer *, complex *, 
 	    complex *);
     integer ioffst, irsign;
@@ -936,7 +937,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1075,8 +1076,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("CLATMT", &i__1);
-	return 0;
+	xerbla_("CLATMT", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1098,7 +1099,7 @@ static logical c_false = FALSE_;
 	    iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1127,7 +1128,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	sscal_(rank, &alpha, &d__[1], &c__1);
@@ -2041,7 +2042,7 @@ static logical c_false = FALSE_;
 
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2216,7 +2217,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of CLATMT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlagge.c b/lapack-netlib/TESTING/MATGEN/dlagge.c
index 448f70a39..dcca84326 100644
--- a/lapack-netlib/TESTING/MATGEN/dlagge.c
+++ b/lapack-netlib/TESTING/MATGEN/dlagge.c
@@ -626,7 +626,7 @@ static doublereal c_b13 = 0.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlagge_(integer *m, integer *n, integer *kl, integer *ku,
+/* Subroutine */ void dlagge_(integer *m, integer *n, integer *kl, integer *ku,
 	 doublereal *d__, doublereal *a, integer *lda, integer *iseed, 
 	doublereal *work, integer *info)
 {
@@ -635,17 +635,18 @@ static doublereal c_b13 = 0.;
     doublereal d__1;
 
     /* Local variables */
-    extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dger_(integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     extern doublereal dnrm2_(integer *, doublereal *, integer *);
     integer i__, j;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *), dgemv_(char *, integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    doublereal *, integer *);
     doublereal wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), dlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(
 	    integer *, integer *, integer *, doublereal *);
     doublereal tau;
 
@@ -684,8 +685,8 @@ static doublereal c_b13 = 0.;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("DLAGGE", &i__1);
-	return 0;
+	xerbla_("DLAGGE", &i__1, 6);
+	return;
     }
 
 /*     initialize A to diagonal matrix */
@@ -708,7 +709,7 @@ static doublereal c_b13 = 0.;
 /*     Quick exit if the user wants a diagonal matrix */
 
     if (*kl == 0 && *ku == 0) {
-	return 0;
+	return;
     }
 
 /*     pre- and post-multiply A by random orthogonal matrices */
@@ -963,7 +964,7 @@ static doublereal c_b13 = 0.;
 	}
 /* L70: */
     }
-    return 0;
+    return;
 
 /*     End of DLAGGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlagsy.c b/lapack-netlib/TESTING/MATGEN/dlagsy.c
index 802ce683c..2812db86d 100644
--- a/lapack-netlib/TESTING/MATGEN/dlagsy.c
+++ b/lapack-netlib/TESTING/MATGEN/dlagsy.c
@@ -615,7 +615,7 @@ static doublereal c_b26 = 1.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlagsy_(integer *n, integer *k, doublereal *d__, 
+/* Subroutine */ void dlagsy_(integer *n, integer *k, doublereal *d__, 
 	doublereal *a, integer *lda, integer *iseed, doublereal *work, 
 	integer *info)
 {
@@ -624,17 +624,17 @@ static doublereal c_b26 = 1.;
     doublereal d__1;
 
     /* Local variables */
-    extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dger_(integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     extern doublereal ddot_(integer *, doublereal *, integer *, doublereal *, 
 	    integer *), dnrm2_(integer *, doublereal *, integer *);
-    extern /* Subroutine */ int dsyr2_(char *, integer *, doublereal *, 
+    extern /* Subroutine */ void dsyr2_(char *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     integer i__, j;
     doublereal alpha;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *), dgemv_(char *, integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    doublereal *, integer *), daxpy_(integer *, doublereal *, 
@@ -642,7 +642,8 @@ static doublereal c_b26 = 1.;
 	    integer *, doublereal *, doublereal *, integer *, doublereal *, 
 	    integer *, doublereal *, doublereal *, integer *);
     doublereal wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), dlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(
 	    integer *, integer *, integer *, doublereal *);
     doublereal tau;
 
@@ -677,8 +678,8 @@ static doublereal c_b26 = 1.;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("DLAGSY", &i__1);
-	return 0;
+	xerbla_("DLAGSY", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -822,7 +823,7 @@ static doublereal c_b26 = 1.;
 	}
 /* L80: */
     }
-    return 0;
+    return;
 
 /*     End of DLAGSY */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlahilb.c b/lapack-netlib/TESTING/MATGEN/dlahilb.c
index 0f79547cf..dccfcd8c7 100644
--- a/lapack-netlib/TESTING/MATGEN/dlahilb.c
+++ b/lapack-netlib/TESTING/MATGEN/dlahilb.c
@@ -635,7 +635,7 @@ static doublereal c_b4 = 0.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlahilb_(integer *n, integer *nrhs, doublereal *a, 
+/* Subroutine */ void dlahilb_(integer *n, integer *nrhs, doublereal *a, 
 	integer *lda, doublereal *x, integer *ldx, doublereal *b, integer *
 	ldb, doublereal *work, integer *info)
 {
@@ -645,9 +645,9 @@ static doublereal c_b4 = 0.;
 
     /* Local variables */
     integer i__, j, m, r__, ti, tm;
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
-	    doublereal *, doublereal *, doublereal *, integer *), 
-	    xerbla_(char *, integer *);
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
+	    doublereal *, doublereal *, doublereal *, integer *); 
+    extern int xerbla_(char *, integer *, ftnlen);
 
 
 /*  -- LAPACK test routine (version 3.8.0) -- */
@@ -691,8 +691,8 @@ static doublereal c_b4 = 0.;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("DLAHILB", &i__1);
-	return 0;
+	xerbla_("DLAHILB", &i__1, 7);
+	return;
     }
     if (*n > 6) {
 	*info = 1;
@@ -745,6 +745,6 @@ static doublereal c_b4 = 0.;
 	}
     }
 
-    return 0;
+    return;
 } /* dlahilb_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlakf2.c b/lapack-netlib/TESTING/MATGEN/dlakf2.c
index 5d838ef35..0b403540a 100644
--- a/lapack-netlib/TESTING/MATGEN/dlakf2.c
+++ b/lapack-netlib/TESTING/MATGEN/dlakf2.c
@@ -615,7 +615,7 @@ static doublereal c_b3 = 0.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlakf2_(integer *m, integer *n, doublereal *a, integer *
+/* Subroutine */ void dlakf2_(integer *m, integer *n, doublereal *a, integer *
 	lda, doublereal *b, doublereal *d__, doublereal *e, doublereal *z__, 
 	integer *ldz)
 {
@@ -625,7 +625,7 @@ static doublereal c_b3 = 0.;
 
     /* Local variables */
     integer i__, j, l, ik, jk, mn;
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
 	    doublereal *, doublereal *, doublereal *, integer *);
     integer mn2;
 
@@ -731,7 +731,7 @@ static doublereal c_b3 = 0.;
 /* L90: */
     }
 
-    return 0;
+    return;
 
 /*     End of DLAKF2 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlarge.c b/lapack-netlib/TESTING/MATGEN/dlarge.c
index 580c50234..5d8a81387 100644
--- a/lapack-netlib/TESTING/MATGEN/dlarge.c
+++ b/lapack-netlib/TESTING/MATGEN/dlarge.c
@@ -600,7 +600,7 @@ static doublereal c_b10 = 0.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlarge_(integer *n, doublereal *a, integer *lda, integer 
+/* Subroutine */ void dlarge_(integer *n, doublereal *a, integer *lda, integer 
 	*iseed, doublereal *work, integer *info)
 {
     /* System generated locals */
@@ -608,17 +608,18 @@ static doublereal c_b10 = 0.;
     doublereal d__1;
 
     /* Local variables */
-    extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dger_(integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     extern doublereal dnrm2_(integer *, doublereal *, integer *);
     integer i__;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *), dgemv_(char *, integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    doublereal *, integer *);
     doublereal wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), dlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(
 	    integer *, integer *, integer *, doublereal *);
     doublereal tau;
 
@@ -650,8 +651,8 @@ static doublereal c_b10 = 0.;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("DLARGE", &i__1);
-	return 0;
+	xerbla_("DLARGE", &i__1, 6);
+	return;
     }
 
 /*     pre- and post-multiply A by random orthogonal matrix */
@@ -697,7 +698,7 @@ static doublereal c_b10 = 0.;
 		 a_dim1 + 1], lda);
 /* L10: */
     }
-    return 0;
+    return;
 
 /*     End of DLARGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlaror.c b/lapack-netlib/TESTING/MATGEN/dlaror.c
index 8e66e5511..d9e2e46ae 100644
--- a/lapack-netlib/TESTING/MATGEN/dlaror.c
+++ b/lapack-netlib/TESTING/MATGEN/dlaror.c
@@ -659,7 +659,7 @@ static integer c__1 = 1;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlaror_(char *side, char *init, integer *m, integer *n, 
+/* Subroutine */ void dlaror_(char *side, char *init, integer *m, integer *n, 
 	doublereal *a, integer *lda, integer *iseed, doublereal *x, integer *
 	info)
 {
@@ -669,24 +669,24 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer kbeg;
-    extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dger_(integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     integer jcol, irow;
     extern doublereal dnrm2_(integer *, doublereal *, integer *);
     integer j;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     extern logical lsame_(char *, char *);
-    extern /* Subroutine */ int dgemv_(char *, integer *, integer *, 
+    extern /* Subroutine */ void dgemv_(char *, integer *, integer *, 
 	    doublereal *, doublereal *, integer *, doublereal *, integer *, 
 	    doublereal *, doublereal *, integer *);
     integer ixfrm, itype, nxfrm;
     doublereal xnorm;
     extern doublereal dlarnd_(integer *, integer *);
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
-	    doublereal *, doublereal *, doublereal *, integer *), 
-	    xerbla_(char *, integer *);
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
+	    doublereal *, doublereal *, doublereal *, integer *); 
+    extern int xerbla_(char *, integer *, ftnlen);
     doublereal factor, xnorms;
 
 
@@ -709,7 +709,7 @@ static integer c__1 = 1;
     /* Function Body */
     *info = 0;
     if (*n == 0 || *m == 0) {
-	return 0;
+	return;
     }
 
     itype = 0;
@@ -734,8 +734,8 @@ static integer c__1 = 1;
     }
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLAROR", &i__1);
-	return 0;
+	xerbla_("DLAROR", &i__1, 6);
+	return;
     }
 
     if (itype == 1) {
@@ -782,8 +782,8 @@ static integer c__1 = 1;
 	factor = xnorms * (xnorms + x[kbeg]);
 	if (abs(factor) < 1e-20) {
 	    *info = 1;
-	    xerbla_("DLAROR", info);
-	    return 0;
+	    xerbla_("DLAROR", info, 6);
+	    return;
 	} else {
 	    factor = 1. / factor;
 	}
@@ -837,7 +837,7 @@ static integer c__1 = 1;
 /* L50: */
 	}
     }
-    return 0;
+    return;
 
 /*     End of DLAROR */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlarot.c b/lapack-netlib/TESTING/MATGEN/dlarot.c
index 16ff2eeb0..1d7a4939d 100644
--- a/lapack-netlib/TESTING/MATGEN/dlarot.c
+++ b/lapack-netlib/TESTING/MATGEN/dlarot.c
@@ -737,7 +737,7 @@ static integer c__1 = 1;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlarot_(logical *lrows, logical *lleft, logical *lright, 
+/* Subroutine */ void dlarot_(logical *lrows, logical *lleft, logical *lright, 
 	integer *nl, doublereal *c__, doublereal *s, doublereal *a, integer *
 	lda, doublereal *xleft, doublereal *xright)
 {
@@ -746,11 +746,11 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer iinc;
-    extern /* Subroutine */ int drot_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void drot_(integer *, doublereal *, integer *, 
 	    doublereal *, integer *, doublereal *, doublereal *);
     integer inext, ix, iy, nt;
     doublereal xt[2], yt[2];
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     integer iyt;
 
 
@@ -799,12 +799,12 @@ static integer c__1 = 1;
 /*     Check for errors */
 
     if (*nl < nt) {
-	xerbla_("DLAROT", &c__4);
-	return 0;
+	xerbla_("DLAROT", &c__4, 6);
+	return;
     }
     if (*lda <= 0 || ! (*lrows) && *lda < *nl - nt) {
-	xerbla_("DLAROT", &c__8);
-	return 0;
+	xerbla_("DLAROT", &c__8, 6);
+	return;
     }
 
 /*     Rotate */
@@ -825,7 +825,7 @@ static integer c__1 = 1;
 	a[iyt] = yt[nt - 1];
     }
 
-    return 0;
+    return;
 
 /*     End of DLAROT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatm1.c b/lapack-netlib/TESTING/MATGEN/dlatm1.c
index 5c2837e22..aa81469e7 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatm1.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatm1.c
@@ -641,7 +641,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatm1_(integer *mode, doublereal *cond, integer *irsign,
+/* Subroutine */ void dlatm1_(integer *mode, doublereal *cond, integer *irsign,
 	 integer *idist, integer *iseed, doublereal *d__, integer *n, integer 
 	*info)
 {
@@ -654,7 +654,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
     integer i__;
     doublereal alpha;
     extern doublereal dlaran_(integer *);
-    extern /* Subroutine */ int xerbla_(char *, integer *), dlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(
 	    integer *, integer *, integer *, doublereal *);
 
 
@@ -679,7 +680,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -699,8 +700,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATM1", &i__1);
-	return 0;
+	xerbla_("DLATM1", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -814,7 +815,7 @@ L120:
 
     }
 
-    return 0;
+    return;
 
 /*     End of DLATM1 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatm5.c b/lapack-netlib/TESTING/MATGEN/dlatm5.c
index 5043da398..94b49d6e3 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatm5.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatm5.c
@@ -778,7 +778,7 @@ static doublereal c_b33 = -1.;
 /* > \endverbatim */
 /* > */
 /*  ===================================================================== */
-/* Subroutine */ int dlatm5_(integer *prtype, integer *m, integer *n, 
+/* Subroutine */ void dlatm5_(integer *prtype, integer *m, integer *n, 
 	doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal *
 	c__, integer *ldc, doublereal *d__, integer *ldd, doublereal *e, 
 	integer *lde, doublereal *f, integer *ldf, doublereal *r__, integer *
@@ -792,7 +792,7 @@ static doublereal c_b33 = -1.;
 
     /* Local variables */
     integer i__, j, k;
-    extern /* Subroutine */ int dgemm_(char *, char *, integer *, integer *, 
+    extern /* Subroutine */ void dgemm_(char *, char *, integer *, integer *, 
 	    integer *, doublereal *, doublereal *, integer *, doublereal *, 
 	    integer *, doublereal *, doublereal *, integer *);
     doublereal imeps, reeps;
@@ -1100,6 +1100,6 @@ static doublereal c_b33 = -1.;
 
 /*     End of DLATM5 */
 
-    return 0;
+    return;
 } /* dlatm5_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatm6.c b/lapack-netlib/TESTING/MATGEN/dlatm6.c
index a1c995d24..fee10724a 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatm6.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatm6.c
@@ -692,7 +692,7 @@ static integer c__60 = 60;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatm6_(integer *type__, integer *n, doublereal *a, 
+/* Subroutine */ void dlatm6_(integer *type__, integer *n, doublereal *a, 
 	integer *lda, doublereal *b, doublereal *x, integer *ldx, doublereal *
 	y, integer *ldy, doublereal *alpha, doublereal *beta, doublereal *wx, 
 	doublereal *wy, doublereal *s, doublereal *dif)
@@ -706,7 +706,7 @@ static integer c__60 = 60;
     doublereal work[100];
     integer i__, j;
     doublereal z__[144]	/* was [12][12] */;
-    extern /* Subroutine */ int dlakf2_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dlakf2_(integer *, integer *, doublereal *, 
 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *,
 	     integer *), dgesvd_(char *, char *, integer *, integer *, 
 	    doublereal *, integer *, doublereal *, doublereal *, integer *, 
@@ -866,7 +866,7 @@ static integer c__60 = 60;
 
     }
 
-    return 0;
+    return;
 
 /*     End of DLATM6 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatm7.c b/lapack-netlib/TESTING/MATGEN/dlatm7.c
index df5eccf84..eb42e2b8b 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatm7.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatm7.c
@@ -627,7 +627,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatm7_(integer *mode, doublereal *cond, integer *irsign,
+/* Subroutine */ void dlatm7_(integer *mode, doublereal *cond, integer *irsign,
 	 integer *idist, integer *iseed, doublereal *d__, integer *n, integer 
 	*rank, integer *info)
 {
@@ -640,7 +640,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
     integer i__;
     doublereal alpha;
     extern doublereal dlaran_(integer *);
-    extern /* Subroutine */ int xerbla_(char *, integer *), dlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(
 	    integer *, integer *, integer *, doublereal *);
 
 
@@ -665,7 +666,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -685,8 +686,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATM7", &i__1);
-	return 0;
+	xerbla_("DLATM7", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -815,7 +816,7 @@ L240:
 
     }
 
-    return 0;
+    return;
 
 /*     End of DLATM7 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatme.c b/lapack-netlib/TESTING/MATGEN/dlatme.c
index c67900789..a92c70ef2 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatme.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatme.c
@@ -841,7 +841,7 @@ static doublereal c_b39 = 1.;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatme_(integer *n, char *dist, integer *iseed, 
+/* Subroutine */ void dlatme_(integer *n, char *dist, integer *iseed, 
 	doublereal *d__, integer *mode, doublereal *cond, doublereal *dmax__, 
 	char *ei, char *rsign, char *upper, char *sim, doublereal *ds, 
 	integer *modes, doublereal *conds, integer *kl, integer *ku, 
@@ -854,7 +854,7 @@ static doublereal c_b39 = 1.;
 
     /* Local variables */
     logical bads;
-    extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void dger_(integer *, integer *, doublereal *, 
 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
 	    integer *);
     integer isim;
@@ -862,10 +862,10 @@ static doublereal c_b39 = 1.;
     logical badei;
     integer i__, j;
     doublereal alpha;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     extern logical lsame_(char *, char *);
-    extern /* Subroutine */ int dgemv_(char *, integer *, integer *, 
+    extern /* Subroutine */ void dgemv_(char *, integer *, integer *, 
 	    doublereal *, doublereal *, integer *, doublereal *, integer *, 
 	    doublereal *, doublereal *, integer *);
     integer iinfo;
@@ -873,22 +873,23 @@ static doublereal c_b39 = 1.;
     integer icols;
     logical useei;
     integer idist;
-    extern /* Subroutine */ int dcopy_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dcopy_(integer *, doublereal *, integer *, 
 	    doublereal *, integer *);
     integer irows;
-    extern /* Subroutine */ int dlatm1_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm1_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *);
     integer ic, jc;
     extern doublereal dlange_(char *, integer *, integer *, doublereal *, 
 	    integer *, doublereal *);
     integer ir, jr;
-    extern /* Subroutine */ int dlarge_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlarge_(integer *, doublereal *, integer *, 
 	    integer *, doublereal *, integer *), dlarfg_(integer *, 
 	    doublereal *, doublereal *, integer *, doublereal *);
     extern doublereal dlaran_(integer *);
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
-	    doublereal *, doublereal *, doublereal *, integer *), 
-	    xerbla_(char *, integer *), dlarnv_(integer *, integer *, 
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
+	    doublereal *, doublereal *, doublereal *, integer *); 
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void dlarnv_(integer *, integer *, 
 	    integer *, doublereal *);
     integer irsign, iupper;
     doublereal xnorms;
@@ -924,7 +925,7 @@ static doublereal c_b39 = 1.;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1042,8 +1043,8 @@ static doublereal c_b39 = 1.;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATME", &i__1);
-	return 0;
+	xerbla_("DLATME", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1064,7 +1065,7 @@ static doublereal c_b39 = 1.;
     dlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], n, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && abs(*mode) != 6) {
 
@@ -1083,7 +1084,7 @@ static doublereal c_b39 = 1.;
 	    alpha = *dmax__ / temp;
 	} else if (*dmax__ != 0.) {
 	    *info = 2;
-	    return 0;
+	    return;
 	} else {
 	    alpha = 0.;
 	}
@@ -1155,7 +1156,7 @@ static doublereal c_b39 = 1.;
 	dlatm1_(modes, conds, &c__0, &c__0, &iseed[1], &ds[1], n, &iinfo);
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by V and V' */
@@ -1163,7 +1164,7 @@ static doublereal c_b39 = 1.;
 	dlarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by S and (1/S) */
@@ -1176,7 +1177,7 @@ static doublereal c_b39 = 1.;
 		dscal_(n, &d__1, &a[j * a_dim1 + 1], &c__1);
 	    } else {
 		*info = 5;
-		return 0;
+		return;
 	    }
 /* L80: */
 	}
@@ -1186,7 +1187,7 @@ static doublereal c_b39 = 1.;
 	dlarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1274,7 +1275,7 @@ static doublereal c_b39 = 1.;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of DLATME */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatmr.c b/lapack-netlib/TESTING/MATGEN/dlatmr.c
index c5d88ad29..e6ce5456d 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatmr.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatmr.c
@@ -979,7 +979,7 @@ static integer c__1 = 1;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatmr_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void dlatmr_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublereal *d__, integer *mode, doublereal *cond, 
 	doublereal *dmax__, char *rsign, char *grade, doublereal *dl, integer 
 	*model, doublereal *condl, doublereal *dr, integer *moder, doublereal 
@@ -996,7 +996,7 @@ static integer c__1 = 1;
     doublereal temp;
     integer isym, i__, j, k;
     doublereal alpha;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     integer ipack;
     extern logical lsame_(char *, char *);
@@ -1006,7 +1006,7 @@ static integer c__1 = 1;
     integer mnsub;
     doublereal onorm;
     integer mxsub, npvts;
-    extern /* Subroutine */ int dlatm1_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm1_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *);
     extern doublereal dlatm2_(integer *, integer *, integer *, integer *, 
 	    integer *, integer *, integer *, integer *, doublereal *, integer 
@@ -1021,7 +1021,7 @@ static integer c__1 = 1;
     extern doublereal dlansb_(char *, char *, integer *, integer *, 
 	    doublereal *, integer *, doublereal *);
     logical fulbnd;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     logical badpvt;
     extern doublereal dlansp_(char *, char *, integer *, doublereal *, 
 	    doublereal *), dlansy_(char *, char *, integer *, 
@@ -1058,7 +1058,7 @@ static integer c__1 = 1;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1252,8 +1252,8 @@ static integer c__1 = 1;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATMR", &i__1);
-	return 0;
+	xerbla_("DLATMR", &i__1, 6);
+	return;
     }
 
 /*     Decide if we can pivot consistently */
@@ -1279,7 +1279,7 @@ static integer c__1 = 1;
     dlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, info);
     if (*info != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && *mode != -6 && *mode != 6) {
 
@@ -1295,7 +1295,7 @@ static integer c__1 = 1;
 	}
 	if (temp == 0. && *dmax__ != 0.) {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 	if (temp != 0.) {
 	    alpha = *dmax__ / temp;
@@ -1316,7 +1316,7 @@ static integer c__1 = 1;
 	dlatm1_(model, condl, &c__0, &idist, &iseed[1], &dl[1], m, info);
 	if (*info != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1326,7 +1326,7 @@ static integer c__1 = 1;
 	dlatm1_(moder, condr, &c__0, &idist, &iseed[1], &dr[1], n, info);
 	if (*info != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1815,7 +1815,7 @@ static integer c__1 = 1;
 /*           Desired scaling impossible */
 
 	    *info = 5;
-	    return 0;
+	    return;
 
 	} else if (*anorm > 1. && onorm < 1. || *anorm < 1. && onorm > 1.) {
 
@@ -1887,6 +1887,6 @@ static integer c__1 = 1;
 
 /*     End of DLATMR */
 
-    return 0;
+    return;
 } /* dlatmr_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatms.c b/lapack-netlib/TESTING/MATGEN/dlatms.c
index 330bd279a..d56d26a50 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatms.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatms.c
@@ -833,7 +833,7 @@ static logical c_false = FALSE_;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatms_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void dlatms_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublereal *d__, integer *mode, doublereal *cond, 
 	doublereal *dmax__, integer *kl, integer *ku, char *pack, doublereal *
 	a, integer *lda, doublereal *work, integer *info)
@@ -851,27 +851,29 @@ static logical c_false = FALSE_;
     integer i__, j, k;
     doublereal s, alpha, angle;
     integer ipack;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     integer ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo, idist, mnmin;
-    extern /* Subroutine */ int dcopy_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dcopy_(integer *, doublereal *, integer *, 
 	    doublereal *, integer *);
     integer iskew;
     doublereal extra, dummy;
-    extern /* Subroutine */ int dlatm1_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm1_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *);
     integer ic, jc, nc;
-    extern /* Subroutine */ int dlagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void dlagge_(integer *, integer *, integer *, 
 	    integer *, doublereal *, doublereal *, integer *, integer *, 
 	    doublereal *, integer *);
     integer il, iendch, ir, jr, ipackg, mr, minlda;
     extern doublereal dlarnd_(integer *, integer *);
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
 	    doublereal *, doublereal *, doublereal *, integer *), 
 	    dlartg_(doublereal *, doublereal *, doublereal *, doublereal *, 
-	    doublereal *), xerbla_(char *, integer *), dlagsy_(
+	    doublereal *);
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void dlagsy_(
 	    integer *, integer *, doublereal *, doublereal *, integer *, 
 	    integer *, doublereal *, integer *), dlarot_(logical *, logical *,
 	     logical *, integer *, doublereal *, doublereal *, doublereal *, 
@@ -908,7 +910,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1043,8 +1045,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATMS", &i__1);
-	return 0;
+	xerbla_("DLATMS", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1065,7 +1067,7 @@ static logical c_false = FALSE_;
     dlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1094,7 +1096,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	dscal_(&mnmin, &alpha, &d__[1], &c__1);
@@ -1723,7 +1725,7 @@ static logical c_false = FALSE_;
 	}
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1885,7 +1887,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of DLATMS */
 
diff --git a/lapack-netlib/TESTING/MATGEN/dlatmt.c b/lapack-netlib/TESTING/MATGEN/dlatmt.c
index 5037405e3..b61e8b779 100644
--- a/lapack-netlib/TESTING/MATGEN/dlatmt.c
+++ b/lapack-netlib/TESTING/MATGEN/dlatmt.c
@@ -843,7 +843,7 @@ static logical c_false = FALSE_;
 /* > \ingroup double_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int dlatmt_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void dlatmt_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublereal *d__, integer *mode, doublereal *cond, 
 	doublereal *dmax__, integer *rank, integer *kl, integer *ku, char *
 	pack, doublereal *a, integer *lda, doublereal *work, integer *info)
@@ -861,28 +861,30 @@ static logical c_false = FALSE_;
     integer i__, j, k;
     doublereal s, alpha, angle;
     integer ipack;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     integer ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo, idist, mnmin;
-    extern /* Subroutine */ int dcopy_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dcopy_(integer *, doublereal *, integer *, 
 	    doublereal *, integer *);
     integer iskew;
     doublereal extra, dummy;
-    extern /* Subroutine */ int dlatm7_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm7_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *, integer 
 	    *);
     integer ic, jc, nc;
-    extern /* Subroutine */ int dlagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void dlagge_(integer *, integer *, integer *, 
 	    integer *, doublereal *, doublereal *, integer *, integer *, 
 	    doublereal *, integer *);
     integer il, iendch, ir, jr, ipackg, mr, minlda;
     extern doublereal dlarnd_(integer *, integer *);
-    extern /* Subroutine */ int dlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void dlaset_(char *, integer *, integer *, 
 	    doublereal *, doublereal *, doublereal *, integer *), 
 	    dlartg_(doublereal *, doublereal *, doublereal *, doublereal *, 
-	    doublereal *), xerbla_(char *, integer *), dlagsy_(
+	    doublereal *);
+    extern int xerbla_(char *, integer *, ftnlen);
+    extern void dlagsy_(
 	    integer *, integer *, doublereal *, doublereal *, integer *, 
 	    integer *, doublereal *, integer *), dlarot_(logical *, logical *,
 	     logical *, integer *, doublereal *, doublereal *, doublereal *, 
@@ -918,7 +920,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1053,8 +1055,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("DLATMT", &i__1);
-	return 0;
+	xerbla_("DLATMT", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1076,7 +1078,7 @@ static logical c_false = FALSE_;
 	    iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1105,7 +1107,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	dscal_(rank, &alpha, &d__[1], &c__1);
@@ -1734,7 +1736,7 @@ static logical c_false = FALSE_;
 	}
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1896,7 +1898,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of DLATMT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slagge.c b/lapack-netlib/TESTING/MATGEN/slagge.c
index 9e3e47712..c3ec8f2b8 100644
--- a/lapack-netlib/TESTING/MATGEN/slagge.c
+++ b/lapack-netlib/TESTING/MATGEN/slagge.c
@@ -626,7 +626,7 @@ static real c_b13 = 0.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slagge_(integer *m, integer *n, integer *kl, integer *ku,
+/* Subroutine */ void slagge_(integer *m, integer *n, integer *kl, integer *ku,
 	 real *d__, real *a, integer *lda, integer *iseed, real *work, 
 	integer *info)
 {
@@ -635,15 +635,16 @@ static real c_b13 = 0.f;
     real r__1;
 
     /* Local variables */
-    extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, 
+    extern /* Subroutine */ void sger_(integer *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     extern real snrm2_(integer *, real *, integer *);
     integer i__, j;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), 
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *), 
 	    sgemv_(char *, integer *, integer *, real *, real *, integer *, 
 	    real *, integer *, real *, real *, integer *);
     real wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), slarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void slarnv_(
 	    integer *, integer *, integer *, real *);
     real tau;
 
@@ -682,8 +683,8 @@ static real c_b13 = 0.f;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("SLAGGE", &i__1);
-	return 0;
+	xerbla_("SLAGGE", &i__1, 6);
+	return;
     }
 
 /*     initialize A to diagonal matrix */
@@ -706,7 +707,7 @@ static real c_b13 = 0.f;
 /*     Quick exit if the user wants a diagonal matrix */
 
     if (*kl == 0 && *ku == 0) {
-	return 0;
+	return;
     }
 
 /*     pre- and post-multiply A by random orthogonal matrices */
@@ -961,7 +962,7 @@ static real c_b13 = 0.f;
 	}
 /* L70: */
     }
-    return 0;
+    return;
 
 /*     End of SLAGGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slagsy.c b/lapack-netlib/TESTING/MATGEN/slagsy.c
index 5038e1c65..62614593d 100644
--- a/lapack-netlib/TESTING/MATGEN/slagsy.c
+++ b/lapack-netlib/TESTING/MATGEN/slagsy.c
@@ -615,7 +615,7 @@ static real c_b26 = 1.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slagsy_(integer *n, integer *k, real *d__, real *a, 
+/* Subroutine */ void slagsy_(integer *n, integer *k, real *d__, real *a, 
 	integer *lda, integer *iseed, real *work, integer *info)
 {
     /* System generated locals */
@@ -623,22 +623,23 @@ static real c_b26 = 1.f;
     real r__1;
 
     /* Local variables */
-    extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, 
+    extern /* Subroutine */ void sger_(integer *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     extern real sdot_(integer *, real *, integer *, real *, integer *), 
 	    snrm2_(integer *, real *, integer *);
     integer i__, j;
-    extern /* Subroutine */ int ssyr2_(char *, integer *, real *, real *, 
+    extern /* Subroutine */ void ssyr2_(char *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     real alpha;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), 
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *), 
 	    sgemv_(char *, integer *, integer *, real *, real *, integer *, 
 	    real *, integer *, real *, real *, integer *), saxpy_(
 	    integer *, real *, real *, integer *, real *, integer *), ssymv_(
 	    char *, integer *, real *, real *, integer *, real *, integer *, 
 	    real *, real *, integer *);
     real wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), slarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void slarnv_(
 	    integer *, integer *, integer *, real *);
     real tau;
 
@@ -673,8 +674,8 @@ static real c_b26 = 1.f;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("SLAGSY", &i__1);
-	return 0;
+	xerbla_("SLAGSY", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -818,7 +819,7 @@ static real c_b26 = 1.f;
 	}
 /* L80: */
     }
-    return 0;
+    return;
 
 /*     End of SLAGSY */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slahilb.c b/lapack-netlib/TESTING/MATGEN/slahilb.c
index 786804ca3..ae633d3c9 100644
--- a/lapack-netlib/TESTING/MATGEN/slahilb.c
+++ b/lapack-netlib/TESTING/MATGEN/slahilb.c
@@ -635,7 +635,7 @@ static real c_b4 = 0.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slahilb_(integer *n, integer *nrhs, real *a, integer *
+/* Subroutine */ void slahilb_(integer *n, integer *nrhs, real *a, integer *
 	lda, real *x, integer *ldx, real *b, integer *ldb, real *work, 
 	integer *info)
 {
@@ -645,7 +645,8 @@ static real c_b4 = 0.f;
 
     /* Local variables */
     integer i__, j, m, r__, ti, tm;
-    extern /* Subroutine */ int xerbla_(char *, integer *), slaset_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void slaset_(
 	    char *, integer *, integer *, real *, real *, real *, integer *);
 
 
@@ -690,8 +691,8 @@ static real c_b4 = 0.f;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("SLAHILB", &i__1);
-	return 0;
+	xerbla_("SLAHILB", &i__1, 7);
+	return;
     }
     if (*n > 6) {
 	*info = 1;
@@ -745,6 +746,6 @@ static real c_b4 = 0.f;
 	}
     }
 
-    return 0;
+    return;
 } /* slahilb_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slakf2.c b/lapack-netlib/TESTING/MATGEN/slakf2.c
index 27efd82a8..67be11662 100644
--- a/lapack-netlib/TESTING/MATGEN/slakf2.c
+++ b/lapack-netlib/TESTING/MATGEN/slakf2.c
@@ -615,7 +615,7 @@ static real c_b3 = 0.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slakf2_(integer *m, integer *n, real *a, integer *lda, 
+/* Subroutine */ void slakf2_(integer *m, integer *n, real *a, integer *lda, 
 	real *b, real *d__, real *e, real *z__, integer *ldz)
 {
     /* System generated locals */
@@ -624,7 +624,7 @@ static real c_b3 = 0.f;
 
     /* Local variables */
     integer i__, j, l, ik, jk, mn;
-    extern /* Subroutine */ int slaset_(char *, integer *, integer *, real *, 
+    extern /* Subroutine */ void slaset_(char *, integer *, integer *, real *, 
 	    real *, real *, integer *);
     integer mn2;
 
@@ -730,7 +730,7 @@ static real c_b3 = 0.f;
 /* L90: */
     }
 
-    return 0;
+    return;
 
 /*     End of SLAKF2 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slarge.c b/lapack-netlib/TESTING/MATGEN/slarge.c
index 20453f08b..6b37e9400 100644
--- a/lapack-netlib/TESTING/MATGEN/slarge.c
+++ b/lapack-netlib/TESTING/MATGEN/slarge.c
@@ -600,7 +600,7 @@ static real c_b10 = 0.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slarge_(integer *n, real *a, integer *lda, integer *
+/* Subroutine */ void slarge_(integer *n, real *a, integer *lda, integer *
 	iseed, real *work, integer *info)
 {
     /* System generated locals */
@@ -608,15 +608,16 @@ static real c_b10 = 0.f;
     real r__1;
 
     /* Local variables */
-    extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, 
+    extern /* Subroutine */ void sger_(integer *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     extern real snrm2_(integer *, real *, integer *);
     integer i__;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), 
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *), 
 	    sgemv_(char *, integer *, integer *, real *, real *, integer *, 
 	    real *, integer *, real *, real *, integer *);
     real wa, wb, wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), slarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void slarnv_(
 	    integer *, integer *, integer *, real *);
     real tau;
 
@@ -648,8 +649,8 @@ static real c_b10 = 0.f;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("SLARGE", &i__1);
-	return 0;
+	xerbla_("SLARGE", &i__1, 6);
+	return;
     }
 
 /*     pre- and post-multiply A by random orthogonal matrix */
@@ -695,7 +696,7 @@ static real c_b10 = 0.f;
 		 a_dim1 + 1], lda);
 /* L10: */
     }
-    return 0;
+    return;
 
 /*     End of SLARGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slaror.c b/lapack-netlib/TESTING/MATGEN/slaror.c
index 26db2e645..48b532dfd 100644
--- a/lapack-netlib/TESTING/MATGEN/slaror.c
+++ b/lapack-netlib/TESTING/MATGEN/slaror.c
@@ -659,7 +659,7 @@ static integer c__1 = 1;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slaror_(char *side, char *init, integer *m, integer *n, 
+/* Subroutine */ void slaror_(char *side, char *init, integer *m, integer *n, 
 	real *a, integer *lda, integer *iseed, real *x, integer *info)
 {
     /* System generated locals */
@@ -668,21 +668,21 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer kbeg, jcol;
-    extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, 
+    extern /* Subroutine */ void sger_(integer *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     integer irow;
     extern real snrm2_(integer *, real *, integer *);
     integer j;
     extern logical lsame_(char *, char *);
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), 
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *), 
 	    sgemv_(char *, integer *, integer *, real *, real *, integer *, 
 	    real *, integer *, real *, real *, integer *);
     integer ixfrm, itype, nxfrm;
     real xnorm;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     real factor;
     extern real slarnd_(integer *, integer *);
-    extern /* Subroutine */ int slaset_(char *, integer *, integer *, real *, 
+    extern /* Subroutine */ void slaset_(char *, integer *, integer *, real *, 
 	    real *, real *, integer *);
     real xnorms;
 
@@ -706,7 +706,7 @@ static integer c__1 = 1;
     /* Function Body */
     *info = 0;
     if (*n == 0 || *m == 0) {
-	return 0;
+	return;
     }
 
     itype = 0;
@@ -731,8 +731,8 @@ static integer c__1 = 1;
     }
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLAROR", &i__1);
-	return 0;
+	xerbla_("SLAROR", &i__1, 6);
+	return;
     }
 
     if (itype == 1) {
@@ -779,8 +779,8 @@ static integer c__1 = 1;
 	factor = xnorms * (xnorms + x[kbeg]);
 	if (abs(factor) < 1e-20f) {
 	    *info = 1;
-	    xerbla_("SLAROR", info);
-	    return 0;
+	    xerbla_("SLAROR", info, 6);
+	    return;
 	} else {
 	    factor = 1.f / factor;
 	}
@@ -834,7 +834,7 @@ static integer c__1 = 1;
 /* L50: */
 	}
     }
-    return 0;
+    return;
 
 /*     End of SLAROR */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slarot.c b/lapack-netlib/TESTING/MATGEN/slarot.c
index 205f7a622..a23e186e1 100644
--- a/lapack-netlib/TESTING/MATGEN/slarot.c
+++ b/lapack-netlib/TESTING/MATGEN/slarot.c
@@ -737,7 +737,7 @@ static integer c__1 = 1;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slarot_(logical *lrows, logical *lleft, logical *lright, 
+/* Subroutine */ void slarot_(logical *lrows, logical *lleft, logical *lright, 
 	integer *nl, real *c__, real *s, real *a, integer *lda, real *xleft, 
 	real *xright)
 {
@@ -746,11 +746,11 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer iinc;
-    extern /* Subroutine */ int srot_(integer *, real *, integer *, real *, 
+    extern /* Subroutine */ void srot_(integer *, real *, integer *, real *, 
 	    integer *, real *, real *);
     integer inext, ix, iy, nt;
     real xt[2], yt[2];
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     integer iyt;
 
 
@@ -799,12 +799,12 @@ static integer c__1 = 1;
 /*     Check for errors */
 
     if (*nl < nt) {
-	xerbla_("SLAROT", &c__4);
-	return 0;
+	xerbla_("SLAROT", &c__4, 6);
+	return;
     }
     if (*lda <= 0 || ! (*lrows) && *lda < *nl - nt) {
-	xerbla_("SLAROT", &c__8);
-	return 0;
+	xerbla_("SLAROT", &c__8, 6);
+	return;
     }
 
 /*     Rotate */
@@ -825,7 +825,7 @@ static integer c__1 = 1;
 	a[iyt] = yt[nt - 1];
     }
 
-    return 0;
+    return;
 
 /*     End of SLAROT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatm1.c b/lapack-netlib/TESTING/MATGEN/slatm1.c
index c5bb56d30..5c4cdc067 100644
--- a/lapack-netlib/TESTING/MATGEN/slatm1.c
+++ b/lapack-netlib/TESTING/MATGEN/slatm1.c
@@ -641,7 +641,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatm1_(integer *mode, real *cond, integer *irsign, 
+/* Subroutine */ void slatm1_(integer *mode, real *cond, integer *irsign, 
 	integer *idist, integer *iseed, real *d__, integer *n, integer *info)
 {
     /* System generated locals */
@@ -652,9 +652,9 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
     real temp;
     integer i__;
     real alpha;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern real slaran_(integer *);
-    extern /* Subroutine */ int slarnv_(integer *, integer *, integer *, real 
+    extern /* Subroutine */ void slarnv_(integer *, integer *, integer *, real 
 	    *);
 
 
@@ -679,7 +679,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -699,8 +699,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATM1", &i__1);
-	return 0;
+	xerbla_("SLATM1", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -815,7 +815,7 @@ L120:
 
     }
 
-    return 0;
+    return;
 
 /*     End of SLATM1 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatm5.c b/lapack-netlib/TESTING/MATGEN/slatm5.c
index 238bf7e97..24ee0915d 100644
--- a/lapack-netlib/TESTING/MATGEN/slatm5.c
+++ b/lapack-netlib/TESTING/MATGEN/slatm5.c
@@ -778,7 +778,7 @@ static real c_b33 = -1.f;
 /* > \endverbatim */
 /* > */
 /*  ===================================================================== */
-/* Subroutine */ int slatm5_(integer *prtype, integer *m, integer *n, real *a,
+/* Subroutine */ void slatm5_(integer *prtype, integer *m, integer *n, real *a,
 	 integer *lda, real *b, integer *ldb, real *c__, integer *ldc, real *
 	d__, integer *ldd, real *e, integer *lde, real *f, integer *ldf, real 
 	*r__, integer *ldr, real *l, integer *ldl, real *alpha, integer *
@@ -791,7 +791,7 @@ static real c_b33 = -1.f;
 
     /* Local variables */
     integer i__, j, k;
-    extern /* Subroutine */ int sgemm_(char *, char *, integer *, integer *, 
+    extern /* Subroutine */ void sgemm_(char *, char *, integer *, integer *, 
 	    integer *, real *, real *, integer *, real *, integer *, real *, 
 	    real *, integer *);
     real imeps, reeps;
@@ -1091,6 +1091,6 @@ static real c_b33 = -1.f;
 
 /*     End of SLATM5 */
 
-    return 0;
+    return;
 } /* slatm5_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatm6.c b/lapack-netlib/TESTING/MATGEN/slatm6.c
index d69028e6c..c50fcd29f 100644
--- a/lapack-netlib/TESTING/MATGEN/slatm6.c
+++ b/lapack-netlib/TESTING/MATGEN/slatm6.c
@@ -692,7 +692,7 @@ static integer c__60 = 60;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatm6_(integer *type__, integer *n, real *a, integer *
+/* Subroutine */ void slatm6_(integer *type__, integer *n, real *a, integer *
 	lda, real *b, real *x, integer *ldx, real *y, integer *ldy, real *
 	alpha, real *beta, real *wx, real *wy, real *s, real *dif)
 {
@@ -705,7 +705,7 @@ static integer c__60 = 60;
     real work[100];
     integer i__, j;
     real z__[144]	/* was [12][12] */;
-    extern /* Subroutine */ int slakf2_(integer *, integer *, real *, integer 
+    extern /* Subroutine */ void slakf2_(integer *, integer *, real *, integer 
 	    *, real *, real *, real *, real *, integer *), sgesvd_(char *, 
 	    char *, integer *, integer *, real *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *, integer *), slacpy_(char *, integer *, integer *, real *, 
@@ -864,7 +864,7 @@ static integer c__60 = 60;
 
     }
 
-    return 0;
+    return;
 
 /*     End of SLATM6 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatm7.c b/lapack-netlib/TESTING/MATGEN/slatm7.c
index 21976e322..e1b32ac9f 100644
--- a/lapack-netlib/TESTING/MATGEN/slatm7.c
+++ b/lapack-netlib/TESTING/MATGEN/slatm7.c
@@ -627,7 +627,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatm7_(integer *mode, real *cond, integer *irsign, 
+/* Subroutine */ void slatm7_(integer *mode, real *cond, integer *irsign, 
 	integer *idist, integer *iseed, real *d__, integer *n, integer *rank, 
 	integer *info)
 {
@@ -639,9 +639,9 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
     real temp;
     integer i__;
     real alpha;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern real slaran_(integer *);
-    extern /* Subroutine */ int slarnv_(integer *, integer *, integer *, real 
+    extern /* Subroutine */ void slarnv_(integer *, integer *, integer *, real 
 	    *);
 
 
@@ -666,7 +666,7 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -686,8 +686,8 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATM7", &i__1);
-	return 0;
+	xerbla_("SLATM7", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -817,7 +817,7 @@ L240:
 
     }
 
-    return 0;
+    return;
 
 /*     End of SLATM7 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatme.c b/lapack-netlib/TESTING/MATGEN/slatme.c
index 96178bee6..a8a6b39a3 100644
--- a/lapack-netlib/TESTING/MATGEN/slatme.c
+++ b/lapack-netlib/TESTING/MATGEN/slatme.c
@@ -841,7 +841,7 @@ static real c_b39 = 1.f;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatme_(integer *n, char *dist, integer *iseed, real *
+/* Subroutine */ void slatme_(integer *n, char *dist, integer *iseed, real *
 	d__, integer *mode, real *cond, real *dmax__, char *ei, char *rsign, 
 	char *upper, char *sim, real *ds, integer *modes, real *conds, 
 	integer *kl, integer *ku, real *anorm, real *a, integer *lda, real *
@@ -853,7 +853,7 @@ static real c_b39 = 1.f;
 
     /* Local variables */
     logical bads;
-    extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, 
+    extern /* Subroutine */ void sger_(integer *, integer *, real *, real *, 
 	    integer *, real *, integer *, real *, integer *);
     integer isim;
     real temp;
@@ -862,28 +862,29 @@ static real c_b39 = 1.f;
     real alpha;
     extern logical lsame_(char *, char *);
     integer iinfo;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     real tempa[1];
     integer icols;
     logical useei;
     integer idist;
-    extern /* Subroutine */ int sgemv_(char *, integer *, integer *, real *, 
+    extern /* Subroutine */ void sgemv_(char *, integer *, integer *, real *, 
 	    real *, integer *, real *, integer *, real *, real *, integer *), scopy_(integer *, real *, integer *, real *, integer *);
     integer irows;
-    extern /* Subroutine */ int slatm1_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void slatm1_(integer *, real *, integer *, integer 
 	    *, integer *, real *, integer *, integer *);
     integer ic, jc, ir, jr;
     extern real slange_(char *, integer *, integer *, real *, integer *, real 
 	    *);
-    extern /* Subroutine */ int slarge_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void slarge_(integer *, real *, integer *, integer 
 	    *, real *, integer *), slarfg_(integer *, real *, real *, integer 
-	    *, real *), xerbla_(char *, integer *);
+	    *, real *);
+    extern int xerbla_(char *, integer *, ftnlen);
     extern real slaran_(integer *);
     integer irsign;
-    extern /* Subroutine */ int slaset_(char *, integer *, integer *, real *, 
+    extern /* Subroutine */ void slaset_(char *, integer *, integer *, real *, 
 	    real *, real *, integer *);
     integer iupper;
-    extern /* Subroutine */ int slarnv_(integer *, integer *, integer *, real 
+    extern /* Subroutine */ void slarnv_(integer *, integer *, integer *, real 
 	    *);
     real xnorms;
     integer jcr;
@@ -918,7 +919,7 @@ static real c_b39 = 1.f;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1036,8 +1037,8 @@ static real c_b39 = 1.f;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATME", &i__1);
-	return 0;
+	xerbla_("SLATME", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1058,7 +1059,7 @@ static real c_b39 = 1.f;
     slatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], n, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && abs(*mode) != 6) {
 
@@ -1077,7 +1078,7 @@ static real c_b39 = 1.f;
 	    alpha = *dmax__ / temp;
 	} else if (*dmax__ != 0.f) {
 	    *info = 2;
-	    return 0;
+	    return;
 	} else {
 	    alpha = 0.f;
 	}
@@ -1149,7 +1150,7 @@ static real c_b39 = 1.f;
 	slatm1_(modes, conds, &c__0, &c__0, &iseed[1], &ds[1], n, &iinfo);
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by V and V' */
@@ -1157,7 +1158,7 @@ static real c_b39 = 1.f;
 	slarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by S and (1/S) */
@@ -1170,7 +1171,7 @@ static real c_b39 = 1.f;
 		sscal_(n, &r__1, &a[j * a_dim1 + 1], &c__1);
 	    } else {
 		*info = 5;
-		return 0;
+		return;
 	    }
 /* L80: */
 	}
@@ -1180,7 +1181,7 @@ static real c_b39 = 1.f;
 	slarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1268,7 +1269,7 @@ static real c_b39 = 1.f;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of SLATME */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatmr.c b/lapack-netlib/TESTING/MATGEN/slatmr.c
index 9d416d23f..01de66f54 100644
--- a/lapack-netlib/TESTING/MATGEN/slatmr.c
+++ b/lapack-netlib/TESTING/MATGEN/slatmr.c
@@ -979,7 +979,7 @@ static integer c__1 = 1;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatmr_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void slatmr_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, real *d__, integer *mode, real *cond, real *dmax__, 
 	char *rsign, char *grade, real *dl, integer *model, real *condl, real 
 	*dr, integer *moder, real *condr, char *pivtng, integer *ipivot, 
@@ -998,13 +998,13 @@ static integer c__1 = 1;
     integer ipack;
     extern logical lsame_(char *, char *);
     real tempa[1];
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     integer iisub, idist, jjsub, mnmin;
     logical dzero;
     integer mnsub;
     real onorm;
     integer mxsub, npvts;
-    extern /* Subroutine */ int slatm1_(integer *, real *, integer *, integer 
+    extern /* Subroutine */ void slatm1_(integer *, real *, integer *, integer 
 	    *, integer *, real *, integer *, integer *);
     extern real slatm2_(integer *, integer *, integer *, integer *, integer *,
 	     integer *, integer *, integer *, real *, integer *, real *, real 
@@ -1017,7 +1017,7 @@ static integer c__1 = 1;
 	    integer *, real *), slange_(char *, integer *, integer *, 
 	    real *, integer *, real *);
     logical fulbnd;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     logical badpvt;
     extern real slansb_(char *, char *, integer *, integer *, real *, integer 
 	    *, real *);
@@ -1057,7 +1057,7 @@ static integer c__1 = 1;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1251,8 +1251,8 @@ static integer c__1 = 1;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATMR", &i__1);
-	return 0;
+	xerbla_("SLATMR", &i__1, 6);
+	return;
     }
 
 /*     Decide if we can pivot consistently */
@@ -1278,7 +1278,7 @@ static integer c__1 = 1;
     slatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, info);
     if (*info != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && *mode != -6 && *mode != 6) {
 
@@ -1294,7 +1294,7 @@ static integer c__1 = 1;
 	}
 	if (temp == 0.f && *dmax__ != 0.f) {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 	if (temp != 0.f) {
 	    alpha = *dmax__ / temp;
@@ -1315,7 +1315,7 @@ static integer c__1 = 1;
 	slatm1_(model, condl, &c__0, &idist, &iseed[1], &dl[1], m, info);
 	if (*info != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1325,7 +1325,7 @@ static integer c__1 = 1;
 	slatm1_(moder, condr, &c__0, &idist, &iseed[1], &dr[1], n, info);
 	if (*info != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1814,7 +1814,7 @@ static integer c__1 = 1;
 /*           Desired scaling impossible */
 
 	    *info = 5;
-	    return 0;
+	    return;
 
 	} else if (*anorm > 1.f && onorm < 1.f || *anorm < 1.f && onorm > 1.f)
 		 {
@@ -1887,6 +1887,6 @@ static integer c__1 = 1;
 
 /*     End of SLATMR */
 
-    return 0;
+    return;
 } /* slatmr_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatms.c b/lapack-netlib/TESTING/MATGEN/slatms.c
index ade2391cc..8e85f5db8 100644
--- a/lapack-netlib/TESTING/MATGEN/slatms.c
+++ b/lapack-netlib/TESTING/MATGEN/slatms.c
@@ -833,7 +833,7 @@ static logical c_false = FALSE_;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatms_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void slatms_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, real *d__, integer *mode, real *cond, real *dmax__, 
 	integer *kl, integer *ku, char *pack, real *a, integer *lda, real *
 	work, integer *info)
@@ -853,22 +853,22 @@ static logical c_false = FALSE_;
     integer ipack, ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     integer idist, mnmin, iskew;
     real extra, dummy;
-    extern /* Subroutine */ int scopy_(integer *, real *, integer *, real *, 
+    extern /* Subroutine */ void scopy_(integer *, real *, integer *, real *, 
 	    integer *), slatm1_(integer *, real *, integer *, integer *, 
 	    integer *, real *, integer *, integer *);
     integer ic, jc, nc, il, iendch, ir, jr, ipackg, mr;
-    extern /* Subroutine */ int slagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void slagge_(integer *, integer *, integer *, 
 	    integer *, real *, real *, integer *, integer *, real *, integer *
 	    );
     integer minlda;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern real slarnd_(integer *, integer *);
     logical iltemp, givens;
     integer ioffst, irsign;
-    extern /* Subroutine */ int slartg_(real *, real *, real *, real *, real *
+    extern /* Subroutine */ void slartg_(real *, real *, real *, real *, real *
 	    ), slaset_(char *, integer *, integer *, real *, real *, real *, 
 	    integer *), slagsy_(integer *, integer *, real *, real *, 
 	    integer *, integer *, real *, integer *), slarot_(logical *, 
@@ -904,7 +904,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1039,8 +1039,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATMS", &i__1);
-	return 0;
+	xerbla_("SLATMS", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1061,7 +1061,7 @@ static logical c_false = FALSE_;
     slatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1090,7 +1090,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	sscal_(&mnmin, &alpha, &d__[1], &c__1);
@@ -1719,7 +1719,7 @@ static logical c_false = FALSE_;
 	}
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1881,7 +1881,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of SLATMS */
 
diff --git a/lapack-netlib/TESTING/MATGEN/slatmt.c b/lapack-netlib/TESTING/MATGEN/slatmt.c
index b3185363b..e34bb8c03 100644
--- a/lapack-netlib/TESTING/MATGEN/slatmt.c
+++ b/lapack-netlib/TESTING/MATGEN/slatmt.c
@@ -843,7 +843,7 @@ static logical c_false = FALSE_;
 /* > \ingroup real_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int slatmt_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void slatmt_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, real *d__, integer *mode, real *cond, real *dmax__, 
 	integer *rank, integer *kl, integer *ku, char *pack, real *a, integer 
 	*lda, real *work, integer *info)
@@ -863,22 +863,22 @@ static logical c_false = FALSE_;
     integer ipack, ioffg;
     extern logical lsame_(char *, char *);
     integer iinfo;
-    extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *);
+    extern /* Subroutine */ void sscal_(integer *, real *, real *, integer *);
     integer idist, mnmin, iskew;
     real extra, dummy;
-    extern /* Subroutine */ int scopy_(integer *, real *, integer *, real *, 
+    extern /* Subroutine */ void scopy_(integer *, real *, integer *, real *, 
 	    integer *), slatm7_(integer *, real *, integer *, integer *, 
 	    integer *, real *, integer *, integer *, integer *);
     integer ic, jc, nc, il, iendch, ir, jr, ipackg, mr;
-    extern /* Subroutine */ int slagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void slagge_(integer *, integer *, integer *, 
 	    integer *, real *, real *, integer *, integer *, real *, integer *
 	    );
     integer minlda;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern real slarnd_(integer *, integer *);
     integer ioffst, irsign;
     logical givens, iltemp;
-    extern /* Subroutine */ int slartg_(real *, real *, real *, real *, real *
+    extern /* Subroutine */ void slartg_(real *, real *, real *, real *, real *
 	    ), slaset_(char *, integer *, integer *, real *, real *, real *, 
 	    integer *), slagsy_(integer *, integer *, real *, real *, 
 	    integer *, integer *, real *, integer *), slarot_(logical *, 
@@ -914,7 +914,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1049,8 +1049,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("SLATMT", &i__1);
-	return 0;
+	xerbla_("SLATMT", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1072,7 +1072,7 @@ static logical c_false = FALSE_;
 	    iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1101,7 +1101,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	sscal_(rank, &alpha, &d__[1], &c__1);
@@ -1730,7 +1730,7 @@ static logical c_false = FALSE_;
 	}
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1892,7 +1892,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of SLATMT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlagge.c b/lapack-netlib/TESTING/MATGEN/zlagge.c
index 70651f850..bfa33cf34 100644
--- a/lapack-netlib/TESTING/MATGEN/zlagge.c
+++ b/lapack-netlib/TESTING/MATGEN/zlagge.c
@@ -628,7 +628,7 @@ static integer c__1 = 1;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlagge_(integer *m, integer *n, integer *kl, integer *ku,
+/* Subroutine */ void zlagge_(integer *m, integer *n, integer *kl, integer *ku,
 	 doublereal *d__, doublecomplex *a, integer *lda, integer *iseed, 
 	doublecomplex *work, integer *info)
 {
@@ -639,7 +639,7 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer i__, j;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *), zgemv_(char *, integer *, integer *, 
@@ -648,7 +648,8 @@ static integer c__1 = 1;
     extern doublereal dznrm2_(integer *, doublecomplex *, integer *);
     doublecomplex wa, wb;
     doublereal wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), zlacgv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void zlacgv_(
 	    integer *, doublecomplex *, integer *), zlarnv_(integer *, 
 	    integer *, integer *, doublecomplex *);
     doublecomplex tau;
@@ -688,8 +689,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("ZLAGGE", &i__1);
-	return 0;
+	xerbla_("ZLAGGE", &i__1, 6);
+	return;
     }
 
 /*     initialize A to diagonal matrix */
@@ -715,7 +716,7 @@ static integer c__1 = 1;
 /*     Quick exit if the user wants a diagonal matrix */
 
     if (*kl == 0 && *ku == 0) {
-	return 0;
+	return;
     }
 
 /*     pre- and post-multiply A by random unitary matrices */
@@ -1026,7 +1027,7 @@ static integer c__1 = 1;
 	}
 /* L70: */
     }
-    return 0;
+    return;
 
 /*     End of ZLAGGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlaghe.c b/lapack-netlib/TESTING/MATGEN/zlaghe.c
index b781dc994..06a8e31a3 100644
--- a/lapack-netlib/TESTING/MATGEN/zlaghe.c
+++ b/lapack-netlib/TESTING/MATGEN/zlaghe.c
@@ -616,7 +616,7 @@ static integer c__1 = 1;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlaghe_(integer *n, integer *k, doublereal *d__, 
+/* Subroutine */ void zlaghe_(integer *n, integer *k, doublereal *d__, 
 	doublecomplex *a, integer *lda, integer *iseed, doublecomplex *work, 
 	integer *info)
 {
@@ -626,18 +626,18 @@ static integer c__1 = 1;
     doublecomplex z__1, z__2, z__3, z__4;
 
     /* Local variables */
-    extern /* Subroutine */ int zher2_(char *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zher2_(char *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *);
     integer i__, j;
     doublecomplex alpha;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *);
     extern /* Double Complex */ VOID zdotc_(doublecomplex *, integer *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *);
-    extern /* Subroutine */ int zgemv_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zgemv_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, integer *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *, integer *), 
 	    zhemv_(char *, integer *, doublecomplex *, doublecomplex *, 
@@ -648,7 +648,8 @@ static integer c__1 = 1;
     extern doublereal dznrm2_(integer *, doublecomplex *, integer *);
     doublecomplex wa, wb;
     doublereal wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), zlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void zlarnv_(
 	    integer *, integer *, integer *, doublecomplex *);
     doublecomplex tau;
 
@@ -683,8 +684,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("ZLAGHE", &i__1);
-	return 0;
+	xerbla_("ZLAGHE", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -862,7 +863,7 @@ static integer c__1 = 1;
 	}
 /* L80: */
     }
-    return 0;
+    return;
 
 /*     End of ZLAGHE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlagsy.c b/lapack-netlib/TESTING/MATGEN/zlagsy.c
index 7f29344a9..431522ebe 100644
--- a/lapack-netlib/TESTING/MATGEN/zlagsy.c
+++ b/lapack-netlib/TESTING/MATGEN/zlagsy.c
@@ -616,7 +616,7 @@ static integer c__1 = 1;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlagsy_(integer *n, integer *k, doublereal *d__, 
+/* Subroutine */ void zlagsy_(integer *n, integer *k, doublereal *d__, 
 	doublecomplex *a, integer *lda, integer *iseed, doublecomplex *work, 
 	integer *info)
 {
@@ -629,13 +629,13 @@ static integer c__1 = 1;
     /* Local variables */
     integer i__, j;
     doublecomplex alpha;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *);
     extern /* Double Complex */ VOID zdotc_(doublecomplex *, integer *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *);
-    extern /* Subroutine */ int zgemv_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zgemv_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, integer *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *, integer *), 
 	    zaxpy_(integer *, doublecomplex *, doublecomplex *, integer *, 
@@ -646,7 +646,8 @@ static integer c__1 = 1;
     integer ii, jj;
     doublecomplex wa, wb;
     doublereal wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), zlacgv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void zlacgv_(
 	    integer *, doublecomplex *, integer *), zlarnv_(integer *, 
 	    integer *, integer *, doublecomplex *);
     doublecomplex tau;
@@ -682,8 +683,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("ZLAGSY", &i__1);
-	return 0;
+	xerbla_("ZLAGSY", &i__1, 6);
+	return;
     }
 
 /*     initialize lower triangle of A to diagonal matrix */
@@ -913,7 +914,7 @@ static integer c__1 = 1;
 	}
 /* L120: */
     }
-    return 0;
+    return;
 
 /*     End of ZLAGSY */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlahilb.c b/lapack-netlib/TESTING/MATGEN/zlahilb.c
index c08c75196..0e9cc1c4a 100644
--- a/lapack-netlib/TESTING/MATGEN/zlahilb.c
+++ b/lapack-netlib/TESTING/MATGEN/zlahilb.c
@@ -646,7 +646,7 @@ static doublecomplex c_b6 = {0.,0.};
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlahilb_(integer *n, integer *nrhs, doublecomplex *a, 
+/* Subroutine */ void zlahilb_(integer *n, integer *nrhs, doublecomplex *a, 
 	integer *lda, doublecomplex *x, integer *ldx, doublecomplex *b, 
 	integer *ldb, doublereal *work, integer *info, char *path)
 {
@@ -671,9 +671,9 @@ static doublecomplex c_b6 = {0.,0.};
     integer i__, j, m, r__;
     char c2[2];
     integer ti, tm;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern logical lsamen_(integer *, char *, char *);
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *);
     doublecomplex tmp;
 
@@ -723,8 +723,8 @@ static doublecomplex c_b6 = {0.,0.};
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("ZLAHILB", &i__1);
-	return 0;
+	xerbla_("ZLAHILB", &i__1, 7);
+	return;
     }
     if (*n > 6) {
 	*info = 1;
@@ -831,6 +831,6 @@ static doublecomplex c_b6 = {0.,0.};
 	    }
 	}
     }
-    return 0;
+    return;
 } /* zlahilb_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlakf2.c b/lapack-netlib/TESTING/MATGEN/zlakf2.c
index ea2390753..f7a90eef7 100644
--- a/lapack-netlib/TESTING/MATGEN/zlakf2.c
+++ b/lapack-netlib/TESTING/MATGEN/zlakf2.c
@@ -616,7 +616,7 @@ static doublecomplex c_b1 = {0.,0.};
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlakf2_(integer *m, integer *n, doublecomplex *a, 
+/* Subroutine */ void zlakf2_(integer *m, integer *n, doublecomplex *a, 
 	integer *lda, doublecomplex *b, doublecomplex *d__, doublecomplex *e, 
 	doublecomplex *z__, integer *ldz)
 {
@@ -627,7 +627,7 @@ static doublecomplex c_b1 = {0.,0.};
 
     /* Local variables */
     integer i__, j, l, ik, jk, mn;
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *);
     integer mn2;
 
@@ -739,7 +739,7 @@ static doublecomplex c_b1 = {0.,0.};
 /* L90: */
     }
 
-    return 0;
+    return;
 
 /*     End of ZLAKF2 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlarge.c b/lapack-netlib/TESTING/MATGEN/zlarge.c
index 9084df6ac..4b06a818d 100644
--- a/lapack-netlib/TESTING/MATGEN/zlarge.c
+++ b/lapack-netlib/TESTING/MATGEN/zlarge.c
@@ -601,7 +601,7 @@ static integer c__1 = 1;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlarge_(integer *n, doublecomplex *a, integer *lda, 
+/* Subroutine */ void zlarge_(integer *n, doublecomplex *a, integer *lda, 
 	integer *iseed, doublecomplex *work, integer *info)
 {
     /* System generated locals */
@@ -611,7 +611,7 @@ static integer c__1 = 1;
 
     /* Local variables */
     integer i__;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *), zgemv_(char *, integer *, integer *, 
@@ -620,7 +620,8 @@ static integer c__1 = 1;
     extern doublereal dznrm2_(integer *, doublecomplex *, integer *);
     doublecomplex wa, wb;
     doublereal wn;
-    extern /* Subroutine */ int xerbla_(char *, integer *), zlarnv_(
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
+    extern void zlarnv_(
 	    integer *, integer *, integer *, doublecomplex *);
     doublecomplex tau;
 
@@ -652,8 +653,8 @@ static integer c__1 = 1;
     }
     if (*info < 0) {
 	i__1 = -(*info);
-	xerbla_("ZLARGE", &i__1);
-	return 0;
+	xerbla_("ZLARGE", &i__1, 6);
+	return;
     }
 
 /*     pre- and post-multiply A by random unitary matrix */
@@ -704,7 +705,7 @@ static integer c__1 = 1;
 		* a_dim1 + 1], lda);
 /* L10: */
     }
-    return 0;
+    return;
 
 /*     End of ZLARGE */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlaror.c b/lapack-netlib/TESTING/MATGEN/zlaror.c
index eeb9f4fcf..6ada57b8a 100644
--- a/lapack-netlib/TESTING/MATGEN/zlaror.c
+++ b/lapack-netlib/TESTING/MATGEN/zlaror.c
@@ -673,7 +673,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlaror_(char *side, char *init, integer *m, integer *n, 
+/* Subroutine */ void zlaror_(char *side, char *init, integer *m, integer *n, 
 	doublecomplex *a, integer *lda, integer *iseed, doublecomplex *x, 
 	integer *info)
 {
@@ -687,25 +687,25 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     integer irow, j;
     extern logical lsame_(char *, char *);
     doublecomplex csign;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *);
     integer ixfrm;
-    extern /* Subroutine */ int zgemv_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zgemv_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, integer *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *, integer *);
     integer itype, nxfrm;
     doublereal xnorm;
     extern doublereal dznrm2_(integer *, doublecomplex *, integer *);
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     doublereal factor;
-    extern /* Subroutine */ int zlacgv_(integer *, doublecomplex *, integer *)
+    extern /* Subroutine */ void zlacgv_(integer *, doublecomplex *, integer *)
 	    ;
     //extern /* Double Complex */ VOID zlarnd_(doublecomplex *, integer *, 
     extern doublecomplex zlarnd_(integer *, 
 	    integer *);
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *);
     doublecomplex xnorms;
 
@@ -729,7 +729,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     /* Function Body */
     *info = 0;
     if (*n == 0 || *m == 0) {
-	return 0;
+	return;
     }
 
     itype = 0;
@@ -756,8 +756,8 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
     }
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLAROR", &i__1);
-	return 0;
+	xerbla_("ZLAROR", &i__1, 6);
+	return;
     }
 
     if (itype == 1) {
@@ -821,8 +821,8 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 	if (abs(factor) < 1e-20) {
 	    *info = 1;
 	    i__2 = -(*info);
-	    xerbla_("ZLAROR", &i__2);
-	    return 0;
+	    xerbla_("ZLAROR", &i__2, 6);
+	    return;
 	} else {
 	    factor = 1. / factor;
 	}
@@ -905,7 +905,7 @@ t by U and the right by UC>       SIDE = 'T'   Multiply A on the left by U and t
 /* L60: */
 	}
     }
-    return 0;
+    return;
 
 /*     End of ZLAROR */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlarot.c b/lapack-netlib/TESTING/MATGEN/zlarot.c
index da7207df2..6b2e88b32 100644
--- a/lapack-netlib/TESTING/MATGEN/zlarot.c
+++ b/lapack-netlib/TESTING/MATGEN/zlarot.c
@@ -740,7 +740,7 @@ static integer c__8 = 8;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlarot_(logical *lrows, logical *lleft, logical *lright, 
+/* Subroutine */ void zlarot_(logical *lrows, logical *lleft, logical *lright, 
 	integer *nl, doublecomplex *c__, doublecomplex *s, doublecomplex *a, 
 	integer *lda, doublecomplex *xleft, doublecomplex *xright)
 {
@@ -753,7 +753,7 @@ static integer c__8 = 8;
     doublecomplex tempx;
     integer ix, iy, nt;
     doublecomplex xt[2], yt[2];
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     integer iyt;
 
 
@@ -805,12 +805,12 @@ static integer c__8 = 8;
 /*     Check for errors */
 
     if (*nl < nt) {
-	xerbla_("ZLAROT", &c__4);
-	return 0;
+	xerbla_("ZLAROT", &c__4, 6);
+	return;
     }
     if (*lda <= 0 || ! (*lrows) && *lda < *nl - nt) {
-	xerbla_("ZLAROT", &c__8);
-	return 0;
+	xerbla_("ZLAROT", &c__8, 6);
+	return;
     }
 
 /*     Rotate */
@@ -888,7 +888,7 @@ static integer c__8 = 8;
 	a[i__1].r = yt[i__2].r, a[i__1].i = yt[i__2].i;
     }
 
-    return 0;
+    return;
 
 /*     End of ZLAROT */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatm1.c b/lapack-netlib/TESTING/MATGEN/zlatm1.c
index 4b84fcc59..c75787730 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatm1.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatm1.c
@@ -648,7 +648,7 @@ static integer c__3 = 3;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatm1_(integer *mode, doublereal *cond, integer *irsign,
+/* Subroutine */ void zlatm1_(integer *mode, doublereal *cond, integer *irsign,
 	 integer *idist, integer *iseed, doublecomplex *d__, integer *n, 
 	integer *info)
 {
@@ -663,11 +663,11 @@ static integer c__3 = 3;
     doublereal alpha;
     doublecomplex ctemp;
     extern doublereal dlaran_(integer *);
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     //extern /* Double Complex */ VOID zlarnd_(doublecomplex *, integer *, 
     extern doublecomplex zlarnd_(integer *, 
 	    integer *);
-    extern /* Subroutine */ int zlarnv_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void zlarnv_(integer *, integer *, integer *, 
 	    doublecomplex *);
 
 
@@ -692,7 +692,7 @@ static integer c__3 = 3;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Set INFO if an error */
@@ -712,8 +712,8 @@ static integer c__3 = 3;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLATM1", &i__1);
-	return 0;
+	xerbla_("ZLATM1", &i__1, 6);
+	return;
     }
 
 /*     Compute D according to COND and MODE */
@@ -848,7 +848,7 @@ L120:
 
     }
 
-    return 0;
+    return;
 
 /*     End of ZLATM1 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatm5.c b/lapack-netlib/TESTING/MATGEN/zlatm5.c
index 047a633c4..753ee0ce6 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatm5.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatm5.c
@@ -780,7 +780,7 @@ static doublecomplex c_b5 = {20.,0.};
 /* > \endverbatim */
 /* > */
 /*  ===================================================================== */
-/* Subroutine */ int zlatm5_(integer *prtype, integer *m, integer *n, 
+/* Subroutine */ void zlatm5_(integer *prtype, integer *m, integer *n, 
 	doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, 
 	doublecomplex *c__, integer *ldc, doublecomplex *d__, integer *ldd, 
 	doublecomplex *e, integer *lde, doublecomplex *f, integer *ldf, 
@@ -797,7 +797,7 @@ static doublecomplex c_b5 = {20.,0.};
     /* Local variables */
     integer i__, j, k;
     doublecomplex imeps, reeps;
-    extern /* Subroutine */ int zgemm_(char *, char *, integer *, integer *, 
+    extern /* Subroutine */ void zgemm_(char *, char *, integer *, integer *, 
 	    integer *, doublecomplex *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *, doublecomplex *, doublecomplex *, 
 	    integer *);
@@ -1281,6 +1281,6 @@ static doublecomplex c_b5 = {20.,0.};
 
 /*     End of ZLATM5 */
 
-    return 0;
+    return;
 } /* zlatm5_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatm6.c b/lapack-netlib/TESTING/MATGEN/zlatm6.c
index 6e25e5023..e20ec7092 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatm6.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatm6.c
@@ -687,7 +687,7 @@ static integer c__24 = 24;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatm6_(integer *type__, integer *n, doublecomplex *a, 
+/* Subroutine */ void zlatm6_(integer *type__, integer *n, doublecomplex *a, 
 	integer *lda, doublecomplex *b, doublecomplex *x, integer *ldx, 
 	doublecomplex *y, integer *ldy, doublecomplex *alpha, doublecomplex *
 	beta, doublecomplex *wx, doublecomplex *wy, doublereal *s, doublereal 
@@ -705,7 +705,7 @@ static integer c__24 = 24;
     integer i__, j;
     doublecomplex z__[64]	/* was [8][8] */;
     doublereal rwork[50];
-    extern /* Subroutine */ int zlakf2_(integer *, integer *, doublecomplex *,
+    extern /* Subroutine */ void zlakf2_(integer *, integer *, doublecomplex *,
 	     integer *, doublecomplex *, doublecomplex *, doublecomplex *, 
 	    doublecomplex *, integer *), zgesvd_(char *, char *, integer *, 
 	    integer *, doublecomplex *, integer *, doublereal *, 
@@ -934,7 +934,7 @@ static integer c__24 = 24;
 	    &c__1, &work[2], &c__24, &rwork[8], &info);
     dif[5] = rwork[7];
 
-    return 0;
+    return;
 
 /*     End of ZLATM6 */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatme.c b/lapack-netlib/TESTING/MATGEN/zlatme.c
index 3e4c1f3a2..ded6ca214 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatme.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatme.c
@@ -812,7 +812,7 @@ static integer c__5 = 5;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatme_(integer *n, char *dist, integer *iseed, 
+/* Subroutine */ void zlatme_(integer *n, char *dist, integer *iseed, 
 	doublecomplex *d__, integer *mode, doublereal *cond, doublecomplex *
 	dmax__, char *rsign, char *upper, char *sim, doublereal *ds, integer *
 	modes, doublereal *conds, integer *kl, integer *ku, doublereal *anorm,
@@ -833,26 +833,26 @@ static integer c__5 = 5;
     integer iinfo;
     doublereal tempa[1];
     integer icols;
-    extern /* Subroutine */ int zgerc_(integer *, integer *, doublecomplex *, 
+    extern /* Subroutine */ void zgerc_(integer *, integer *, doublecomplex *, 
 	    doublecomplex *, integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *);
     integer idist;
-    extern /* Subroutine */ int zscal_(integer *, doublecomplex *, 
+    extern /* Subroutine */ void zscal_(integer *, doublecomplex *, 
 	    doublecomplex *, integer *), zgemv_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, integer *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *, integer *);
     integer irows;
-    extern /* Subroutine */ int zcopy_(integer *, doublecomplex *, integer *, 
+    extern /* Subroutine */ void zcopy_(integer *, doublecomplex *, integer *, 
 	    doublecomplex *, integer *), dlatm1_(integer *, doublereal *, 
 	    integer *, integer *, integer *, doublereal *, integer *, integer 
 	    *), zlatm1_(integer *, doublereal *, integer *, integer *, 
 	    integer *, doublecomplex *, integer *, integer *);
     integer ic, jc, ir;
     doublereal ralpha;
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     extern doublereal zlange_(char *, integer *, integer *, doublecomplex *, 
 	    integer *, doublereal *);
-    extern /* Subroutine */ int zdscal_(integer *, doublereal *, 
+    extern /* Subroutine */ void zdscal_(integer *, doublereal *, 
 	    doublecomplex *, integer *), zlarge_(integer *, doublecomplex *, 
 	    integer *, integer *, doublecomplex *, integer *), zlarfg_(
 	    integer *, doublecomplex *, doublecomplex *, integer *, 
@@ -861,10 +861,10 @@ static integer c__5 = 5;
     extern doublecomplex zlarnd_(integer *, 
 	    integer *);
     integer irsign;
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *);
     integer iupper;
-    extern /* Subroutine */ int zlarnv_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void zlarnv_(integer *, integer *, integer *, 
 	    doublecomplex *);
     doublecomplex xnorms;
     integer jcr;
@@ -898,7 +898,7 @@ static integer c__5 = 5;
 /*     Quick return if possible */
 
     if (*n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -990,8 +990,8 @@ static integer c__5 = 5;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLATME", &i__1);
-	return 0;
+	xerbla_("ZLATME", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1012,7 +1012,7 @@ static integer c__5 = 5;
     zlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], n, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && abs(*mode) != 6) {
 
@@ -1032,7 +1032,7 @@ static integer c__5 = 5;
 	    alpha.r = z__1.r, alpha.i = z__1.i;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	zscal_(n, &alpha, &d__[1], &c__1);
@@ -1069,7 +1069,7 @@ static integer c__5 = 5;
 	dlatm1_(modes, conds, &c__0, &c__0, &iseed[1], &ds[1], n, &iinfo);
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by V and V' */
@@ -1077,7 +1077,7 @@ static integer c__5 = 5;
 	zlarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
 
 /*        Multiply by S and (1/S) */
@@ -1090,7 +1090,7 @@ static integer c__5 = 5;
 		zdscal_(n, &d__1, &a[j * a_dim1 + 1], &c__1);
 	    } else {
 		*info = 5;
-		return 0;
+		return;
 	    }
 /* L50: */
 	}
@@ -1100,7 +1100,7 @@ static integer c__5 = 5;
 	zlarge_(n, &a[a_offset], lda, &iseed[1], &work[1], &iinfo);
 	if (iinfo != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1214,7 +1214,7 @@ static integer c__5 = 5;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of ZLATME */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatmr.c b/lapack-netlib/TESTING/MATGEN/zlatmr.c
index c0c2d0b3b..c42a25f2b 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatmr.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatmr.c
@@ -999,7 +999,7 @@ static integer c__1 = 1;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatmr_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void zlatmr_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublecomplex *d__, integer *mode, doublereal *cond,
 	 doublecomplex *dmax__, char *rsign, char *grade, doublecomplex *dl, 
 	integer *model, doublereal *condl, doublecomplex *dr, integer *moder, 
@@ -1024,7 +1024,7 @@ static integer c__1 = 1;
     integer mnsub;
     doublereal onorm;
     integer mxsub, npvts;
-    extern /* Subroutine */ int zlatm1_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void zlatm1_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublecomplex *, integer *, integer *);
     extern /* Double Complex */ VOID zlatm2_(doublecomplex *, integer *, 
 	    integer *, integer *, integer *, integer *, integer *, integer *, 
@@ -1039,11 +1039,11 @@ static integer c__1 = 1;
     logical fulbnd;
     extern doublereal zlangb_(char *, integer *, integer *, integer *, 
 	    doublecomplex *, integer *, doublereal *);
-    extern /* Subroutine */ int xerbla_(char *, integer *);
+    extern /* Subroutine */ int xerbla_(char *, integer *, ftnlen);
     logical badpvt;
     extern doublereal zlange_(char *, integer *, integer *, doublecomplex *, 
 	    integer *, doublereal *);
-    extern /* Subroutine */ int zdscal_(integer *, doublereal *, 
+    extern /* Subroutine */ void zdscal_(integer *, doublereal *, 
 	    doublecomplex *, integer *);
     extern doublereal zlansb_(char *, char *, integer *, integer *, 
 	    doublecomplex *, integer *, doublereal *);
@@ -1083,7 +1083,7 @@ static integer c__1 = 1;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1284,8 +1284,8 @@ static integer c__1 = 1;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLATMR", &i__1);
-	return 0;
+	xerbla_("ZLATMR", &i__1, 6);
+	return;
     }
 
 /*     Decide if we can pivot consistently */
@@ -1311,7 +1311,7 @@ static integer c__1 = 1;
     zlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, info);
     if (*info != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
     if (*mode != 0 && *mode != -6 && *mode != 6) {
 
@@ -1327,7 +1327,7 @@ static integer c__1 = 1;
 	}
 	if (temp == 0. && (dmax__->r != 0. || dmax__->i != 0.)) {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 	if (temp != 0.) {
 	    z__1.r = dmax__->r / temp, z__1.i = dmax__->i / temp;
@@ -1367,7 +1367,7 @@ static integer c__1 = 1;
 	zlatm1_(model, condl, &c__0, &idist, &iseed[1], &dl[1], m, info);
 	if (*info != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -1377,7 +1377,7 @@ static integer c__1 = 1;
 	zlatm1_(moder, condr, &c__0, &idist, &iseed[1], &dr[1], n, info);
 	if (*info != 0) {
 	    *info = 4;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2032,7 +2032,7 @@ static integer c__1 = 1;
 /*           Desired scaling impossible */
 
 	    *info = 5;
-	    return 0;
+	    return;
 
 	} else if (*anorm > 1. && onorm < 1. || *anorm < 1. && onorm > 1.) {
 
@@ -2104,6 +2104,6 @@ static integer c__1 = 1;
 
 /*     End of ZLATMR */
 
-    return 0;
+    return;
 } /* zlatmr_ */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatms.c b/lapack-netlib/TESTING/MATGEN/zlatms.c
index ca0191a13..f0bb66a44 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatms.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatms.c
@@ -846,7 +846,7 @@ static logical c_false = FALSE_;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatms_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void zlatms_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublereal *d__, integer *mode, doublereal *cond, 
 	doublereal *dmax__, integer *kl, integer *ku, char *pack, 
 	doublecomplex *a, integer *lda, doublecomplex *work, integer *info)
@@ -869,40 +869,41 @@ static logical c_false = FALSE_;
     integer ipack;
     doublereal realc;
     integer ioffg;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     extern logical lsame_(char *, char *);
     integer iinfo;
     doublecomplex ctemp;
     integer idist, mnmin, iskew;
     doublecomplex extra, dummy;
-    extern /* Subroutine */ int dlatm1_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm1_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *);
     integer ic, jc, nc, il;
     doublecomplex ct;
     integer iendch, ir, jr, ipackg, mr, minlda;
     extern doublereal dlarnd_(integer *, integer *);
     doublecomplex st;
-    extern /* Subroutine */ int zlagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void zlagge_(integer *, integer *, integer *, 
 	    integer *, doublereal *, doublecomplex *, integer *, integer *, 
 	    doublecomplex *, integer *), zlaghe_(integer *, integer *, 
 	    doublereal *, doublecomplex *, integer *, integer *, 
-	    doublecomplex *, integer *), xerbla_(char *, integer *);
+	    doublecomplex *, integer *);
+    extern int xerbla_(char *, integer *, ftnlen);
     logical iltemp, givens;
     integer ioffst, irsign;
     //extern /* Double Complex */ VOID zlarnd_(doublecomplex *, integer *, 
     extern doublecomplex zlarnd_(integer *, 
 	    integer *);
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *), zlartg_(doublecomplex *, doublecomplex *, doublereal *, 
 	    doublecomplex *, doublecomplex *);
     logical ilextr;
-    extern /* Subroutine */ int zlagsy_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void zlagsy_(integer *, integer *, doublereal *, 
 	    doublecomplex *, integer *, integer *, doublecomplex *, integer *)
 	    ;
     logical topdwn;
     integer ir1, ir2, isympk;
-    extern /* Subroutine */ int zlarot_(logical *, logical *, logical *, 
+    extern /* Subroutine */ void zlarot_(logical *, logical *, logical *, 
 	    integer *, doublecomplex *, doublecomplex *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *);
     integer jch, llb, jkl, jku, uub;
@@ -934,7 +935,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1073,8 +1074,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLATMS", &i__1);
-	return 0;
+	xerbla_("ZLATMS", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1095,7 +1096,7 @@ static logical c_false = FALSE_;
     dlatm1_(mode, cond, &irsign, &idist, &iseed[1], &d__[1], &mnmin, &iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1124,7 +1125,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	dscal_(&mnmin, &alpha, &d__[1], &c__1);
@@ -2038,7 +2039,7 @@ static logical c_false = FALSE_;
 
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2213,7 +2214,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of ZLATMS */
 
diff --git a/lapack-netlib/TESTING/MATGEN/zlatmt.c b/lapack-netlib/TESTING/MATGEN/zlatmt.c
index 47279858b..551b17628 100644
--- a/lapack-netlib/TESTING/MATGEN/zlatmt.c
+++ b/lapack-netlib/TESTING/MATGEN/zlatmt.c
@@ -854,7 +854,7 @@ static logical c_false = FALSE_;
 /* > \ingroup complex16_matgen */
 
 /*  ===================================================================== */
-/* Subroutine */ int zlatmt_(integer *m, integer *n, char *dist, integer *
+/* Subroutine */ void zlatmt_(integer *m, integer *n, char *dist, integer *
 	iseed, char *sym, doublereal *d__, integer *mode, doublereal *cond, 
 	doublereal *dmax__, integer *rank, integer *kl, integer *ku, char *
 	pack, doublecomplex *a, integer *lda, doublecomplex *work, integer *
@@ -876,14 +876,14 @@ static logical c_false = FALSE_;
     doublecomplex s;
     doublereal alpha, angle, realc;
     integer ipack, ioffg;
-    extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, 
+    extern /* Subroutine */ void dscal_(integer *, doublereal *, doublereal *, 
 	    integer *);
     extern logical lsame_(char *, char *);
     integer iinfo, idist, mnmin;
     doublecomplex extra;
     integer iskew;
     doublecomplex dummy, ztemp;
-    extern /* Subroutine */ int dlatm7_(integer *, doublereal *, integer *, 
+    extern /* Subroutine */ void dlatm7_(integer *, doublereal *, integer *, 
 	    integer *, integer *, doublereal *, integer *, integer *, integer 
 	    *);
     integer ic, jc, nc, il;
@@ -891,26 +891,27 @@ static logical c_false = FALSE_;
     integer iendch, ir, jr, ipackg, mr, minlda;
     extern doublereal dlarnd_(integer *, integer *);
     doublecomplex st;
-    extern /* Subroutine */ int zlagge_(integer *, integer *, integer *, 
+    extern /* Subroutine */ void zlagge_(integer *, integer *, integer *, 
 	    integer *, doublereal *, doublecomplex *, integer *, integer *, 
 	    doublecomplex *, integer *), zlaghe_(integer *, integer *, 
 	    doublereal *, doublecomplex *, integer *, integer *, 
-	    doublecomplex *, integer *), xerbla_(char *, integer *);
+	    doublecomplex *, integer *);
+    extern int xerbla_(char *, integer *, ftnlen);
     integer ioffst, irsign;
     logical givens, iltemp;
     //extern /* Double Complex */ VOID zlarnd_(doublecomplex *, integer *, 
     extern doublecomplex zlarnd_(integer *, 
 	    integer *);
-    extern /* Subroutine */ int zlaset_(char *, integer *, integer *, 
+    extern /* Subroutine */ void zlaset_(char *, integer *, integer *, 
 	    doublecomplex *, doublecomplex *, doublecomplex *, integer *), zlartg_(doublecomplex *, doublecomplex *, doublereal *, 
 	    doublecomplex *, doublecomplex *);
     logical ilextr;
-    extern /* Subroutine */ int zlagsy_(integer *, integer *, doublereal *, 
+    extern /* Subroutine */ void zlagsy_(integer *, integer *, doublereal *, 
 	    doublecomplex *, integer *, integer *, doublecomplex *, integer *)
 	    ;
     integer ir1, ir2, isympk;
     logical topdwn;
-    extern /* Subroutine */ int zlarot_(logical *, logical *, logical *, 
+    extern /* Subroutine */ void zlarot_(logical *, logical *, logical *, 
 	    integer *, doublecomplex *, doublecomplex *, doublecomplex *, 
 	    integer *, doublecomplex *, doublecomplex *);
     integer jch, llb, jkl, jku, uub;
@@ -942,7 +943,7 @@ static logical c_false = FALSE_;
 /*     Quick return if possible */
 
     if (*m == 0 || *n == 0) {
-	return 0;
+	return;
     }
 
 /*     Decode DIST */
@@ -1081,8 +1082,8 @@ static logical c_false = FALSE_;
 
     if (*info != 0) {
 	i__1 = -(*info);
-	xerbla_("ZLATMT", &i__1);
-	return 0;
+	xerbla_("ZLATMT", &i__1, 6);
+	return;
     }
 
 /*     Initialize random number generator */
@@ -1104,7 +1105,7 @@ static logical c_false = FALSE_;
 	    iinfo);
     if (iinfo != 0) {
 	*info = 1;
-	return 0;
+	return;
     }
 
 /*     Choose Top-Down if D is (apparently) increasing, */
@@ -1133,7 +1134,7 @@ static logical c_false = FALSE_;
 	    alpha = *dmax__ / temp;
 	} else {
 	    *info = 2;
-	    return 0;
+	    return;
 	}
 
 	dscal_(rank, &alpha, &d__[1], &c__1);
@@ -2046,7 +2047,7 @@ static logical c_false = FALSE_;
 
 	if (iinfo != 0) {
 	    *info = 3;
-	    return 0;
+	    return;
 	}
     }
 
@@ -2221,7 +2222,7 @@ static logical c_false = FALSE_;
 	}
     }
 
-    return 0;
+    return;
 
 /*     End of ZLATMT */