fehlerstelle gefunden: DGEMM in glsfit. details noch unklar.

src/glsfit.f

@@ -74,6 +74,7 @@ c     debug options
 
       dbglvl=0
 
+
 c     check for existence of solution:
       IF (N.LE.NTREND) THEN
          IERR=1
@@ -108,6 +109,8 @@ c     determine CMAT**1/2
             CWORK(J,I)=CMAT(I,J)
  10      CONTINUE
  20   CONTINUE
+
+
 c     eigenvalue decomposition CMAT = V * D * V':
       call DSYEV( 'V', 'U', N, CWORK, LDCWRK, CWRK2(1,1), WORK, LWORK, 
      $            IERR )
@@ -194,9 +197,15 @@ c     calculate covariance for estimator BETA
 c     intermediate step: COVBTA=FMAT^T * CMAT^-1 * FMAT
       CALL DGEMM('T','N',NTREND,N,N,ONE,
      .            FMAT,LDF,CMINV,LDCINV,ZERO,CWORK,LDCWRK)
+ccc
+c      goto 123
 
+ccccccc FEHLER ccccccccc:
       CALL DGEMM('N','N',NTREND,NTREND,N,ONE,
      .            CWORK,LDCWRK,FMAT,LDF,ZERO,COVBTA,LDCVBT)
+
+ccc
+      goto 123
       DO 130 I=1,NTREND
          DO 120 J=1,NTREND
             IF( I.EQ.J) THEN
@@ -343,5 +352,8 @@ c     final step
             COVBTA(I,J)=COVBTA(I,J)*SGSQR
  180     CONTINUE
  190  CONTINUE
+
+ 123  continue
+
       RETURN
       END