anchor.diff:        Y(3*i -1)= D(9*i-5)*X1 + D(9*i-4)*X2 + D(9*i-3)*X3
         Y(3*i   )= D(9*i-2)*X1 + D(9*i-1)*X2 + D(9*i  )*X3
       enddo
  +    !$OMP END DO
  +    !$OMP END PARALLEL
  +
       call MATJAD(hecMAT%N, MJAD, IAJAD, JAJAD, AJAD, JADORD, X, Y, WP1, WP2, WP3)
     end subroutine hecmw_JAD_MATVEC_33
   
  
positive.code:         W2(IXX)=W2(IXX) + AJAD(K*9-5)*X1 + AJAD(K*9-4)*X2 + AJAD(K*9-3)*X3
          W3(IXX)=W3(IXX) + AJAD(K*9-2)*X1 + AJAD(K*9-1)*X2 + AJAD(K*9-0)*X3
        enddo
      enddo
  
      do I=1,N
        Y(3*I-2)=Y(3*I-2)+W1(JADORD(I))
        Y(3*I-1)=Y(3*I-1)+W2(JADORD(I))
        Y(3*I  )=Y(3*I  )+W3(JADORD(I))
      enddo
    end subroutine MATJAD
  
  end module hecmw_JAD_TYPE_33
  
negative.code: 
      D => hecMAT%D
  
      !$OMP PARALLEL PRIVATE(i)
      !$OMP DO
      do i= 1, hecMAT%N
        X1= X(3*i-2)
        X2= X(3*i-1)
  