anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer nstate,nprop,ierr
  +      integer nstate,nstate0,nprop,ierr
         integer iddmat(nstr,nstr)
         character errstrng*(*)
         logical matchg
  -      double precision state(nstr,nstate),dstate(nstr,nstate),ee(nstr)
  +      double precision state(nstr,nstate),dstate(nstr,nstate)
  +      double precision state0(nstate0),ee(nstr)
         double precision dmat(nddmat),prop(nprop),tmax
   c
   c...  included dimension and type statements
  
positive.code: c       
  
  c version
  c $Id: mat_1.f,v 1.14 2005/01/18 19:32:56 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  
negative.code: c
  c...  subroutine arguments
  c
        integer nstate,nstate0,ierr
        double precision state(nstr,nstate),state0(nstate0),ee(nstr)
        double precision dmat(nddmat)
        character errstrng*(*)
  c
        call dcopy(nstr,ee,ione,state(1,2),ione)
        call dcopy(nstr,state0,ione,state(1,1),ione)
        call dspmv("u",nstr,one,dmat,state(1,2),ione,one,state(1,1),ione)
        return
        end
  c
  c
        subroutine td_matinit_1(state,dstate,state0,dmat,prop,rtimdat,
       & rgiter,ntimdat,iddmat,tmax,nstate,nstate0,nprop,matchg,ierr,
       & errstrng)
  c
  c...  subroutine to form the material matrix for an integration point
  c     for the time-dependent solution.  This routine is meant to be
  