Merge pull request dftbplus#1432 from bhourahine/trivialHybTypeErr

Correct real to int and make fn pure

src/dftbp/dftb/hybridxc.F90

@@ -1068,22 +1068,25 @@ subroutine updateCoords_kpts(this, env, symNeighbourList, nNeighbourCamSym, skOv
 
   contains
 
-    !> Returns the number of non-zero elements in a descending array of reals.
+    !> Returns the number of non-zero elements in a descending array of non-negative reals.
     function getNumberOfNonZeroElements(array) result(nNonZeroEntries)
 
       !> Descending, one-dimensional, real-valued array to search
       real(dp), intent(in) :: array(:)
 
       !> Number of non-zero entries
-      real(dp) :: nNonZeroEntries
+      integer :: nNonZeroEntries
 
       !! iterates over all array elements
       integer :: ii
 
       nNonZeroEntries = 0
 
       do ii = 1, size(array)
-        if (array(ii) < 1e-16_dp) return
+        if (array(ii) < epsilon(1.0_dp)) then
+          @:ASSERT(all(array(ii:) >= 0.0_dp))
+          return
+        end if
         nNonZeroEntries = ii
       end do