Comment échantillonner des points dans l'espace hyperbolique?

L'espace hyperbolique dans le modèle du demi-espace supérieur de Poincaré ressemble à $\Bbb R^n$ ordinaire mais avec la notion d'angle et de distance déformée d'une manière relativement simple. Dans l'espace euclidien, je peux échantillonner un point aléatoire uniformément dans une balle de plusieurs manières, par exemple en générant $n$ échantillons gaussiens indépendants pour obtenir une direction, et échantillonner séparément une coordonnée radiale $r$ en échantillonnant uniformément $s$ partir de $\left[0, \frac1{n+1}R^{n+1}\right]$, où$R$est le rayon et avec$r = \left((n+1)s\right)^{\frac1{n+1}}$ . Dans le demi-plan supérieur hyperbolique, une sphère est toujours une sphère, seul son centre ne sera pas le centre de la métrique euclidienne, nous pourrions donc faire de même.

Si nous voulons échantillonner selon une distribution non uniforme, mais toujours de manière isotrope, par exemple une distribution gaussienne, cela ne semble pas si facile. Dans l'espace euclidien, nous pourrions simplement générer un échantillon gaussien pour chaque coordonnée (cela ne fonctionne que pour la distribution gaussienne), ou générer de manière équivalente un échantillon gaussien multidimensionnel. Existe-t-il un moyen direct de convertir cet échantillon en échantillon dans l'espace hyperbolique?

Une approche alternative pourrait être de générer d'abord une direction distribuée uniformément (par exemple à partir de $n$ échantillons gaussiens) puis un échantillon gaussien pour la composante radiale, et enfin de générer l'image sous la carte exponentielle dans la direction spécifiée pour la longueur spécifiée. Une variante serait de simplement prendre l'échantillon gaussien euclidien et de le cartographier sous la carte exponentielle.

Mes questions:

• quelle serait une bonne et efficace façon d'obtenir un échantillon gaussien avec une moyenne et un écart-type donnés dans l'espace hyperbolique?
• les méthodes que je décris ci-dessus fournissent-elles l'échantillonnage souhaité?
• quelqu'un a-t-il déjà élaboré la formule?
• comment cela se généralise-t-il à d'autres métriques et à d'autres distributions de probabilités?

Merci d'avance.

ÉDITER

Je viens de réaliser que même dans le cas d'un échantillonnage uniforme, ces questions demeurent; même si une sphère est une sphère, une distribution uniforme ne serait pas décrite par une fonction constante sur une balle.

Je suis en train de faire ça pour moi. Je pense que l'analogue le plus approprié au gaussien serait le noyau de chaleur dans l'espace hyperbolique. Heureusement, cela a déjà été compris: https://www.math.uni-bielefeld.de/~grigor/nog.pdf (également disponible dans un Bulletin de la London Mathematical Society ).

$e^{-dist^2/constant}$

Voici un échantillon uniforme pour la boule de rayon 3 centrée à l'origine:

Si vous le souhaitez, je serais heureux d'en dire plus. Je pensais juste que je mettrais cela en place, car il y avait clairement un certain intérêt pour cela, au moins dans le passé.

La constante pi n'est qu'une constante dans l'espace euclidien. La valeur de pi est différente dans d'autres géométries. Le paramètre pi modifie la masse de probabilité sous la gaussienne. Le paramètre pi est utilisé pour normaliser les probabilités. Je commence tout juste à étudier cela.

J'ai conclu il y a quelque temps que l'espace passe d'hyperbolique à euclidien à sphérique à mesure que le nombre de sigmas augmente. J'étais heureux de parcourir une discussion des cercles dans chaque espace et pi en fonction des espaces Lp via le paramètre p.

$n$