anchor.diff:          else
             call solve_LINEQ_mumps_contact(hecMESH,hecMAT,fstrMAT)
           endif
  +      endif
  +
  +      t2 = hecmw_wtime()
  +      if (hecmw_mat_get_timelog(hecMAT) .eq. 1) then
  +        write(*,*) ' solve time :', t2 - t1
         endif
   
         if(paraContactFlag.and.present(conMAT)) then
  
positive.code:       else
          call hecmw_mat_ass_equation( hecMESH, hecMAT )
        endif
  
        if( hecMAT%Iarray(99)==1 )then
          call solve_LINEQ_iter_contact(hecMESH,hecMAT,fstrMAT)
  
negative.code: 
        real(kind=kreal)                         :: factor
        real(kind=kreal) :: resid
        real(kind=kreal) :: t1, t2
  
        factor = 1.0d0
        if( present(rf) )factor = rf
  