anchor.diff:         qC%Norm=SQRT(norm)
       ENDIF
     END FUNCTION SpAMM_NormReduce
  +
  +  FUNCTION NormReduce_dns(A) RESULT(norm)
  +
  +    REAL(DOUBLE), DIMENSION(:,:) :: A
  +    REAL(DOUBLE) :: norm
  +    INTEGER :: I,J
  +    norm=zero
  +    DO I=1,SIZE(A,1)
  +       DO J=1,SIZE(A,2)
  +          norm=norm+A(I,J)**2
  +       ENDDO
  +    ENDDO
  +    norm=SQRT(norm)
  +  END FUNCTION NormReduce_dns
   
     RECURSIVE SUBROUTINE SpAMM_Multiply(qC,qA,qB,count)
       IMPLICIT NONE
  
positive.code:   A=Zero
    A(1:N_OLD,1:N_OLD)=A_NOPADDING(1:N_OLD,1:N_OLD)
    qP=>Dense2Quad(A)
    !--------------------------------------------------
    DEALLOCATE(A)
    DEALLOCATE(A_NOPADDING)
  
negative.code:   CALL NewQuNode(qTmp1,init=.TRUE.)
    CALL NewQuNode(qTmp2,init=.TRUE.)
    !
  !! BETA CAROTENE EXAMPLE
  !!$ ./a.out bc  bc 256  296 1D-7
  !!$BLOCK_SIZE=1
  !!$1D-4, Tr(P)=148.00410675, O(N)/N^3 =  0.03703, TrE =   -453.64326697675500
  