anchor.diff:    INTEGER, PARAMETER :: INT8=SELECTED_INT_KIND(18) !--Integer*8
     !
     INTEGER, PARAMETER :: SINGLE=KIND(0.0)           !--Real*4
  -  INTEGER, PARAMETER :: DOUBLE=KIND(0.D0)          !--Real*8
  +  INTEGER, PARAMETER :: DOUBLE=KIND(0.00)          !--Real*8
     REAL(DOUBLE),PARAMETER :: Zero=0D0,Half=5D-1,One=1D0,Two=2D0
     !
     TYPE QuTree
  
positive.code:   REAL(DOUBLE)       :: SpAMM_tolerance, SpAMM_multiplies, SpAMM_dimension
    INTEGER            :: SpAMM_tiles, SpAMM_levels, SpAMM_quadcount
    TYPE(SpAMM_cubes),POINTER :: SpAMM_stream
    INTEGER, PARAMETER :: SpAMM_BLOCK_SIZE=256
    !----------------------------------------------------------------
    !===============================================================================
    !  Interface blocks for generic linear algebra routines
  
negative.code:   Occ2 = 0.D0
    Occ3 = 0.D0
    DO I=1,40
      CALL TC2(qP,qTmp1,qTmp2,Half*FLOAT(NEl),Occ0,I)
      Occ3 = Occ2
      Occ2 = Occ1
      Occ1 = Occ0
  