anchor.diff:         P=>Copy(Tmp1)
       ENDIF
       TrP=Trace(P)
  +!    CALL Filter_qutree(P,SpAMM_tolerance*1D-2)
  +    !
     END SUBROUTINE TC2
   
   
  
positive.code:                         SpAMM_dimension
    INTEGER            :: SpAMM_tiles, SpAMM_levels, SpAMM_quadcount
    TYPE(SpAMM_cubes),POINTER :: SpAMM_stream
    INTEGER, PARAMETER :: SpAMM_BLOCK_SIZE=4
    !----------------------------------------------------------------
    !===============================================================================
    !  Interface blocks for generic linear algebra routines
  
negative.code:     !
      CALL NewQuNode(qA,init=.TRUE.)
      CALL SpAMM_Mat2Quad(A,qA)
      qA%Norm=SQRT(SpAMM_NormReduce(qA))
    END FUNCTION Dense2Quad
  
    RECURSIVE SUBROUTINE SpAMM_Mat2Quad(A,qA)
  