anchor.diff:  !!$1D-7, Tr(P)=148.00003758, O(N)/N^3 =  0.34967, TrE =   -453.48904088853374
   !!$1D-8  Tr(P)=148.00000463, O(N)/N^3 =  0.52270, TrE =   -453.48864186234346
   !!$1D-10 Tr(P)=148.00000006, O(N)/N^3 =  0.81716, TrE =   -453.48859793980563
  -
   !!$
   !!$BLOCK_SIZE=2 (DOUBLE)
   !!$1D-4, Tr(P)=147.99964372, O(N)/N^3 =  0.11211, TrE =   -453.53013549760510
  
positive.code:     CR2 = ABS(2.D0*TrP - TrP2 - Norm) ! CR2 = Occupation error criteria
  
      WRITE(*,33)I,NORM,TrP,SpAMM_multiplies/SpAMM_dimension**3
  33  FORMAT(I4,", N =",F8.5,", Tr(P)=",F8.5,", O(N)/N^3 = ",F8.5)
  
      IF (CR1 < CR2) THEN               ! Too many states
  
  
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
  