reldiff_m.f90 Source File

module reldiff_m
    use, intrinsic :: iso_c_binding
    implicit none

    contains

    pure real function relative_difference(a, b)
      real(C_DOUBLE), intent(in) :: a
      real(C_DOUBLE), intent(in) :: b

      relative_difference = abs(a - b) / max(abs(a), abs(b))
    end function relative_difference
end module reldiff_m