anchor.diff:        end
   c
   c version
  -c $Id: pquadtet.f,v 1.2 2004/07/07 19:53:01 willic3 Exp $
  +c $Id: pquadtet.f,v 1.3 2004/08/02 21:21:54 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code:       infetype(4)=nee
  c
        do l=1,ngauss
          do i=1,nen
            rr=r(i)*gauss(1,l)
            ss=s(i)*gauss(2,l)
            tt=t(i)*gauss(3,l)
            uu=one-rr-ss-tt
            if(i.le.4) then
              sh(4,i,l)=rr*(two*rr-one)+ss*(two*ss-one)+tt*(two*tt-one)+
       &       uu*(two*uu-one)
              sh(1,i,l)=r(i)*(four*rr-one)-u(i)*(four*uu-one)
              sh(2,i,l)=s(i)*(four*ss-one)-u(i)*(four*uu-one)
              sh(3,i,l)=t(i)*(four*tt-one)-u(i)*(four*uu-one)
            else
              sh(4,i,l)=max(rr,one)*max(ss,one)*max(tt,one)*max(uu,one)
              drr=eight*r(i)*max(ss,one)*max(tt,one)*max(uu,one)
              dss=eight*s(i)*max(rr,one)*max(tt,one)*max(uu,one)
              dtt=eight*t(i)*max(rr,one)*max(ss,one)*max(uu,one)
              duu=(-eight)*u(i)*max(rr,one)*max(ss,one)*max(tt,one)
              sh(1,i,l)=drr-duu
              sh(2,i,l)=dss-duu
              sh(3,i,l)=dtt-duu
            end if
          end do
        end do
        call dcopy(nshsize,sh,ione,shj,ione)
  c
  
negative.code: c
  c...  local constants
  c
  c*      double precision r(10),s(10),t(10),u(10)
  c*      data r/ 1d0, 0d0, 0d0, 0d0,0.5d0,0.0d0,0.5d0,0.5d0,0.0d0,0.0d0/
  c*      data s/ 0d0, 1d0, 0d0, 0d0,0.5d0,0.5d0,0.0d0,0.0d0,0.5d0,0.0d0/
  c*      data t/ 0d0, 0d0, 1d0, 0d0,0.0d0,0.5d0,0.5d0,0.0d0,0.0d0,0.5d0/
  c*      data u/ 0d0, 0d0, 0d0, 1d0,0.0d0,0.0d0,0.0d0,0.5d0,0.5d0,0.5d0/
  c
  c...  intrinsic functions
  c
  