anchor.diff:        end
   c
   c version
  -c $Id: td_strs_mat_cmp_ss.f,v 1.3 2004/07/09 01:38:46 willic3 Exp $
  +c $Id: td_strs_mat_cmp_ss.f,v 1.4 2004/07/21 18:21:19 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code:       double precision dx(ndof,numnp),tfault(ndof,numfn)
        double precision s(neemax*neemax),stemp(neemax*neemax)
        double precision state(nstr,nstatesz),dstate(nstr,nstatesz)
        double precision dmat(nddmat,ndmatsz),prop(nprop)
        double precision gauss(nsd+1,ngaussmax,netypes)
        double precision sh(nsd+1,nenmax,ngaussmax,netypes)
        double precision shj(nsd+1,nenmax,ngaussmax,netypes)
  
negative.code:         do l=1,ngauss
            call td_strs_mat(state(1,indstateg),dstate(1,indstateg),
       &     ee(nstr*(l-1)),dmat(1,inddmatg),prop,rtimdat,rgiter,ntimdat,
       &     iddmat,tmax,nstate,nprop,matchg,ierr,errstrng)
            if(ierr.ne.izero) return
            tminmax=min(tmax,tminmax)
            call dcopy(nstr,dstate(1,indstateg),ione,scur(nstr*(l-1)),
       &     ione)
            indstateg=indstateg+incstate
  