anchor.diff:        end
   c
   c version
  -c $Id: update_state.f,v 1.4 2004/07/12 21:07:20 willic3 Exp $
  +c $Id: update_state.f,v 1.5 2004/08/02 21:32:09 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code:             do l=1,ngauss
                call dcopy(nstr,dstate(1,indstateg),ione,
       &         state(1,indstateg),ione)
                call daxpy(nstr,sub,state(1,indstateg+nstr),ione,
       &         dstate(1,indstateg+nstr),ione)
                call daxpy(nstr,one,dstate(1,indstateg+nstr),ione,
       &         state(1,indstateg+nstr),ione)
              end do
              indstateg=indstateg+incstate
            end do
          else if(nstate.eq.ithree) then
            do ielg=matgpt,matgpt+nmatel-1
  
negative.code:             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
          else if(nstate.eq.ifour) then
            do ielg=matgpt,matgpt+nmatel-1
  