]> git.donarmstrong.com Git - lilypond.git/blob - buildscripts/output-distance.py
remove debugging gob.
[lilypond.git] / buildscripts / output-distance.py
1 #!@TARGET_PYTHON@
2 import sys
3 import optparse
4 import os
5 import math
6
7 ## so we can call directly as buildscripts/output-distance.py
8 me_path = os.path.abspath (os.path.split (sys.argv[0])[0])
9 sys.path.insert (0, me_path + '/../python/')
10
11
12 import safeeval
13
14
15 X_AXIS = 0
16 Y_AXIS = 1
17 INFTY = 1e6
18
19 OUTPUT_EXPRESSION_PENALTY = 1
20 ORPHAN_GROB_PENALTY = 1
21 options = None
22
23 ################################################################
24 # system interface.
25 temp_dir = None
26 class TempDirectory:
27     def __init__ (self):
28         import tempfile
29         self.dir = tempfile.mkdtemp ()
30         print 'dir is', self.dir
31     def __del__ (self):
32         print 'rm -rf %s' % self.dir 
33         os.system ('rm -rf %s' % self.dir )
34     def __call__ (self):
35         return self.dir
36
37
38 def get_temp_dir  ():
39     global temp_dir
40     if not temp_dir:
41         temp_dir = TempDirectory ()
42     return temp_dir ()
43
44 def read_pipe (c):
45     print 'pipe' , c
46     return os.popen (c).read ()
47
48 def system (c):
49     print 'system' , c
50     s = os.system (c)
51     if s :
52         raise Exception ("failed")
53     return
54
55 def shorten_string (s):
56     threshold = 15 
57     if len (s) > 2*threshold:
58         s = s[:threshold] + '..' + s[-threshold:]
59     return s
60
61 def max_distance (x1, x2):
62     dist = 0.0
63
64     for (p,q) in zip (x1, x2):
65         dist = max (abs (p-q), dist)
66         
67     return dist
68
69
70 def compare_png_images (old, new, dest_dir):
71     def png_dims (f):
72         m = re.search ('([0-9]+) x ([0-9]+)', read_pipe ('file %s' % f))
73         
74         return tuple (map (int, m.groups ()))
75
76     dest = os.path.join (dest_dir, new.replace ('.png', '.compare.jpeg'))
77     try:
78         dims1 = png_dims (old)
79         dims2 = png_dims (new)
80     except AttributeError:
81         ## hmmm. what to do?
82         system ('touch %(dest)s' % locals ())
83         return
84     
85     dims = (min (dims1[0], dims2[0]),
86             min (dims1[1], dims2[1]))
87
88     dir = get_temp_dir ()
89     system ('convert -depth 8 -crop %dx%d+0+0 %s %s/crop1.png' % (dims + (old, dir)))
90     system ('convert -depth 8 -crop %dx%d+0+0 %s %s/crop2.png' % (dims + (new, dir)))
91
92     system ('compare -depth 8 %(dir)s/crop1.png %(dir)s/crop2.png %(dir)s/diff.png' % locals ())
93
94     system ("convert  -depth 8 %(dir)s/diff.png -blur 0x3 -negate -channel alpha,blue -type TrueColorMatte -fx 'intensity'    %(dir)s/matte.png" % locals ())
95
96     system ("composite -quality 65 %(dir)s/matte.png %(new)s %(dest)s" % locals ())
97
98
99 ################################################################
100 # interval/bbox arithmetic.
101
102 empty_interval = (INFTY, -INFTY)
103 empty_bbox = (empty_interval, empty_interval)
104
105 def interval_is_empty (i):
106     return i[0] > i[1]
107
108 def interval_length (i):
109     return max (i[1]-i[0], 0) 
110     
111 def interval_union (i1, i2):
112     return (min (i1[0], i2[0]),
113             max (i1[1], i2[1]))
114
115 def interval_intersect (i1, i2):
116     return (max (i1[0], i2[0]),
117             min (i1[1], i2[1]))
118
119 def bbox_is_empty (b):
120     return (interval_is_empty (b[0])
121             or interval_is_empty (b[1]))
122
123 def bbox_union (b1, b2):
124     return (interval_union (b1[X_AXIS], b2[X_AXIS]),
125             interval_union (b2[Y_AXIS], b2[Y_AXIS]))
126             
127 def bbox_intersection (b1, b2):
128     return (interval_intersect (b1[X_AXIS], b2[X_AXIS]),
129             interval_intersect (b2[Y_AXIS], b2[Y_AXIS]))
130
131 def bbox_area (b):
132     return interval_length (b[X_AXIS]) * interval_length (b[Y_AXIS])
133
134 def bbox_diameter (b):
135     return max (interval_length (b[X_AXIS]),
136                 interval_length (b[Y_AXIS]))
137                 
138
139 def difference_area (a, b):
140     return bbox_area (a) - bbox_area (bbox_intersection (a,b))
141
142 class GrobSignature:
143     def __init__ (self, exp_list):
144         (self.name, self.origin, bbox_x,
145          bbox_y, self.output_expression) = tuple (exp_list)
146         
147         self.bbox = (bbox_x, bbox_y)
148         self.centroid = (bbox_x[0] + bbox_x[1], bbox_y[0] + bbox_y[1])
149
150     def __repr__ (self):
151         return '%s: (%.2f,%.2f), (%.2f,%.2f)\n' % (self.name,
152                                                    self.bbox[0][0],
153                                                    self.bbox[0][1],
154                                                    self.bbox[1][0],
155                                                    self.bbox[1][1])
156                                                  
157     def axis_centroid (self, axis):
158         return apply (sum, self.bbox[axis])  / 2 
159     
160     def centroid_distance (self, other, scale):
161         return max_distance (self.centroid, other.centroid) / scale 
162         
163     def bbox_distance (self, other):
164         divisor = bbox_area (self.bbox) + bbox_area (other.bbox)
165
166         if divisor:
167             return (difference_area (self.bbox, other.bbox) +
168                     difference_area (other.bbox, self.bbox)) / divisor
169         else:
170             return 0.0
171         
172     def expression_distance (self, other):
173         if self.output_expression == other.output_expression:
174             return 0
175         else:
176             return 1
177
178 ################################################################
179 # single System.
180
181 class SystemSignature:
182     def __init__ (self, grob_sigs):
183         d = {}
184         for g in grob_sigs:
185             val = d.setdefault (g.name, [])
186             val += [g]
187
188         self.grob_dict = d
189         self.set_all_bbox (grob_sigs)
190
191     def set_all_bbox (self, grobs):
192         self.bbox = empty_bbox
193         for g in grobs:
194             self.bbox = bbox_union (g.bbox, self.bbox)
195
196     def closest (self, grob_name, centroid):
197         min_d = INFTY
198         min_g = None
199         try:
200             grobs = self.grob_dict[grob_name]
201
202             for g in grobs:
203                 d = max_distance (g.centroid, centroid)
204                 if d < min_d:
205                     min_d = d
206                     min_g = g
207
208
209             return min_g
210
211         except KeyError:
212             return None
213     def grobs (self):
214         return reduce (lambda x,y: x+y, self.grob_dict.values(), [])
215
216 ################################################################
217 ## comparison of systems.
218
219 class SystemLink:
220     def __init__ (self, system1, system2):
221         self.system1 = system1
222         self.system2 = system2
223         
224         self.link_list_dict = {}
225         self.back_link_dict = {}
226
227
228         ## pairs
229         self.orphans = []
230
231         ## pair -> distance
232         self.geo_distances = {}
233
234         ## pairs
235         self.expression_changed = []
236
237         self._geometric_distance = None
238         self._expression_change_count = None
239         self._orphan_count = None
240         
241         for g in system1.grobs ():
242
243             ## skip empty bboxes.
244             if bbox_is_empty (g.bbox):
245                 continue
246             
247             closest = system2.closest (g.name, g.centroid)
248             
249             self.link_list_dict.setdefault (closest, [])
250             self.link_list_dict[closest].append (g)
251             self.back_link_dict[g] = closest
252
253
254     def calc_geometric_distance (self):
255         total = 0.0
256         for (g1,g2) in self.back_link_dict.items ():
257             if g2:
258                 d = g1.bbox_distance (g2)
259                 if d:
260                     self.geo_distances[(g1,g2)] = d
261
262                 total += d
263
264         self._geometric_distance = total
265     
266     def calc_orphan_count (self):
267         count = 0
268         for (g1, g2) in self.back_link_dict.items ():
269             if g2 == None:
270                 self.orphans.append ((g1, None))
271                 
272                 count += 1
273
274         self._orphan_count = count
275     
276     def calc_output_exp_distance (self):
277         d = 0
278         for (g1,g2) in self.back_link_dict.items ():
279             if g2:
280                 d += g1.expression_distance (g2)
281
282         self._expression_change_count = d
283
284     def output_expression_details_string (self):
285         return ', '.join ([g1.name for g1 in self.expression_changed])
286     
287     def geo_details_string (self):
288         results = [(d, g1,g2) for ((g1, g2), d) in self.geo_distances.items()]
289         results.sort ()
290         results.reverse ()
291         
292         return ', '.join (['%s: %f' % (g1.name, d) for (d, g1, g2) in results])
293
294     def orphan_details_string (self):
295         return ', '.join (['%s-None' % g1.name for (g1,g2) in self.orphans if g2==None])
296
297     def geometric_distance (self):
298         if self._geometric_distance == None:
299             self.calc_geometric_distance ()
300         return self._geometric_distance
301     
302     def orphan_count (self):
303         if self._orphan_count == None:
304             self.calc_orphan_count ()
305             
306         return self._orphan_count
307     
308     def output_expression_change_count (self):
309         if self._expression_change_count == None:
310             self.calc_output_exp_distance ()
311         return self._expression_change_count
312         
313     def distance (self):
314         return (self.output_expression_change_count (),
315                 self.orphan_count (),
316                 self.geometric_distance ())
317     
318 def read_signature_file (name):
319     print 'reading', name
320     
321     entries = open (name).read ().split ('\n')
322     def string_to_tup (s):
323         return tuple (map (float, s.split (' '))) 
324
325     def string_to_entry (s):
326         fields = s.split('@')
327         fields[2] = string_to_tup (fields[2])
328         fields[3] = string_to_tup (fields[3])
329
330         return tuple (fields)
331     
332     entries = [string_to_entry (e) for e in entries
333                if e and not e.startswith ('#')]
334
335     grob_sigs = [GrobSignature (e) for e in entries]
336     sig = SystemSignature (grob_sigs)
337     return sig
338
339
340 ################################################################
341 # different systems of a .ly file.
342
343 hash_to_original_name = {}
344
345 class FileLink:
346     def __init__ (self, f1, f2):
347         self._distance = None
348         self.file_names = (f1, f2)
349         
350     def text_record_string (self):
351         return '%-30f %-20s\n' % (self.distance (),
352                                   self.name ())
353     def calc_distance (self):
354         return 0.0
355
356     def distance (self):
357         if self._distance == None:
358            self._distance = self.calc_distance ()
359
360         return self._distance
361     
362         
363     def name (self):
364         base = os.path.basename (self.file_names[1])
365         base = os.path.splitext (base)[0]
366         
367         base = hash_to_original_name.get (base, base)
368         base = os.path.splitext (base)[0]
369         return base
370     
371     def extension (self):
372         return os.path.splitext (self.file_names[1])[1]
373
374     def link_files_for_html (self, dest_dir):
375         for f in self.file_names:
376             link_file (f, os.path.join (dest_dir, f))
377
378     def get_distance_details (self):
379         return ''
380
381     def get_cell (self, oldnew):
382         return ''
383     
384     def get_file (self, oldnew):
385         return self.file_names[oldnew]
386     
387     def html_record_string (self, dest_dir):
388         self.link_files_for_html (dest_dir)
389         
390         dist = self.distance()
391         
392         details = self.get_distance_details ()
393         if details:
394             details_base = os.path.splitext (self.file_names[1])[0]
395             details_base += '.details.html'
396             fn = dest_dir + '/'  + details_base
397             open (fn, 'w').write (details)
398
399             details = '<br>(<a href="%(details_base)s">details</a>)' % locals ()
400
401         cell1 = self.get_cell (0)
402         cell2 = self.get_cell (1)
403
404         name = self.name () + self.extension ()
405         file1 = self.get_file (0)
406         file2 = self.get_file (1)
407         
408         return '''<tr>
409 <td>
410 %(dist)f
411 %(details)s
412 </td>
413 <td>%(cell1)s<br><font size=-2><a href="%(file1)s"><tt>%(name)s</tt></font></td>
414 <td>%(cell2)s<br><font size=-2><a href="%(file2)s"><tt>%(name)s</tt></font></td>
415 </tr>''' % locals ()
416
417
418 class FileCompareLink (FileLink):
419     def __init__ (self, f1, f2):
420         FileLink.__init__ (self, f1, f2)
421         self.contents = (self.get_content (self.file_names[0]),
422                          self.get_content (self.file_names[1]))
423         
424
425     def calc_distance (self):
426         ## todo: could use import MIDI to pinpoint
427         ## what & where changed.
428
429         if self.contents[0] == self.contents[1]:
430             return 0.0
431         else:
432             return 100.0;
433         
434     def get_content (self, f):
435         print 'reading', f
436         s = open (f).read ()
437         return s
438
439
440 class GitFileCompareLink (FileCompareLink):
441     def get_cell (self, oldnew):
442         return self.contents[oldnew]
443     
444     def calc_distance (self):
445         if self.contents[0] == self.contents[1]:
446             d = 0.0
447         else:
448             d = 1.0001 *options.threshold
449
450         return d
451         
452 class TextFileCompareLink (FileCompareLink):
453     def calc_distance (self):
454         import difflib
455         diff = difflib.unified_diff (self.contents[0].strip().split ('\n'),
456                                      self.contents[1].strip().split ('\n'),
457                                      fromfiledate = self.file_names[0],
458                                      tofiledate = self.file_names[1]
459                                      )
460         
461         self.diff_lines =  [l for l in diff]
462         self.diff_lines = self.diff_lines[2:]
463         
464         return float (len ([l for l in self.diff_lines if l[0] in '-+']))
465         
466     def get_cell (self, oldnew):
467         str = ''
468         if oldnew == 1:
469             str = '\n'.join ([d.replace ('\n','') for d in self.diff_lines])
470         str = '<font size="-2"><pre>%s</pre></font>' % str
471         return str
472
473         
474 class ProfileFileLink (FileCompareLink):
475     def __init__ (self, f1, f2):
476         FileCompareLink.__init__ (self, f1, f2)
477         self.results = [{}, {}]
478     
479     def get_cell (self, oldnew):
480         str = ''
481         for k in ('time', 'cells'):
482             if oldnew==0:
483                 str += '%-8s: %d\n' %  (k, int (self.results[oldnew][k]))
484             else:
485                 str += '%-8s: %8d (%5.3f)\n' % (k, int (self.results[oldnew][k]),
486                                          self.get_ratio (k))
487
488         return '<pre>%s</pre>' % str
489             
490     def get_ratio (self, key):
491         (v1,v2) = (self.results[0].get (key, -1),
492                    self.results[1].get (key, -1))
493
494         if v1 <= 0 or v2 <= 0:
495             return 0.0
496
497         return (v1 - v2) / float (v1+v2)
498     
499     def calc_distance (self):
500         for oldnew in (0,1):
501             def note_info (m):
502                 self.results[oldnew][m.group(1)] = float (m.group (2))
503             
504             re.sub ('([a-z]+): ([-0-9.]+)\n',
505                     note_info, self.contents[oldnew])
506
507         dist = 0.0
508         factor = {'time': 1.0 ,
509                   'cells': 10.0,
510                   }
511         
512         for k in ('time', 'cells'):
513             dist += math.exp (self.get_ratio (k) * factor[k]) - 1
514
515         dist = min (dist, 100)
516         return dist
517
518     
519 class MidiFileLink (FileCompareLink):
520     def get_content (self, f):
521         s = FileCompareLink.get_content (self, f)
522         s = re.sub ('LilyPond [0-9.]+', '', s)
523         return s
524     
525     def get_cell (self, oldnew):
526         str = ''
527         if oldnew == 1 and self.distance () > 0:
528             str = 'changed' 
529         return str
530     
531
532 class SignatureFileLink (FileLink):
533     def __init__ (self, f1, f2 ):
534         FileLink.__init__ (self, f1, f2)
535         self.system_links = {}
536
537     def add_system_link (self, link, number):
538         self.system_links[number] = link
539
540     def calc_distance (self):
541         d = 0.0
542
543         orphan_distance = 0.0
544         for l in self.system_links.values ():
545             d = max (d, l.geometric_distance ())
546             orphan_distance += l.orphan_count ()
547             
548         return d + orphan_distance
549
550     def source_file (self):
551         for ext in ('.ly', '.ly.txt'):
552             if os.path.exists (self.base_names[1] + ext):
553                 return self.base_names[1] + ext
554         return ''
555     
556     def add_file_compare (self, f1, f2):
557         system_index = [] 
558
559         def note_system_index (m):
560             system_index.append (int (m.group (1)))
561             return ''
562         
563         base1 = re.sub ("-([0-9]+).signature", note_system_index, f1)
564         base2 = re.sub ("-([0-9]+).signature", note_system_index, f2)
565
566         self.base_names = (os.path.normpath (base1),
567                            os.path.normpath (base2))
568
569         def note_original (match):
570             hash_to_original_name[os.path.basename (self.base_names[1])] = match.group (1)
571             return ''
572         
573         ## ugh: drop the .ly.txt
574         for ext in ('.ly', '.ly.txt'):
575             try:
576                 re.sub (r'\\sourcefilename "([^"]+)"',
577                         note_original, open (base1 + ext).read ())
578             except IOError:
579                 pass
580                 
581         s1 = read_signature_file (f1)
582         s2 = read_signature_file (f2)
583
584         link = SystemLink (s1, s2)
585
586         self.add_system_link (link, system_index[0])
587
588     
589     def create_images (self, dest_dir):
590
591         files_created = [[], []]
592         for oldnew in (0, 1):
593             pat = self.base_names[oldnew] + '.eps'
594
595             for f in glob.glob (pat):
596                 infile = f
597                 outfile = (dest_dir + '/' + f).replace ('.eps', '.png')
598
599                 mkdir (os.path.split (outfile)[0])
600                 cmd = ('gs -sDEVICE=png16m -dGraphicsAlphaBits=4 -dTextAlphaBits=4 '
601                        ' -r101 '
602                        ' -sOutputFile=%(outfile)s -dNOSAFER -dEPSCrop -q -dNOPAUSE '
603                        ' %(infile)s  -c quit '  % locals ())
604
605                 files_created[oldnew].append (outfile)
606                 system (cmd)
607
608         return files_created
609     
610     def link_files_for_html (self, dest_dir):
611         FileLink.link_files_for_html (self, dest_dir)
612         to_compare = [[], []]
613
614         exts = []
615         if options.create_images:
616             to_compare = self.create_images (dest_dir)
617         else:
618             exts += ['.png', '-page*png']
619         
620         for ext in exts:            
621             for oldnew in (0,1):
622                 for f in glob.glob (self.base_names[oldnew] + ext):
623                     dst = dest_dir + '/' + f
624                     link_file (f, dst)
625
626                     if f.endswith ('.png'):
627                         to_compare[oldnew].append (f)
628                         
629         if options.compare_images:                
630             for (old, new) in zip (to_compare[0], to_compare[1]):
631                 compare_png_images (old, new, dest_dir)
632
633
634     def get_cell (self, oldnew):
635         def img_cell (ly, img, name):
636             if not name:
637                 name = 'source'
638             else:
639                 name = '<tt>%s</tt>' % name
640                 
641             return '''
642 <a href="%(img)s">
643 <img src="%(img)s" style="border-style: none; max-width: 500px;">
644 </a><br>
645 ''' % locals ()
646         def multi_img_cell (ly, imgs, name):
647             if not name:
648                 name = 'source'
649             else:
650                 name = '<tt>%s</tt>' % name
651
652             imgs_str = '\n'.join (['''<a href="%s">
653 <img src="%s" style="border-style: none; max-width: 500px;">
654 </a><br>''' % (img, img) 
655                                   for img in imgs])
656
657
658             return '''
659 %(imgs_str)s
660 ''' % locals ()
661
662
663
664         def cell (base, name):
665             pat = base + '-page*.png'
666             pages = glob.glob (pat)
667
668             if pages:
669                 return multi_img_cell (base + '.ly', sorted (pages), name)
670             else:
671                 return img_cell (base + '.ly', base + '.png', name)
672
673
674
675         str = cell (os.path.splitext (self.file_names[oldnew])[0], self.name ())  
676         if options.compare_images and oldnew == 1:
677             str = str.replace ('.png', '.compare.jpeg')
678             
679         return str
680
681
682     def get_distance_details (self):
683         systems = self.system_links.items ()
684         systems.sort ()
685
686         html = ""
687         for (c, link) in systems:
688             e = '<td>%d</td>' % c
689             for d in link.distance ():
690                 e += '<td>%f</td>' % d
691             
692             e = '<tr>%s</tr>' % e
693
694             html += e
695
696             e = '<td>%d</td>' % c
697             for s in (link.output_expression_details_string (),
698                       link.orphan_details_string (),
699                       link.geo_details_string ()):
700                 e += "<td>%s</td>" % s
701
702             
703             e = '<tr>%s</tr>' % e
704             html += e
705             
706         original = self.name ()
707         html = '''<html>
708 <head>
709 <title>comparison details for %(original)s</title>
710 </head>
711 <body>
712 <table border=1>
713 <tr>
714 <th>system</th>
715 <th>output</th>
716 <th>orphan</th>
717 <th>geo</th>
718 </tr>
719
720 %(html)s
721 </table>
722
723 </body>
724 </html>
725 ''' % locals ()
726         return html
727
728
729 ################################################################
730 # Files/directories
731
732 import glob
733 import re
734
735 def compare_signature_files (f1, f2):
736     s1 = read_signature_file (f1)
737     s2 = read_signature_file (f2)
738     
739     return SystemLink (s1, s2).distance ()
740
741 def paired_files (dir1, dir2, pattern):
742     """
743     Search DIR1 and DIR2 for PATTERN.
744
745     Return (PAIRED, MISSING-FROM-2, MISSING-FROM-1)
746
747     """
748
749     files = []
750     for d in (dir1,dir2):
751         found = [os.path.split (f)[1] for f in glob.glob (d + '/' + pattern)]
752         found = dict ((f, 1) for f in found)
753         files.append (found)
754         
755     pairs = []
756     missing = []
757     for f in files[0].keys ():
758         try:
759             files[1].pop (f)
760             pairs.append (f)
761         except KeyError:
762             missing.append (f)
763
764     return (pairs, files[1].keys (), missing)
765     
766 class ComparisonData:
767     def __init__ (self):
768         self.result_dict = {}
769         self.missing = []
770         self.added = []
771         self.file_links = {}
772
773     def compare_trees (self, dir1, dir2):
774         self.compare_directories (dir1, dir2)
775         
776         (root, dirs, files) = os.walk (dir1).next ()
777         for d in dirs:
778             d1 = os.path.join (dir1, d)
779             d2 = os.path.join (dir2, d)
780
781             if os.path.islink (d1) or os.path.islink (d2):
782                 continue
783             
784             if os.path.isdir (d2):
785                 self.compare_trees (d1, d2)
786     
787     def compare_directories (self, dir1, dir2):
788         for ext in ['signature', 'midi', 'log', 'profile', 'gittxt']:
789             (paired, m1, m2) = paired_files (dir1, dir2, '*.' + ext)
790
791             self.missing += [(dir1, m) for m in m1] 
792             self.added += [(dir2, m) for m in m2] 
793
794             for p in paired:
795                 if (options.max_count
796                     and len (self.file_links) > options.max_count):
797                     continue
798                 
799                 f2 = dir2 +  '/' + p
800                 f1 = dir1 +  '/' + p
801                 self.compare_files (f1, f2)
802
803     def compare_files (self, f1, f2):
804         if f1.endswith ('signature'):
805             self.compare_signature_files (f1, f2)
806         else:
807             ext = os.path.splitext (f1)[1]
808             klasses = {
809                 '.midi': MidiFileLink,
810                 '.log' : TextFileCompareLink,
811                 '.profile': ProfileFileLink,
812                 '.gittxt': GitFileCompareLink, 
813                 }
814             
815             if klasses.has_key (ext):
816                 self.compare_general_files (klasses[ext], f1, f2)
817
818     def compare_general_files (self, klass, f1, f2):
819         name = os.path.split (f1)[1]
820
821         file_link = klass (f1, f2)
822         self.file_links[name] = file_link
823         
824     def compare_signature_files (self, f1, f2):
825         name = os.path.split (f1)[1]
826         name = re.sub ('-[0-9]+.signature', '', name)
827         
828         file_link = None
829         try:
830             file_link = self.file_links[name]
831         except KeyError:
832             generic_f1 = re.sub ('-[0-9]+.signature', '.ly', f1)
833             generic_f2 = re.sub ('-[0-9]+.signature', '.ly', f2)
834             file_link = SignatureFileLink (generic_f1, generic_f2)
835             self.file_links[name] = file_link
836
837         file_link.add_file_compare (f1, f2)
838
839     def write_changed (self, dest_dir, threshold):
840         (changed, below, unchanged) = self.thresholded_results (threshold)
841
842         str = '\n'.join ([os.path.splitext (link.file_names[1])[0]
843                         for link in changed])
844         fn = dest_dir + '/changed.txt'
845         
846         open_write_file (fn).write (str)
847                 
848     def thresholded_results (self, threshold):
849         ## todo: support more scores.
850         results = [(link.distance(), link)
851                    for link in self.file_links.values ()]
852         results.sort ()
853         results.reverse ()
854
855         unchanged = [r for (d,r) in results if d == 0.0]
856         below = [r for (d,r) in results if threshold >= d > 0.0]
857         changed = [r for (d,r) in results if d > threshold]
858
859         return (changed, below, unchanged)
860                 
861     def write_text_result_page (self, filename, threshold):
862         out = None
863         if filename == '':
864             out = sys.stdout
865         else:
866             print 'writing "%s"' % filename
867             out = open_write_file (filename)
868
869         (changed, below, unchanged) = self.thresholded_results (threshold)
870
871         
872         for link in changed:
873             out.write (link.text_record_string ())
874
875         out.write ('\n\n')
876         out.write ('%d below threshold\n' % len (below))
877         out.write ('%d unchanged\n' % len (unchanged))
878         
879     def create_text_result_page (self, dir1, dir2, dest_dir, threshold):
880         self.write_text_result_page (dest_dir + '/index.txt', threshold)
881         
882     def create_html_result_page (self, dir1, dir2, dest_dir, threshold):
883         dir1 = dir1.replace ('//', '/')
884         dir2 = dir2.replace ('//', '/')
885
886         (changed, below, unchanged) = self.thresholded_results (threshold)
887
888
889         html = ''
890         old_prefix = os.path.split (dir1)[1]
891         for link in changed:
892             html += link.html_record_string (dest_dir)
893
894
895         short_dir1 = shorten_string (dir1)
896         short_dir2 = shorten_string (dir2)
897         html = '''<html>
898 <table rules="rows" border bordercolor="blue">
899 <tr>
900 <th>distance</th>
901 <th>%(short_dir1)s</th>
902 <th>%(short_dir2)s</th>
903 </tr>
904 %(html)s
905 </table>
906 </html>''' % locals()
907
908         html += ('<p>')
909         below_count = len (below)
910
911         if below_count:
912             html += ('<p>%d below threshold</p>' % below_count)
913             
914         html += ('<p>%d unchanged</p>' % len (unchanged))
915
916         dest_file = dest_dir + '/index.html'
917         open_write_file (dest_file).write (html)
918         
919     def print_results (self, threshold):
920         self.write_text_result_page ('', threshold)
921
922 def compare_trees (dir1, dir2, dest_dir, threshold):
923     data = ComparisonData ()
924     data.compare_trees (dir1, dir2)
925     data.print_results (threshold)
926
927     if os.path.isdir (dest_dir):
928         system ('rm -rf %s '% dest_dir)
929
930     data.write_changed (dest_dir, threshold)
931     data.create_html_result_page (dir1, dir2, dest_dir, threshold)
932     data.create_text_result_page (dir1, dir2, dest_dir, threshold)
933     
934 ################################################################
935 # TESTING
936
937 def mkdir (x):
938     if not os.path.isdir (x):
939         print 'mkdir', x
940         os.makedirs (x)
941
942 def link_file (x, y):
943     mkdir (os.path.split (y)[0])
944     try:
945         print x, '->', y
946         os.link (x, y)
947     except OSError, z:
948         print 'OSError', x, y, z
949         raise OSError
950     
951 def open_write_file (x):
952     d = os.path.split (x)[0]
953     mkdir (d)
954     return open (x, 'w')
955
956
957 def system (x):
958     
959     print 'invoking', x
960     stat = os.system (x)
961     assert stat == 0
962
963
964 def test_paired_files ():
965     print paired_files (os.environ["HOME"] + "/src/lilypond/scripts/",
966                         os.environ["HOME"] + "/src/lilypond-stable/buildscripts/", '*.py')
967                   
968     
969 def test_compare_trees ():
970     system ('rm -rf dir1 dir2')
971     system ('mkdir dir1 dir2')
972     system ('cp 20{-*.signature,.ly,.png,.eps,.log,.profile} dir1')
973     system ('cp 20{-*.signature,.ly,.png,.eps,.log,.profile} dir2')
974     system ('cp 20expr{-*.signature,.ly,.png,.eps,.log,.profile} dir1')
975     system ('cp 19{-*.signature,.ly,.png,.eps,.log,.profile} dir2/')
976     system ('cp 19{-*.signature,.ly,.png,.eps,.log,.profile} dir1/')
977     system ('cp 19-1.signature 19.sub-1.signature')
978     system ('cp 19.ly 19.sub.ly')
979     system ('cp 19.profile 19.sub.profile')
980     system ('cp 19.log 19.sub.log')
981     system ('cp 19.png 19.sub.png')
982     system ('cp 19.eps 19.sub.eps')
983
984     system ('cp 20multipage* dir1')
985     system ('cp 20multipage* dir2')
986     system ('cp 19multipage-1.signature dir2/20multipage-1.signature')
987
988     
989     system ('mkdir -p dir1/subdir/ dir2/subdir/')
990     system ('cp 19.sub{-*.signature,.ly,.png,.eps,.log,.profile} dir1/subdir/')
991     system ('cp 19.sub{-*.signature,.ly,.png,.eps,.log,.profile} dir2/subdir/')
992     system ('cp 20grob{-*.signature,.ly,.png,.eps,.log,.profile} dir2/')
993     system ('cp 20grob{-*.signature,.ly,.png,.eps,.log,.profile} dir1/')
994     system ('echo HEAD is 1 > dir1/tree.gittxt')
995     system ('echo HEAD is 2 > dir2/tree.gittxt')
996
997     ## introduce differences
998     system ('cp 19-1.signature dir2/20-1.signature')
999     system ('cp 19.profile dir2/20.profile')
1000     system ('cp 19.png dir2/20.png')
1001     system ('cp 19multipage-page1.png dir2/20multipage-page1.png')
1002     system ('cp 20-1.signature dir2/subdir/19.sub-1.signature')
1003     system ('cp 20.png dir2/subdir/19.sub.png')
1004     system ("sed 's/: /: 1/g'  20.profile > dir2/subdir/19.sub.profile")
1005
1006     ## radical diffs.
1007     system ('cp 19-1.signature dir2/20grob-1.signature')
1008     system ('cp 19-1.signature dir2/20grob-2.signature')
1009     system ('cp 19multipage.midi dir1/midi-differ.midi')
1010     system ('cp 20multipage.midi dir2/midi-differ.midi')
1011     system ('cp 19multipage.log dir1/log-differ.log')
1012     system ('cp 19multipage.log dir2/log-differ.log &&  echo different >> dir2/log-differ.log &&  echo different >> dir2/log-differ.log')
1013
1014     compare_trees ('dir1', 'dir2', 'compare-dir1dir2', options.threshold)
1015
1016
1017 def test_basic_compare ():
1018     ly_template = r"""
1019
1020 \version "2.10.0"
1021 #(define default-toplevel-book-handler
1022   print-book-with-defaults-as-systems )
1023
1024 #(ly:set-option (quote no-point-and-click))
1025
1026 \sourcefilename "my-source.ly"
1027  
1028 %(papermod)s
1029 \header { tagline = ##f }
1030 \score {
1031 <<
1032 \new Staff \relative c {
1033   c4^"%(userstring)s" %(extragrob)s
1034   }
1035 \new Staff \relative c {
1036   c4^"%(userstring)s" %(extragrob)s
1037   }
1038 >>
1039 \layout{}
1040 }
1041
1042 """
1043
1044     dicts = [{ 'papermod' : '',
1045                'name' : '20',
1046                'extragrob': '',
1047                'userstring': 'test' },
1048              { 'papermod' : '#(set-global-staff-size 19.5)',
1049                'name' : '19',
1050                'extragrob': '',
1051                'userstring': 'test' },
1052              { 'papermod' : '',
1053                'name' : '20expr',
1054                'extragrob': '',
1055                'userstring': 'blabla' },
1056              { 'papermod' : '',
1057                'name' : '20grob',
1058                'extragrob': 'r2. \\break c1',
1059                'userstring': 'test' },
1060              ]
1061
1062     for d in dicts:
1063         open (d['name'] + '.ly','w').write (ly_template % d)
1064         
1065     names = [d['name'] for d in dicts]
1066
1067     system ('lilypond -ddump-profile -dseparate-log-files -ddump-signatures --png -b eps ' + ' '.join (names))
1068     
1069
1070     multipage_str = r'''
1071     #(set-default-paper-size "a6")
1072     \score {
1073       \relative {c1 \pageBreak c1 }
1074       \layout {}
1075       \midi {}
1076     }
1077     '''
1078
1079     open ('20multipage.ly', 'w').write (multipage_str.replace ('c1', 'd1'))
1080     open ('19multipage.ly', 'w').write ('#(set-global-staff-size 19.5)\n' + multipage_str)
1081     system ('lilypond -dseparate-log-files -ddump-signatures --png 19multipage 20multipage ')
1082  
1083     test_compare_signatures (names)
1084     
1085 def test_compare_signatures (names, timing=False):
1086
1087     import time
1088
1089     times = 1
1090     if timing:
1091         times = 100
1092
1093     t0 = time.clock ()
1094
1095     count = 0
1096     for t in range (0, times):
1097         sigs = dict ((n, read_signature_file ('%s-1.signature' % n)) for n in names)
1098         count += 1
1099
1100     if timing:
1101         print 'elapsed', (time.clock() - t0)/count
1102
1103
1104     t0 = time.clock ()
1105     count = 0
1106     combinations = {}
1107     for (n1, s1) in sigs.items():
1108         for (n2, s2) in sigs.items():
1109             combinations['%s-%s' % (n1, n2)] = SystemLink (s1,s2).distance ()
1110             count += 1
1111
1112     if timing:
1113         print 'elapsed', (time.clock() - t0)/count
1114
1115     results = combinations.items ()
1116     results.sort ()
1117     for k,v in results:
1118         print '%-20s' % k, v
1119
1120     assert combinations['20-20'] == (0.0,0.0,0.0)
1121     assert combinations['20-20expr'][0] > 0.0
1122     assert combinations['20-19'][2] < 10.0
1123     assert combinations['20-19'][2] > 0.0
1124
1125
1126 def run_tests ():
1127     dir = 'test-output-distance'
1128
1129     do_clean = not os.path.exists (dir)
1130
1131     print 'test results in ', dir
1132     if do_clean:
1133         system ('rm -rf ' + dir)
1134         system ('mkdir ' + dir)
1135         
1136     os.chdir (dir)
1137     if do_clean:
1138         test_basic_compare ()
1139         
1140     test_compare_trees ()
1141     
1142 ################################################################
1143 #
1144
1145 def main ():
1146     p = optparse.OptionParser ("output-distance - compare LilyPond formatting runs")
1147     p.usage = 'output-distance.py [options] tree1 tree2'
1148     
1149     p.add_option ('', '--test-self',
1150                   dest="run_test",
1151                   action="store_true",
1152                   help='run test method')
1153     
1154     p.add_option ('--max-count',
1155                   dest="max_count",
1156                   metavar="COUNT",
1157                   type="int",
1158                   default=0, 
1159                   action="store",
1160                   help='only analyze COUNT signature pairs')
1161
1162     p.add_option ('', '--threshold',
1163                   dest="threshold",
1164                   default=0.3,
1165                   action="store",
1166                   type="float",
1167                   help='threshold for geometric distance')
1168
1169     p.add_option ('--no-compare-images',
1170                   dest="compare_images",
1171                   default=True,
1172                   action="store_false",
1173                   help="Don't run graphical comparisons")
1174
1175     p.add_option ('--create-images',
1176                   dest="create_images",
1177                   default=False,
1178                   action="store_true",
1179                   help="Create PNGs from EPSes")
1180
1181     p.add_option ('-o', '--output-dir',
1182                   dest="output_dir",
1183                   default=None,
1184                   action="store",
1185                   type="string",
1186                   help='where to put the test results [tree2/compare-tree1tree2]')
1187
1188     global options
1189     (options, a) = p.parse_args ()
1190
1191     if options.run_test:
1192         run_tests ()
1193         sys.exit (0)
1194
1195     if len (a) != 2:
1196         p.print_usage()
1197         sys.exit (2)
1198
1199     name = options.output_dir
1200     if not name:
1201         name = a[0].replace ('/', '')
1202         name = os.path.join (a[1], 'compare-' + shorten_string (name))
1203     
1204     compare_trees (a[0], a[1], name, options.threshold)
1205
1206 if __name__ == '__main__':
1207     main()
1208