anchor.diff:      NDOF2=NDOF*NDOF
   
       W=0.0d0
  +
       do NZ=1,MJAD
  +      !$OMP PARALLEL PRIVATE(K,IXX,idof,jdof)
  +      !$OMP DO
         do K=IAJAD(NZ),IAJAD(NZ+1)-1
           IXX = K-IAJAD(NZ)+1
           do idof = 1, NDOF
  
positive.code:           end do
          end do
        enddo
      enddo
  
      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
    end subroutine MATJAD
  
  end module hecmw_JAD_TYPE_nn
  
negative.code:         end do
        end do
      enddo
      !$OMP END DO
      !$OMP END PARALLEL
      call MATJAD(hecMAT%N,hecMAT%NDOF, MJAD, IAJAD, JAJAD, AJAD, JADORD, X, Y, WP)
    end subroutine hecmw_JAD_MATVEC_nn
  
  