anchor.diff:        BTmat%A => new_A
       endif
     end subroutine sort_and_uniq_rows
  -
  -
   
     subroutine hecmw_localmat_add(Amat, Bmat, Cmat)
       implicit none
  
positive.code:   integer(kind=kint), parameter :: cGID = 3       !< index for global ID (used only when cNCOL_ITEM==3)
  
    integer(kind=kint), parameter :: DEBUG = 0
  
  contains
  
  
negative.code:     enddo
      !$omp end do
      !$omp end parallel
      t1 = hecmw_wtime()
      if (TIMER >= 3) write(0, '(A,f10.4)') "###### multiply_mat_mat (2) : ",t1-t0
      t0 = hecmw_wtime()
      call sort_and_uniq_rows(Cmat)
      t1 = hecmw_wtime()
      if (TIMER >= 3) write(0, '(A,f10.4)') "###### multiply_mat_mat (3) : ",t1-t0
    end subroutine multiply_mat_mat
  
    subroutine blk_matmul_add(ndof, A, B, AB)
  