anchor.diff:              indstateg=indstate
               ngauss=infetype(1,ietype)
               do l=1,ngauss
  -              call dcopy(nstr,dstate(1,indstateg),ione,
  +              call daxpy(nstr2,sub,state(1,indstateg),ione,
  +     &         dstate(1,indstateg),ione)
  +              call daxpy(nstr4,one,dstate(1,indstateg),ione,
        &         state(1,indstateg),ione)
  -              call daxpy(nstr,sub,state(1,indstateg+ione),ione,
  -     &         dstate(1,indstateg+ione),ione)
  -              call daxpy(nstr,one,dstate(1,indstateg+ione),ione,
  -     &         state(1,indstateg+ione),ione)
  -              call daxpy(nstr,one,dstate(1,indstateg+itwo),ione,
  -     &         state(1,indstateg+itwo),ione)
  -              call daxpy(nstr,one,dstate(1,indstateg+ithree),ione,
  -     &         state(1,indstateg+ithree),ione)
                 indstateg=indstateg+nstate
               end do
             end do
  
positive.code:             indstateg=indstate
              ngauss=infetype(1,ietype)
              do l=1,ngauss
                call dcopy(nstr,dstate(1,indstateg),ione,
       &         state(1,indstateg),ione)
                call daxpy(nstr,sub,state(1,indstateg+ione),ione,
       &         dstate(1,indstateg+ione),ione)
                call daxpy(nstr,one,dstate(1,indstateg+ione),ione,
       &         state(1,indstateg+ione),ione)
                call daxpy(nstr,one,dstate(1,indstateg+itwo),ione,
       &         state(1,indstateg+itwo),ione)
                indstateg=indstateg+nstate
              end do
            end do
  
negative.code:       end
  c
  c version
  c $Id: update_state.f,v 1.6 2004/08/12 02:33:46 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  