this is musing from an electrical engineer that knows something about discrete modeling of analog or continuous-time systems.
systems that are continuous are often described by continuous differential equations. if the diff eqs. are linear, there is a way (using Laplace Transform) to describe the system exactly and solve for a closed-form solution. but if the system has little non-linear components in them, sometimes the only way to understand the system is to simulate it with a discrete-time approximation. probably the simplest discrete-time approximation is
Euler's forward method . now when programming a computer to simulate a physical system will involve turning those differential equations into difference equations, using, say the Euler method. in doing so you will compute dimensionless constants that will serve as coefficients in the discrete-time simulation. the sampling period and other time parameters will play a role in these dimensionless coefficients.
similar to
cellular automata, you can take equations of physical interaction, like Maxwell's equations, and Schrodinger, and Einstein's field equation, and continuity equations for conserved quantities of physical stuff, and you can turn them into difference equations that will involve the sampling period (in the time dimension) of the simulation and cell size (in the 3 length dimensions).
now if you were to select a quantization unit in time and length that were essentially the
Planck units, but defined so that these constants are removed from the above equations of physical law (which means these constants are set to 1): 4 \pi G = c = \hbar = \epsilon_0 = 1, then when you turn the differential equations into simple difference equations (Eulers method is plenty adequate because the discrete time and cell width are virtually infinitesimally small). but because these constants of nature disappear from expressions of physical law, there are no arbitrary scaling constants in the discrete-time simulator. we don't have Nature taking this quantity of flux density and pulling this number (\frac{1}{\epsilon_0}) out of her butt and scaling it to convert it into field strength. their quantities are equal so then maybe
they, the two physical quantities, are the same thing. maybe flux density is the very same thing as field strength in Maxwells equations as well as the discrete-time simulation of it.
but that happens (ditch the arbitrary scaling constants) only if you run your discrete-time, discrete-space simulator with discrete units of the rationalized Planck time and Planck length (where 4 \pi G = c = \hbar = \epsilon_0 = 1 ).
because the Planck scale is soooo... ridiculously tiny, there is no way human beings can ever hope to measure anything around that scale to observe any of this discrete phenomena. way smaller than anything in the atomic or subatomic scale. so i like to pretend it could be true. i don't think there is any way to find out one way or another.