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:     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:        ENDIF
      ENDIF
  
      WRITE(*,*)' SpAMM SpAMM SpAMM SpAMM SpAMM SpAMM '
      CALL Multiply(P,P,P2)             ! The only multiplication is a square
      WRITE(*,*)' SpAMM SpAMM SpAMM SpAMM SpAMM SpAMM '
  
      dP2=MATMUL(dP,dP)
      dP2b=Quad2Dense(P2)
  