anchor.diff:              HTMPJ = ZERO 
               HTMPK = ZERO 
               DO 270 L = 1,K
  -C             Hamilton Product
  -              HTMPS = HTMPS + A(1,L,I)*B(1,J,L) + A(2,L,I)*B(2,J,L) +
  -     $                        A(3,L,I)*B(3,J,L) + A(4,L,I)*B(4,J,L)
  -
  -              HTMPI = HTMPI - A(1,L,I)*B(2,J,L) + A(2,L,I)*B(1,J,L) -
  -     $                        A(3,L,I)*B(4,J,L) + A(4,L,I)*B(3,J,L)
  -
  -              HTMPJ = HTMPJ - A(1,L,I)*B(3,J,L) + A(2,L,I)*B(4,J,L) +
  -     $                        A(3,L,I)*B(1,J,L) - A(4,L,I)*B(2,J,L)
  -
  -              HTMPK = HTMPK - A(1,L,I)*B(4,J,L) - A(2,L,I)*B(3,J,L) +
  -     $                        A(3,L,I)*B(2,J,L) + A(4,L,I)*B(1,J,L)
  +C             Hamilton Product (IMPLIED CONJ(B))
  +              HTMPS = HTMPS + 
  +     $          A(1,L,I)*CONJG(B(1,J,L)) + A(2,L,I)*CONJG(B(2,J,L))
  +              HTMPJ = HTMPJ - 
  +     $          A(1,L,I)*B(2,J,L) + A(2,L,I)*B(1,J,L)
     270       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
     260     CONTINUE
           ELSE
  
positive.code:           DO 80 J = 1,N
            DO 80 I = 1,M
              HTMPS = ZERO 
              HTMPI = ZERO 
              HTMPJ = ZERO 
              HTMPK = ZERO 
              DO 90 L = 1,K
  C             Hamilton Product (IMPLIED CONJ(A))
                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)
    90        CONTINUE
  C
              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)
              ENDIF
    80      CONTINUE
          ELSE
  
negative.code:           DO 10 I = 1,M
              C(1,I,J) = ZERO
              C(2,I,J) = ZERO
    10      CONTINUE
          ELSE
            DO 20 J = 1,N
            DO 20 I = 1,M
              C(1,I,J) = BETA * C(1,I,J) 
              C(2,I,J) = BETA * C(2,I,J) 
    20      CONTINUE
          ENDIF
          RETURN ! Nothing more to do
  