pub fn log_tensor_approx(f: &Tensor2) -> Tensor2
Approximate matrix logarithm of a tensor close to identity: ln(I + X) ≈ X - X²/2 + X³/3.
Only accurate for ||X|| < 0.5 (small deformation regime).