Line data    Source code

       1             : !--------------------------------------------------------------------------------
       2             : ! Copyright (c) 2021 Peter Grünberg Institut, Forschungszentrum Jülich, Germany
       3             : ! This file is part of FLEUR and available as free software under the conditions
       4             : ! of the MIT license as expressed in the LICENSE file in more detail.
       5             : !--------------------------------------------------------------------------------
       6             : 
       7             : module m_fermifct
       8             : CONTAINS
       9           0 :   real elemental function fermifct(e,e_f,tkb)
      10             :     implicit NONE
      11             :     REAL, INTENT(IN) :: e,e_f,tkb
      12             : 
      13             :     REAL:: expo
      14           0 :     expo= EXP(-ABS(e-e_f)/tkb)
      15           0 :     IF (e<e_f) THEN
      16           0 :        fermifct = 1./ (expo+1.)
      17             :     ELSE
      18           0 :        fermifct= expo/ (expo+1.)
      19             :     ENDIF
      20           0 :   end function
      21             : end module