X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=lily%2Ftie-formatting-problem.cc;h=f7682c904c56db5af97bcd4c9c3f0c240e6a7b84;hb=3f8485925e8c879fe4c9ae86acef9804126c3b91;hp=0ea0f8c966dfcf37c4258d1adda5ccdc83d5e708;hpb=d3ce419548f0159fd704c777e2a8435d2be055ed;p=lilypond.git diff --git a/lily/tie-formatting-problem.cc b/lily/tie-formatting-problem.cc index 0ea0f8c966..f7682c904c 100644 --- a/lily/tie-formatting-problem.cc +++ b/lily/tie-formatting-problem.cc @@ -1,41 +1,64 @@ /* - tie-formatting-problem.cc -- implement Tie_formatting_problem6 + tie-formatting-problem.cc -- implement Tie_formatting_problem source file of the GNU LilyPond music typesetter - (c) 2005 Han-Wen Nienhuys + (c) 2005--2007 Han-Wen Nienhuys */ #include "tie-formatting-problem.hh" -#include "tie-configuration.hh" -#include "directional-element-interface.hh" -#include "staff-symbol-referencer.hh" -#include "tie.hh" -#include "item.hh" -#include "spanner.hh" +#include "paper-column.hh" #include "bezier.hh" -#include "stem.hh" +#include "directional-element-interface.hh" +#include "libc-extension.hh" +#include "misc.hh" #include "note-head.hh" #include "rhythmic-head.hh" +#include "spanner.hh" +#include "staff-symbol-referencer.hh" +#include "stem.hh" +#include "tie-configuration.hh" +#include "tie.hh" #include "warn.hh" +#include "pointer-group-interface.hh" +#include "output-def.hh" + +void +Tie_formatting_problem::print_ties_configuration (Ties_configuration const *ties) +{ + for (vsize i = 0; i < ties->size (); i++) + { + char const *man_pos = (specifications_[i].has_manual_position_) ? "(M)" : ""; + char const *man_dir = (specifications_[i].has_manual_dir_) ? "(M)" : ""; + char const *dir = (ties->at (i).dir_ == UP) ? "up" : "dn"; + + printf ("(P%d%s, %s%s) ", ties->at (i).position_, man_pos, dir, man_dir); + } + printf ("\n"); +} Interval -Tie_formatting_problem::get_attachment (Real y) const +Tie_formatting_problem::get_attachment (Real y, Drul_array columns) const { - Interval attachments; + Interval attachments (0,0); Direction d = LEFT; do { - attachments[d] = skyline_height (chord_outlines_[d], y, -d); + Tuple2 key (columns[d], int (d)); + Chord_outline_map::const_iterator i (chord_outlines_.find (key)); + if (i == chord_outlines_.end ()) + programming_error ("Cannot find chord outline"); + else + attachments[d] = i->second.height (y); } while (flip (&d) != LEFT); return attachments; } -Tie_formatting_problem::Tie_formatting_problem() +Tie_formatting_problem::Tie_formatting_problem () { x_refpoint_ = 0; } @@ -48,20 +71,22 @@ Tie_formatting_problem::~Tie_formatting_problem () } void -Tie_formatting_problem::set_chord_outline (Link_array bounds, - Direction d) +Tie_formatting_problem::set_column_chord_outline (vector bounds, + Direction dir, + int column_rank) { Real staff_space = Staff_symbol_referencer::staff_space (bounds[0]); - Array boxes; + vector boxes; + vector head_boxes; Grob *stem = 0; - for (int i = 0; i < bounds.size (); i++) + for (vsize i = 0; i < bounds.size (); i++) { Grob *head = bounds[i]; if (!Note_head::has_interface (head)) continue; - + if (!stem) stem = unsmob_grob (head->get_object ("stem")); @@ -70,91 +95,184 @@ Tie_formatting_problem::set_chord_outline (Link_array bounds, (p+1) * 0.5 * staff_space); Interval x = head->extent (x_refpoint_, X_AXIS); - boxes.push (Box (x, y)); + head_boxes.push_back (Box (x, y)); + boxes.push_back (Box (x, y)); Grob *dots = Rhythmic_head::get_dots (head); - if (d == LEFT && dots) + if (dir == LEFT && dots) { Interval x = dots->extent (x_refpoint_, X_AXIS); - Interval y (-0.5, 0.5); int p = int (Staff_symbol_referencer::get_position (dots)); - y.translate (p); + /* + TODO: shouldn't this use column-rank dependent key? + */ dot_positions_.insert (p); dot_x_.unite (x); + + Interval y (dots->extent (dots, Y_AXIS)); + y.translate (p * staff_space * 0.5); - y *= staff_space * 0.5; - // boxes.push (Box (x, y)); + boxes.push_back (Box (x, y)); } } - chord_outlines_[d] = empty_skyline (-d); + Tuple2 key (column_rank, int (dir)); - if (bounds[0]->break_status_dir ()) + if (stem) { - Real x = robust_relative_extent (bounds[0], x_refpoint_, X_AXIS)[-d]; - chord_outlines_[d].elem_ref (0).height_ = x; - } + if (Stem::is_normal_stem (stem)) + { - for (int i = 0; i < boxes.size (); i++) - insert_extent_into_skyline (&chord_outlines_[d] , - boxes[i], Y_AXIS, -d); + Interval x; + x.add_point (stem->relative_coordinate (x_refpoint_, X_AXIS)); + x.widen (staff_space / 20); // ugh. + Interval y; + Real stem_end_position = + Stem::is_cross_staff (stem) + ? get_grob_direction (stem) * infinity_f + : Stem::stem_end_position (stem) * staff_space * .5; - if (stem - && !Stem::is_invisible (stem)) - { - Interval x; - x.add_point (stem->relative_coordinate (x_refpoint_, X_AXIS)); - x.widen (staff_space / 20); // ugh. - Interval y; - y.add_point (Stem::stem_end_position (stem) * staff_space * .5); + y.add_point (stem_end_position); - Direction stemdir = get_grob_direction (stem); - y.add_point (Stem::head_positions (stem)[-stemdir] - * staff_space * .5); - - insert_extent_into_skyline (&chord_outlines_[d], Box (x,y), Y_AXIS, -d); + Direction stemdir = get_grob_direction (stem); + y.add_point (Stem::head_positions (stem)[-stemdir] + * staff_space * .5); + /* + add extents of stem. + */ + boxes.push_back (Box (x, y)); + stem_extents_[key].unite (Box (x, y)); - if (d == LEFT) + if (dir == LEFT) + { + Box flag_box = Stem::get_translated_flag (stem).extent_box (); + flag_box.translate ( Offset (x[RIGHT], X_AXIS)); + boxes.push_back (flag_box); + } + } + else { - Box flag_box = Stem::get_translated_flag (stem).extent_box (); - flag_box.translate( Offset (x[RIGHT], X_AXIS)); - insert_extent_into_skyline (&chord_outlines_[d], flag_box, - Y_AXIS, -d); + Grob *head = Stem::support_head (stem); + + /* + In case of invisible stem, don't pass x-center of heads. + */ + Real x_center = head->extent (x_refpoint_, X_AXIS).center (); + Interval x_ext; + x_ext[-dir] = x_center; + Interval y_ext; + for (vsize j = 0; j < head_boxes.size (); j++) + y_ext.unite (head_boxes[j][Y_AXIS]); + + boxes.push_back (Box (x_ext, y_ext)); + } + + extract_grob_set (stem, "note-heads", heads); + for (vsize i = 0; i < heads.size (); i ++) + { + if (find (bounds.begin (), bounds.end (), dynamic_cast (heads[i])) == bounds.end ()) + { + /* + other untied notes in the same chord. + */ + + Interval y = Staff_symbol_referencer::extent_in_staff (heads[i]); + Interval x = heads[i]->extent (x_refpoint_, X_AXIS); + boxes.push_back (Box (x, y)); + } + + Grob *acc = unsmob_grob (heads[i]->get_object ("accidental-grob")); + if (acc && acc->is_live () && dir == RIGHT) + { + boxes.push_back (Box (acc->extent (x_refpoint_, X_AXIS), + Staff_symbol_referencer::extent_in_staff (acc))); + } + + head_positions_[column_rank].add_point (int (Staff_symbol_referencer::get_position (heads[i]))); } + } Direction updowndir = DOWN; do { - Interval x ; + Interval x; Interval y; - if (boxes.size()) + if (head_boxes.size ()) { - Box b = boxes.boundary (updowndir, 0); + Box b = boundary (head_boxes, updowndir, 0); x = b[X_AXIS]; - x[-d] = b[X_AXIS].linear_combination (-d / 2); + x[-dir] = b[X_AXIS].linear_combination (-dir / 2); y[-updowndir] = b[Y_AXIS][updowndir]; y[updowndir] = updowndir * infinity_f; } if (!x.is_empty ()) - insert_extent_into_skyline (&chord_outlines_[d], - Box (x,y), - Y_AXIS, -d); + boxes.push_back (Box (x, y)); } while (flip (&updowndir) != DOWN); + + /* todo: the horizon_padding is somewhat arbitrary */ + chord_outlines_[key] = Skyline (boxes, details_.skyline_padding_, Y_AXIS, -dir); + if (bounds[0]->break_status_dir ()) + { + Real x = robust_relative_extent (bounds[0], x_refpoint_, X_AXIS)[-dir]; + + chord_outlines_[key].set_minimum_height (x); + } + else + { + Interval x; + for (vsize j = 0; j < head_boxes.size (); j++) + { + x.unite (head_boxes[j][X_AXIS]); + } + + chord_outlines_[key].set_minimum_height (x[dir]); + } + + head_extents_[key].set_empty (); + for (vsize i = 0; i < head_boxes.size (); i++) + { + head_extents_[key].unite (head_boxes[i]); + } } +void +Tie_formatting_problem::set_chord_outline (vector bounds, + Direction dir) + +{ + vector ranks; + for (vsize i = 0; i < bounds.size (); i++) + ranks.push_back (bounds[i]->get_column ()->get_rank ()); + + vector_sort (ranks, less ()); + uniq (ranks); + + for (vsize i = 0; i < ranks.size (); i++) + { + vector col_items; + for (vsize j = 0; j < bounds.size (); j ++) + { + if (bounds[j]->get_column ()->get_rank () == ranks[i]) + col_items.push_back (bounds[j]); + } + + set_column_chord_outline (col_items, dir, ranks[i]); + } +} + + void Tie_formatting_problem::from_tie (Grob *tie) { - Link_array ties; - ties.push (tie); - + vector ties; + ties.push_back (tie); from_ties (ties); details_.from_grob (tie); @@ -167,13 +285,13 @@ Tie_formatting_problem::common_x_refpoint () const } void -Tie_formatting_problem::from_ties (Link_array const &ties) +Tie_formatting_problem::from_ties (vector const &ties) { - if (ties.is_empty ()) + if (ties.empty ()) return; x_refpoint_ = ties[0]; - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties.size (); i++) { x_refpoint_ = dynamic_cast (ties[i])->get_bound (LEFT)->common_refpoint (x_refpoint_, X_AXIS); x_refpoint_ = dynamic_cast (ties[i])->get_bound (RIGHT)->common_refpoint (x_refpoint_, X_AXIS); @@ -184,13 +302,24 @@ Tie_formatting_problem::from_ties (Link_array const &ties) Direction d = LEFT; do { - Link_array bounds; + vector bounds; - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties.size (); i++) { Item *it = dynamic_cast (ties[i])->get_bound (d); - - bounds.push (it); + if (it->break_status_dir ()) + { + Item *sep + = dynamic_cast (unsmob_grob (ties[i]->get_object ("separation-item"))); + if (sep && sep->get_column () == it->get_column ()) + it = sep; + + bounds.push_back (it); + } + else + { + bounds.push_back (it); + } } set_chord_outline (bounds, d); @@ -198,35 +327,35 @@ Tie_formatting_problem::from_ties (Link_array const &ties) while (flip (&d) != LEFT); - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties.size (); i++) { Tie_specification spec; + spec.from_grob (ties[i]); - if (scm_is_number (ties[i]->get_property_data (ly_symbol2scm ("direction")))) + do { - spec.manual_dir_ = to_dir (ties[i]->get_property ("direction")); - spec.has_manual_dir_ = true; + spec.note_head_drul_[d] = Tie::head (ties[i], d); + spec.column_ranks_[d] = Tie::get_column_rank (ties[i], d); } - - spec.position_ = Tie::get_position (ties[i]); - - specifications_.push (spec); + while (flip (&d) != LEFT); + specifications_.push_back (spec); } } void -Tie_formatting_problem::from_lv_ties (Link_array const &lv_ties) +Tie_formatting_problem::from_semi_ties (vector const &semi_ties, Direction head_dir) { - if (lv_ties.is_empty ()) - return ; - - details_.from_grob (lv_ties[0]); - Link_array heads; + if (semi_ties.empty ()) + return; - for (int i = 0; i < lv_ties.size (); i++) + details_.from_grob (semi_ties[0]); + vector heads; + + int column_rank = -1; + for (vsize i = 0; i < semi_ties.size (); i++) { Tie_specification spec; - Item *head = unsmob_item (lv_ties[i]->get_object ("note-head")); + Item *head = unsmob_item (semi_ties[i]->get_object ("note-head")); if (!head) programming_error ("LV tie without head?!"); @@ -235,141 +364,360 @@ Tie_formatting_problem::from_lv_ties (Link_array const &lv_ties) { spec.position_ = int (Staff_symbol_referencer::get_position (head)); } + + spec.from_grob (semi_ties[i]); - heads.push (head); - specifications_.push (spec); + spec.note_head_drul_[head_dir] = head; + column_rank = Tie::get_column_rank (semi_ties[i], head_dir); + spec.column_ranks_ = Drul_array (column_rank, column_rank); + heads.push_back (head); + specifications_.push_back (spec); } - x_refpoint_ = lv_ties [0]; - for (int i = 0; i < lv_ties.size (); i++) - { - x_refpoint_ = lv_ties[i]->common_refpoint (x_refpoint_, X_AXIS); - } + x_refpoint_ = semi_ties [0]; + for (vsize i = 0; i < semi_ties.size (); i++) + x_refpoint_ = semi_ties[i]->common_refpoint (x_refpoint_, X_AXIS); + for (vsize i = 0; i < heads.size (); i++) + x_refpoint_ = heads[i]->common_refpoint (x_refpoint_, X_AXIS); - set_chord_outline (heads, LEFT); + set_chord_outline (heads, head_dir); - Real right_most = - infinity_f; + Tuple2 head_key (column_rank, head_dir); + Tuple2 open_key (column_rank, -head_dir); + Real extremal = chord_outlines_[head_key].max_height (); - for (int i = 0; i < chord_outlines_[LEFT].size (); i++) - { - right_most = max (right_most, chord_outlines_[LEFT][i].height_); - } + chord_outlines_[open_key] = Skyline (head_dir); + chord_outlines_[open_key].set_minimum_height (extremal - head_dir * 1.5); +} - Skyline_entry right_entry; - right_entry.width_.set_full (); - right_entry.height_ = right_most + 1.5; - - chord_outlines_[RIGHT].push (right_entry); + +Tie_specification +Tie_formatting_problem::get_tie_specification (int i) const +{ + return specifications_[i]; } + +/* + Return configuration, create it if necessary. +*/ Tie_configuration* -Tie_formatting_problem::get_configuration (int pos, Direction dir) +Tie_formatting_problem::get_configuration (int pos, Direction dir, Drul_array columns, + bool tune_dy) const { - pair key (pos, dir); + int key_components[] = { + pos, dir, columns[LEFT], columns[RIGHT] + }; + Tuple key (key_components); + Tie_configuration_map::const_iterator f = possibilities_.find (key); - if (f != possibilities_.end ()) { return (*f).second; } - Tie_configuration *conf = generate_configuration (pos,dir); - possibilities_[key] = conf; + Tie_configuration *conf = generate_configuration (pos, dir, columns, tune_dy); + ((Tie_formatting_problem*) this)->possibilities_[key] = conf; return conf; } Tie_configuration* -Tie_formatting_problem::generate_configuration (int pos, Direction dir) const +Tie_formatting_problem::generate_configuration (int pos, Direction dir, + Drul_array columns, bool y_tune) const { Tie_configuration *conf = new Tie_configuration; conf->position_ = pos; conf->dir_ = dir; + + conf->column_ranks_ = columns; + Real y = conf->position_ * 0.5 * details_.staff_space_; if (dot_positions_.find (pos) != dot_positions_.end ()) { - conf->delta_y_ += 0.25 * details_.staff_space_; + conf->delta_y_ += dir * 0.25 * details_.staff_space_; + y_tune = false; } - conf->attachment_x_ = get_attachment (y + conf->delta_y_); + if (y_tune + && max (fabs (get_head_extent (columns[LEFT], LEFT, Y_AXIS)[dir] - y), + fabs (get_head_extent (columns[RIGHT], RIGHT, Y_AXIS)[dir] - y)) < 0.25 + && !Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos)) + { + conf->delta_y_ = + (get_head_extent (columns[LEFT], LEFT, Y_AXIS)[dir] - y) + + dir * details_.outer_tie_vertical_gap_; + } - Real h = conf->height (details_); - if (!conf->delta_y_) + if (y_tune) { - if (h < 0.5 * details_.staff_space_ - && !Staff_symbol_referencer::on_staffline (details_.staff_symbol_referencer_, pos)) - { - conf->center_tie_vertically (details_); - } - else if (h < 0.5 * details_.staff_space_ - && Staff_symbol_referencer::on_staffline (details_.staff_symbol_referencer_, pos)) + conf->attachment_x_ = get_attachment (y + conf->delta_y_, conf->column_ranks_); + Real h = conf->height (details_); + + /* + TODO: + + - should make sliding criterion, should flatten ties if + + - they're just the wrong (ie. touching line at top & bottom) + size. + + */ + if (head_positions_slice (columns[LEFT]).contains (pos) + || head_positions_slice (columns[RIGHT]).contains (pos) + || abs (pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_)) { - conf->delta_y_ += dir * 0.2 * details_.staff_space_; + if (h < details_.intra_space_threshold_ * 0.5 * details_.staff_space_) + { + if (!Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos) + && abs (pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_)) + { + conf->center_tie_vertically (details_); + } + else if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, pos)) + { + conf->delta_y_ += dir * + details_.tip_staff_line_clearance_ * 0.5 * details_.staff_space_; + } + } + else + { + Real top_y = y + conf->delta_y_ + conf->dir_ * h; + Real top_pos = top_y / (0.5*details_.staff_space_); + int round_pos = int (my_round (top_pos)); + + /* TODO: should use other variable? */ + Real clearance = details_.center_staff_line_clearance_; + if (fabs (top_pos - round_pos) < clearance + && Staff_symbol_referencer::on_staff_line (details_.staff_symbol_referencer_, + round_pos)) + { + Real new_y = (round_pos + clearance * conf->dir_) * 0.5 * details_.staff_space_; + conf->delta_y_ = (new_y - top_y); + } + } } + } + conf->attachment_x_ = get_attachment (y + conf->delta_y_, conf->column_ranks_); + if (conf->height (details_) < details_.intra_space_threshold_ * 0.5 * details_.staff_space_) + { + /* + This is less sensible for long ties, since those are more + horizontal. + */ + Interval close_by = get_attachment (y + + conf->delta_y_ + + (dir * details_.intra_space_threshold_ * 0.25 + * details_.staff_space_), + conf->column_ranks_); + + conf->attachment_x_.intersect (close_by); } - + conf->attachment_x_.widen ( - details_.x_gap_); + + if (conf->column_span_length ()) + { + /* + avoid the stems that we attach to as well. We don't do this + for semities (span length = 0) + + It would be better to check D against HEAD-DIRECTION if + applicable. + */ + Direction d = LEFT; + do + { + Real y = conf->position_ * details_.staff_space_ * 0.5 + conf->delta_y_; + if (get_stem_extent (conf->column_ranks_[d], d, X_AXIS).is_empty () + || !get_stem_extent (conf->column_ranks_[d], d, Y_AXIS).contains (y)) + continue; + + conf->attachment_x_[d] = + d * min (d * conf->attachment_x_[d], + d * (get_stem_extent (conf->column_ranks_[d], d, X_AXIS)[-d] - d * details_.stem_gap_)); + } + while (flip (&d) != LEFT); + } return conf; } +Interval +Tie_formatting_problem::get_head_extent (int col, Direction d, Axis a) const +{ + Column_extent_map::const_iterator i = head_extents_.find (Tuple2 (col, int (d))); + if (i != head_extents_.end ()) + return (*i).second[a]; + else + return Interval (); +} + +Interval +Tie_formatting_problem::get_stem_extent (int col, Direction d, Axis a) const +{ + Column_extent_map::const_iterator i = stem_extents_.find (Tuple2 (col, int (d))); + if (i != stem_extents_.end ()) + return (*i).second[a]; + else + return Interval (); +} + +/** + TIE_IDX and TIES_CONF are optional. + */ Real -Tie_formatting_problem::score_aptitude (Tie_configuration const &conf, - int tie_position) const -{ - Real wrong_direction_offset_penalty_; - Real distance_penalty_factor_; - - wrong_direction_offset_penalty_ = 10; - distance_penalty_factor_ = 5; - +Tie_formatting_problem::score_aptitude (Tie_configuration *conf, + Tie_specification const &spec, + Ties_configuration *ties_conf, int tie_idx) const +{ Real penalty = 0.0; - Real curve_y = conf.position_ * details_.staff_space_ * 0.5 + conf.delta_y_; - Real tie_y = tie_position * details_.staff_space_ * 0.5; - if (sign (curve_y - tie_y) != conf.dir_) - penalty += wrong_direction_offset_penalty_; + Real curve_y = conf->position_ * details_.staff_space_ * 0.5 + conf->delta_y_; + Real tie_y = spec.position_ * details_.staff_space_ * 0.5; + if (sign (curve_y - tie_y) != conf->dir_) + { + Real p = details_.wrong_direction_offset_penalty_; + if (ties_conf) + ties_conf->add_tie_score (p, tie_idx, "wrong dir"); + else + penalty += p; + } + + { + Real relevant_dist = max (fabs (curve_y - tie_y) - 0.5, 0.0); + Real p = details_.vertical_distance_penalty_factor_ * convex_amplifier (1.0, 0.9, relevant_dist); + if (ties_conf) + ties_conf->add_tie_score (p, tie_idx, "vdist"); + else + penalty += p; + } + + Direction d = LEFT; + do + { + if (!spec.note_head_drul_[d]) + continue; + + Interval head_x = spec.note_head_drul_[d]->extent (x_refpoint_, X_AXIS); + Real dist = head_x.distance (conf->attachment_x_[d]); + + + /* + TODO: flatten with log or sqrt. + */ + Real p = details_.horizontal_distance_penalty_factor_ + * convex_amplifier (1.25, 1.0, dist); + if (ties_conf) + ties_conf->add_tie_score (p, tie_idx, + (d == LEFT) ? "lhdist" : "rhdist"); + else + penalty += p; + + } + while (flip (&d) != LEFT); - penalty += distance_penalty_factor_ * fabs (curve_y - tie_y); + if (ties_conf + && ties_conf->size () == 1) + { + Direction d = LEFT; + Drul_array stems (0, 0); + do + { + if (!spec.note_head_drul_[d]) + continue; + + Grob *stem = unsmob_grob (spec.note_head_drul_[d]->get_object ("stem")); + if (stem + && Stem::is_normal_stem (stem)) + stems[d] = stem; + } + while (flip (&d) != LEFT); + + bool tie_stem_dir_ok = true; + bool tie_position_dir_ok = true; + if (stems[LEFT] && !stems[RIGHT]) + tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[LEFT]); + else if (!stems[LEFT] && stems[RIGHT]) + tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[RIGHT]); + else if (stems[LEFT] && stems[RIGHT] + && get_grob_direction (stems[LEFT]) == get_grob_direction (stems[RIGHT])) + tie_stem_dir_ok = conf->dir_ != get_grob_direction (stems[LEFT]); + else if (spec.position_) + tie_position_dir_ok = conf->dir_ == sign (spec.position_); + + if (!tie_stem_dir_ok) + ties_conf->add_score (details_.same_dir_as_stem_penalty_, "tie/stem dir"); + if (!tie_position_dir_ok) + ties_conf->add_score (details_.same_dir_as_stem_penalty_, "tie/pos dir"); + } + while (flip (&d) != LEFT); + + return penalty; } -Real -Tie_formatting_problem::score_configuration (Tie_configuration const &conf) const -{ - Real length_penalty_factor = 1.0; - Real min_length = 0.333; - Real staff_line_clearance = 0.1; - Real staff_line_collision_penalty = 5; - Real dot_collision_clearance = 0.25; - Real dot_collision_penalty = 10; - +Slice +Tie_formatting_problem::head_positions_slice (int rank) const +{ + Position_extent_map::const_iterator i (head_positions_.find (rank)); + if (i != head_positions_.end ()) + { + return (*i).second; + } + Slice empty; + return empty; +} - Real penalty = 0.0; - Real length = conf.attachment_x_.length (); - if (length < min_length) - penalty += length_penalty_factor / max (0.01, length); +/* + Score a configuration, ie. how well these ties looks without regard + to the note heads that they should connect to. + */ +void +Tie_formatting_problem::score_configuration (Tie_configuration *conf) const +{ + if (conf->scored_) + { + return ; + } + + Real length = conf->attachment_x_.length (); - Real tip_pos = conf.position_ + conf.delta_y_ / 0.5 * details_.staff_space_; + conf->add_score (details_.min_length_penalty_factor_ + * peak_around (0.33 * details_.min_length_, details_.min_length_, length), + "minlength"); + + Real tip_pos = conf->position_ + conf->delta_y_ / 0.5 * details_.staff_space_; Real tip_y = tip_pos * details_.staff_space_ * 0.5; - Real height = conf.height (details_); + Real height = conf->height (details_); - Real top_y = tip_y + conf.dir_ * height; + Real top_y = tip_y + conf->dir_ * height; Real top_pos = 2 * top_y / details_.staff_space_; Real round_top_pos = rint (top_pos); - if (fabs (top_pos - round_top_pos) < staff_line_clearance - && Staff_symbol_referencer::on_staffline (details_.staff_symbol_referencer_, + if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, int (round_top_pos)) && Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_) > top_y) { - penalty += staff_line_collision_penalty; + conf->add_score ( + details_.staff_line_collision_penalty_ + * peak_around (0.1 * details_.center_staff_line_clearance_, + details_.center_staff_line_clearance_, + fabs (top_pos - round_top_pos)), + "line center"); } - - if (fabs (tip_pos - rint (tip_pos)) < staff_line_clearance - && Staff_symbol_referencer::on_staffline (details_.staff_symbol_referencer_, - int (rint (tip_pos)))) + + int rounded_tip_pos = int (rint (tip_pos)); + if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, rounded_tip_pos) + && (head_positions_slice (conf->column_ranks_[LEFT]).contains (rounded_tip_pos) + || head_positions_slice (conf->column_ranks_[RIGHT]).contains (rounded_tip_pos) + || abs (rounded_tip_pos) < 2 * Staff_symbol_referencer::staff_radius (details_.staff_symbol_referencer_)) + ) { - penalty += staff_line_collision_penalty; + conf->add_score (details_.staff_line_collision_penalty_ + * peak_around (0.1 * details_.tip_staff_line_clearance_, + details_.tip_staff_line_clearance_, + fabs (tip_pos - rint (tip_pos))), + "tipline"); } if (!dot_x_.is_empty ()) @@ -377,7 +725,7 @@ Tie_formatting_problem::score_configuration (Tie_configuration const &conf) cons /* use left edge? */ Real x = dot_x_.center (); - Bezier b = conf.get_transformed_bezier (details_); + Bezier b = conf->get_transformed_bezier (details_); if (b.control_point_extent (X_AXIS).contains (x)) { Real y = b.get_other_coordinate (X_AXIS, x); @@ -386,104 +734,57 @@ Tie_formatting_problem::score_configuration (Tie_configuration const &conf) cons i != dot_positions_.end (); i ++) { int dot_pos = (*i); - if (fabs (dot_pos * details_.staff_space_ * 0.5 - y) < dot_collision_clearance) - { - penalty += dot_collision_penalty; - } + conf->add_score (details_.dot_collision_penalty_ + * peak_around (.1 * details_.dot_collision_clearance_, + details_.dot_collision_clearance_, + fabs (dot_pos * details_.staff_space_ * 0.5 - y)), + "dot collision"); } - } - } - - return penalty; -} - -Tie_configuration -Tie_formatting_problem::find_optimal_tie_configuration (int pos, Direction dir) const -{ - Link_array confs; - - int region_size = 3; - for (int i = 0; i < region_size; i ++) - { - confs.push (generate_configuration (pos + i * dir, dir)); - } - - Array scores; - - int best_idx = -1; - Real best_score = 1e6; - for (int i = 0; i < confs.size (); i ++) - { - Real score = 0.0; - score += score_configuration (*confs[i]); - score += score_aptitude (*confs[i], pos); - - if (score < best_score) - { - best_score = score; - best_idx = i; } } - Tie_configuration best = *confs[best_idx]; - for (int i = 0; i < confs.size (); i++) - delete confs[i]; - - return best; -} - -Tie_specification::Tie_specification () -{ - has_manual_position_ = false; - has_manual_dir_ = false; - position_ = 0; - manual_position_ = 0; - manual_dir_ = CENTER; + conf->scored_ = true; } - -Real -Tie_formatting_problem::score_ties_aptitude (Ties_configuration const &ties) const +void +Tie_formatting_problem::score_ties_aptitude (Ties_configuration *ties) const { - Real score = 0.0; - if (ties.size () != specifications_.size ()) + if (ties->size () != specifications_.size ()) { programming_error ("Huh? Mismatch between sizes."); - return infinity_f; + return; } - for (int i = 0; i < ties.size (); i++) - score += score_aptitude (ties[i], specifications_[i].position_); - - return score; + for (vsize i = 0; i < ties->size (); i++) + score_aptitude (&ties->at (i), specifications_[i], + ties, i); } -Real -Tie_formatting_problem::score_ties (Ties_configuration const &ties) const +void +Tie_formatting_problem::score_ties (Ties_configuration *ties) const { - return score_ties_configuration (ties) - + score_ties_aptitude (ties); + if (ties->scored_) + return; + + score_ties_configuration (ties); + score_ties_aptitude (ties); + ties->scored_ = true; } -Real -Tie_formatting_problem::score_ties_configuration (Ties_configuration const &ties) const +void +Tie_formatting_problem::score_ties_configuration (Ties_configuration *ties) const { - const Real tie_monotonicity_penalty = 100; - const Real tie_collision_penalty = 30; - const Real tie_collision_distance = 0.25; - - Real score = 0.0; - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties->size (); i++) { - score += score_configuration (ties[i]); + score_configuration (&ties->at (i)); + ties->add_tie_score (ties->at (i).score (), i, "conf"); } - - + Real last_edge = 0.0; Real last_center = 0.0; - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties->size (); i++) { - Bezier b (ties[i].get_transformed_bezier (details_)); + Bezier b (ties->at (i).get_transformed_bezier (details_)); Real center = b.curve_point (0.5)[Y_AXIS]; Real edge = b.curve_point (0.0)[Y_AXIS]; @@ -491,25 +792,42 @@ Tie_formatting_problem::score_ties_configuration (Ties_configuration const &ties if (i) { if (edge <= last_edge) - score += tie_monotonicity_penalty; + ties->add_score (details_.tie_column_monotonicity_penalty_, "monoton edge"); if (center <= last_center) - score += tie_monotonicity_penalty; - - score += - tie_collision_penalty - * max (tie_collision_distance - fabs (center - last_center), 0.0) / tie_collision_distance; - score += - tie_collision_penalty - * max (tie_collision_distance - fabs (edge - last_edge), 0.0) / tie_collision_distance; + ties->add_score (details_.tie_column_monotonicity_penalty_, "monoton cent"); + + ties->add_score (details_.tie_tie_collision_penalty_ * + peak_around (0.1 * details_.tie_tie_collision_distance_, + details_.tie_tie_collision_distance_, + fabs (center - last_center)), + "tietie center"); + ties->add_score (details_.tie_tie_collision_penalty_ * + peak_around (0.1 * details_.tie_tie_collision_distance_, + details_.tie_tie_collision_distance_, + fabs (edge - last_edge)), "tietie edge"); } last_edge = edge; last_center = center; } - return score; + if (ties->size () > 1) + { + ties->add_score (details_.outer_tie_length_symmetry_penalty_factor_ + * fabs (ties->at (0).attachment_x_.length () - ties->back ().attachment_x_.length ()), + "length symm"); + + ties->add_score (details_.outer_tie_vertical_distance_symmetry_penalty_factor_ + * fabs (fabs (specifications_[0].position_ * 0.5 * details_.staff_space_ + - (ties->at (0).position_ * 0.5 * details_.staff_space_ + + ties->at (0).delta_y_)) + - + fabs (specifications_.back ().position_ * 0.5 * details_.staff_space_ + - (ties->back ().position_ * 0.5 * details_.staff_space_ + + ties->back ().delta_y_))), + "pos symmetry"); + } } - /* Generate with correct X-attachments and beziers, copying delta_y_ from TIES_CONFIG if necessary. @@ -518,16 +836,18 @@ Ties_configuration Tie_formatting_problem::generate_ties_configuration (Ties_configuration const &ties_config) { Ties_configuration copy; - for (int i = 0; i < ties_config.size (); i++) + for (vsize i = 0; i < ties_config.size (); i++) { - Tie_configuration * ptr = get_configuration (ties_config[i].position_, ties_config[i].dir_); - if (specifications_[i].has_manual_position_) + Tie_configuration * ptr = get_configuration (ties_config[i].position_, ties_config[i].dir_, + ties_config[i].column_ranks_, + !specifications_[i].has_manual_delta_y_); + if (specifications_[i].has_manual_delta_y_) { ptr->delta_y_ = (specifications_[i].manual_position_ - ties_config[i].position_) * 0.5 * details_.staff_space_; } - copy.push (*ptr); + copy.push_back (*ptr); } return copy; @@ -537,26 +857,31 @@ Ties_configuration Tie_formatting_problem::generate_base_chord_configuration () { Ties_configuration ties_config; - - - for (int i = 0; i < specifications_.size (); i ++) + for (vsize i = 0; i < specifications_.size (); i ++) { Tie_configuration conf; if (specifications_[i].has_manual_dir_) conf.dir_ = specifications_[i].manual_dir_; if (specifications_[i].has_manual_position_) { - conf.position_ = (int) round (specifications_[i].manual_position_); - conf.delta_y_ = (specifications_[i].manual_position_ - conf.position_) - * 0.5 * details_.staff_space_; + conf.position_ = (int) my_round (specifications_[i].manual_position_); + if (specifications_[i].has_manual_delta_y_) + conf.delta_y_ = (specifications_[i].manual_position_ - conf.position_) + * 0.5 * details_.staff_space_; } else - conf.position_ = specifications_[i].position_; - - ties_config.push (conf); + { + conf.position_ = specifications_[i].position_; + } + + conf.column_ranks_ = specifications_[i].column_ranks_; + ties_config.push_back (conf); } set_ties_config_standard_directions (&ties_config); + for (vsize i = 0; i < ties_config.size (); i++) + if (!specifications_[i].manual_position_) + ties_config[i].position_ += ties_config[i].dir_; ties_config = generate_ties_configuration (ties_config); @@ -564,66 +889,103 @@ Tie_formatting_problem::generate_base_chord_configuration () } Ties_configuration -Tie_formatting_problem::generate_optimal_chord_configuration () +Tie_formatting_problem::find_best_variation (Ties_configuration const &base, + vector const &vars) { - Ties_configuration base = generate_base_chord_configuration (); - Array vars = get_variations (base); - Ties_configuration best = base; - Real best_score = score_ties_configuration (best); - + /* This simply is 1-opt: we have K substitions, and we try applying exactly every one for each. */ - for (int i = 0; i < vars.size (); i++) + for (vsize i = 0; i < vars.size (); i++) { - Ties_configuration variant = base; + Ties_configuration variant (base); variant[vars[i].index_] = *vars[i].suggestion_; - Real score = (score_ties_configuration (variant) - + score_ties_aptitude (variant)); - if (score < best_score) + variant.reset_score (); + score_ties (&variant); + + if (variant.score () < best.score ()) { best = variant; - best_score = score; } } return best; } + + + +Ties_configuration +Tie_formatting_problem::generate_optimal_configuration () +{ + Ties_configuration base = generate_base_chord_configuration (); + score_ties (&base); + + vector vars; + if (specifications_.size () > 1) + vars = generate_collision_variations (base); + else + vars = generate_single_tie_variations (base); + + Ties_configuration best = find_best_variation (base, vars); + + if (specifications_.size () > 1) + { + vars = generate_extremal_tie_variations (best); + best = find_best_variation (best, vars); + } + return best; +} void Tie_formatting_problem::set_ties_config_standard_directions (Ties_configuration *tie_configs) { - if (tie_configs->is_empty ()) + if (tie_configs->empty ()) return ; + + if (!tie_configs->at (0).dir_) + { + if (tie_configs->size () == 1) + tie_configs->at (0).dir_ = Direction (sign (tie_configs->at (0).position_)); + + if (!tie_configs->at (0).dir_) + tie_configs->at (0).dir_ = DOWN; + } - if (!tie_configs->elem (0).dir_) - tie_configs->elem_ref (0).dir_ = DOWN; - if (!tie_configs->top().dir_) - tie_configs->top().dir_ = UP; + if (!tie_configs->back ().dir_) + tie_configs->back ().dir_ = UP; /* Seconds */ - for (int i = 1; i < tie_configs->size (); i++) + for (vsize i = 1; i < tie_configs->size (); i++) { - Real diff = (tie_configs->elem (i-1).position_ - - tie_configs->elem (i).position_); - - if (fabs (diff) <= 1) + Real diff = (tie_configs->at (i).position_ + -tie_configs->at (i-1).position_); + + Real span_diff + = specifications_[i].column_span () - specifications_[i-1].column_span (); + if (span_diff && fabs (diff) <= 2) + { + if (span_diff > 0) + tie_configs->at (i).dir_ = UP; + else if (span_diff < 0) + tie_configs->at (i-1).dir_ = DOWN; + } + else if (fabs (diff) <= 1) { - if (!tie_configs->elem (i-1).dir_) - tie_configs->elem_ref (i-1).dir_ = DOWN; - if (!tie_configs->elem (i).dir_) - tie_configs->elem_ref (i).dir_ = UP; + if (!tie_configs->at (i-1).dir_) + tie_configs->at (i-1).dir_ = DOWN; + if (!tie_configs->at (i).dir_) + tie_configs->at (i).dir_ = UP; } } - for (int i = 1; i < tie_configs->size() - 1; i++) + for (vsize i = 1; i + 1 < tie_configs->size (); i++) { - Tie_configuration &conf = tie_configs->elem_ref (i); + Tie_configuration &conf = tie_configs->at (i); if (conf.dir_) continue; @@ -642,14 +1004,75 @@ Tie_configuration_variation::Tie_configuration_variation () suggestion_ = 0; } -Array -Tie_formatting_problem::get_variations (Ties_configuration const &ties) +vector +Tie_formatting_problem::generate_extremal_tie_variations (Ties_configuration const &ties) const +{ + vector vars; + Direction d = DOWN; + do + { + if (boundary (ties, d, 0).dir_ == d + && !boundary (specifications_, d, 0).has_manual_position_) + for (int i = 1; i <= details_.multi_tie_region_size_; i++) + { + Tie_configuration_variation var; + var.index_ = (d == DOWN) ? 0 : ties.size () - 1; + var.suggestion_ = get_configuration (boundary (ties, d, 0).position_ + + d * i, d, + boundary (ties, d, 0).column_ranks_, + true); + vars.push_back (var); + } + } + while (flip (&d) != DOWN); + + return vars; +} + +vector +Tie_formatting_problem::generate_single_tie_variations (Ties_configuration const &ties) const +{ + vector vars; + + int sz = details_.single_tie_region_size_; + if (specifications_[0].has_manual_position_) + sz = 1; + for (int i = 0; i < sz; i ++) + { + Direction d = LEFT; + do + { + if (i == 0 + && ties[0].dir_ == d) + continue; + + int p = ties[0].position_ + i * d; + + if (!specifications_[0].has_manual_dir_ + || d == specifications_[0].manual_dir_) + { + Tie_configuration_variation var; + var.index_ = 0; + var.suggestion_ = get_configuration (p, + d, specifications_[0].column_ranks_, + !specifications_[0].has_manual_delta_y_); + vars.push_back (var); + } + } + while (flip (&d) != LEFT); + } + return vars; +} + + +vector +Tie_formatting_problem::generate_collision_variations (Ties_configuration const &ties) const { Real center_distance_tolerance = 0.25; - Array vars; + vector vars; Real last_center = 0.0; - for (int i = 0; i < ties.size (); i++) + for (vsize i = 0; i < ties.size (); i++) { Bezier b (ties[i].get_transformed_bezier (details_)); @@ -663,56 +1086,121 @@ Tie_formatting_problem::get_variations (Ties_configuration const &ties) { Tie_configuration_variation var; var.index_ = i; - var.suggestion_ = get_configuration (ties[i].position_, - -ties[i].dir_); + var.suggestion_ = get_configuration (specifications_[i].position_ + - ties[i].dir_, + - ties[i].dir_, + + ties[i].column_ranks_, + !specifications_[i].has_manual_delta_y_ + ); - vars.push (var); + vars.push_back (var); } if (!specifications_[i-1].has_manual_dir_) { Tie_configuration_variation var; var.index_ = i-1; - var.suggestion_ = get_configuration (ties[i-1].position_, - -ties[i-1].dir_); + var.suggestion_ = get_configuration (specifications_[i-1].position_ + - ties[i-1].dir_, + - ties[i-1].dir_, + specifications_[i-1].column_ranks_, + !specifications_[i-1].has_manual_delta_y_ + ); + + vars.push_back (var); + } + + if (i == 1 && !specifications_[i-1].has_manual_position_ + && ties[i-1].dir_ == DOWN) + { + Tie_configuration_variation var; + var.index_ = i-1; + var.suggestion_ = get_configuration (specifications_[i-1].position_ - 1, DOWN, + specifications_[i-1].column_ranks_, + !specifications_[i-1].has_manual_delta_y_ - vars.push (var); + ); + vars.push_back (var); + } + if (i == ties.size () && !specifications_[i].has_manual_position_ + && ties[i].dir_ == UP) + { + Tie_configuration_variation var; + var.index_ = i; + var.suggestion_ = get_configuration (specifications_[i].position_ + + 1, UP, + specifications_[i].column_ranks_, + !specifications_[i].has_manual_delta_y_ + ); + vars.push_back (var); } } + else if (dot_positions_.find (ties[i].position_) != dot_positions_.end () + && !specifications_[i].has_manual_position_) + { + Tie_configuration_variation var; + var.index_ = i; + var.suggestion_ = get_configuration (ties[i].position_ + ties[i].dir_, + ties[i].dir_, + ties[i].column_ranks_, + !specifications_[i].has_manual_delta_y_ + ); + vars.push_back (var); + } + } last_center = center; } + return vars; - } void Tie_formatting_problem::set_manual_tie_configuration (SCM manual_configs) { - int k = 0; + vsize k = 0; for (SCM s = manual_configs; - scm_is_pair (s) && k < specifications_.size(); s = scm_cdr (s)) + scm_is_pair (s) && k < specifications_.size (); s = scm_cdr (s)) { SCM entry = scm_car (s); - if (!scm_is_pair (entry)) - continue; - - Tie_specification &spec = specifications_[k]; - - if (scm_is_number (scm_cdr (entry))) - { - spec.has_manual_dir_ = true; - spec.manual_dir_ = Direction (scm_to_int (scm_cdr (entry))); - } - if (scm_is_number (scm_car (entry))) + if (scm_is_pair (entry)) { - spec.has_manual_position_ = true; - spec.manual_position_ = scm_to_double (scm_car (entry)); - } + Tie_specification &spec = specifications_[k]; + + if (scm_is_number (scm_car (entry))) + { + spec.has_manual_position_ = true; + spec.manual_position_ = scm_to_double (scm_car (entry)); + spec.has_manual_delta_y_ = (scm_inexact_p (scm_car (entry)) == SCM_BOOL_T); + } + if (scm_is_number (scm_cdr (entry))) + { + spec.has_manual_dir_ = true; + spec.manual_dir_ = Direction (scm_to_int (scm_cdr (entry))); + } + } k ++; } } + +void +Tie_formatting_problem::set_debug_scoring (Ties_configuration const &base) +{ +#if DEBUG_TIE_SCORING + if (to_boolean (x_refpoint_->layout () + ->lookup_variable (ly_symbol2scm ("debug-tie-scoring")))) + { + for (vsize i = 0; i < base.size (); i++) + { + string card = base.complete_tie_card (i); + specifications_[i].tie_grob_->set_property ("quant-score", + ly_string2scm (card)); + } + } +#endif +}