anchor.diff:      ! choose slave DOFs to be eliminated with Lag. DOFs
       call choose_slaves(hecMAT, fstrMAT, iw2, iwS, wSL)
       call make_wSU(fstrMAT, hecMAT%N, ndof, iw2, wSU)
  -    write(0,*) 'INFO: Slave DOFs successfully chosen'
  +    write(0,*) 'INFO: Slave DOFs successfully chosen', hecmw_wtime()-t1
   
       ! make transformation matrix and its transpose
       call make_Tmat(hecMAT, fstrMAT, iw2, wSL, Tmat)
       !call write_Tmat(Tmat, 0)
       call make_Ttmat(hecMAT, fstrMAT, iw2, iwS, wSU, Ttmat)
       !call write_Tmat(Ttmat, 0)
  -    write(0,*) 'INFO: Making Tmat and Ttmat done'
  +    write(0,*) 'INFO: Making Tmat and Ttmat done', hecmw_wtime()-t1
   
       ! make 3x3-block version of Tmat
       call blocking_Tmat(Tmat, ndof, BTmat)
  
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:     hecMATLag%Iarray=hecMAT%Iarray
      hecMATLag%Rarray=hecMAT%Rarray
  
      write(0,*) 'INFO: made hecMATLag', hecmw_wtime()-t1
  
      call hecmw_solve_33(hecMESH, hecMATLag)
  
      do i = 1, hecMAT%N*ndof+fstrMAT%num_lagrange
  