anchor.diff:        integer l
         double precision shd(4,nenmax,ngaussmax),b(6,3*nenmax)
         double precision shbar(4,nenmax),db(6,3*nenmax),det(ngaussmax)
  +      double precision dtmp(6,6)
         double precision vol
  +c
  +cdebug      integer idb,jdb
   c
   c...form shape functions for each integration point
   c
  
positive.code: c
        do l=1,ngauss
          call bmatrix(b,shd(1,1,l),shbar,nen)
          call dspmv("u",nstr,det(l),dmat(1,l),b,ione,zero,db,ione)
          call dgemm("t","n",nee,nee,nstr,one,b,nstr,db,nstr,one,s,nee)
        end do
        return
        end
  c
  c version
  c $Id: stiff_ss.f,v 1.5 2005/03/19 01:49:49 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c...  subroutine arguments
  c
        integer iel,nen,ngauss,nee,ierr
        integer iddmat(nstr,nstr),ien(nen)
        character errstrng*(*)
        double precision xl(nsd,nen),dmat(nddmat,ngauss)
        double precision gauss(nsd+1,ngauss),sh(nsd+1,nen,ngauss)
  