anchor.diff:    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
  
positive.code: 
  !!$BLOCK_SIZE=256
  !!$      Tr(P)=148.00000000, O(N)/N^3 =  1.00000, TrE =   -453.48859746697627
  
  
    Occ0 = 0.D0
  
negative.code:        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
  