anchor.diff:  c
   c...  user-defined functions
   c
  -      double precision gquad,dgquad
  +      double precision gquad,dgquad,hquad,dhquadr,dhquads,dhquadt
   c
   c...  local variables
   c
         integer i,l,l1,l2,l3,nshsize,ngssize
  -      double precision g1,w1,w2,rr,ss,tt,rrw,ssw,ttw,drr,dss,dtt,beta
  -      double precision betai
  +      double precision g1,w1,w2,rr,ss,tt,rrw,ssw,ttw,drr,dss,dtt
  +      double precision uu,dur,dus,dut,v,vi,rp,ri,sp,si,tp,ti
   c
   c...  function definitions
   c
  -      gquad(beta,betai)=(half*(one-abs(betai))+half)*(one+beta*betai+
  -     & (abs(betai)-one)*beta*beta)
  -      dgquad(beta,betai)=half*betai-two*beta*(one-abs(betai))
  +      gquad(v,vi)=(half*(one-abs(vi))+half)*(one+v*vi+(abs(vi)-one)*v*v)
  +      dgquad(v,vi)=(half*(one-abs(vi))+half)*(vi+two*v*(abs(vi)-one))
  +      hquad(rp,sp,tp,ri,si,ti)=abs(ri*si*ti)*(ri*rp+si*sp+ti*tp-two)+
  +     & one-abs(ri*si*ti)
  +      dhquadr(ri,si,ti)=abs(ri*si*ti)*ri
  +      dhquads(ri,si,ti)=abs(ri*si*ti)*si
  +      dhquadt(ri,si,ti)=abs(ri*si*ti)*ti
   c
   c...  definitions
   c
  
positive.code: c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nsd,nenmax,ngaussmax,intord
        integer ngauss,nen
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
        double precision r(20),s(20),t(20)
        data r/-1d0, 1d0, 1d0,-1d0,-1d0, 1d0, 1d0,-1d0,
  
negative.code:       end
  c
  c version
  c $Id: pquadhex.f,v 1.2 2004/07/06 19:14:36 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  