The Eikonal equation arises naturally in the limit of the second order Aviles-Giga functional whose $Gamma$-convergence is a long standing challenging problem. The theory of entropy solutions of the Eikonal equation plays a central role in the variational analysis of this problem. Establishing fine structures of entropy solutions of the Eikonal equation, e.g. concentration of entropy measures on $mathcal{H}^1$-rectifiable sets in $2$D, is arguably the key missing part for a proof of the full $Gamma$-convergence of the Aviles-Giga functional. In the first part of this work, for $pin left(1,frac{4}{3}right]$ we establish an $L^p$ version of the main theorem of Ghiraldin and Lamy [Comm. Pure Appl. Math. 73 (2020), no. 2, 317-349]. Specifically we show that if $m$ is a solution to the Eikonal equation, then $min B^{frac{1}{3}}_{3p,infty,loc}$ is equivalent to all entropy productions of $m$ being in $L^p_{loc}$. This result also shows that as a consequence of a weak form of the Aviles-Giga conjecture (namely the conjecture that all solutions to the Eikonal equation whose entropy productions are in $L^p_{loc}$ are rigid) - the rigidity/flexibility threshold of the Eikonal equation is exactly the space $ B^{frac{1}{3}}_{3,infty,loc}$. In the second part of this paper, under the assumption that all entropy productions are in $L^p_{loc}$, we establish a factorization formula for entropy productions of solutions of the Eikonal equation in terms of the two Jin-Kohn entropies. A consequence of this formula is control of all entropy productions by the Jin-Kohn entropies in the $L^p$ setting - this is a strong extension of an earlier result of the authors [Annales de lInstitut Henri Poincar{e}. Analyse Non Lin{e}aire 35 (2018), no. 2, 481-516].