anchor.diff:  c
   c
         subroutine formdf_ss(
  -     & bdisp,neq,                                                       ! force
  +     & bintern,neq,                                                     ! force
        & x,d,deld,numnp,                                                  ! global
        & s,stemp,                                                         ! stiff
        & dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,             ! elemnt
  
positive.code: c
            call fill(p,zero,nee)
            call dsymv("u",nee,one,s,nee,dld,ione,zero,p,ione)
            call addfor(bdisp,p,lm(1,indien),lmx(1,indien),neq,nee)
  cdebug          write(6,*) "P vector:"
  cdebug          write(6,*) (p(idebug),idebug=1,nee)
  cdebug          write(6,*) "B vector:"
  cdebug          write(6,*) (bdisp(idebug),idebug=1,neq)
   150      continue
          end do
          matgpt=matgpt+nmatel
        end do
        return
        end
  c
  c version
  c $Id: formdf_ss.f,v 1.5 2004/08/02 21:14:29 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:       integer matgpt,imat,matmodel,nmatel,imatvar,ngtest,ielg,iel,indien
        integer ietype,nen,inddmat,ngauss,nee,ngaussdim,i
        double precision p(60),dld(60)
  cdebug      integer idb
  c
  cdebug      write(6,*) "Hello from formdf_ss_f!"
  c
  