In this paper we analyse the structure of approximate solutions to the compatible two well problem with the constraint that the surface energy of the solution is less than some fixed constant. We prove a quantitative estimate that can be seen as a two well analogue of the Liouville theorem of Friesecke James Müller.
Let $H=\bigl(\begin{smallmatrix} \sigma& 0 0 & \sigma^{-1} \end{smallmatrix}\bigr)$ for $\sigma>0$. Let $0<\zeta_1<1<\zeta_2<\infty$. Let $K:=SO\left(2\right)\cup SO\left(2\right)H$. Let $u\in W^{2,1}\left(Q_{1}\left(0\right)\right)$ be a $\xCone$ invertible bilipschitz function with $\mathrm{Lip}\left(u\right)<\zeta_2$, $\mathrm{Lip}\left(u^{-1}\right)<\zeta_1^{-1}$. 
There exists positive constants $\mathfrak{c}_1<1$ and $\mathfrak{c}_2>1$ depending only on σ, $\zeta_1$, $\zeta_2$ such that if $\epsilon\in\left(0,\mathfrak{c}_1\right)$ and u satisfies the following inequalities \[ \int_{Q_{1}\left(0\right)} {\rm d}\left(Du\left(z\right),K\right) {\rm d}L^2 z\leq \epsilon \]\[ \int_{Q_{1}\left(0\right)} \left|D^2 u\left(z\right)\right| {\rm d}L^2 z\leq \mathfrak{c}_1, \] then there exists $J\in\left\{Id,H\right\}$ and $R\in SO\left(2\right)$ such that \[ \int_{Q_{\mathfrak{c}_1}\left(0\right)} \left|Du\left(z\right)-RJ\right| {\rm d}L^2 z\leq \mathfrak{c}_2\epsilon^{\frac{1}{800}}. \]