anchor.diff:  
   !!$BLOCK_SIZE=256
   !!$      Tr(P)=148.00000000, O(N)/N^3 =  1.00000, TrE =   -453.48859746697627
  +!!
  +!! 4,3 160 ATOM NANOTUBE EXAMPLE (HF/STO-2G)
  +!! ./a.out tube_4_3__160_Z_18363_Geom#1_Base#2_Clone#1 744 890 1D-4
  +!!$BLOCK_SIZE=2 (DOUBLE)
  +!! 1D-4, Tr(P)=444.99287617, O(N)/N^3 =  0.02251, TrE =   -1738.7744315310249
  +!! 1D-5, Tr(P)=444.99706812, O(N)/N^3 =  0.07218, TrE =   -1738.6751392952985
  +!! 1D-6  Tr(P)=445.00047176, O(N)/N^3 =  0.14222, TrE =   -1738.6525685266417
  +!! 1D-7, Tr(P)=444.99997010, O(N)/N^3 =  0.22885, TrE =   -1738.6515660877203
  +!! 1D-8, Tr(P)=445.00000412, O(N)/N^3 =  0.29912, TrE =   -1738.6515677903926
  +!!$BLOCK_SIZE=2 (SINGLE)
  +!! 1D-8, Tr(P)=445.0000610,    O(/N^3 =  0.12500, TrE =   -1738.6519   << .125 vs .299 diff in computing the norm?
  +!!$BLOCK_SIZE=1024 (DOUBLE)
  +!!      Tr(P)=445.00000000, O(N)/N^3 =  1.00000,  TrE =   -1738.6515695910239
   
   
     Occ0 = 0.D0
  
positive.code:   CALL NewQuNode(qTmp1,init=.TRUE.)
    CALL NewQuNode(qTmp2,init=.TRUE.)
    !
  
  !!$ ./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
  
negative.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 =",F12.5,", Tr(P)=",F12.5,", O(N)/N^3 = ",F8.5)
  
      IF (CR1 < CR2) THEN               ! Too many states
  
  