anchor.diff:    INTEGER           ::  N,I,J,M,II,JJ
     REAL(DOUBLE),ALLOCATABLE, DIMENSION(:,:) ::  A, B, C1,C2,CDiff
     TYPE(QuTree),POINTER  :: qA,qB,qC
  -  REAL(DOUBLE)          :: CDiffNorm,CErrBound
  +  REAL(DOUBLE)          :: CNorm,CDiffNorm,CErrBound,CElementErr
     CHARACTER *100        :: Buffer
   
     CALL GETARG(1,BUFFER)
  
positive.code:        CALL NewQuNode(qC%Quad11,SpAMM_quadcount)
      ! SpAMM=Sparse Approximate Matrix-Matrix
  
  
      IF(qA%Siz<=2*BLOCK_SIZE)THEN
  
      do00x00=qA%Quad00%Norm*qB%Quad00%Norm>SpAMM_tolerance
      do01x10=qA%Quad01%Norm*qB%Quad10%Norm>SpAMM_tolerance
  
negative.code:        CTmp=MATMUL(qA%Blok,qB%Blok)
         ! Increment
         qC%Blok=qC%Blok+CTmp
  !!$
  !!$    IF(qA%Box(2,1)==7.AND.qB%Box(1,2)==6)THEN
  !!$!    IF(qC%Box(1,1)==7.AND.qC%Box(2,2)==6)THEN
  !!$!    WRITE(*,33)qA%Box,qB%Box,qC%Box,CTMP
  !!$    WRITE(*,*)qC%Box(1,1),qC%Box(2,2),CTMP,qC%Blok,qA%Norm*qB%Norm
  !!$ ENDIF
  !!$33  format(4(I2,', '), ' || ',4(I2,', '),' || ',4(I2,', '),F20.10)
         ! Bounds
         qC%Norm=qC%Norm+SUM(CTmp(1:BLOCK_SIZE,1:BLOCK_SIZE)**2)
         RETURN
  