cwhmisc: div.prot – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

div.prot

Protected division

num/den, but num/0 -> .Machine$double.xmax^(3/4)

div.prot(num,den)

den,num

real, numerator and denominator

num/den, if is.infinite(num/den) then .Machine$double.xmax^(3/4), the ^(3/4) for getting something well below Inf.

Christian W. Hoffmann <christian@echoffmann.ch>

de <- .Machine$double.eps
  v<-c(0,de/c(1,2,4,8))
  div.prot(1,v)
 #  1.55252e+231  4.50360e+15  9.00720e+15  1.80144e+16  3.60288e+16

Miscellaneous Functions for Math, Plotting, Printing, Statistics, Strings, and Tools

v6.6

GPL (>= 2)

Authors

Christian W. Hoffmann

Initial release

2018-08-24, 10:40:10