anchor.diff:          !       this symmetry index.  Ignored if less than sym0 or greater than
           !       sym_max.
   
  -        use const, only: i0, p
  +        use const, only: i0, p, depsilon
           use checking, only: check_allocate
           use errors, only: stop_all
   
  
positive.code:                                                             sp_eigv_sum = sp_eigv_sum + &
                                                                  basis_fns(tmp_occ_list(iel))%sp_eigv
                                                          end do
                                                          if (sp_eigv_sum < eigv_sum) then
                                                              curr_occ_list = tmp_occ_list
                                                              eigv_sum = sp_eigv_sum
                                                          end if
  
negative.code:                                             do iel = 1, sys%nel
                                                  sp_eigv_sum = sp_eigv_sum + basis_fns(tmp_occ_list(iel))%sp_eigv
                                              end do
                                              if (sp_eigv_sum+depsilon < eigv_sum) then
                                                  curr_occ_list = tmp_occ_list
                                                  eigv_sum = sp_eigv_sum
                                              end if
  