public static double LogDiff(double lna, double lnc) { if (lna > lnc) return lna + Math.Exp(1.0 - Math.Exp(lnc - lna)); return Double.NegativeInfinity; }