anchor.diff:        Y(4*i-1)= D(16*i- 7)*X1 + D(16*i- 6)*X2 + D(16*i- 5)*X3 + D(16*i- 4)*X4
         Y(4*i  )= D(16*i- 3)*X1 + D(16*i- 2)*X2 + D(16*i- 1)*X3 + D(16*i- 0)*X4
       enddo
  +    !$OMP END DO
       call MATJAD(hecMAT%N, MJAD, IAJAD, JAJAD, AJAD, JADORD, X, Y, WP1, WP2, WP3, WP4)
     end subroutine hecmw_JAD_MATVEC_44
   
  
positive.code:         W3(IXX)=W3(IXX) + AJAD(K*16- 7)*X1 + AJAD(K*16- 6)*X2 + AJAD(K*16- 5)*X3 + AJAD(K*16- 4)*X4
          W4(IXX)=W4(IXX) + AJAD(K*16- 3)*X1 + AJAD(K*16- 2)*X2 + AJAD(K*16- 1)*X3 + AJAD(K*16- 0)*X4
        enddo
      enddo
  
      do I=1,N
        Y(4*I-3)=Y(4*I-3)+W1(JADORD(I))
        Y(4*I-2)=Y(4*I-2)+W2(JADORD(I))
        Y(4*I-1)=Y(4*I-1)+W3(JADORD(I))
        Y(4*I  )=Y(4*I  )+W4(JADORD(I))
      enddo
    end subroutine MATJAD
  
  end module hecmw_JAD_TYPE_44
  
negative.code:     integer(kind=kint) :: I, K, NZ, IXX
      real(kind=kreal)   :: X1, X2, X3, X4
  
      !$OMP PARALLEL PRIVATE(I,K,X1,X2,X3,IXX)
      !$OMP DO
      do I=1,N
        W1(I)=0.D0
        W2(I)=0.D0
        W3(I)=0.D0
        W4(I)=0.D0
      enddo
      !$OMP END DO
  
      do NZ=1,MJAD
        !$OMP DO
        do K=IAJAD(NZ),IAJAD(NZ+1)-1
          X1=X(JAJAD(K)*4-3)
          X2=X(JAJAD(K)*4-2)
  