anchor.diff:                DO 180 I = 1,M
                   C(1,I,J) = ZERO
                   C(2,I,J) = ZERO
  -                C(3,I,J) = ZERO
  -                C(4,I,J) = ZERO
     180         CONTINUE
               ELSEIF ( .NOT. BONE ) THEN
                 DO 190 I = 1,M
                   C(1,I,J) = BETA * C(1,I,J) 
                   C(2,I,J) = BETA * C(2,I,J) 
  -                C(3,I,J) = BETA * C(3,I,J) 
  -                C(4,I,J) = BETA * C(4,I,J) 
     190         CONTINUE
               ENDIF
   C
               DO 200 L = 1,K
                 HTMPS = ALPHA * B(1,J,L)
  -              HTMPI = ALPHA * B(2,J,L)
  -              HTMPJ = ALPHA * B(3,J,L)
  -              HTMPK = ALPHA * B(4,J,L)
  +              HTMPJ = ALPHA * B(2,J,L)
                 DO 210 I = 1,M
   C               Hamilton Product
                   C(1,I,J) = C(1,I,J) + 
  -     $            A(1,I,L)*HTMPS - A(2,I,L)*HTMPI -
  -     $            A(3,I,L)*HTMPJ - A(4,I,L)*HTMPK
  -
  +     $            A(1,I,L)*HTMPS - A(2,I,L)*CONJG(HTMPJ)
                   C(2,I,J) = C(2,I,J) + 
  -     $            A(1,I,L)*HTMPI + A(2,I,L)*HTMPS +
  -     $            A(3,I,L)*HTMPK - A(4,I,L)*HTMPJ
  -
  -                C(3,I,J) = C(3,I,J) + 
  -     $            A(1,I,L)*HTMPJ - A(2,I,L)*HTMPK +
  -     $            A(3,I,L)*HTMPS + A(4,I,L)*HTMPI
  -
  -                C(4,I,J) = C(4,I,J) + 
  -     $            A(1,I,L)*HTMPK + A(2,I,L)*HTMPJ -
  -     $            A(3,I,L)*HTMPI + A(4,I,L)*HTMPS
  +     $            A(1,I,L)*HTMPJ + A(2,I,L)*CONJG(HTMPS)
     210         CONTINUE
     200       CONTINUE
     170     CONTINUE
  
positive.code:           DO 220 J = 1,N
            DO 220 I = 1,M
              HTMPS = ZERO 
              HTMPI = ZERO 
              HTMPJ = ZERO 
              HTMPK = ZERO 
              DO 230 L = 1,K
  C             Hamilton Product (IMPLIED CONJ(A) + CONJ(B))
                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)
    230       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
    220     CONTINUE
          ELSE
  
negative.code:             HTMPJ = ZERO 
              HTMPK = ZERO 
              DO 270 L = 1,K
  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 * 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
  