2 This file is part of LilyPond, the GNU music typesetter.
4 Copyright (C) 2005--2011 Han-Wen Nienhuys <hanwen@xs4all.nl>
7 LilyPond is free software: you can redistribute it and/or modify
8 it under the terms of the GNU General Public License as published by
9 the Free Software Foundation, either version 3 of the License, or
10 (at your option) any later version.
12 LilyPond is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
17 You should have received a copy of the GNU General Public License
18 along with LilyPond. If not, see <http://www.gnu.org/licenses/>.
22 #include "tie-formatting-problem.hh"
24 #include "axis-group-interface.hh"
25 #include "paper-column.hh"
27 #include "directional-element-interface.hh"
28 #include "libc-extension.hh"
30 #include "note-head.hh"
31 #include "rhythmic-head.hh"
33 #include "staff-symbol-referencer.hh"
35 #include "tie-configuration.hh"
38 #include "pointer-group-interface.hh"
39 #include "output-def.hh"
42 Tie_formatting_problem::print_ties_configuration (Ties_configuration const *ties)
44 for (vsize i = 0; i < ties->size (); i++)
46 char const *man_pos = (specifications_[i].has_manual_position_) ? "(M)" : "";
47 char const *man_dir = (specifications_[i].has_manual_dir_) ? "(M)" : "";
48 char const *dir = (ties->at (i).dir_ == UP) ? "up" : "dn";
50 printf ("(P%d%s, %s%s) ", ties->at (i).position_, man_pos, dir, man_dir);
56 Tie_formatting_problem::get_attachment (Real y, Drul_array<int> columns) const
58 Interval attachments (0, 0);
62 Tuple2<int> key (columns[d], int (d));
63 Chord_outline_map::const_iterator i (chord_outlines_.find (key));
64 if (i == chord_outlines_.end ())
65 programming_error ("Cannot find chord outline");
67 attachments[d] = i->second.height (y);
69 while (flip (&d) != LEFT);
74 Tie_formatting_problem::Tie_formatting_problem ()
78 use_horizontal_spacing_ = true;
81 Tie_formatting_problem::~Tie_formatting_problem ()
83 for (Tie_configuration_map::const_iterator i (possibilities_.begin ());
84 i != possibilities_.end (); i++)
89 Tie_formatting_problem::set_column_chord_outline (vector<Item *> bounds,
93 Real staff_space = Staff_symbol_referencer::staff_space (bounds[0]);
96 vector<Box> head_boxes;
99 for (vsize i = 0; i < bounds.size (); i++)
101 Grob *head = bounds[i];
102 if (!Note_head::has_interface (head))
106 stem = unsmob_grob (head->get_object ("stem"));
108 Real p = Staff_symbol_referencer::get_position (head);
109 Interval y ((p - 1) * 0.5 * staff_space,
110 (p + 1) * 0.5 * staff_space);
112 Interval x = head->extent (x_refpoint_, X_AXIS);
113 head_boxes.push_back (Box (x, y));
114 boxes.push_back (Box (x, y));
116 Grob *dots = Rhythmic_head::get_dots (head);
117 if (dir == LEFT && dots)
119 Interval x = dots->extent (x_refpoint_, X_AXIS);
120 int p = int (Staff_symbol_referencer::get_position (dots));
123 TODO: shouldn't this use column-rank dependent key?
125 dot_positions_.insert (p);
128 Interval y (dots->extent (dots, Y_AXIS));
129 y.translate (p * staff_space * 0.5);
131 boxes.push_back (Box (x, y));
135 Tuple2<int> key (column_rank, int (dir));
139 if (Stem::is_normal_stem (stem))
142 x.add_point (stem->relative_coordinate (x_refpoint_, X_AXIS));
143 x.widen (staff_space / 20); // ugh.
146 Real stem_end_position = 0.0;
147 if (Stem::is_cross_staff (stem))
148 stem_end_position = get_grob_direction (stem) * infinity_f;
151 if (use_horizontal_spacing_ || !Stem::get_beam (stem))
152 stem_end_position = stem->extent (stem, Y_AXIS)[get_grob_direction (stem)];
154 // May want to change this to the stem's pure height...
155 stem_end_position = Stem::note_head_positions (stem)[get_grob_direction (stem)]
159 y.add_point (stem_end_position);
161 Direction stemdir = get_grob_direction (stem);
162 y.add_point (Stem::head_positions (stem)[-stemdir]
168 boxes.push_back (Box (x, y));
170 stem_extents_[key].unite (Box (x, y));
174 Grob *flag = Stem::flag (stem);
177 Grob* commony = stem->common_refpoint (flag, Y_AXIS);
178 boxes.push_back (Box (flag->extent (x_refpoint_, X_AXIS),
179 flag->extent (commony, Y_AXIS)));
185 Grob *head = Stem::support_head (stem);
188 In case of invisible stem, don't pass x-center of heads.
190 Real x_center = head->extent (x_refpoint_, X_AXIS).center ();
192 x_ext[-dir] = x_center;
193 x_ext[dir] = infinity_f * dir;
195 for (vsize j = 0; j < head_boxes.size (); j++)
196 y_ext.unite (head_boxes[j][Y_AXIS]);
198 boxes.push_back (Box (x_ext, y_ext));
201 extract_grob_set (stem, "note-heads", heads);
202 for (vsize i = 0; i < heads.size (); i++)
204 if (find (bounds.begin (), bounds.end (), dynamic_cast<Item *> (heads[i])) == bounds.end ())
207 other untied notes in the same chord.
210 Interval y = Staff_symbol_referencer::extent_in_staff (heads[i]);
211 Interval x = heads[i]->extent (x_refpoint_, X_AXIS);
212 boxes.push_back (Box (x, y));
215 Grob *acc = unsmob_grob (heads[i]->get_object ("accidental-grob"));
217 acc->get_property ("stencil"); /* trigger tie-related suicide */
219 if (acc && acc->is_live () && dir == RIGHT)
221 boxes.push_back (Box (acc->extent (x_refpoint_, X_AXIS),
222 Staff_symbol_referencer::extent_in_staff (acc)));
225 head_positions_[column_rank].add_point (int (Staff_symbol_referencer::get_position (heads[i])));
230 Direction updowndir = DOWN;
235 if (head_boxes.size ())
237 Box b = boundary (head_boxes, updowndir, 0);
239 x[-dir] = b[X_AXIS].linear_combination (-dir / 2);
240 y[-updowndir] = b[Y_AXIS][updowndir];
241 y[updowndir] = updowndir * infinity_f;
245 boxes.push_back (Box (x, y));
247 while (flip (&updowndir) != DOWN);
249 /* todo: the horizon_padding is somewhat arbitrary */
250 chord_outlines_[key] = Skyline (boxes, details_.skyline_padding_, Y_AXIS, -dir);
251 if (bounds[0]->break_status_dir ())
253 Interval iv (Axis_group_interface::staff_extent (bounds[0], x_refpoint_, X_AXIS, y_refpoint_, Y_AXIS));
255 iv.add_point (bounds[0]->relative_coordinate (x_refpoint_, X_AXIS));
257 chord_outlines_[key].set_minimum_height (iv[-dir]);
262 for (vsize j = 0; j < head_boxes.size (); j++)
264 x.unite (head_boxes[j][X_AXIS]);
267 chord_outlines_[key].set_minimum_height (x[dir]);
270 head_extents_[key].set_empty ();
271 for (vsize i = 0; i < head_boxes.size (); i++)
273 head_extents_[key].unite (head_boxes[i]);
278 Tie_formatting_problem::set_chord_outline (vector<Item *> bounds,
283 for (vsize i = 0; i < bounds.size (); i++)
284 ranks.push_back (bounds[i]->get_column ()->get_rank ());
286 vector_sort (ranks, less<int> ());
289 for (vsize i = 0; i < ranks.size (); i++)
291 vector<Item *> col_items;
292 for (vsize j = 0; j < bounds.size (); j++)
294 if (bounds[j]->get_column ()->get_rank () == ranks[i])
295 col_items.push_back (bounds[j]);
298 set_column_chord_outline (col_items, dir, ranks[i]);
303 Tie_formatting_problem::from_tie (Grob *tie)
306 ties.push_back (tie);
309 details_.from_grob (tie);
313 Tie_formatting_problem::common_x_refpoint () const
319 Tie_formatting_problem::from_ties (vector<Grob *> const &ties)
324 x_refpoint_ = ties[0];
325 y_refpoint_ = ties[0];
326 for (vsize i = 0; i < ties.size (); i++)
328 Spanner *tie = dynamic_cast<Spanner *> (ties[i]);
329 Item *l = tie->get_bound (LEFT);
330 Item *r = tie->get_bound (RIGHT);
332 x_refpoint_ = l->common_refpoint (x_refpoint_, X_AXIS);
333 x_refpoint_ = r->common_refpoint (x_refpoint_, X_AXIS);
335 if (!l->break_status_dir ())
336 y_refpoint_ = l->common_refpoint (y_refpoint_, Y_AXIS);
337 if (!r->break_status_dir ())
338 y_refpoint_ = r->common_refpoint (y_refpoint_, Y_AXIS);
341 details_.from_grob (ties[0]);
346 vector<Item *> bounds;
348 for (vsize i = 0; i < ties.size (); i++)
350 Item *it = dynamic_cast<Spanner *> (ties[i])->get_bound (d);
351 if (it->break_status_dir ())
352 it = it->get_column ();
354 bounds.push_back (it);
357 set_chord_outline (bounds, d);
359 while (flip (&d) != LEFT);
361 for (vsize i = 0; i < ties.size (); i++)
363 Tie_specification spec;
364 spec.from_grob (ties[i]);
368 spec.note_head_drul_[d] = Tie::head (ties[i], d);
369 spec.column_ranks_[d] = Tie::get_column_rank (ties[i], d);
371 while (flip (&d) != LEFT);
372 specifications_.push_back (spec);
377 Tie_formatting_problem::from_semi_ties (vector<Grob *> const &semi_ties, Direction head_dir)
379 if (semi_ties.empty ())
382 use_horizontal_spacing_ = false;
383 details_.from_grob (semi_ties[0]);
384 vector<Item *> heads;
386 int column_rank = -1;
387 for (vsize i = 0; i < semi_ties.size (); i++)
389 Tie_specification spec;
390 Item *head = unsmob_item (semi_ties[i]->get_object ("note-head"));
393 programming_error ("LV tie without head?!");
397 spec.position_ = int (Staff_symbol_referencer::get_position (head));
400 spec.from_grob (semi_ties[i]);
402 spec.note_head_drul_[head_dir] = head;
403 column_rank = Tie::get_column_rank (semi_ties[i], head_dir);
404 spec.column_ranks_ = Drul_array<int> (column_rank, column_rank);
405 heads.push_back (head);
406 specifications_.push_back (spec);
409 x_refpoint_ = semi_ties[0];
410 y_refpoint_ = semi_ties[0];
412 for (vsize i = 0; i < semi_ties.size (); i++)
414 x_refpoint_ = semi_ties[i]->common_refpoint (x_refpoint_, X_AXIS);
415 y_refpoint_ = semi_ties[i]->common_refpoint (y_refpoint_, Y_AXIS);
417 for (vsize i = 0; i < heads.size (); i++)
419 x_refpoint_ = heads[i]->common_refpoint (x_refpoint_, X_AXIS);
420 y_refpoint_ = heads[i]->common_refpoint (y_refpoint_, Y_AXIS);
423 set_chord_outline (heads, head_dir);
425 Tuple2<int> head_key (column_rank, head_dir);
426 Tuple2<int> open_key (column_rank, -head_dir);
427 Real extremal = chord_outlines_[head_key].max_height ();
429 chord_outlines_[open_key] = Skyline (head_dir);
430 chord_outlines_[open_key].set_minimum_height (extremal - head_dir * 1.5);
434 Tie_formatting_problem::get_tie_specification (int i) const
436 return specifications_[i];
440 Return configuration, create it if necessary.
443 Tie_formatting_problem::get_configuration (int pos, Direction dir, Drul_array<int> columns,
446 int key_components[] =
448 pos, dir, columns[LEFT], columns[RIGHT]
450 Tuple<int, 4> key (key_components);
452 Tie_configuration_map::const_iterator f = possibilities_.find (key);
453 if (f != possibilities_.end ())
458 Tie_configuration *conf = generate_configuration (pos, dir, columns, tune_dy);
459 ((Tie_formatting_problem *) this)->possibilities_[key] = conf;
464 Tie_formatting_problem::generate_configuration (int pos, Direction dir,
465 Drul_array<int> columns, bool y_tune) const
467 Tie_configuration *conf = new Tie_configuration;
468 conf->position_ = pos;
471 conf->column_ranks_ = columns;
473 Real y = conf->position_ * 0.5 * details_.staff_space_;
475 if (dot_positions_.find (pos) != dot_positions_.end ())
477 conf->delta_y_ += dir * 0.25 * details_.staff_space_;
482 && max (fabs (get_head_extent (columns[LEFT], LEFT, Y_AXIS)[dir] - y),
483 fabs (get_head_extent (columns[RIGHT], RIGHT, Y_AXIS)[dir] - y)) < 0.25
484 && !Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos))
487 = (get_head_extent (columns[LEFT], LEFT, Y_AXIS)[dir] - y)
488 + dir * details_.outer_tie_vertical_gap_;
493 conf->attachment_x_ = get_attachment (y + conf->delta_y_, conf->column_ranks_);
494 Real h = conf->height (details_);
499 - should make sliding criterion, should flatten ties if
501 - they're just the wrong (ie. touching line at top & bottom)
505 if (head_positions_slice (columns[LEFT]).contains (pos)
506 || head_positions_slice (columns[RIGHT]).contains (pos)
507 || abs (pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_))
509 if (h < details_.intra_space_threshold_ * 0.5 * details_.staff_space_)
511 if (!Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos)
512 && abs (pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_))
514 conf->center_tie_vertically (details_);
516 else if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos))
518 conf->delta_y_ += dir *
519 details_.tip_staff_line_clearance_ * 0.5 * details_.staff_space_;
524 Real top_y = y + conf->delta_y_ + conf->dir_ * h;
525 Real top_pos = top_y / (0.5 * details_.staff_space_);
526 int round_pos = int (my_round (top_pos));
528 /* TODO: should use other variable? */
529 Real clearance = details_.center_staff_line_clearance_;
530 if (fabs (top_pos - round_pos) < clearance
531 && Staff_symbol_referencer::on_staff_line (details_.staff_symbol_referencer_,
534 Real new_y = (round_pos + clearance * conf->dir_) * 0.5 * details_.staff_space_;
535 conf->delta_y_ = (new_y - top_y);
540 conf->attachment_x_ = get_attachment (y + conf->delta_y_, conf->column_ranks_);
541 if (conf->height (details_) < details_.intra_space_threshold_ * 0.5 * details_.staff_space_)
544 This is less sensible for long ties, since those are more
547 Interval close_by = get_attachment (y
549 + (dir * details_.intra_space_threshold_ * 0.25
550 * details_.staff_space_),
551 conf->column_ranks_);
553 conf->attachment_x_.intersect (close_by);
556 conf->attachment_x_.widen ( - details_.x_gap_);
558 if (conf->column_span_length ())
561 avoid the stems that we attach to as well. We don't do this
562 for semities (span length = 0)
564 It would be better to check D against HEAD-DIRECTION if
570 Real y = conf->position_ * details_.staff_space_ * 0.5 + conf->delta_y_;
571 if (get_stem_extent (conf->column_ranks_[d], d, X_AXIS).is_empty ()
572 || !get_stem_extent (conf->column_ranks_[d], d, Y_AXIS).contains (y))
575 conf->attachment_x_[d]
576 = d * min (d * conf->attachment_x_[d],
577 d * (get_stem_extent (conf->column_ranks_[d], d, X_AXIS)[-d] - d * details_.stem_gap_));
579 while (flip (&d) != LEFT);
585 Tie_formatting_problem::get_head_extent (int col, Direction d, Axis a) const
587 Column_extent_map::const_iterator i = head_extents_.find (Tuple2<int> (col, int (d)));
588 if (i != head_extents_.end ())
589 return (*i).second[a];
595 Tie_formatting_problem::get_stem_extent (int col, Direction d, Axis a) const
597 Column_extent_map::const_iterator i = stem_extents_.find (Tuple2<int> (col, int (d)));
598 if (i != stem_extents_.end ())
599 return (*i).second[a];
605 TIE_IDX and TIES_CONF are optional.
608 Tie_formatting_problem::score_aptitude (Tie_configuration *conf,
609 Tie_specification const &spec,
610 Ties_configuration *ties_conf, int tie_idx) const
613 Real curve_y = conf->position_ * details_.staff_space_ * 0.5 + conf->delta_y_;
614 Real tie_y = spec.position_ * details_.staff_space_ * 0.5;
615 if (sign (curve_y - tie_y) != conf->dir_)
617 Real p = details_.wrong_direction_offset_penalty_;
619 ties_conf->add_tie_score (p, tie_idx, "wrong dir");
625 Real relevant_dist = max (fabs (curve_y - tie_y) - 0.5, 0.0);
626 Real p = details_.vertical_distance_penalty_factor_ * convex_amplifier (1.0, 0.9, relevant_dist);
628 ties_conf->add_tie_score (p, tie_idx, "vdist");
636 if (!spec.note_head_drul_[d])
639 Interval head_x = spec.note_head_drul_[d]->extent (x_refpoint_, X_AXIS);
640 Real dist = head_x.distance (conf->attachment_x_[d]);
643 TODO: flatten with log or sqrt.
645 Real p = details_.horizontal_distance_penalty_factor_
646 * convex_amplifier (1.25, 1.0, dist);
648 ties_conf->add_tie_score (p, tie_idx,
649 (d == LEFT) ? "lhdist" : "rhdist");
654 while (flip (&d) != LEFT);
657 && ties_conf->size () == 1)
660 Drul_array<Grob *> stems (0, 0);
663 if (!spec.note_head_drul_[d])
666 Grob *stem = unsmob_grob (spec.note_head_drul_[d]->get_object ("stem"));
668 && Stem::is_normal_stem (stem))
671 while (flip (&d) != LEFT);
673 bool tie_stem_dir_ok = true;
674 bool tie_position_dir_ok = true;
675 if (stems[LEFT] && !stems[RIGHT])
676 tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[LEFT]);
677 else if (!stems[LEFT] && stems[RIGHT])
678 tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[RIGHT]);
679 else if (stems[LEFT] && stems[RIGHT]
680 && get_grob_direction (stems[LEFT]) == get_grob_direction (stems[RIGHT]))
681 tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[LEFT]);
682 else if (spec.position_)
683 tie_position_dir_ok = conf->dir_ == sign (spec.position_);
685 if (!tie_stem_dir_ok)
686 ties_conf->add_score (details_.same_dir_as_stem_penalty_, "tie/stem dir");
687 if (!tie_position_dir_ok)
688 ties_conf->add_score (details_.same_dir_as_stem_penalty_, "tie/pos dir");
695 Tie_formatting_problem::head_positions_slice (int rank) const
697 Position_extent_map::const_iterator i (head_positions_.find (rank));
698 if (i != head_positions_.end ())
707 Score a configuration, ie. how well these ties looks without regard
708 to the note heads that they should connect to.
711 Tie_formatting_problem::score_configuration (Tie_configuration *conf) const
718 Real length = conf->attachment_x_.length ();
721 = peak_around (0.33 * details_.min_length_, details_.min_length_, length);
722 conf->add_score (details_.min_length_penalty_factor_
723 * length_penalty, "minlength");
725 Real tip_pos = conf->position_ + conf->delta_y_ / 0.5 * details_.staff_space_;
726 Real tip_y = tip_pos * details_.staff_space_ * 0.5;
727 Real height = conf->height (details_);
729 Real top_y = tip_y + conf->dir_ * height;
730 Real top_pos = 2 * top_y / details_.staff_space_;
731 Real round_top_pos = rint (top_pos);
732 if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_,
734 && Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_) > top_y)
736 conf->add_score (details_.staff_line_collision_penalty_
737 * peak_around (0.1 * details_.center_staff_line_clearance_,
738 details_.center_staff_line_clearance_,
739 fabs (top_pos - round_top_pos)),
743 int rounded_tip_pos = int (rint (tip_pos));
744 if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, rounded_tip_pos)
745 && (head_positions_slice (conf->column_ranks_[LEFT]).contains (rounded_tip_pos)
746 || head_positions_slice (conf->column_ranks_[RIGHT]).contains (rounded_tip_pos)
747 || abs (rounded_tip_pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_))
750 conf->add_score (details_.staff_line_collision_penalty_
751 * peak_around (0.1 * details_.tip_staff_line_clearance_,
752 details_.tip_staff_line_clearance_,
753 fabs (tip_pos - rint (tip_pos))),
757 if (!dot_x_.is_empty ())
760 Real x = dot_x_.center ();
762 Bezier b = conf->get_transformed_bezier (details_);
763 if (b.control_point_extent (X_AXIS).contains (x))
765 Real y = b.get_other_coordinate (X_AXIS, x);
767 for (set<int>::const_iterator i (dot_positions_.begin ());
768 i != dot_positions_.end (); i++)
771 conf->add_score (details_.dot_collision_penalty_
772 * peak_around (.1 * details_.dot_collision_clearance_,
773 details_.dot_collision_clearance_,
774 fabs (dot_pos * details_.staff_space_ * 0.5 - y)),
780 conf->scored_ = true;
784 Tie_formatting_problem::score_ties_aptitude (Ties_configuration *ties) const
786 if (ties->size () != specifications_.size ())
788 programming_error ("Huh? Mismatch between sizes.");
792 for (vsize i = 0; i < ties->size (); i++)
793 score_aptitude (&ties->at (i), specifications_[i],
798 Tie_formatting_problem::score_ties (Ties_configuration *ties) const
803 score_ties_configuration (ties);
804 score_ties_aptitude (ties);
805 ties->scored_ = true;
809 Tie_formatting_problem::score_ties_configuration (Ties_configuration *ties) const
811 for (vsize i = 0; i < ties->size (); i++)
813 score_configuration (&ties->at (i));
814 ties->add_tie_score (ties->at (i).score (), i, "conf");
817 Real last_edge = 0.0;
818 Real last_center = 0.0;
819 for (vsize i = 0; i < ties->size (); i++)
821 Bezier b (ties->at (i).get_transformed_bezier (details_));
823 Real center = b.curve_point (0.5)[Y_AXIS];
824 Real edge = b.curve_point (0.0)[Y_AXIS];
828 if (edge <= last_edge)
829 ties->add_score (details_.tie_column_monotonicity_penalty_, "monoton edge");
830 if (center <= last_center)
831 ties->add_score (details_.tie_column_monotonicity_penalty_, "monoton cent");
833 ties->add_score (details_.tie_tie_collision_penalty_ *
834 peak_around (0.1 * details_.tie_tie_collision_distance_,
835 details_.tie_tie_collision_distance_,
836 fabs (center - last_center)),
838 ties->add_score (details_.tie_tie_collision_penalty_ *
839 peak_around (0.1 * details_.tie_tie_collision_distance_,
840 details_.tie_tie_collision_distance_,
841 fabs (edge - last_edge)), "tietie edge");
845 last_center = center;
848 if (ties->size () > 1)
850 ties->add_score (details_.outer_tie_length_symmetry_penalty_factor_
851 * fabs (ties->at (0).attachment_x_.length () - ties->back ().attachment_x_.length ()),
854 ties->add_score (details_.outer_tie_vertical_distance_symmetry_penalty_factor_
855 * fabs (fabs (specifications_[0].position_ * 0.5 * details_.staff_space_
856 - (ties->at (0).position_ * 0.5 * details_.staff_space_
857 + ties->at (0).delta_y_))
859 fabs (specifications_.back ().position_ * 0.5 * details_.staff_space_
860 - (ties->back ().position_ * 0.5 * details_.staff_space_
861 + ties->back ().delta_y_))),
867 Generate with correct X-attachments and beziers, copying delta_y_
868 from TIES_CONFIG if necessary.
871 Tie_formatting_problem::generate_ties_configuration (Ties_configuration const &ties_config)
873 Ties_configuration copy;
874 for (vsize i = 0; i < ties_config.size (); i++)
876 Tie_configuration *ptr = get_configuration (ties_config[i].position_, ties_config[i].dir_,
877 ties_config[i].column_ranks_,
878 !specifications_[i].has_manual_delta_y_);
879 if (specifications_[i].has_manual_delta_y_)
882 = (specifications_[i].manual_position_ - ties_config[i].position_)
883 * 0.5 * details_.staff_space_;
885 copy.push_back (*ptr);
892 Tie_formatting_problem::generate_base_chord_configuration ()
894 Ties_configuration ties_config;
895 for (vsize i = 0; i < specifications_.size (); i++)
897 Tie_configuration conf;
898 if (specifications_[i].has_manual_dir_)
899 conf.dir_ = specifications_[i].manual_dir_;
900 if (specifications_[i].has_manual_position_)
902 conf.position_ = (int) my_round (specifications_[i].manual_position_);
903 if (specifications_[i].has_manual_delta_y_)
904 conf.delta_y_ = (specifications_[i].manual_position_ - conf.position_)
905 * 0.5 * details_.staff_space_;
909 conf.position_ = specifications_[i].position_;
912 conf.column_ranks_ = specifications_[i].column_ranks_;
913 ties_config.push_back (conf);
916 set_ties_config_standard_directions (&ties_config);
917 for (vsize i = 0; i < ties_config.size (); i++)
918 if (!specifications_[i].manual_position_)
919 ties_config[i].position_ += ties_config[i].dir_;
921 ties_config = generate_ties_configuration (ties_config);
927 Tie_formatting_problem::find_best_variation (Ties_configuration const &base,
928 vector<Tie_configuration_variation> const &vars)
930 Ties_configuration best = base;
933 This simply is 1-opt: we have K substitions, and we try applying
934 exactly every one for each.
936 for (vsize i = 0; i < vars.size (); i++)
938 Ties_configuration variant (base);
939 for (vsize j = 0; j < vars[i].index_suggestion_pairs_.size (); j++)
940 variant[vars[i].index_suggestion_pairs_[j].first] = *vars[i].index_suggestion_pairs_[j].second;
942 variant.reset_score ();
943 score_ties (&variant);
945 if (variant.score () < best.score ())
955 Tie_formatting_problem::generate_optimal_configuration ()
957 Ties_configuration base = generate_base_chord_configuration ();
960 vector<Tie_configuration_variation> vars;
961 if (specifications_.size () > 1)
962 vars = generate_collision_variations (base);
964 vars = generate_single_tie_variations (base);
966 Ties_configuration best = find_best_variation (base, vars);
968 if (specifications_.size () > 1)
970 vars = generate_extremal_tie_variations (best);
971 best = find_best_variation (best, vars);
977 Tie_formatting_problem::set_ties_config_standard_directions (Ties_configuration *tie_configs)
979 if (tie_configs->empty ())
982 if (!tie_configs->at (0).dir_)
984 if (tie_configs->size () == 1)
985 tie_configs->at (0).dir_ = Direction (sign (tie_configs->at (0).position_));
987 if (!tie_configs->at (0).dir_)
988 tie_configs->at (0).dir_
989 = (tie_configs->size () > 1) ? DOWN : details_.neutral_direction_;
992 if (!tie_configs->back ().dir_)
993 tie_configs->back ().dir_ = UP;
998 for (vsize i = 1; i < tie_configs->size (); i++)
1000 Real diff = (tie_configs->at (i).position_
1001 - tie_configs->at (i - 1).position_);
1004 = specifications_[i].column_span () - specifications_[i - 1].column_span ();
1005 if (span_diff && fabs (diff) <= 2)
1008 tie_configs->at (i).dir_ = UP;
1009 else if (span_diff < 0)
1010 tie_configs->at (i - 1).dir_ = DOWN;
1012 else if (fabs (diff) <= 1)
1014 if (!tie_configs->at (i - 1).dir_)
1015 tie_configs->at (i - 1).dir_ = DOWN;
1016 if (!tie_configs->at (i).dir_)
1017 tie_configs->at (i).dir_ = UP;
1021 for (vsize i = 1; i + 1 < tie_configs->size (); i++)
1023 Tie_configuration &conf = tie_configs->at (i);
1027 Direction position_dir
1028 = Direction (sign (conf.position_));
1030 position_dir = DOWN;
1032 conf.dir_ = position_dir;
1036 vector<Tie_configuration_variation>
1037 Tie_formatting_problem::generate_extremal_tie_variations (Ties_configuration const &ties) const
1039 vector<Tie_configuration_variation> vars;
1041 for (int i = 1; i <= details_.multi_tie_region_size_; i++)
1043 Drul_array<Tie_configuration *> configs (0, 0);
1046 const Tie_configuration &config = boundary (ties, d, 0);
1047 if (config.dir_ == d
1048 && !boundary (specifications_, d, 0).has_manual_position_)
1050 Tie_configuration_variation var;
1051 configs[d] = get_configuration (config.position_ + d * i, d,
1052 config.column_ranks_,
1054 var.add_suggestion ((d == DOWN) ? 0 : ties.size () - 1,
1056 vars.push_back (var);
1059 while (flip (&d) != DOWN);
1060 if (configs[LEFT] && configs[RIGHT])
1062 Tie_configuration_variation var;
1063 var.add_suggestion (0, configs[DOWN]);
1064 var.add_suggestion (ties.size () - 1, configs[UP]);
1065 vars.push_back (var);
1072 vector<Tie_configuration_variation>
1073 Tie_formatting_problem::generate_single_tie_variations (Ties_configuration const &ties) const
1075 vector<Tie_configuration_variation> vars;
1077 int sz = details_.single_tie_region_size_;
1078 if (specifications_[0].has_manual_position_)
1080 for (int i = 0; i < sz; i++)
1086 && ties[0].dir_ == d)
1089 int p = ties[0].position_ + i * d;
1091 if (!specifications_[0].has_manual_dir_
1092 || d == specifications_[0].manual_dir_)
1094 Tie_configuration_variation var;
1095 var.add_suggestion (0,
1096 get_configuration (p,
1097 d, specifications_[0].column_ranks_,
1098 !specifications_[0].has_manual_delta_y_));
1099 vars.push_back (var);
1102 while (flip (&d) != LEFT);
1107 vector<Tie_configuration_variation>
1108 Tie_formatting_problem::generate_collision_variations (Ties_configuration const &ties) const
1110 Real center_distance_tolerance = 0.25;
1112 vector<Tie_configuration_variation> vars;
1113 Real last_center = 0.0;
1114 for (vsize i = 0; i < ties.size (); i++)
1116 Bezier b (ties[i].get_transformed_bezier (details_));
1118 Real center = b.curve_point (0.5)[Y_AXIS];
1122 if (center <= last_center + center_distance_tolerance)
1124 if (!specifications_[i].has_manual_dir_)
1126 Tie_configuration_variation var;
1127 var.add_suggestion (i,
1128 get_configuration (specifications_[i].position_
1132 ties[i].column_ranks_,
1133 !specifications_[i].has_manual_delta_y_
1136 vars.push_back (var);
1139 if (!specifications_[i - 1].has_manual_dir_)
1141 Tie_configuration_variation var;
1142 var.add_suggestion (i - 1,
1143 get_configuration (specifications_[i - 1].position_
1146 specifications_[i - 1].column_ranks_,
1147 !specifications_[i - 1].has_manual_delta_y_));
1149 vars.push_back (var);
1152 if (i == 1 && !specifications_[i - 1].has_manual_position_
1153 && ties[i - 1].dir_ == DOWN)
1155 Tie_configuration_variation var;
1156 var.add_suggestion (i - 1,
1157 get_configuration (specifications_[i - 1].position_ - 1, DOWN,
1158 specifications_[i - 1].column_ranks_,
1159 !specifications_[i - 1].has_manual_delta_y_
1161 vars.push_back (var);
1163 if (i == ties.size () && !specifications_[i].has_manual_position_
1164 && ties[i].dir_ == UP)
1166 Tie_configuration_variation var;
1167 var.add_suggestion (i,
1168 get_configuration (specifications_[i].position_
1170 specifications_[i].column_ranks_,
1171 !specifications_[i].has_manual_delta_y_
1173 vars.push_back (var);
1176 else if (dot_positions_.find (ties[i].position_) != dot_positions_.end ()
1177 && !specifications_[i].has_manual_position_)
1179 Tie_configuration_variation var;
1180 var.add_suggestion (i,
1181 get_configuration (ties[i].position_ + ties[i].dir_,
1183 ties[i].column_ranks_,
1184 !specifications_[i].has_manual_delta_y_
1186 vars.push_back (var);
1191 last_center = center;
1198 Tie_formatting_problem::set_manual_tie_configuration (SCM manual_configs)
1201 for (SCM s = manual_configs;
1202 scm_is_pair (s) && k < specifications_.size (); s = scm_cdr (s))
1204 SCM entry = scm_car (s);
1205 if (scm_is_pair (entry))
1207 Tie_specification &spec = specifications_[k];
1209 if (scm_is_number (scm_car (entry)))
1211 spec.has_manual_position_ = true;
1212 spec.manual_position_ = scm_to_double (scm_car (entry));
1213 spec.has_manual_delta_y_ = (scm_inexact_p (scm_car (entry)) == SCM_BOOL_T);
1216 if (scm_is_number (scm_cdr (entry)))
1218 spec.has_manual_dir_ = true;
1219 spec.manual_dir_ = Direction (scm_to_int (scm_cdr (entry)));
1227 Tie_formatting_problem::set_debug_scoring (Ties_configuration const &base)
1229 #if DEBUG_TIE_SCORING
1230 if (to_boolean (x_refpoint_->layout ()
1231 ->lookup_variable (ly_symbol2scm ("debug-tie-scoring"))))
1233 for (vsize i = 0; i < base.size (); i++)
1235 string card = base.complete_tie_card (i);
1236 specifications_[i].tie_grob_->set_property ("annotation",
1237 ly_string2scm (card));