anchor.diff:        end
   c
   c version
  -c $Id: elas_strs_mat_cmp_ss.f,v 1.7 2004/07/21 19:14:50 willic3 Exp $
  +c $Id: elas_strs_mat_cmp_ss.f,v 1.8 2004/08/02 21:11:35 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
        ngtest=izero
        if(imatvar.eq.izero) ngtest=ngaussmax
        incstate=nstr*nstate
  c
  c...  loop over elements in a material group
  c
  
negative.code: c     scur
  c
          do l=1,ngauss
            call elas_strs_mat(dstate(1,indstateg),ee(nstr*(l-1)+1),
       &     dmat(1,inddmatg),nstate,ierr,errstrng)
            if(ierr.ne.izero) return
            call dcopy(nstr,dstate(1,indstateg),ione,scur(nstr*(l-1)+1),
       &     ione)
            indstateg=indstateg+nstate
            inddmatg=inddmatg+ione
          end do
  c
  c...  compute equivalent nodal loads and add them to global load
  