Given any f, a locally finitely piecewise affine homeomorphism of \(\Omega \subset \mathbb {R}^d\) onto \(\Delta \subset \mathbb {R}^d\) (for \(d=3, 4\) ) such that \(f\in W^{1,p}(\Omega , \mathbb {R}^d)\) and \(f^{-1}\in W^{1,q}(\Delta , \mathbb {R}^d)\) , \(1\le p,q < \infty \) and any \(\varepsilon >0\) we construct a diffeomorphism \(\tilde{f}\) such that \(\begin{aligned} \Vert f-\tilde{f}\Vert _{W^{1,p}(\Omega ,\mathbb {R}^d)} + \Vert f^{-1}-\tilde{f}^{-1}\Vert _{W^{1,q}(\Delta ,\mathbb {R}^d)} < \varepsilon . \end{aligned}\)