anchor.diff:      type(T_matrix) :: Tmat, Ttmat
       type(T_matrix) :: BTmat, BTtmat, BTtKT
       type(hecmwST_matrix) :: hecTKT
  +    real(kind=kreal) :: t1
  +
  +    t1 = hecmw_wtime()
  +    write(0,*) 'INFO: solve_eliminate start', hecmw_wtime()-t1
   
       ndof=hecMAT%NDOF
       allocate(iw2(hecMAT%N*ndof))
  
positive.code:     integer :: i, ls, le, l, j, jb_lag, ib_lag, idof, jdof, ilag, k, ks, ke
      integer, allocatable :: iwUr(:), iwUc(:), iwLr(:), iwLc(:)
      type(hecmwST_matrix) :: hecMATLag
  
      write(0,*) 'INFO: solve_no_eliminate'
  
      ndof = hecMAT%NDOF
      ndof2 = ndof*ndof
  
negative.code: 
      call free_hecmat(hecMATLag)
  
      write(0,*) 'INFO: solve_no_eliminate end', hecmw_wtime()-t1
    end subroutine solve_no_eliminate
  
  end module m_solve_LINEQ_iter_contact
  