To sample points from isotropic power-law kernels
*f_2(s) = f(||s||)* such as `siaf.powerlaw`

on a
bounded domain (i.e., *||s|| < \code{ub}*), it is
convenient to switch to polar coordinates *(r,θ)*,
which have a density proportional to
*r f_2((r \cos(θ), r \sin(θ))) = r f(r)*
(independent of the angle *θ* due to isotropy).
The angle is thus simply drawn uniformly in *[0,2π)*, and
*r* can be sampled by the inversion method, where numeric root
finding is used for the quantiles (since the quantile function is not
available in closed form).

1 | ```
siaf.simulatePC(intrfr)
``` |

`intrfr` |
a function computing the integral of |

a function with arguments `(n, siafpars, type, ub)`

, which
samples `n`

points from the spatial kernel *f_2(s)* within the
disc of radius `ub`

, where `siafpars`

and `type`

are
passed as second and third argument to `intrfr`

.
The environment of the returned function will be the caller's environment.

Sebastian Meyer

1 2 3 4 5 |

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.