pub fn geometric_factor(eta: f64, c: f64, mu: f64) -> f64
Expected one-step geometric factor 1 - 2 eta c(eta) mu from V(x_k) <= (1 - 2 eta c(eta) mu)^k V(x_0).
1 - 2 eta c(eta) mu
V(x_k) <= (1 - 2 eta c(eta) mu)^k V(x_0)