X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=lily%2Ftie-formatting-problem.cc;h=0e29d3ca50a1f028c7c1f87270e333bf31538208;hb=0c289671d105bea58592607cb508aef59aee882b;hp=36c81fe0a2b1a4d25d0b42dc6e5d043814d37509;hpb=94af2fb150edb6889a47bf5b4052a966f115e12a;p=lilypond.git diff --git a/lily/tie-formatting-problem.cc b/lily/tie-formatting-problem.cc index 36c81fe0a2..0e29d3ca50 100644 --- a/lily/tie-formatting-problem.cc +++ b/lily/tie-formatting-problem.cc @@ -1,33 +1,57 @@ /* - 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--2006 Han-Wen Nienhuys */ #include "tie-formatting-problem.hh" +#include "paper-column.hh" +#include "bezier.hh" #include "directional-element-interface.hh" +#include "item.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 "item.hh" -#include "spanner.hh" -#include "stem.hh" -#include "note-head.hh" -#include "rhythmic-head.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 ("Can't find chord outline"); + else + attachments[d] = i->second.height (y); } while (flip (&d) != LEFT); @@ -47,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")); @@ -69,34 +95,26 @@ 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); 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); - - if (bounds[0]->break_status_dir ()) - { - Real x = robust_relative_extent (bounds[0], x_refpoint_, X_AXIS)[-d]; - chord_outlines_[d].elem_ref (0).height_ = x; - } - - for (int i = 0; i < boxes.size (); i++) - insert_extent_into_skyline (&chord_outlines_[d] , - boxes[i], Y_AXIS, -d); + Tuple2 key (column_rank, int (dir)); if (stem && !Stem::is_invisible (stem)) @@ -110,49 +128,110 @@ Tie_formatting_problem::set_chord_outline (Link_array bounds, 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); + 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)); - insert_extent_into_skyline (&chord_outlines_[d], flag_box, - Y_AXIS, -d); + boxes.push_back (flag_box); } } + else if (stem) + { + 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)); + } 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); + + chord_outlines_[key] = Skyline (boxes, 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); @@ -165,13 +244,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); @@ -182,188 +261,409 @@ 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); + bounds.push_back (it); } set_chord_outline (bounds, d); } while (flip (&d) != LEFT); + + + for (vsize i = 0; i < ties.size (); i++) + { + Tie_specification spec; + + spec.get_tie_manual_settings (ties[i]); + + + + do + { + spec.note_head_drul_[d] = Tie::head (ties[i], d); + spec.column_ranks_[d] = + dynamic_cast (ties[i])->get_bound (d)->get_column ()->get_rank (); + } + 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 ; + if (semi_ties.empty ()) + return; - details_.from_grob (lv_ties[0]); - Link_array heads; - 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++) { - Item *head = unsmob_item (lv_ties[i]->get_object ("note-head")); + Tie_specification spec; + Item *head = unsmob_item (semi_ties[i]->get_object ("note-head")); + if (!head) - continue; + programming_error ("LV tie without head?!"); + + if (head) + { + spec.position_ = int (Staff_symbol_referencer::get_position (head)); + } + + spec.get_tie_manual_settings (semi_ties[i]); - heads.push (head); + spec.note_head_drul_[head_dir] = head; + column_rank = dynamic_cast (head)->get_column ()->get_rank (); + 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) 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); + ((Tie_formatting_problem*) this)->possibilities_[key] = conf; return conf; } Tie_configuration* -Tie_formatting_problem::generate_configuration (int pos, Direction dir) +Tie_formatting_problem::generate_configuration (int pos, Direction dir, + Drul_array columns) 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_; + bool y_tune = true; 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; + } + + 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_; } - conf->attachment_x_ = get_attachment (y + conf->delta_y_); - 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->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 (h < details_.intra_space_threshold_ * 0.5 * details_.staff_space_) { - conf->center_tie_vertically (details_); + 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 if (h < 0.5 * details_.staff_space_ - && Staff_symbol_referencer::on_staffline (details_.staff_symbol_referencer_, pos)) + else { - conf->delta_y_ += dir * 0.2 * details_.staff_space_; + 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) -{ - 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 p = details_.vertical_distance_penalty_factor_ * fabs (curve_y - tie_y); + 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_ * 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); return penalty; } - -Real -Tie_formatting_problem::score_configuration (Tie_configuration const &conf) +void +Tie_formatting_problem::score_configuration (Tie_configuration *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 penalty = 0.0; - Real length = conf.attachment_x_.length (); - if (length < min_length) - penalty += length_penalty_factor / max (0.01, length); + 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)))) + if (Staff_symbol_referencer::on_line (details_.staff_symbol_referencer_, + int (rint (tip_pos)))) { - 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"); } - - return penalty; + + if (!dot_x_.is_empty ()) + { + /* use left edge? */ + Real x = dot_x_.center (); + + 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); + + for (set::const_iterator i (dot_positions_.begin ()); + i != dot_positions_.end (); i ++) + { + int dot_pos = (*i); + 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"); + } + } + } + + conf->scored_ = true; } Tie_configuration -Tie_formatting_problem::find_optimal_tie_configuration (int pos, Direction dir) +Tie_formatting_problem::find_optimal_tie_configuration (Tie_specification const &spec) const { - Link_array confs; + vector confs; + + int pos = spec.position_; + Direction dir = spec.manual_dir_; - int region_size = 3; - for (int i = 0; i < region_size; i ++) + for (int i = 0; i < details_.single_tie_region_size_; i ++) { - confs.push (generate_configuration (pos + i * dir, dir)); + confs.push_back (generate_configuration (pos + i * dir, dir, + spec.column_ranks_)); + + if (spec.has_manual_position_) + { + confs.back ()->delta_y_ + = (spec.manual_position_ - spec.position_) + * 0.5 * details_.staff_space_; + + break; + } } - Array scores; + vector scores; int best_idx = -1; Real best_score = 1e6; - for (int i = 0; i < confs.size (); i ++) + for (vsize i = 0; i < confs.size (); i ++) { - Real score = 0.0; - score += score_configuration (*confs[i]); - score += score_aptitude (*confs[i], pos); + score_configuration (confs[i]); + Real score = score_aptitude (confs[i], spec, 0, 0) + + confs[i]->score (); if (score < best_score) { @@ -372,9 +672,434 @@ Tie_formatting_problem::find_optimal_tie_configuration (int pos, Direction dir) } } - Tie_configuration best = *confs[best_idx]; - for (int i = 0; i < confs.size (); i++) + if (best_idx < 0) + programming_error ("No best tie configuration found."); + + Tie_configuration best + = (best_idx >= 0) ? *confs[best_idx] : *confs[0]; + + for (vsize 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; + note_head_drul_[LEFT] = + note_head_drul_[RIGHT] = 0; + column_ranks_[RIGHT] = + column_ranks_[LEFT] = 0; +} + + +void +Tie_specification::get_tie_manual_settings (Grob *tie) +{ + if (scm_is_number (tie->get_property_data ("direction"))) + { + manual_dir_ = to_dir (tie->get_property ("direction")); + has_manual_dir_ = true; + } + + position_ = Tie::get_position (tie); + if (scm_is_number (tie->get_property ("staff-position"))) + { + manual_position_ = scm_to_double (tie->get_property ("staff-position")); + has_manual_position_ = true; + position_ = int (my_round (manual_position_)); + } +} + +int +Tie_specification::column_span () const +{ + return column_ranks_[RIGHT] - column_ranks_[LEFT]; +} + +void +Tie_formatting_problem::score_ties_aptitude (Ties_configuration *ties) const +{ + if (ties->size () != specifications_.size ()) + { + programming_error ("Huh? Mismatch between sizes."); + return; + } + + for (vsize i = 0; i < ties->size (); i++) + score_aptitude (&ties->at (i), specifications_[i], + ties, i); +} + +void +Tie_formatting_problem::score_ties (Ties_configuration *ties) const +{ + if (ties->scored_) + return; + + score_ties_configuration (ties); + score_ties_aptitude (ties); + ties->scored_ = true; +} + +void +Tie_formatting_problem::score_ties_configuration (Ties_configuration *ties) const +{ + for (vsize i = 0; i < ties->size (); 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 (vsize i = 0; i < ties->size (); i++) + { + 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]; + + if (i) + { + if (edge <= last_edge) + ties->add_score (details_.tie_column_monotonicity_penalty_, "monoton edge"); + if (center <= last_center) + 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; + } + + 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. +*/ +Ties_configuration +Tie_formatting_problem::generate_ties_configuration (Ties_configuration const &ties_config) +{ + Ties_configuration copy; + for (vsize i = 0; i < ties_config.size (); i++) + { + Tie_configuration * ptr = get_configuration (ties_config[i].position_, ties_config[i].dir_, + ties_config[i].column_ranks_); + if (specifications_[i].has_manual_position_) + { + ptr->delta_y_ + = (specifications_[i].manual_position_ - ties_config[i].position_) + * 0.5 * details_.staff_space_; + } + copy.push_back (*ptr); + } + + return copy; +} + +Ties_configuration +Tie_formatting_problem::generate_base_chord_configuration () +{ + Ties_configuration ties_config; + 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) my_round (specifications_[i].manual_position_); + conf.delta_y_ = (specifications_[i].manual_position_ - conf.position_) + * 0.5 * details_.staff_space_; + } + else + { + 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); + + return ties_config; +} + +Ties_configuration +Tie_formatting_problem::find_best_variation (Ties_configuration const &base, + vector vars) +{ + Ties_configuration best = base; + + /* + This simply is 1-opt: we have K substitions, and we try applying + exactly every one for each. + */ + for (vsize i = 0; i < vars.size (); i++) + { + Ties_configuration variant (base); + variant[vars[i].index_] = *vars[i].suggestion_; + + variant.reset_score (); + score_ties (&variant); + + if (variant.score () < best.score ()) + { + best = variant; + } + } + + return best; +} + + + +Ties_configuration +Tie_formatting_problem::generate_optimal_chord_configuration () +{ + Ties_configuration base = generate_base_chord_configuration (); + vector vars = generate_collision_variations (base); + + score_ties (&base); + Ties_configuration best = find_best_variation (base, vars); + 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->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->back ().dir_) + tie_configs->back ().dir_ = UP; + + /* + Seconds + */ + for (vsize i = 1; i < tie_configs->size (); i++) + { + 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->at (i-1).dir_) + tie_configs->at (i-1).dir_ = DOWN; + if (!tie_configs->at (i).dir_) + tie_configs->at (i).dir_ = UP; + } + } + + for (vsize i = 1; i < tie_configs->size() - 1; i++) + { + Tie_configuration &conf = tie_configs->at (i); + if (conf.dir_) + continue; + + Direction position_dir = + Direction (sign (conf.position_)); + if (!position_dir) + position_dir = DOWN; + + conf.dir_ = position_dir; + } +} + +Tie_configuration_variation::Tie_configuration_variation () +{ + index_ = 0; + suggestion_ = 0; +} + +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_); + vars.push_back (var); + } + } + while (flip (&d) != DOWN); + + return vars; +} + + +vector +Tie_formatting_problem::generate_collision_variations (Ties_configuration const &ties) const +{ + Real center_distance_tolerance = 0.25; + + vector vars; + Real last_center = 0.0; + for (vsize i = 0; i < ties.size (); i++) + { + Bezier b (ties[i].get_transformed_bezier (details_)); + + Real center = b.curve_point (0.5)[Y_AXIS]; + + if (i) + { + if (center <= last_center + center_distance_tolerance) + { + if (!specifications_[i].has_manual_dir_) + { + Tie_configuration_variation var; + var.index_ = i; + var.suggestion_ = get_configuration (specifications_[i].position_ + - ties[i].dir_, + - ties[i].dir_, + + ties[i].column_ranks_ + ); + + vars.push_back (var); + } + + if (!specifications_[i-1].has_manual_dir_) + { + Tie_configuration_variation var; + var.index_ = i-1; + var.suggestion_ = get_configuration (specifications_[i-1].position_ + - ties[i-1].dir_, + - ties[i-1].dir_, + specifications_[i-1].column_ranks_); + + 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_); + 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_); + 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_); + vars.push_back (var); + } + + } + + last_center = center; + } + + + return vars; +} + +void +Tie_formatting_problem::set_manual_tie_configuration (SCM manual_configs) +{ + vsize k = 0; + for (SCM s = manual_configs; + scm_is_pair (s) && k < specifications_.size (); s = scm_cdr (s)) + { + SCM entry = scm_car (s); + if (scm_is_pair (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)); + } + if (scm_is_number (scm_cdr (entry))) + { + spec.has_manual_dir_ = true; + spec.manual_dir_ = Direction (scm_to_int (scm_cdr (entry))); + } + } + k ++; + } +} +