anchor.diff:      ! else
       !   write(0,*) 'DEBUG: BTtmat2 is correct'
       ! endif
  -    if (DEBUG >= 4) then
  +    if (DEBUG >= 3) then
         write(700+hecmw_comm_get_rank(),*) 'DEBUG: BTtmat(MPC)'
  -      call hecmw_localmat_write(BTtmat,700+hecmw_comm_get_rank())
  +      if (DEBUG == 3) then
  +        call hecmw_localmat_write_ij(BTtmat,700+hecmw_comm_get_rank())
  +      else
  +        call hecmw_localmat_write(BTtmat,700+hecmw_comm_get_rank())
  +      endif
       endif
       call hecmw_trimatmul_TtKT(hecMESH, BTtmat, hecMAT, BTmat, iwS, n_mpc, hecTKT)
   
  
positive.code:     enddo
    end subroutine hecmw_localmat_write
  
    subroutine hecmw_localmat_blocking(Tmat, ndof, BTmat)
      implicit none
      type (hecmwST_local_matrix), intent(in) :: Tmat
  
negative.code:     call hecmw_localmat_init_with_hecmat(W1mat, hecMAT)
      if (DEBUG >= 3) then
        write(700+hecmw_comm_get_rank(),*) 'BKmat (hecMAT)'
        if (DEBUG == 3) then
          call hecmw_localmat_write_ij(W1mat, 700+hecmw_comm_get_rank())
        else
          call hecmw_localmat_write(W1mat, 700+hecmw_comm_get_rank())
        endif
      endif
      call hecmw_localmat_add(BKmat, W1mat, W2mat)
      call hecmw_localmat_free(BKmat)
  