anchor.diff:    INTEGER, PARAMETER :: INT8=SELECTED_INT_KIND(18) !--Integer*8
     !
     INTEGER, PARAMETER :: SINGLE=KIND(0.0)           !--Real*4
  -  INTEGER, PARAMETER :: DOUBLE=KIND(0.00)          !--Real*8
  +  INTEGER, PARAMETER :: DOUBLE=KIND(0.0)           !--Real*8
     REAL(DOUBLE),PARAMETER :: Zero=0D0,Half=5D-1,One=1D0,Two=2D0
     !
     TYPE QuTree
  
positive.code:     REAL(DOUBLE) :: Norm, CR1, CR2, TrP, TrP2,    TrdP,TrdP2
      INTEGER      :: I,  II,J
  
      INTEGER, PARAMETER:: NN=8
      REAL(DOUBLE),DIMENSION(NN,NN) :: dP,dP2,dTmp1, dTmp2,dP2b
      !-------------------------------------------------------------------------------
      dP=Quad2Dense(P)
  
      IF (I.EQ.1) THEN
  
negative.code:        ALLOCATE(qA%Quad01)
         ALLOCATE(qA%Quad10)
         ALLOCATE(qA%Quad11)
  
         Tier=qA%Lev+1
         qA%Quad00%Lev=Tier
         qA%Quad01%Lev=Tier
  