2 slur.cc -- implement Slur
4 source file of the GNU LilyPond music typesetter
6 (c) 1996--2000 Han-Wen Nienhuys <hanwen@cs.uu.nl>
7 Jan Nieuwenhuizen <janneke@gnu.org>
12 * begin and end should be treated as a/acknowledge Scripts.
13 * broken slur should have uniform trend
16 #include "directional-element-interface.hh"
17 #include "group-interface.hh"
20 #include "paper-def.hh"
21 #include "note-column.hh"
23 #include "paper-column.hh"
24 #include "molecule.hh"
27 #include "bezier-bow.hh"
29 #include "cross-staff.hh"
30 #include "group-interface.hh"
31 #include "staff-symbol-referencer.hh"
35 class Slur_bezier_bow : public Bezier_bow
38 Slur_bezier_bow (Array<Offset> encompass, Direction dir);
39 Array<Real> area_x_gradients_array (Real area);
41 Real enclosed_area_f () const;
42 Real fit_factor () const;
43 void minimise_enclosed_area (Paper_def* paper_l, Real default_height);
46 Slur_bezier_bow::Slur_bezier_bow (Array<Offset> encompass, Direction dir)
47 : Bezier_bow (encompass, dir)
52 Slur_bezier_bow::blow_fit ()
54 Real len = curve_.control_[3][X_AXIS];
55 Real h = curve_.control_[1][Y_AXIS] * fit_factor () / len;
56 curve_.control_[1][Y_AXIS] = h * len;
57 curve_.control_[2][Y_AXIS] = h * len;
58 curve_.assert_sanity ();
63 Slur_bezier_bow::enclosed_area_f () const
66 for (int i=0; i < encompass_.size (); i++)
72 x = Interval (0, encompass_[1][X_AXIS] / 2);
74 curve_.get_other_coordinate (X_AXIS,
78 else if (i == encompass_.size () - 1)
80 x = Interval ((encompass_[i-1][X_AXIS] + encompass_[i][X_AXIS])/2,
81 encompass_[i][X_AXIS]);
83 (curve_.get_other_coordinate (X_AXIS,
84 (x[MIN] + x[MAX]) / 2)));
88 x = Interval ((encompass_[i-1][X_AXIS] + encompass_[i][X_AXIS]) / 2,
89 (encompass_[i][X_AXIS] + encompass_[i+1][X_AXIS]) / 2);
90 y = Interval (encompass_[i][Y_AXIS],
91 (curve_.get_other_coordinate (X_AXIS, x[MIN])
92 + curve_.get_other_coordinate (X_AXIS,
93 (x[MIN] + x[MAX]) / 2)
94 + curve_.get_other_coordinate (X_AXIS, x[MAX])) / 3);
97 Real da = x.length () * y.length ();
104 Slur_bezier_bow::area_x_gradients_array (Real area)
106 Real len = curve_.control_[3][X_AXIS];
107 Real grow = len / 10.0;
109 for (int i=0; i < 2; i++)
111 Real r = curve_.control_[i+1][X_AXIS];
112 curve_.control_[i+1][X_AXIS] += grow;
113 da[i] = (enclosed_area_f () - area) / grow;
114 curve_.control_[i+1][X_AXIS] = r;
120 Slur_bezier_bow::minimise_enclosed_area (Paper_def* paper_l,
123 Real length = curve_.control_[3][X_AXIS];
124 Real sb = paper_l->get_var ("slur_beautiful");
125 Real beautiful = length * default_height * sb;
127 DEBUG_OUT << to_str ("Beautiful: %f\n", beautiful);
128 DEBUG_OUT << to_str ("Length: %f\n", length);
129 DEBUG_OUT << to_str ("D-height: %f\n", default_height);
130 DEBUG_OUT << to_str ("FitFac: %f\n", fit_factor ());
132 if (fit_factor () > 1.0)
135 Real pct_c0 = paper_l->get_var ("bezier_pct_c0");
136 Real pct_c3 = paper_l->get_var ("bezier_pct_c3");
137 Real pct_in_max = paper_l->get_var ("bezier_pct_in_max");
138 Real pct_out_max = paper_l->get_var ("bezier_pct_out_max");
139 Real steps = paper_l->get_var ("bezier_area_steps");
141 for (int i=0; i < steps; i++)
143 Real area = enclosed_area_f ();
145 DEBUG_OUT << to_str ("Init area: %f\n", area);
147 if (area <= beautiful)
150 Array<Real> da = area_x_gradients_array (area);
153 Real pct = pct_c0 + pct_c3 * length * length * length;
154 pct *= (steps - i) / steps;
155 if (da[0] > 0 || da[1] < 0)
156 pct = pct <? pct_out_max;
158 pct = pct <? pct_in_max;
160 Real u = (abs (curve_.control_[1][X_AXIS] / da[0])
161 <? abs ((curve_.control_[3][X_AXIS]
162 - curve_.control_[2][X_AXIS]) / da[1]));
164 DEBUG_OUT << to_str ("pct: %f\n", pct);
165 DEBUG_OUT << to_str ("u: %f\n", u);
167 DEBUG_OUT << to_str ("da: (%f, %f)\n", da[0], da[1]);
168 DEBUG_OUT << to_str ("da*u: (%f, %f)\n", da[0]*u*pct, da[1]*u*pct);
169 DEBUG_OUT << to_str ("cx: (%f, %f)\n", curve_.control_[1][X_AXIS],
170 curve_.control_[2][X_AXIS]);
172 curve_.control_[1][X_AXIS] -= da[0] * u * pct;
173 curve_.control_[2][X_AXIS] -= da[1] * u * pct;
176 Real area = enclosed_area_f ();
177 DEBUG_OUT << to_str ("Exarea: %f\n", area);
183 max ( encompass.y / curve.y )
187 Slur_bezier_bow::fit_factor () const
189 Real x1 = encompass_[0][X_AXIS];
190 Real x2 = encompass_.top ()[X_AXIS];
193 for (int i=1; i < encompass_.size ()-1; i++)
195 if (encompass_[i][X_AXIS] > x1 && encompass_[i][X_AXIS] < x2)
197 Real y = curve_.get_other_coordinate (X_AXIS, encompass_[i][X_AXIS]);
200 Real f = encompass_[i][Y_AXIS] / y;
201 factor = factor >? f;
221 dy_f_drul_[LEFT] = dy_f_drul_[RIGHT] = 0.0;
222 dx_f_drul_[LEFT] = dx_f_drul_[RIGHT] = 0.0;
224 set_elt_pointer ("note-columns", SCM_EOL);
225 set_elt_property ("control-points", SCM_EOL);
229 Slur::add_column (Note_column*n)
231 if (!gh_pair_p (n->get_elt_pointer ("note-heads")))
232 warning (_ ("Putting slur over rest. Ignoring."));
235 Pointer_group_interface (this, "note-columns").add_element (n);
241 Slur::de_uglyfy (Slur_bezier_bow* bb, Real default_height)
243 Real length = bb->curve_.control_[3][X_AXIS] ;
244 Real ff = bb->fit_factor ();
245 for (int i = 1; i < 3; i++)
247 Real ind = abs (bb->curve_.control_[(i-1)*3][X_AXIS]
248 - bb->curve_.control_[i][X_AXIS]) / length;
249 Real h = bb->curve_.control_[i][Y_AXIS] * ff / length;
251 Real f = default_height / length;
252 Real c1 = paper_l ()->get_var ("bezier_control1");
253 Real c2 = paper_l ()->get_var ("bezier_control2");
254 Real c3 = paper_l ()->get_var ("bezier_control3");
259 else if (h > c2 + c3 * ind)
264 bb->curve_.control_[i][Y_AXIS] = h * length;
267 bb->curve_.assert_sanity ();
271 Slur::get_default_dir () const
273 Link_array<Note_column> encompass_arr =
274 Pointer_group_interface__extract_elements (this, (Note_column*)0, "note-columns");
277 for (int i=0; i < encompass_arr.size (); i ++)
279 if (encompass_arr[i]->dir () < 0)
289 Slur::do_add_processing ()
291 Link_array<Note_column> encompass_arr =
292 Pointer_group_interface__extract_elements (this, (Note_column*)0, "note-columns");
294 if (encompass_arr.size ())
296 set_bound (LEFT, encompass_arr[0]);
297 if (encompass_arr.size () > 1)
298 set_bound (RIGHT, encompass_arr.top ());
305 Slur::encompass_offset (Note_column const* col) const
308 Stem* stem_l = col->stem_l ();
309 Direction dir = directional_element (this).get ();
313 warning (_ ("Slur over rest?"));
314 o[X_AXIS] = col->relative_coordinate (0, X_AXIS);
315 o[Y_AXIS] = col->extent (Y_AXIS)[dir];
318 Direction stem_dir = directional_element (stem_l).get ();
319 o[X_AXIS] = stem_l->relative_coordinate (0, X_AXIS);
322 Simply set x to middle of notehead
325 o[X_AXIS] -= 0.5 * stem_dir * col->extent (X_AXIS).length ();
327 if ((stem_dir == dir)
328 && !stem_l->extent (Y_AXIS).empty_b ())
330 o[Y_AXIS] = stem_l->extent (Y_AXIS)[dir];
334 o[Y_AXIS] = col->extent (Y_AXIS)[dir];
338 leave a gap: slur mustn't touch head/stem
340 o[Y_AXIS] += dir * paper_l ()->get_var ("slur_y_free");
341 o[Y_AXIS] -= calc_interstaff_dist (stem_l, this);
346 Slur::after_line_breaking ()
349 set_control_points ();
357 Slur::set_extremities ()
359 Link_array<Note_column> encompass_arr =
360 Pointer_group_interface__extract_elements (this, (Note_column*)0, "note-columns");
362 if (!encompass_arr.size ())
364 set_elt_property ("transparent", SCM_BOOL_T);
365 set_extent_callback (0, X_AXIS);
366 set_extent_callback (0, Y_AXIS);
370 if (!directional_element (this).get ())
371 directional_element (this).set (get_default_dir ());
375 Slur and tie placement [OSU]
378 * x = centre of head - d * x_gap_f
381 * y = length < 5ss : horizontal tangent + d * 0.25 ss
382 y = length >= 5ss : y next interline - d * 0.25 ss
385 Real staff_space = paper_l ()->get_var ("interline");
386 Real half_staff_space = staff_space / 2;
388 Real x_gap_f = paper_l ()->get_var ("slur_x_gap");
389 Real y_gap_f = paper_l ()->get_var ("slur_y_gap");
391 Drul_array<Note_column*> note_column_drul;
392 note_column_drul[LEFT] = encompass_arr[0];
393 note_column_drul[RIGHT] = encompass_arr.top ();
395 bool fix_broken_b = false;
397 Direction my_dir = directional_element (this).get ();
405 if ((note_column_drul[d] == get_bound (d))
406 && note_column_drul[d]->first_head ()
407 && (note_column_drul[d]->stem_l ()))
409 Stem* stem_l = note_column_drul[d]->stem_l ();
411 side directly attached to note head;
412 no beam getting in the way
414 if ((stem_l->extent (Y_AXIS).empty_b ()
415 || !((stem_l->get_direction () == my_dir) && (my_dir != d)))
416 && !((my_dir == stem_l->get_direction ())
417 && stem_l->beam_l () && (stem_l->beam_count (-d) >= 1)))
419 dx_f_drul_[d] = get_bound (d)->extent (X_AXIS).length () / 2;
420 dx_f_drul_[d] -= d * x_gap_f;
422 if (stem_l->get_direction () != my_dir)
424 dy_f_drul_[d] = note_column_drul[d]->extent (Y_AXIS)[my_dir];
428 dy_f_drul_[d] = stem_l->chord_start_f ()
429 + my_dir * half_staff_space;
431 dy_f_drul_[d] += my_dir * y_gap_f;
434 side attached to (visible) stem
438 dx_f_drul_[d] = stem_l->relative_coordinate (0, X_AXIS)
439 - get_bound (d)->relative_coordinate (0, X_AXIS);
441 side attached to beamed stem
443 if (stem_l->beam_l () && (stem_l->beam_count (-d) >= 1))
445 dy_f_drul_[d] = stem_l->extent (Y_AXIS)[my_dir];
446 dy_f_drul_[d] += my_dir * 2 * y_gap_f;
449 side attached to notehead, with stem getting in the way
453 dx_f_drul_[d] -= d * x_gap_f;
455 dy_f_drul_[d] = stem_l->chord_start_f ()
456 + my_dir * half_staff_space;
457 dy_f_drul_[d] += my_dir * y_gap_f;
466 dx_f_drul_[d] = get_broken_left_end_align ();
469 broken: should get y from other piece, so that slur
470 continues up/down trend
472 for now: be horizontal..
477 while (flip (&d) != LEFT);
479 int cross_count = cross_staff_count ();
480 bool interstaff_b = (0 < cross_count) && (cross_count < encompass_arr.size ());
482 Drul_array<Offset> info_drul;
483 Drul_array<Real> interstaff_interval;
487 info_drul[d] = encompass_offset (encompass_arr.boundary (d, 0));
488 interstaff_interval[d] = - calc_interstaff_dist (encompass_arr.boundary (d,0),
491 while (flip (&d) != LEFT);
493 Real interstaff_f = interstaff_interval[RIGHT] - interstaff_interval[LEFT];
497 Direction d = (encompass_arr.top () != get_bound (RIGHT)) ?
499 dy_f_drul_[d] = info_drul[d][Y_AXIS];
502 dy_f_drul_[d] -= interstaff_interval[d];
503 if (cross_count) // interstaff_i ?
505 dy_f_drul_[LEFT] += interstaff_interval[d];
506 dy_f_drul_[RIGHT] += interstaff_interval[d];
512 dy_f_drul_[RIGHT] += interstaff_f;
517 Slur::cross_staff_count ()const
519 Link_array<Note_column> encompass_arr =
520 Pointer_group_interface__extract_elements (this, (Note_column*)0, "note-columns");
524 for (int i = 0; i < encompass_arr.size (); i++)
526 if (calc_interstaff_dist (encompass_arr[i], this))
534 Slur::get_encompass_offset_arr () const
536 Link_array<Note_column> encompass_arr =
537 Pointer_group_interface__extract_elements (this, (Note_column*)0, "note-columns");
539 Array<Offset> offset_arr;
543 check non-disturbed slur
544 FIXME: x of ends off by a tiny bit!!
546 offset_arr.push (Offset (0, dy_f_drul_[LEFT]));
547 offset_arr.push (Offset (0, dy_f_drul_[RIGHT]));
551 Offset origin (relative_coordinate (0, X_AXIS), 0);
554 int last = encompass_arr.size () - 2;
556 offset_arr.push (Offset (dx_f_drul_[LEFT], dy_f_drul_[LEFT]));
562 int cross_count = cross_staff_count ();
563 bool cross_b = cross_count && cross_count < encompass_arr.size ();
564 if (encompass_arr[0] != get_bound (LEFT))
567 Real is = calc_interstaff_dist (encompass_arr[0], this);
569 offset_arr[0][Y_AXIS] += is;
575 if (encompass_arr.top () != get_bound (RIGHT))
580 for (int i = first; i <= last; i++)
582 Offset o (encompass_offset (encompass_arr[i]));
583 offset_arr.push (o - origin);
586 offset_arr.push (Offset (spanner_length ()+ dx_f_drul_[RIGHT],
594 Slur::get_rods () const
599 r.item_l_drul_[LEFT] = get_bound (LEFT);
600 r.item_l_drul_[RIGHT] = get_bound (RIGHT);
601 r.distance_f_ = paper_l ()->get_var ("slur_x_minimum");
610 Ugh should have dash-length + dash-period
613 Slur::do_brew_molecule () const
615 Real thick = paper_l ()->get_var ("slur_thickness");
616 Bezier one = get_curve ();
619 SCM d = get_elt_property ("dashed");
621 a = lookup_l ()->dashed_slur (one, thick, thick * gh_scm2double (d));
623 a = lookup_l ()->slur (one, directional_element (this).get () * thick, thick);
629 Slur::set_control_points ()
631 Slur_bezier_bow bb (get_encompass_offset_arr (),
632 directional_element (this).get ());
634 Real staff_space = Staff_symbol_referencer_interface (this).staff_space ();
635 Real h_inf = paper_l ()->get_var ("slur_height_limit_factor") * staff_space;
636 Real r_0 = paper_l ()->get_var ("slur_ratio");
638 bb.set_default_bezier (h_inf, r_0);
640 if (bb.fit_factor () > 1.0)
642 Real length = bb.curve_.control_[3][X_AXIS];
643 Real default_height = bb.get_default_height (h_inf, r_0, length);
644 bb.minimise_enclosed_area (paper_l(), default_height);
646 Real bff = paper_l ()->get_var ("slur_force_blowfit");
647 bb.curve_.control_[1][Y_AXIS] *= bff;
648 bb.curve_.control_[2][Y_AXIS] *= bff;
651 Real sb = paper_l ()->get_var ("slur_beautiful");
652 Real beautiful = length * default_height * sb;
653 Real area = bb.enclosed_area_f ();
656 Slurs that fit beautifully are not ugly
658 if (area > beautiful)
659 de_uglyfy (&bb, default_height);
662 Bezier b = bb.get_bezier ();
665 SCM controls = SCM_EOL;
667 controls = gh_cons ( ly_offset2scm (b.control_[i]), controls);
669 set_elt_property ("control-points", controls);
674 Slur::get_curve () const
679 if (!directional_element (this).get ())
680 ((Slur*)this)->set_extremities ();
682 if (!gh_pair_p (get_elt_property ("control-points")))
683 ((Slur*)this)->set_control_points ();
686 for (SCM s= get_elt_property ("control-points"); s != SCM_EOL; s = gh_cdr (s))
688 b.control_[i] = ly_scm2offset (gh_car (s));
692 Array<Offset> enc (get_encompass_offset_arr ());
693 Direction dir = directional_element (this).get ();
695 Real x1 = enc[0][X_AXIS];
696 Real x2 = enc.top ()[X_AXIS];
699 for (int i=1; i < enc.size ()-1; i++)
701 Real x = enc[i][X_AXIS];
704 Real y = b.get_other_coordinate (X_AXIS, x);
705 off = off >? dir * (enc[i][Y_AXIS] - y);
708 b.translate (Offset (0, dir * off));