anchor.diff:          ttc = DOT_PRODUCT(TT, spfunc)
           tt0 = DOT_PRODUCT(T0, spfunc)
           ttn = DOT_PRODUCT(TN, spfunc)
  -        CALL MatlMatrix( gausses(LX), D3, D, tincr, coordsys, ttc )
  +        CALL MatlMatrix( gausses(LX), D3, D, tincr, coordsys, ttc, isEp )
   
           ina(1) = ttc
           IF( matlaniso ) THEN
  
positive.code: 
        ELSE
  
          CALL MatlMatrix( gausses(LX), D3, D, tincr , coordsys)
  
        END IF
  
  
negative.code:       ! UPDATE STRAIN and STRESS
        ! ========================================================
  
        if( isElastoplastic(mtype) .OR. mtype == NORTON )then
          isEp = 1
        else
          isEp = 0
        endif
        !gausses(LX)%pMaterial%mtype = ELASTIC
  
        EPSTH = 0.0D0
        IF( PRESENT(tt) .AND. PRESENT(t0) ) THEN
  