anchor.diff:  c
   c...  local variables
   c
  -      integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
  -      double precision g1,g2,rr,ss,tt,uu,drr,dss,dtt,duu
  +c*      integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
  +      integer nen,ngauss,nec,nee,l,nshsize,ngssize
  +c*      double precision g1,g2,rr,ss,tt,uu,drr,dss,dtt,duu
  +      double precision g1,g2,rr,ss,tt,uu
         double precision tetvol
  +cdebug      integer idb,jdb
  +cdebug      double precision dbsum
  +c
  +cdebug      write(6,*) "Hello from pquadtet_f!"
   c
   c...  definitions
   c
  
positive.code:       end
  c
  c version
  c $Id: pquadtet.f,v 1.2 2004/07/07 19:53:01 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  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
  