Is there an easy way to carry out multivariate power series reversion in Sage?
I have an automorphism $f$ of a power series ring $\mathbb{Z}[[x_1, \dots ,x_n]]$ specified by the n-tuple of power series $f(x_1), f(x_2), \dots , f(x_n)$. I would like to calculate the inverse of this automorphism. This is the multivariate case of power series reversion. Univariate power series reversion is implemented in Sage, but to my knowledge, the multivariate case is not implemented. The answer to this question from 7 years ago indicates that multivariate power series reversion was not yet implemented in Sage at that time.
As far as I know, my options are to either implement the multivariate version of Lagrange inversion in Sage, or try to find some other computer algebra software that might already have this feature. I hope I am wrong, though. Are there any plans to add multivariate power series reversion to Sage? Is there any test code somewhere that might already do this?
Thanks!