anchor.diff:        end
   c
   c version
  -c $Id: formf_ss.f,v 1.1 2004/06/16 21:14:06 willic3 Exp $
  +c $Id: formf_ss.f,v 1.2 2004/06/21 21:03:39 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c...  form element stiffness matrix
  c
          call formes_ss(
       &   x,nsd,ndof,numnp,                                              ! global
       &   s,stemp,                                                       ! stiff
       &   dmat(1,inddmat),ien(indien),lm(1,indien),iddmat,nstr,nddmat,   ! elemnt
       &   iel,                                                           ! elemnt
       &   gauss(1,1,ietype),sh(1,1,1,ietype),shj(1,1,1,ietype),          ! eltype
       &   ngauss,ngaussdim,nen,nee,                                      ! eltype
       &   skew,nskdim,numrot,                                            ! skew
       &   getshape,bmatrix,                                              ! bbar
       &   ierr)                                                          ! errcode
  c
          if(ierr.ne.izero) return
  c
          call lflteq(dl,dfault(1,i),nfault(1,i),ien(indien),nen,ndof)
  	call dsymv("u",nee,one,s,nee,dl,ione,zero,p,ione)
          call addfor(bdisp,p,lm(1,indien),lmx(1,indien),neq,nee)
        end do
  
negative.code:       double precision sh(nsd+1,nenmax,ngaussmax,netypes)
        double precision shj(nsd+1,nenmax,ngaussmax,netypes)
        double precision skew(nskdim,numnp)
  c
  c...  external routines
  c
  