  We present a numerical method to evaluate worldline (WL) path integrals
defined on a curved Euclidean space, sampled with Monte Carlo (MC) techniques.
In particular, we adopt an algorithm known as YLOOPS with a slight modification
due to the introduction of a quadratic term which has the function of
stabilizing and speeding up the convergence. Our method, as the perturbative
counterparts, treats the non-trivial measure and deviation of the kinetic term
from flat, as interaction terms. Moreover, the numerical discretization adopted
in the present WLMC is realized with respect to the proper time of the
associated bosonic point-particle, hence such procedure may be seen as an
analogue of the time-slicing (TS) discretization already introduced to
construct quantum path integrals in curved space. As a result, a TS
counter-term is taken into account during the computation. The method is tested
against existing analytic calculations of the heat kernel for a free bosonic
point-particle in a D-dimensional maximally symmetric space.
