cwhmisc (version 3.0)

div.prot: Protected division

Description

1/x, but 1/0 -> 0

Usage

div.prot(x, eps = .Machine$double.eps)

Arguments

x

array to invert.

eps

limit: if abs(x) < eps then 1/x -> 0.

Examples

div.prot(c(0,.Machine$double.eps,1.0e-10,1,1.0e100))
  div.prot(c(0,.Machine$double.eps,1.0e-10,1,1.0e100), 100*.Machine$double.eps)