pub fn exp_tensor_pade(a: &Tensor2) -> Option<Tensor2>
Approximate matrix exponential via Padé (2,2) approximation.
exp(A) ≈ (I - A/2 + A²/12)^{-1} (I + A/2 + A²/12).