This optimizes the complex dilogarithm Li2 using targeted mappings
(reflection/inversion), a real-axis fast path, and a Bernoulli-in-x
expansion (x = -log(1 - z)), inspired by this Rust implementation[1]
(and references therein). The exported API remains unchanged, and
runtime improvements are measured by an added BenchmarkLi2.
Performance improvements range from a few percent to large factors,
depending on the region of the complex plane.
[1]https://github.com/Expander/polylog.rs