anchor.diff:    REAL(DOUBLE)       :: SpAMM_tolerance, SpAMM_multiplies, SpAMM_dimension
     INTEGER            :: SpAMM_tiles, SpAMM_levels, SpAMM_quadcount
     TYPE(SpAMM_cubes),POINTER :: SpAMM_stream
  -  INTEGER, PARAMETER :: SpAMM_BLOCK_SIZE=256
  +  INTEGER, PARAMETER :: SpAMM_BLOCK_SIZE=2
     !----------------------------------------------------------------
     !===============================================================================
     !  Interface blocks for generic linear algebra routines
  
positive.code:   Occ2 = 0.D0
    Occ3 = 0.D0
    DO I=1,40
      CALL TC2(qP,qTmp1,qTmp2,Half*DBLE(NEl),Occ0,I)
      Occ3 = Occ2
      Occ2 = Occ1
      Occ1 = Occ0
  
negative.code: !!$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
  !!$1D-5, Tr(P)=148.00000991, O(N)/N^3 =  0.24183, TrE =   -453.49417884535342
  !!$1D-6, Tr(P)=148.00000677, O(N)/N^3 =  0.42200, TrE =   -453.48909106857525
  !!$1D-7, Tr(P)=148.00000486, O(N)/N^3 =  0.61572, TrE =   -453.48864613262896
  !!$1D-8, Tr(P)=148.00000090, O(N)/N^3 =  0.77632, TrE =   -453.48860167043756
  !!$1D-10,Tr(P)=148.00000000, O(N)/N^3 =  0.95261, TrE =   -453.48859749814483
  !!$BLOCK_SIZE=2 (SINGLE)
  !!$1D-4, Tr(P)=147.99964905, O(N)/N^3 =  0.11212, TrE =   -453.53015
  !!$1D-5, Tr(P)=148.00001526, O(N)/N^3 =  0.24183, TrE =   -453.49417
  !!$1D-6, Tr(P)=148.00000000, O(N)/N^3 =  0.42211, TrE =   -453.48907
  !!$1D-7, Tr(P)=148.00001526, O(N)/N^3 =  0.61580, TrE =   -453.48868
  !!$1D-8, Tr(P)=148.00001526, O(N)/N^3 =  0.77637, TrE =   -453.48865
  !!$1D-10,Tr(P)=148.00000000, O(N)/N^3 =  0.95261, TrE =   -453.48865
  
  !!$BLOCK_SIZE=256
  !!$      Tr(P)=148.00000000, O(N)/N^3 =  1.00000, TrE =   -453.48859746697627
  
  