Smoothness of interpolatory multivariate subdivision schemes in Lie groups