anchor.diff:  
     END SUBROUTINE Print_Quad
   
  +  SUBROUTINE NewQuNode(qA,count)
  +    INTEGER :: count
  +    TYPE(QuTree), POINTER :: qA
  +    ALLOCATE(qA);  count=count+1
  +    qA%Num=count
  +    qA%Norm=Zero
  +    NULLIFY(qA%Quad00)
  +    NULLIFY(qA%Quad01)
  +    NULLIFY(qA%Quad10)
  +    NULLIFY(qA%Quad11)
  +  END SUBROUTINE NewQuNode
  +
  +
     FUNCTION SpAMM(qA,qB,tolerance) RESULT(qC)
       IMPLICIT NONE
       TYPE(QuTree), POINTER                :: qA,qB
  
positive.code: PROGRAM SpAMM_TEST
    USE SpAMM_TYPES
    IMPLICIT NONE
    INTEGER           ::  N,I,J,M,II,JJ
    REAL(DOUBLE),ALLOCATABLE, DIMENSION(:,:) ::  A, B, C1,C2,CDiff
    TYPE(QuTree),POINTER  :: qA,qB,qC
    REAL(DOUBLE)          :: CNorm,CDiffNorm,CErrBound,CElementErr
    CHARACTER *100        :: Buffer
  
  
negative.code:   INTEGER, PARAMETER :: DOUBLE=KIND(0.D0)          !--Real*8
    REAL(DOUBLE),PARAMETER :: Zero=0D0
    !
    TYPE QuTree
       INTEGER                      :: Siz
       INTEGER                      :: Lev     ! Level of this box
  