| 577 |
|
where $\hat{f}_\textrm{M}$ is the unit vector pointing along the force |
| 578 |
|
vector computed using method M. Each of these $\theta$ values was |
| 579 |
|
accumulated in a distribution function and weighted by the area on the |
| 580 |
< |
unit sphere. Non-linear Gaussian fits were used to measure the width |
| 581 |
< |
of the resulting distributions. |
| 580 |
> |
unit sphere. Since this distribution is a measure of angular error |
| 581 |
> |
between two different electrostatic summation methods, there is no |
| 582 |
> |
{\it a priori} reason for the profile to adhere to any specific |
| 583 |
> |
shape. Thus, gaussian fits were used to measure the width of the |
| 584 |
> |
resulting distributions. |
| 585 |
|
% |
| 586 |
|
%\begin{figure} |
| 587 |
|
%\centering |
