anchor.diff:      integer :: nsize
   
       if(present(chunk_size)) then
  -      nsize = chunk_size
  +      nsize = chunk_size &
  +            * (NXY_e - NXY_s + 1)
       else
         nsize = (NK_e  - NK_s  + 1) * NBoccmax &
               * (NXY_e - NXY_s + 1)
  
positive.code:     FLOP = FLOP_scalar + (FLOP_reduction + FLOP_scatter) * sum(Mps(a_tbl(:)))
  
      if(present(chunk_size)) then
        nsize = chunk_size
      else
        nsize = (NK_e  - NK_s  + 1) * NBoccmax &
              * (NXY_e - NXY_s + 1)
  
negative.code:     integer :: nsize
  
      if(present(chunk_size)) then
        nsize = chunk_size &
              * (NXY_e - NXY_s + 1)
      else
        nsize = (NK_e  - NK_s  + 1) * NBoccmax &
              * (NXY_e - NXY_s + 1)
  