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