anchor.diff:  
         real(kind=kreal)                         :: factor
         real(kind=kreal) :: resid
  +      real(kind=kreal) :: t1, t2
   
         factor = 1.0d0
         if( present(rf) )factor = rf
  
positive.code:         else
            call solve_LINEQ_mumps_contact(hecMESH,hecMAT,fstrMAT)
          endif
        endif
  
        if(paraContactFlag.and.present(conMAT)) then
  
negative.code:       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)
  