also return after initializing x and y

Author: thijssteel

SRC/cggglm.f

@@ -278,6 +278,7 @@ SUBROUTINE CGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK,
          DO I = 1, P
             Y(I) = CZERO
          END DO
+         RETURN
       END IF
 *
 *     Compute the GQR factorization of matrices A and B:

SRC/dggglm.f

@@ -277,6 +277,7 @@ SUBROUTINE DGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK,
          DO I = 1, P
             Y(I) = ZERO
          END DO
+         RETURN
       END IF
 *
 *     Compute the GQR factorization of matrices A and B:

SRC/sggglm.f

@@ -277,6 +277,7 @@ SUBROUTINE SGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK,
          DO I = 1, P
             Y(I) = ZERO
          END DO
+         RETURN
       END IF
 *
 *     Compute the GQR factorization of matrices A and B:

SRC/zggglm.f

@@ -278,6 +278,7 @@ SUBROUTINE ZGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK,
          DO I = 1, P
             Y(I) = CZERO
          END DO
+         RETURN
       END IF
 *
 *     Compute the GQR factorization of matrices A and B: