anchor.diff:            DO 100 J = 1,N
             DO 100 I = 1,M
               HTMPS = ZERO 
  -            HTMPI = ZERO 
  -            HTMPJ = ZERO 
  -            HTMPK = ZERO 
  +            HTMPJ = ZERO 
               DO 110 L = 1,K
   C             Hamilton Product
  -              HTMPS = HTMPS + A(1,L,I)*B(1,L,J) - A(2,L,I)*B(2,L,J) -
  -     $                        A(3,L,I)*B(3,L,J) - A(4,L,I)*B(4,L,J)
  -
  -              HTMPI = HTMPI + A(1,L,I)*B(2,L,J) + A(2,L,I)*B(1,L,J) +
  -     $                        A(3,L,I)*B(4,L,J) - A(4,L,I)*B(3,L,J)
  -
  -              HTMPJ = HTMPJ + A(1,L,I)*B(3,L,J) - A(2,L,I)*B(4,L,J) +
  -     $                        A(3,L,I)*B(1,L,J) + A(4,L,I)*B(2,L,J)
  -
  -              HTMPK = HTMPK + A(1,L,I)*B(4,L,J) + A(2,L,I)*B(3,L,J) -
  -     $                        A(3,L,I)*B(2,L,J) + A(4,L,I)*B(1,L,J)
  +              HTMPS = HTMPS + 
  +     $          A(1,L,I)*B(1,L,J) - A(2,L,I)*CONJG(B(2,L,J))
  +              HTMPJ = HTMPJ + 
  +     $          A(1,L,I)*B(2,L,J) + A(2,L,I)*CONJG(B(1,L,J))
     110       CONTINUE
   
               IF( BZERO ) THEN
                 C(1,I,J) = ALPHA * HTMPS
  -              C(2,I,J) = ALPHA * HTMPI
  -              C(3,I,J) = ALPHA * HTMPJ
  -              C(4,I,J) = ALPHA * HTMPK
  -            ELSE
  -              C(1,I,J) = ALPHA * HTMPS + BETA * C(1,I,J)
  -              C(2,I,J) = ALPHA * HTMPI + BETA * C(2,I,J)
  -              C(3,I,J) = ALPHA * HTMPJ + BETA * C(3,I,J)
  -              C(4,I,J) = ALPHA * HTMPK + BETA * C(4,I,J)
  +              C(2,I,J) = ALPHA * HTMPJ
  +            ELSE
  +              C(1,I,J) = ALPHA * HTMPS + BETA * C(1,I,J)
  +              C(2,I,J) = ALPHA * HTMPJ + BETA * C(2,I,J)
               ENDIF
     100     CONTINUE
           ENDIF
  
positive.code:      $          BZERO, BONE
        INTEGER*4 M,N,K, LDA,LDB,LDC, I,J,L
        REAL*8    ALPHA, BETA
        REAL*8    A(4,LDA,*), B(4,LDB,*), C(4,LDC,*)
  C
        REAL*8    HTMPS, HTMPI, HTMPJ, HTMPK
  C
        REAL*8    ONE, ZERO
        PARAMETER (ONE=1.0D+0, ZERO=0.0D+0)
  C
  
        NOTA  = LSAME(TRANSA,'N')
  
negative.code:           DO 240 J = 1,N
            DO 240 I = 1,M
              HTMPS = ZERO 
              HTMPJ = ZERO 
              DO 250 L = 1,K
  C             Hamilton Product (IMPLIED CONJ(A))
                HTMPS = HTMPS + 
       $          CONJG(A(1,L,I))*B(1,J,L) + A(2,L,I)*CONJG(B(2,J,L))
                HTMPJ = HTMPJ + 
       $          CONJG(A(1,L,I))*B(2,J,L) - A(2,L,I)*CONJG(B(1,J,L))
    250       CONTINUE
  C
              IF( BZERO ) THEN
                C(1,I,J) = ALPHA * HTMPS
                C(2,I,J) = ALPHA * HTMPJ
              ELSE
                C(1,I,J) = ALPHA * HTMPS + BETA * C(1,I,J)
                C(2,I,J) = ALPHA * HTMPJ + BETA * C(2,I,J)
              ENDIF
    240     CONTINUE
          ENDIF
  