##// END OF EJS Templates
ancestors: simplify symmetric difference...
Matt Mackall -
r6427:6b704ef9 default
parent child Browse files
Show More
@@ -1,134 +1,130 b''
1 # ancestor.py - generic DAG ancestor algorithm for mercurial
1 # ancestor.py - generic DAG ancestor algorithm for mercurial
2 #
2 #
3 # Copyright 2006 Matt Mackall <mpm@selenic.com>
3 # Copyright 2006 Matt Mackall <mpm@selenic.com>
4 #
4 #
5 # This software may be used and distributed according to the terms
5 # This software may be used and distributed according to the terms
6 # of the GNU General Public License, incorporated herein by reference.
6 # of the GNU General Public License, incorporated herein by reference.
7
7
8 import heapq
8 import heapq
9
9
10 def ancestor(a, b, pfunc):
10 def ancestor(a, b, pfunc):
11 """
11 """
12 return the least common ancestor of nodes a and b or None if there
12 return the least common ancestor of nodes a and b or None if there
13 is no such ancestor.
13 is no such ancestor.
14
14
15 pfunc must return a list of parent vertices
15 pfunc must return a list of parent vertices
16 """
16 """
17
17
18 if a == b:
18 if a == b:
19 return a
19 return a
20
20
21 # find depth from root of all ancestors
21 # find depth from root of all ancestors
22 visit = [a, b]
22 visit = [a, b]
23 depth = {}
23 depth = {}
24 while visit:
24 while visit:
25 vertex = visit[-1]
25 vertex = visit[-1]
26 pl = pfunc(vertex)
26 pl = pfunc(vertex)
27 if not pl:
27 if not pl:
28 depth[vertex] = 0
28 depth[vertex] = 0
29 visit.pop()
29 visit.pop()
30 else:
30 else:
31 for p in pl:
31 for p in pl:
32 if p == a or p == b: # did we find a or b as a parent?
32 if p == a or p == b: # did we find a or b as a parent?
33 return p # we're done
33 return p # we're done
34 if p not in depth:
34 if p not in depth:
35 visit.append(p)
35 visit.append(p)
36 if visit[-1] == vertex:
36 if visit[-1] == vertex:
37 depth[vertex] = min([depth[p] for p in pl]) - 1
37 depth[vertex] = min([depth[p] for p in pl]) - 1
38 visit.pop()
38 visit.pop()
39
39
40 # traverse ancestors in order of decreasing distance from root
40 # traverse ancestors in order of decreasing distance from root
41 def ancestors(vertex):
41 def ancestors(vertex):
42 h = [(depth[vertex], vertex)]
42 h = [(depth[vertex], vertex)]
43 seen = {}
43 seen = {}
44 while h:
44 while h:
45 d, n = heapq.heappop(h)
45 d, n = heapq.heappop(h)
46 if n not in seen:
46 if n not in seen:
47 seen[n] = 1
47 seen[n] = 1
48 yield (d, n)
48 yield (d, n)
49 for p in pfunc(n):
49 for p in pfunc(n):
50 heapq.heappush(h, (depth[p], p))
50 heapq.heappush(h, (depth[p], p))
51
51
52 def generations(vertex):
52 def generations(vertex):
53 sg, s = None, {}
53 sg, s = None, {}
54 for g, v in ancestors(vertex):
54 for g, v in ancestors(vertex):
55 if g != sg:
55 if g != sg:
56 if sg:
56 if sg:
57 yield sg, s
57 yield sg, s
58 sg, s = g, {v:1}
58 sg, s = g, {v:1}
59 else:
59 else:
60 s[v] = 1
60 s[v] = 1
61 yield sg, s
61 yield sg, s
62
62
63 x = generations(a)
63 x = generations(a)
64 y = generations(b)
64 y = generations(b)
65 gx = x.next()
65 gx = x.next()
66 gy = y.next()
66 gy = y.next()
67
67
68 # increment each ancestor list until it is closer to root than
68 # increment each ancestor list until it is closer to root than
69 # the other, or they match
69 # the other, or they match
70 try:
70 try:
71 while 1:
71 while 1:
72 if gx[0] == gy[0]:
72 if gx[0] == gy[0]:
73 for v in gx[1]:
73 for v in gx[1]:
74 if v in gy[1]:
74 if v in gy[1]:
75 return v
75 return v
76 gy = y.next()
76 gy = y.next()
77 gx = x.next()
77 gx = x.next()
78 elif gx[0] > gy[0]:
78 elif gx[0] > gy[0]:
79 gy = y.next()
79 gy = y.next()
80 else:
80 else:
81 gx = x.next()
81 gx = x.next()
82 except StopIteration:
82 except StopIteration:
83 return None
83 return None
84
84
85 def symmetricdifference(a, b, pfunc):
85 def symmetricdifference(a, b, pfunc):
86 """symmetric difference of the sets of ancestors of a and b
86 """symmetric difference of the sets of ancestors of a and b
87
87
88 I.e. revisions that are ancestors of a or b, but not both.
88 I.e. revisions that are ancestors of a or b, but not both.
89 """
89 """
90 # basic idea:
90 # basic idea:
91 # - mark a and b with different colors
91 # - mark a and b with different colors
92 # - walk the graph in topological order with the help of a heap;
92 # - walk the graph in topological order with the help of a heap;
93 # for each revision r:
93 # for each revision r:
94 # - if r has only one color, we want to return it
94 # - if r has only one color, we want to return it
95 # - add colors[r] to its parents
95 # - add colors[r] to its parents
96 #
96 #
97 # We keep track of the number of revisions in the heap that
97 # We keep track of the number of revisions in the heap that
98 # we may be interested in. We stop walking the graph as soon
98 # we may be interested in. We stop walking the graph as soon
99 # as this number reaches 0.
99 # as this number reaches 0.
100 if a == b:
100 if a == b:
101 return [a]
101 return [a]
102
102
103 WHITE = 1
103 WHITE = 1
104 BLACK = 2
104 BLACK = 2
105 ALLCOLORS = WHITE | BLACK
105 ALLCOLORS = WHITE | BLACK
106 colors = {a: WHITE, b: BLACK}
106 colors = {a: WHITE, b: BLACK}
107
107
108 visit = [-a, -b]
108 visit = [-a, -b]
109 heapq.heapify(visit)
109 heapq.heapify(visit)
110 n_wanted = len(visit)
110 interesting = len(visit)
111 ret = []
112
111
113 while n_wanted:
112 while interesting:
114 r = -heapq.heappop(visit)
113 r = -heapq.heappop(visit)
115 wanted = colors[r] != ALLCOLORS
114 if colors[r] != ALLCOLORS:
116 n_wanted -= wanted
115 interesting -= 1
117 if wanted:
118 ret.append(r)
119
116
120 for p in pfunc(r):
117 for p in pfunc(r):
121 if p not in colors:
118 if p not in colors:
122 # first time we see p; add it to visit
119 # first time we see p; add it to visit
123 n_wanted += wanted
124 colors[p] = colors[r]
120 colors[p] = colors[r]
121 if colors[p] != ALLCOLORS:
122 interesting += 1
125 heapq.heappush(visit, -p)
123 heapq.heappush(visit, -p)
126 elif colors[p] != ALLCOLORS and colors[p] != colors[r]:
124 elif colors[p] != ALLCOLORS and colors[p] != colors[r]:
127 # at first we thought we wanted p, but now
125 # at first we thought we wanted p, but now
128 # we know we don't really want it
126 # we know we don't really want it
129 n_wanted -= 1
130 colors[p] |= colors[r]
127 colors[p] |= colors[r]
128 interesting -= 1
131
129
132 del colors[r]
130 return [r for r in colors if colors[r] != ALLCOLORS]
133
134 return ret
General Comments 0
You need to be logged in to leave comments. Login now