Consider a proper $C^1$-function $f:\mathbb{R}^n\rightarrow \mathbb{R}^m$ where $n\geq m$.
If $c\in \mathbb{R}^m$ is a regular value of $f$, then we know that $f^{-1}(c)$ becomes a compact submanifold of $\mathbb{R}^n$.
For a sufficiently small $\epsilon$, consider a $C^1$-perturbation $g:\mathbb{R}^n\rightarrow \mathbb{R}^m$ of $f$; that is, $$ \|f-g\|<\epsilon \quad\text{ and }\quad \|Df-Dg\|<\epsilon.$$
My question is, is $f^{-1}(c)$ is diffeomorphic to $g^{-1}(c)$ for sufficienly small $\epsilon$?
If not, what additional conditions are needed for the diffeomorphism?
In fact, my interests lie in functions like $(x_1,x_2,…,x_d)\mapsto x^2_1+⋯+x^2_d$ or $(x_1,x_2)\mapsto(x^2_1−x^2_2,x_1x_2)$. So, maybe we can add some ad hoc assumptions to avoid pathological examples.