anchor.diff:  200 nump1 = num + 1
       IF ( nump1>search ) search = nump1
       DO j = search, Neqns
  -       node = Perm(j)
  -       IF ( Marker(node)>=0 ) THEN
  -          ndeg = Deg(node)
  -          IF ( ndeg<=thresh ) GOTO 300
  -          IF ( ndeg<mindeg ) mindeg = ndeg
  -       ENDIF
  +      node = Perm(j)
  +      IF ( Marker(node)>=0 ) THEN
  +        ndeg = Deg(node)
  +        IF ( ndeg<=thresh ) GOTO 300
  +        IF ( ndeg<mindeg ) mindeg = ndeg
  +      ENDIF
       ENDDO
       GOTO 100
   
  
positive.code:     INTEGER:: i
  
      DO i = 1, N
         C(3,i) = A(3,i) - A(1,i)*B(2,i)
         C(1,i) = A(1,i)*B(1,i)
         C(3,i) = C(3,i)*B(3,i)
         C(1,i) = C(1,i) - C(3,i)*B(2,i)
  
         C(4,i) = A(4,i) - A(2,i)*B(2,i)
         C(2,i) = A(2,i)*B(1,i)
         C(4,i) = C(4,i)*B(3,i)
         C(2,i) = C(2,i) - C(4,i)*B(2,i)
      ENDDO
    END SUBROUTINE V2PROD
  
  
negative.code:     INTEGER:: N
  
      DO i = 1, N
        CALL INVXX(Zz(1,i),Zln(1,i),Diag(1,i),Ndeg)
      ENDDO
    END SUBROUTINE VXPROD
  END MODULE HECMW_SOLVER_DIRECT
  