Zangwill says (Zangwill, Modern Electrodynamics, 2012 edition, page 40, undersection 2.3.1, "Lorentz Averaging")
...An example is the force on the charge density $\rho$ and current density $j$ in a volume $V$ due to electromagnetic fields $E$ and $B$ \begin{equation}F=\iiint (\rho E + j\times B)d\tau \end{equation}(2.40) Interpreted as a microscopic formula, the direct substitution of $\rho$ and $j$ confirms that (2.40) reproduces the Lorentz force law (2.1) for each microscopic particle. However, it is generally the case that $〈\rho E〉 \neq〈\rho 〉〈E〉$ and $〈j × B〉 \neq 〈j〉 × 〈B〉$. How, then, shall we compute the force on a macroscopic body? The answer, not often stated explicitly, is that we simply assume that (2.40) remains valid when all the variables are interpreted macroscopically. No ambiguities arise as long as F is the total force on an isolated sample of macroscopic matter in vacuum.
Let's say we have a point charge $q$ , and I do the spatial averaging of it, and let the resulting macroscopic charge density be $\rho_{macro}$. If I understand Zangwill correctly, the force on this point charge must be equal to $\iiint \rho_{macro}E_{macro}d\tau$ (where $E_{macro}$ is the macroscopic averaged electric field). Is this the correct interpretation of what he has said?
I understand that the force on the point particle is $qE_{micro}$ where $E_{micro}$ is the "true" microscopic electric field. I wonder whether this is equal to the above mentioned integral.
(We are only dealing with electrostatics so we are only bothered with the electric fields)