1 Olivier Devillers
Walking in Poisson Delaunay triangulations
[D. & Hemsley, 2016] [Chenavier & D.,2018] [de Castro & D., 2018] [D. & Noizet, 2018]
Walking in Poisson Delaunay triangulations Olivier Devillers [D. - - PowerPoint PPT Presentation
Walking in Poisson Delaunay triangulations Olivier Devillers [D. & Hemsley, 2016] [Chenavier & D.,2018] [de Castro & D., 2018] [D. & Noizet, 2018] 1 Walking in Delaunay triangulations Straight walk 2 - 1 Walking in Delaunay
1 Olivier Devillers
[D. & Hemsley, 2016] [Chenavier & D.,2018] [de Castro & D., 2018] [D. & Noizet, 2018]
2 - 1 Straight walk
2 - 2 Straight walk
2 - 3 Straight walk
2 - 4 Straight walk
2 - 5 Straight walk
2 - 6 Straight walk
2 - 7 Straight walk Exit edge ? One orientation predicate
2 - 8 Straight walk End of walk ? A second orientation predicate
2 - 9 Straight walk Two orientation predicates per edge
3 - 1 Visibility walk
3 - 2 Visibility walk
3 - 3 Visibility walk
3 - 4 Visibility walk
3 - 5 Visibility walk
3 - 6 Visibility walk
3 - 7 Visibility walk Triangle with two exits One orientation predicate
3 - 8 Visibility walk Triangle with one exit 1.5 orientation predicate
One predicate Two predicates if this neighbor tried first if this neighbor tried first
3 - 9 Visibility walk 1.25 orientation predicate per edge ?
4 - 1
Straight walk Visibility walk How many edges crossed ?
4 - 2
Straight walk Visibility walk How many edges crossed ? Worst case in a triangulation (non Delaunay) 2n
4 - 3
Straight walk Visibility walk How many edges crossed ? Worst case in a triangulation (non Delaunay) 2n ∞
4 - 4
4 - 5 May loop
4 - 6
Straight walk Visibility walk How many edges crossed ? Worst case in a triangulation (non Delaunay) 2n ∞ ≥ 2
3
√n randomized
4 - 7 p p random choice
[D., Pion, & Teillaud 2002]
4 - 8 p p random choice
[D., Pion, & Teillaud 2002]
4 - 9
Straight walk Visibility walk How many edges crossed ? Worst case in a Delaunay triangulation 2n 2n
4 - 10 Not Delaunay May loop
4 - 11 Green power Red power < Power decreases
4 - 12
Straight walk Visibility walk How many edges crossed ? Worst case in a Delaunay triangulation 2n 2n O(√n) random
[Devroye, Lemaire, & Moreau, 2004]
64 3π2
√n + O( 1
√n) ≃ 2.16√n
O(√n)
uniform in domain Stretch in infinite Poisson distribution Stretch in infinite Poisson distribution
[D. & Hemsley, 2016]
5 - 1 s t Walk between vertices
5 - 2 s t Walk between vertices
5 - 3
Shortest path
s t Walk between vertices
5 - 4
Upper path
s t Walk between vertices
5 - 5
Compass walk
s t Walk between vertices
5 - 6
Voronoi path
s t Walk between vertices
5 - 7
Voronoi path with shortcuts
s t Walk between vertices
5 - 8
Shortest path Upper path Voronoi path with shortcuts Compass walk
s t Walk between vertices
6 - 1 Walk between vertices, worst case
Shortest path
6 - 2 Walk between vertices, worst case
Shortest path
Search this ”subgraph”
6 - 3 Walk between vertices, worst case
Shortest path
Search this ”subgraph”
6 - 4 Walk between vertices, worst case
Shortest path
Search this ”subgraph”
6 - 5 Walk between vertices, worst case
Shortest path
Search this ”subgraph” Dobkin, Friedman, and Supowit 1987 Keil and Gutwin 1989 Xia 2011 5.08
1+ √ 5 2
π 2.42 1.998
2π 3 cos(π/6)
Upper bound
6 - 6 Walk between vertices, worst case
Shortest path
Search this ”subgraph” Chew 1989
Bose, Devroye, L¨
Xia & Zhang 2011 1.5708
π 2
1.5846 1.5932 Lower bound
7 - 1 Walk between vertices, worst case
Voronoi path
Unbounded
7 - 2 Walk between vertices, worst case
Voronoi path
Unbounded
7 - 3 Walk between vertices, worst case
Voronoi path
Unbounded
8 - 1 Walk between vertices, worst case
Upper path
Unbounded
8 - 2 Walk between vertices, worst case
Upper path
Unbounded
8 - 3 Walk between vertices, worst case
Upper path
Unbounded
9 - 1 Walk between vertices, worst case
Compass walk
Unbounded
[Bose & Morin 2004]
9 - 2 Walk between vertices, worst case
Compass walk
Unbounded
[Bose & Morin 2004]
α α
9 - 3 Walk between vertices, worst case
Compass walk
Unbounded
[Bose & Morin 2004]
α α −ǫ −2ǫ
9 - 4 Walk between vertices, worst case
Compass walk
Unbounded
[Bose & Morin 2004]
α α −ǫ −2ǫ
10 Walk between vertices
Upper path Voronoi path shortcuts Shortest path Compass walk
11 - 1
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27
11 - 2
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
numerical integration
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ] [ D . & N
z e t , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
11 - 3
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
numerical integration
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ] [ D . & N
z e t , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
12 - 1
E [length] = E
n
1[triangle is Delaunay]1[first edge above st]length(first edge)
Poisson Delaunay triangulation, rate n
12 - 2
E [length] = E
n
1[triangle is Delaunay]1[first edge above st]length(first edge)
Poisson Delaunay triangulation, rate n
= n3
Slivnyak-Mecke
s t
12 - 3
E [length] = E
n
1[triangle is Delaunay]1[first edge above st]length(first edge)
Poisson Delaunay triangulation, rate n
= n3
= n3 ∞
r=0
1
xz=0
r
yz=−r
2 r32A(triangle)dα1:3dyzdxzdr
Blaschke-Petkantschin
s t α1 α2 α3 r z
12 - 4
E [length] = E
n
1[triangle is Delaunay]1[first edge above st]length(first edge)
Poisson Delaunay triangulation, rate n
= n3
= n3 ∞
r=0
1
xz=0
r
yz=−r
2 r32A(triangle)dα1:3dyzdxzdr = 4n3 ∞
r=0
e−nπr2r5dr
1
h= yz
r =−1
2 A(triangle)dα1:3dh
12 - 5
E [length] = E
n
1[triangle is Delaunay]1[first edge above st]length(first edge)
Poisson Delaunay triangulation, rate n
= n3
= n3 ∞
r=0
1
xz=0
r
yz=−r
2 r32A(triangle)dα1:3dyzdxzdr = 4n3 · 1 π3n3 · 35π 12 = 35 3π2 = 4n3 ∞
r=0
e−nπr2r5dr
1
h= yz
r =−1
2 A(triangle)dα1:3dh
t 1 h α1 α2 α3 arcsin h
13
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
numerical integration
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ D . & N
z e t , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
14 - 1
Shortest path
E [length(shortest path)]
14 - 2
Shortest path
E [length(shortest path)] limit density → ∞
14 - 3
Shortest path
E [length(shortest path)] limit density → ∞ by subadditivity remove start and target from point set look at shortest path between closest neighbor = between path from start to target negligible
15 Bad edge = almost horizontal edge P [bad] = small constant difficult dependencies to handle Many bad edges ⇐ length close to 1
16 - 1 Bad cell = cell with enough bad edges Many bad cells ⇐ length close to 1 Still dependencies Make a grid
16 - 2 Split cells in 4 families Many bad cells in 1 family ⇐ length close to 1 No more dependencies Make a grid
16 - 3 Make a grid If E [♯ ∈ cell] = constant ♯ possible paths = 4n (too big)
16 - 4 Make a grid If E [♯ ∈ cell] = constant ♯ possible paths = 4n (too big) big cells √n × √n 4
√n × n
17 - 1 A good cell ?
No Delaunay circles go outside
17 - 2 A good cell ?
No Delaunay circles go outside No short path from left to right
17 - 3 A good cell ?
No Delaunay circles go outside No short path from left to right Not enough edges to make a short path
17 - 4 A good cell ?
No Delaunay circles go outside No short path from left to right Not enough edges to make a short path Choose ♯ points in cell, Choose what ”short path” means
P
≤ O
√n
1 + 10−10
18
Shortest path Upper path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
numerical integration
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ D . & N
z e t , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
Voronoi path in higher dimension
[de Castro & D., 2018]
19 - 1
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
s t
19 - 2
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
Voronoi path in higher dimension
19 - 3
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
E [ℓ(V PX)] = 1
2
∞
−∞
∞
·ru1u2| det(JΦ)|dα1,1 . . . dα1,d−1dα2,1 . . . dα2,d−1drdx
Integral form
19 - 4
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
E [ℓ(V PX)] = 1
2
∞
−∞
∞
·ru1u2| det(JΦ)|dα1,1 . . . dα1,d−1dα2,1 . . . dα2,d−1drdx
Integral form Use Taylor expansion to be able to integrate
19 - 5
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
Γ d
2
4 24d−5d π2(2d − 2)!
d − 1 4d2 − 1 √ 2 ≤ E [ℓ(V PX)] ≤ Γ d
2
4 24d−5d π2(2d − 2)!
1 4d − 2 √ 2
19 - 6
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
asymptotic behavior
− 1 4 √ 2dπ + O(d
3 2 )
+ 3 4 √ 2dπ + O(d
3 2 )
between
19 - 7
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
asymptotic behavior
19 - 8
Voronoi path
1.27
4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
d k lower bound correct upper bound ≃ value ≃ 3 41
788984278470257640690697143 745000536337515228912680960
√ 2 1.49770 1.500
4523370364712510658076963509 4264485828690604413776035840
√ 2 1.50007 4 7
102494570 8729721 √ 2 π2
1.6823 1.698
121774997 10270260 √ 2 π2
1.6990 5 3
135 104
√ 2 1.8357 1.875
21305 16016
√ 2 1.8812 6 1
3014656 225225 √ 2 π2
1.9179 2.04
753664 51975 √ 2 π2
2.0778 7 1
210 143
√ 2 2.0768 2.2
225 143
√ 2 2.2252 8 1
2080374784 134008875 √ 2 π2
2.2244 2.3
130023424 7882875 √ 2 π2
2.3635
numerical integration
20
[ D . & N
z e t , 2 1 8 ]
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
numerical integration
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ] [ O n g
n g w
k ]
21 - 1
Compass walk
Complicated dependencies Stretch factor
21 - 2
Compass walk
Complicated dependencies Stretch factor = ℓ1 + ℓ2 + . . . ℓk x1 + x2 + . . . xk edge lengths length of horizontal projections
21 - 3
Compass walk
Complicated dependencies Stretch factor = ℓ1 + ℓ2 + . . . ℓk x1 + x2 + . . . xk If independent: E [ℓ] E [x] = 175π 512 ≃ 1.0738
21 - 4
Compass walk
Complicated dependencies Stretch factor = ℓ1 + ℓ2 + . . . ℓk x1 + x2 + . . . xk If independent: E [ℓ] E [x] = 175π 512 ≃ 1.0738 Experimental: There is a positive bias 1.06777
21 - 5
Compass walk
Complicated dependencies Stretch factor = ℓ1 + ℓ2 + . . . ℓk x1 + x2 + . . . xk If independent: E [ℓ] E [x] = 175π 512 ≃ 1.0738 Experimental: There is a positive bias 1.06777 Experimental:
E [ℓ] E [x] ≃ 1.061 0.988 ≃ 1.0738
21 - 6
Compass walk
Complicated dependencies Stretch factor = ℓ1 + ℓ2 + . . . ℓk x1 + x2 + . . . xk If independent: E [ℓ] E [x] = 175π 512 ≃ 1.0738 Experimental: There is a positive bias 1.06777 Experimental:
E [ℓ] E [x] ≃ 1.061 0.988 ≃ 1.0738 E [ℓ] E [x] ≃ 1.116 1.045 ≃ 1.0681 E [ℓ] E [x] ≃ 1.117 1.046 ≃ 1.0678
at depth 1 at depth 2
22 - 1
[ D . & N
z e t , 2 1 8 ]
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
22 - 2
[ D . & N
z e t , 2 1 8 ]
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
Locally defined path ”easy” to analyze
(computation may be difficult)
22 - 3
[ D . & N
z e t , 2 1 8 ]
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
Locally defined path ”easy” to analyze
(computation may be difficult)
Incrementally defined path dependency issues the tree can be analyzed
22 - 4
[ D . & N
z e t , 2 1 8 ]
Shortest path Upper path Voronoi path Shortened V. path Compass walk
Expected length (experiments) Euclidean length 1 1.04 1.07 1.16 1.18 1.27 theory ≥ 1 + 10−11 1.16
35 3π2 ≃ 1.18 4 π ≃ 1.27
[ B a c c e l l i e t a l . , 2 ]
[ C h e n a v i e r & D . , 2 1 8 ]
[ C h e n a v i e r & D . , 2 1 8 ]
Locally defined path ”easy” to analyze
(computation may be difficult)
Incrementally defined path dependency issues the tree can be analyzed Incrementally defined path also dependency issues no idea about tight bounds
23 Thank you