Function max_star 

pub fn max_star(a: f64, b: f64, jacobian: Jacobian) -> f64

The max* operation ln(e^a + e^b) of the log-MAP algorithms; max(a, b) with Jacobian::MaxLog or when the maximum is not finite.

Examples

use comnum_math::{Jacobian, max_star};
assert!((max_star(0.0, 0.0, Jacobian::Exact) - 2f64.ln()).abs() < 1e-15);
assert_eq!(max_star(0.0, -1.0, Jacobian::MaxLog), 0.0);