anchor.diff:      NDOF = hecMAT%NDOF
       NDOF2 = NDOF*NDOF
       Y=0.0d0
  +
  +    !$OMP PARALLEL PRIVATE(i, idof, jdof)
  +    !$OMP DO
       do i= 1, hecMAT%N
         do idof=1,hecMAT%NDOF
           do jdof=1,hecMAT%NDOF
  
positive.code:     NDOF2=NDOF*NDOF
  
      W=0.0d0
      do NZ=1,MJAD
        do K=IAJAD(NZ),IAJAD(NZ+1)-1
          IXX = K-IAJAD(NZ)+1
          do idof = 1, NDOF
  
negative.code:           end do
          end do
        enddo
        !$OMP END DO
        !$OMP END PARALLEL
      enddo
  
      !$OMP PARALLEL PRIVATE(I,idof)
      !$OMP DO
      do I=1,N
        do idof = 1, NDOF
          Y(NDOF*(I-1)+idof)=Y(NDOF*(I-1)+idof)+W(idof,JADORD(I))
        end do
      enddo
      !$OMP END DO
      !$OMP END PARALLEL
    end subroutine MATJAD
  
  end module hecmw_JAD_TYPE_nn
  