Function dinv 

fn dinv(d: f64) -> f64

The initial guess’s inverse kernel: 1/D where |D| ≥ 0.065, else ±1/0.065 by the sign of D (+ at D = 0).