cleaning f2c mess / trivial stuff

1 files changed, 17 insertions, 26 deletions

diff --git a/unsupported/Eigen/src/NonLinear/qform.h b/unsupported/Eigen/src/NonLinear/qform.h
index a9920e7b3..844c8c2f6 100644
--- a/unsupported/Eigen/src/NonLinear/qform.h
+++ b/unsupported/Eigen/src/NonLinear/qform.h
@@ -4,10 +4,10 @@ void ei_qform(int m, int n, Scalar *q, int
 	ldq, Scalar *wa)
 {
     /* System generated locals */
-    int q_dim1, q_offset, i__1, i__2, i__3;
+    int q_dim1, q_offset;
 
     /* Local variables */
-    int i__, j, k, l, jm1, np1;
+    int i, j, k, l, jm1, np1;
     Scalar sum, temp;
     int minmn;
 
@@ -25,12 +25,10 @@ void ei_qform(int m, int n, Scalar *q, int
     if (minmn < 2) {
 	goto L30;
     }
-    i__1 = minmn;
-    for (j = 2; j <= i__1; ++j) {
+    for (j = 2; j <= minmn; ++j) {
 	jm1 = j - 1;
-	i__2 = jm1;
-	for (i__ = 1; i__ <= i__2; ++i__) {
-	    q[i__ + j * q_dim1] = 0.;
+	for (i = 1; i <= jm1; ++i) {
+	    q[i + j * q_dim1] = 0.;
 /* L10: */
 	}
 /* L20: */
@@ -43,11 +41,9 @@ L30:
     if (m < np1) {
 	goto L60;
     }
-    i__1 = m;
-    for (j = np1; j <= i__1; ++j) {
-	i__2 = m;
-	for (i__ = 1; i__ <= i__2; ++i__) {
-	    q[i__ + j * q_dim1] = 0.;
+    for (j = np1; j <= m; ++j) {
+	for (i = 1; i <= m; ++i) {
+	    q[i + j * q_dim1] = 0.;
 /* L40: */
 	}
 	q[j + j * q_dim1] = 1.;
@@ -57,31 +53,26 @@ L60:
 
 /*     accumulate q from its factored form. */
 
-    i__1 = minmn;
-    for (l = 1; l <= i__1; ++l) {
+    for (l = 1; l <= minmn; ++l) {
 	k = minmn - l + 1;
-	i__2 = m;
-	for (i__ = k; i__ <= i__2; ++i__) {
-	    wa[i__] = q[i__ + k * q_dim1];
-	    q[i__ + k * q_dim1] = 0.;
+	for (i = k; i <= m; ++i) {
+	    wa[i] = q[i + k * q_dim1];
+	    q[i + k * q_dim1] = 0.;
 /* L70: */
 	}
 	q[k + k * q_dim1] = 1.;
 	if (wa[k] == 0.) {
 	    goto L110;
 	}
-	i__2 = m;
-	for (j = k; j <= i__2; ++j) {
+	for (j = k; j <= m; ++j) {
 	    sum = 0.;
-	    i__3 = m;
-	    for (i__ = k; i__ <= i__3; ++i__) {
-		sum += q[i__ + j * q_dim1] * wa[i__];
+	    for (i = k; i <= m; ++i) {
+		sum += q[i + j * q_dim1] * wa[i];
 /* L80: */
 	    }
 	    temp = sum / wa[k];
-	    i__3 = m;
-	    for (i__ = k; i__ <= i__3; ++i__) {
-		q[i__ + j * q_dim1] -= temp * wa[i__];
+	    for (i = k; i <= m; ++i) {
+		q[i + j * q_dim1] -= temp * wa[i];
 /* L90: */
 	    }
 /* L100: */