      ; log(), with rxode2's floor of sqrt(DBL_EPSILON)
      IPRED = RX_PRED_F_
      IF (IPRED .LE. 1.4901161E-8) IPRED = 1.4901161E-8
      IPRED = DLOG(IPRED)
