anchor.diff:      ENDIF
       add_quadcount=1
       CALL SpAMM_Add(qA,qB,qC,count=add_quadcount)
  -    qC%Norm=SQRT(qC%Norm)
  +    qC%Norm=SQRT(SpAMM_NormReduce(qC))
     END SUBROUTINE Add_qutree_pls_qutree
   
     SUBROUTINE NodeMssg(Mssg,Node,count)
  
positive.code:   SUBROUTINE Add_qutree_pls_qutree(qA,qB,qC)
      IMPLICIT NONE
      TYPE(QuTree), POINTER       :: qA,qB,qC
      INTEGER :: add_quadcount
  
      IF(.NOT.ASSOCIATED(qC))THEN
         CALL NewQuNode(qC,init=.TRUE.)
      ELSE
  
negative.code: 
    RECURSIVE SUBROUTINE Print_Quad(qA)
      TYPE(QuTree),POINTER :: qA
  
      IF(ASSOCIATED(qA))THEN
         WRITE(*,111)qA%Box(:,1),qA%Box(:,2),qA%Num,qA%Siz,qA%Norm
      ENDIF
  
      IF(ASSOCIATED(qA%Quad00))&
           WRITE(*,111)qA%Quad00%Box(:,1),qA%Quad00%Box(:,2),qA%Quad00%Num,qA%Quad00%Siz,qA%Quad00%Norm
      IF(ASSOCIATED(qA%Quad01))&
           WRITE(*,111)qA%Quad01%Box(:,1),qA%Quad01%Box(:,2),qA%Quad01%Num,qA%Quad01%Siz,qA%Quad01%Norm
      IF(ASSOCIATED(qA%Quad10))&
           WRITE(*,111)qA%Quad10%Box(:,1),qA%Quad10%Box(:,2),qA%Quad10%Num,qA%Quad10%Siz,qA%Quad10%Norm
      IF(ASSOCIATED(qA%Quad11))&
           WRITE(*,111)qA%Quad11%Box(:,1),qA%Quad11%Box(:,2),qA%Quad11%Num,qA%Quad11%Siz,qA%Quad11%Norm
      WRITE(*,*)' ========================================'
  