  use, intrinsic :: ieee_arithmetic
  real :: x
  x = ieee_value(0., ieee_positive_inf)
  print *, ieee_rem(x, 1.)
  print *, ieee_rem(ieee_value(0., ieee_positive_inf), 1.)
  end
