  In the study of femtoscopic correlations in high-energy physics, besides
Bose--Einstein correlations, one has to take final-state interactions into
account. Amongst them, Coulomb interactions play a prominent role in the case
of charged particles. Recent measurements have shown that in heavy-ion
collisions, Bose--Einstein correlations can be best described by L\'evy-type
sources instead of the more common Gaussian assumption. Furthermore,
three-dimensional measurements have indicated that, depending on the choice of
frame, a deviation from spherical symmetry observed under the assumption of
Gaussian source functions persists in the case of L\'evy-type sources. To
clarify such three-dimensional L\'evy-type correlation measurements, it is thus
important to study the effect of Coulomb interactions in the case of
non-spherical L\'evy sources. We calculated the Coulomb correction factor
numerically in the case of such a source function for assorted kinematic
domains and parameter values using the Metropolis--Hastings algorithm and
compared our results with previous methods to treat Coulomb interactions in the
presence of L\'evy sources.
