]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/beam.cc
Fix some bugs in the dynamic engraver and PostScript backend
[lilypond.git] / lily / beam.cc
index 8a98a3e96535d662df112e3e6aca7c0645dd7dc2..9dd66cef9a4124658bf447b14d81dc1f9ecac1d1 100644 (file)
@@ -3,7 +3,7 @@
 
   source file of the GNU LilyPond music typesetter
 
-  (c) 1997--2005 Han-Wen Nienhuys <hanwen@xs4all.nl>
+  (c) 1997--2006 Han-Wen Nienhuys <hanwen@xs4all.nl>
   Jan Nieuwenhuizen <janneke@gnu.org>
 */
 
 */
 
 #include "beam.hh"
-#include "interval-set.hh"
+
+#include "beaming-pattern.hh"
 #include "directional-element-interface.hh"
-#include "beaming.hh"
-#include "misc.hh"
+#include "main.hh"
+#include "international.hh"
+#include "interval-set.hh"
+#include "item.hh"
 #include "least-squares.hh"
-#include "stem.hh"
-#include "output-def.hh"
 #include "lookup.hh"
+#include "misc.hh"
+#include "output-def.hh"
 #include "pointer-group-interface.hh"
-#include "staff-symbol-referencer.hh"
-#include "item.hh"
 #include "spanner.hh"
+#include "staff-symbol-referencer.hh"
+#include "stem.hh"
 #include "warn.hh"
 
-#if DEBUG_QUANTING
+#if DEBUG_BEAM_SCORING
 #include "text-interface.hh" // debug output.
 #include "font-interface.hh" // debug output.
 #endif
 
+#include <map>
+
+
+Beam_stem_segment::Beam_stem_segment ()
+{
+  max_connect_ = 1000;         // infinity
+  stem_ = 0;
+  width_ = 0.0;
+  stem_x_ = 0.0;
+  rank_ = 0;
+  stem_index_ = 0;
+  dir_ = CENTER;
+}
+
+Beam_segment::Beam_segment ()
+{
+  vertical_count_ = 0;
+}
+
 void
 Beam::add_stem (Grob *me, Grob *s)
 {
@@ -89,7 +111,7 @@ Beam::get_beam_count (Grob *me)
   int m = 0;
 
   extract_grob_set (me, "stems", stems);
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *stem = stems[i];
       m = max (m, (Stem::beam_multiplicity (stem).length () + 1));
@@ -98,12 +120,6 @@ Beam::get_beam_count (Grob *me)
 }
 
 
-/* After pre-processing all directions should be set.
-   Several post-processing routines (stem, slur, script) need stem/beam
-   direction.
-   Currenly, this means that beam has set all stem's directions.
-   [Alternatively, stems could set its own directions, according to
-   their beam, during 'final-pre-processing'.] */
 MAKE_SCHEME_CALLBACK (Beam, calc_direction, 1);
 SCM
 Beam::calc_direction (SCM smob)
@@ -113,56 +129,48 @@ Beam::calc_direction (SCM smob)
   /* Beams with less than 2 two stems don't make much sense, but could happen
      when you do
 
-     [r8 c8 r8].
+     r8[ c8 r8]
 
-     For a beam that  only has one stem, we try to do some disappearance magic:
-     we revert the flag, and move on to The Eternal Engraving Fields. */
+  */
 
-  Direction d = CENTER;
+  Direction dir = CENTER;
 
   int count = visible_stem_count (me);
   if (count < 2)
     {
       extract_grob_set (me, "stems", stems);
-      if (stems.size () == 1)
+      if (stems.size () == 0)
        {
-         me->warning (_ ("removing beam with less than two stems"));
-
-         stems[0]->set_object ("beam", SCM_EOL);
+         me->warning (_ ("removing beam with no stems"));
          me->suicide ();
 
          return SCM_UNSPECIFIED;
        }
-      else if (stems.size () == 0)
+      else 
        {
-         me->suicide ();
-         return SCM_UNSPECIFIED;
-       }
-      else
-       {
-         /*
-           ugh.
+         Grob *stem = first_visible_stem (me);
 
-           can happen in stem-tremolo case.
-           TODO: fixme.
-          */
-         d = Stem::get_default_dir (stems[0]);
+         /*
+           ugh: stems[0] case happens for chord tremolo.
+         */
+         dir = to_dir ((stem ? stem : stems[0])->get_property ("default-direction"));
        }
     }
 
-  
   if (count >= 1)
     {
-      d = get_default_dir (me);
+      if (!dir)
+       dir = get_default_dir (me);
+      
       consider_auto_knees (me);
     }
 
-  if (d)
+  if (dir)
     {
-      set_stem_directions (me, d);
+      set_stem_directions (me, dir);
     }
   
-  return scm_from_int (d);
+  return scm_from_int (dir);
 }
 
 
@@ -222,7 +230,7 @@ Beam::calc_beaming (SCM smob)
   
   SCM last_beaming = scm_cons (SCM_EOL, scm_list_1 (scm_from_int (0)));
   Direction last_dir = CENTER;
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *this_stem = stems[i];
       SCM this_beaming = this_stem->get_property ("beaming");
@@ -276,39 +284,186 @@ Beam::calc_beaming (SCM smob)
   return SCM_EOL;
 }
 
-/*
-  I really enjoy spaghetti, but spaghetti should be kept on a plate
-  with a little garlic and olive oil. This is too much.
+bool
+operator <(Beam_stem_segment const &a,
+          Beam_stem_segment const &b)
+{
+  return a.rank_ < b.rank_;
+}
 
-  rewrite-me
-*/
-MAKE_SCHEME_CALLBACK (Beam, print, 1);
-SCM
-Beam::print (SCM grob)
+typedef map<int, vector<Beam_stem_segment> >  Position_stem_segments_map; 
+
+vector<Beam_segment>
+Beam::get_beam_segments (Grob *me_grob, Grob **common)
 {
-  Spanner *me = unsmob_spanner (grob);
+  Spanner *me = dynamic_cast<Spanner*> (me_grob);
 
   extract_grob_set (me, "stems", stems);
-  Grob *xcommon = common_refpoint_of_array (stems, me, X_AXIS);
+  Grob *commonx = common_refpoint_of_array (stems, me, X_AXIS);
+
+  commonx = me->get_bound (LEFT)->common_refpoint (commonx, X_AXIS);
+  commonx = me->get_bound (RIGHT)->common_refpoint (commonx, X_AXIS);
+
+  *common = commonx;
+  
+  int gap_count = robust_scm2int (me->get_property ("gap-count"), 0);
+  Real gap_length = robust_scm2double (me->get_property ("gap"), 0.0);
+
+  Position_stem_segments_map stem_segments;
+  Real lt = me->layout ()->get_dimension (ly_symbol2scm ("line-thickness"));
+
+  Slice ranks;
+  
+  for (vsize i = 0; i < stems.size (); i++)
+    {
+      Grob *stem = stems[i];
+      Real stem_width = robust_scm2double (stem->get_property ("thickness"), 1.0) * lt;
+      Real stem_x = stem->relative_coordinate (commonx, X_AXIS);
+      SCM beaming = stem->get_property ("beaming");
+      Direction d = LEFT;
+      do
+       {
+         for (SCM s = index_get_cell (beaming, d);
+              scm_is_pair (s); s = scm_cdr (s))
+           {
+             if (!scm_is_integer (scm_car (s)))
+               continue;
+
+             int beam_rank = scm_to_int (scm_car (s));
+             ranks.add_point (beam_rank);
+           }
+         
+         for (SCM s = index_get_cell (beaming, d);
+              scm_is_pair (s); s = scm_cdr (s))
+           {
+             if (!scm_is_integer (scm_car (s)))
+               continue;
+         
+             int beam_rank = scm_to_int (scm_car (s));
+             Beam_stem_segment seg;
+             seg.stem_ = stem;
+             seg.stem_x_ = stem_x;
+             seg.rank_ = 2 * i  + (d+1)/2;
+             seg.width_ = stem_width;
+             seg.stem_index_ = i;
+             seg.dir_ = d;
+             seg.max_connect_ = robust_scm2int (stem->get_property ("max-beam-connect"), 1000);
+             
+             Direction stem_dir = get_grob_direction (stem);
+             
+             seg.gapped_
+               = (stem_dir * beam_rank < (stem_dir * ranks[-stem_dir] + gap_count));
+             stem_segments[beam_rank].push_back (seg);
+           }
+       }
+      while (flip (&d) != LEFT);
+    }
+
+  Drul_array<Real> break_overshoot
+    = robust_scm2drul (me->get_property ("break-overshoot"),
+                      Drul_array<Real> (-0.5, 0.0));
+
+  vector<Beam_segment> segments;
+  for (Position_stem_segments_map::const_iterator i (stem_segments.begin ());
+       i != stem_segments.end (); i++)
+    {
+      vector<Beam_stem_segment> segs = (*i).second;
+      vector_sort (segs, default_compare);
+
+      Beam_segment current;
+
+      int vertical_count =  (*i).first;
+      for (vsize j = 0; j < segs.size (); j++)
+       {
+         /*
+           event_dir == LEFT: left edge of a beamsegment.
+          */
+         Direction event_dir = LEFT;
+         do
+           {
+             Drul_array<bool> on_bound (j == 0 && event_dir==LEFT,
+                                        j == segs.size() - 1 && event_dir==RIGHT);
+             Drul_array<bool> inside (j > 0, j < segs.size()-1);
+             bool event = on_bound[event_dir]
+               || abs (segs[j].rank_ - segs[j+event_dir].rank_) > 1
+               || (abs (vertical_count) >= segs[j].max_connect_
+                   || abs (vertical_count) >= segs[j + event_dir].max_connect_);
+             
+             if (!event)
+               continue;
+
+             current.vertical_count_ = vertical_count;
+             current.horizontal_[event_dir] = segs[j].stem_x_;
+             if (segs[j].dir_ == event_dir)
+               {
+                 if (on_bound[event_dir]
+                     && me->get_bound (event_dir)->break_status_dir ())
+                   {
+                     current.horizontal_[event_dir]
+                       = (me->get_bound (event_dir)->extent (commonx, X_AXIS)[RIGHT]
+                          + event_dir * break_overshoot[event_dir]);
+                   }
+                 else
+                   {
+                     Real notehead_width = 
+                       Stem::duration_log (segs[j].stem_) == 1
+                       ? 1.98
+                       : 1.32; // URG.
+                     
+                     if (inside[event_dir])
+                       notehead_width = min (notehead_width,
+                                             fabs (segs[j+ event_dir].stem_x_
+                                                   - segs[j].stem_x_)/2);
+                     
+                     current.horizontal_[event_dir] += event_dir * notehead_width;
+                   }
+               }
+             else
+               {
+                 current.horizontal_[event_dir] += event_dir * segs[j].width_/2;
+                 if (segs[j].gapped_)
+                   current.horizontal_[event_dir] -= event_dir * gap_length;  
+               }
+
+             if (event_dir == RIGHT)
+               {
+                 segments.push_back (current);
+                 current = Beam_segment();
+               }
+           }
+         while (flip (&event_dir) != LEFT);
+       }
+      
+    }
 
-  xcommon = me->get_bound (LEFT)->common_refpoint (xcommon, X_AXIS);
-  xcommon = me->get_bound (RIGHT)->common_refpoint (xcommon, X_AXIS);
+  return segments;
+}
 
-  Real x0, dx;
+MAKE_SCHEME_CALLBACK(Beam, print, 1);
+SCM
+Beam::print (SCM grob)
+{
+  Spanner *me = unsmob_spanner (grob);
+  Grob *commonx = 0;
+  vector<Beam_segment> segments = get_beam_segments (me, &commonx);
+
+  Interval span;
   if (visible_stem_count (me))
     {
-      // ugh -> use commonx
-      x0 = first_visible_stem (me)->relative_coordinate (xcommon, X_AXIS);
-      dx = last_visible_stem (me)->relative_coordinate (xcommon, X_AXIS) - x0;
+      span[LEFT] = first_visible_stem (me)->relative_coordinate (commonx, X_AXIS);
+      span[RIGHT] = last_visible_stem (me)->relative_coordinate (commonx, X_AXIS);
     }
   else
     {
-      x0 = stems[0]->relative_coordinate (xcommon, X_AXIS);
-      dx = stems.top ()->relative_coordinate (xcommon, X_AXIS) - x0;
+      extract_grob_set (me, "stems", stems);      
+      span[LEFT] = stems[0]->relative_coordinate (commonx, X_AXIS);
+      span[RIGHT] = stems.back ()->relative_coordinate (commonx, X_AXIS);
     }
 
+  Real blot = me->layout ()->get_dimension (ly_symbol2scm ("blot-diameter"));
+
   SCM posns = me->get_property ("quantized-positions");
-  Drul_array<Real> pos;
+  Interval pos;
   if (!is_number_pair (posns))
     {
       programming_error ("no beam positions?");
@@ -320,191 +475,46 @@ Beam::print (SCM grob)
   scale_drul (&pos, Staff_symbol_referencer::staff_space (me));
 
   Real dy = pos[RIGHT] - pos[LEFT];
-  Real slope = (dy && dx) ? dy / dx : 0;
+  Real slope = (dy && span.length ()) ? dy / span.length ()  : 0;
 
   Real thick = get_thickness (me);
-  Real bdy = get_beam_translation (me);
-
-  SCM last_beaming = SCM_EOL;
-  Real last_xposn = -1;
-  Real last_stem_width = -1;
-
-  Real gap_length = robust_scm2double (me->get_property ("gap"), 0.0);
+  Real beam_dy = get_beam_translation (me);
 
+  Direction feather_dir = to_dir (me->get_property ("grow-direction"));
+  
   Stencil the_beam;
-  Real lt = me->layout ()->get_dimension (ly_symbol2scm ("linethickness"));
-
-  for (int i = 0; i <= stems.size (); i++)
+  for (vsize i = 0; i < segments.size (); i ++)
     {
-      Grob *stem = (i < stems.size ()) ? stems[i] : 0;
-
-      SCM this_beaming = stem ? stem->get_property ("beaming") : SCM_EOL;
-      Real xposn = stem ? stem->relative_coordinate (xcommon, X_AXIS) : 0.0;
-      Real stem_width = stem ? robust_scm2double (stem->get_property ("thickness"), 1.0) * lt : 0;
-      Direction stem_dir = stem ? to_dir (stem->get_property ("direction")) : CENTER;
-      /*
-       We do the space left of ST, with lfliebertjes pointing to the
-       right from the left stem, and rfliebertjes pointing left from
-       right stem.
-      */
-      SCM left = (i > 0) ? scm_cdr (last_beaming) : SCM_EOL;
-      SCM right = stem ? scm_car (this_beaming) : SCM_EOL;
-
-      Array<int> full_beams;
-      Array<int> lfliebertjes;
-      Array<int> rfliebertjes;
-
-      for (SCM s = left;
-          scm_is_pair (s); s = scm_cdr (s))
+      Real local_slope = slope;
+      if (feather_dir)
        {
-         int b = scm_to_int (scm_car (s));
-         if (scm_c_memq (scm_car (s), right) != SCM_BOOL_F)
-           full_beams.push (b);
-         else
-           lfliebertjes.push (b);
-       }
-      for (SCM s = right;
-          scm_is_pair (s); s = scm_cdr (s))
-       {
-         int b = scm_to_int (scm_car (s));
-         if (scm_c_memq (scm_car (s), left) == SCM_BOOL_F)
-           rfliebertjes.push (b);
-       }
-
-      Drul_array<Real> break_overshoot
-       = robust_scm2drul (me->get_property ("break-overshoot"),
-                          Drul_array<Real> (-0.5, 0.0));
-
-      Real w = (i > 0 && stem)
-       ? (xposn - last_xposn)
-       : break_overshoot[ (i == 0) ? LEFT : RIGHT];
-
-      Real stem_offset = 0.0;
-      if (i > 0)
-       {
-         w += last_stem_width / 2;
-         stem_offset = -last_stem_width / 2;
-       }
-
-      if (stem)
-       w += stem_width / 2;
-
-      Real blot = me->layout ()->get_dimension (ly_symbol2scm ("blotdiameter"));
-      Stencil whole = Lookup::beam (slope, w, thick, blot);
-      Stencil gapped;
-
-      int gap_count = 0;
-      if (scm_is_number (me->get_property ("gap-count")))
-       {
-         gap_count = scm_to_int (me->get_property ("gap-count"));
-         gapped = Lookup::beam (slope, w - 2 * gap_length, thick, blot);
-
-         full_beams.sort (default_compare);
-         if (stem_dir == UP)
-           full_beams.reverse ();
-       }
-
-      int k = 0;
-      for (int j = full_beams.size (); j--;)
-       {
-         Stencil b (whole);
-
-         if (k++ < gap_count)
-           {
-             b = gapped;
-             b.translate_axis (gap_length, X_AXIS);
-           }
-         b.translate_axis (last_xposn - x0 + stem_offset, X_AXIS);
-         b.translate_axis (slope * (last_xposn - x0) + bdy * full_beams[j], Y_AXIS);
-
-         the_beam.add_stencil (b);
-       }
-
-      if (lfliebertjes.size () || rfliebertjes.size ())
-       {
-         Real nw_f;
-
-         if (stem)
-           {
-             int t = Stem::duration_log (stem);
-             // ugh. hardcoded.
-             if (t == 1)
-               nw_f = 1.98;
-             else
-               nw_f = 1.32;
-           }
-         else
-           nw_f = break_overshoot[RIGHT] / 2;
-
-         /* Half beam should be one note-width,
-            but let's make sure two half-beams never touch */
-         Real lw = nw_f;
-         Real rw = nw_f;
-         if (i > 0)
-           rw = min (nw_f, ((xposn - last_xposn) / 2));
-         else
-           {
-             if (me->get_bound (LEFT)->break_status_dir ())
-               rw = xposn - me->get_bound (LEFT)->extent (xcommon, X_AXIS)[RIGHT]
-                 + break_overshoot[LEFT];
-             else
-               rw = 1.0;       // ugh.
-           }
-         
-         if (stem)
-           lw = min (nw_f, ((xposn - last_xposn) / 2));
-         else
-           {
-             lw = me->get_bound (RIGHT)->relative_coordinate (xcommon, X_AXIS)
-               - last_xposn
-               + break_overshoot[RIGHT];
-           }
-         rw += stem_width / 2;
-         lw += last_stem_width / 2;
-
-         Stencil rhalf = Lookup::beam (slope, rw, thick, blot);
-         Stencil lhalf = Lookup::beam (slope, lw, thick, blot);
-         for (int j = lfliebertjes.size (); j--;)
-           {
-             Stencil b (lhalf);
-             b.translate_axis (last_xposn - x0 - last_stem_width /2,
-                               X_AXIS);
-             b.translate_axis (slope * (last_xposn - x0)
-                               + bdy * lfliebertjes[j],
-                               Y_AXIS);
-             the_beam.add_stencil (b);
-           }
-         for (int j = rfliebertjes.size (); j--;)
-           {
-             Stencil b (rhalf);
-             b.translate_axis (xposn - x0 - rw + stem_width / 2, X_AXIS);
-             b.translate_axis (slope * (xposn - x0 - rw)
-                               + bdy * rfliebertjes[j], Y_AXIS);
-             the_beam.add_stencil (b);
-           }
+         local_slope += feather_dir * segments[i].vertical_count_ * beam_dy / span.length ();
        }
+      
+      Stencil b = Lookup::beam (local_slope, segments[i].horizontal_.length (), thick, blot);
 
-      last_xposn = xposn;
-      last_stem_width = stem_width;
-      last_beaming = this_beaming;
+      b.translate_axis (segments[i].horizontal_[LEFT], X_AXIS);
+      
+      b.translate_axis (local_slope
+                       * (segments[i].horizontal_[LEFT] - span.linear_combination (feather_dir))
+                       + pos.linear_combination (feather_dir)
+                       + beam_dy * segments[i].vertical_count_, Y_AXIS);
+      the_beam.add_stencil (b);      
     }
-
-  the_beam.translate_axis (x0 - me->relative_coordinate (xcommon, X_AXIS),
-                          X_AXIS);
-  the_beam.translate_axis (pos[LEFT], Y_AXIS);
-
-#if (DEBUG_QUANTING)
+        
+#if (DEBUG_BEAM_SCORING)
   SCM quant_score = me->get_property ("quant-score");
-  SCM debug = me->layout ()->lookup_variable (ly_symbol2scm ("debug-beam-quanting"));
+  SCM debug = me->layout ()->lookup_variable (ly_symbol2scm ("debug-beam-scoring"));
   if (to_boolean (debug) && scm_is_string (quant_score))
     {
+      extract_grob_set (me, "stems", stems);      
 
       /*
        This code prints the demerits for each beam. Perhaps this
        should be switchable for those who want to twiddle with the
        parameters.
       */
-      String str;
+      string str;
       SCM properties = Font_interface::text_font_alist_chain (me);
 
       Direction stem_dir = stems.size () ? to_dir (stems[0]->get_property ("direction")) : UP;
@@ -517,28 +527,47 @@ Beam::print (SCM grob)
     }
 #endif
 
+  the_beam.translate_axis (-me->relative_coordinate (commonx, X_AXIS), X_AXIS);
   return the_beam.smobbed_copy ();
 }
-
 Direction
 Beam::get_default_dir (Grob *me)
 {
-  Drul_array<int> total;
-  total[UP] = total[DOWN] = 0;
-  Drul_array<int> count;
-  count[UP] = count[DOWN] = 0;
-
   extract_grob_set (me, "stems", stems);
 
-  for (int i = 0; i < stems.size (); i++)
+  Drul_array<Real> extremes (0.0, 0.0);
+  for (iterof (s, stems); s != stems.end (); s++)
+    {
+      Interval positions = Stem::head_positions (*s);
+      Direction d = DOWN;
+      do
+       {
+         if (sign (positions[d]) == d)
+           extremes[d] = d * max (d * positions[d], d * extremes[d]);
+       }
+      while (flip (&d) != DOWN);
+    }
+
+  Drul_array<int> total (0, 0);
+  Drul_array<int> count (0, 0);
+
+  bool force_dir = false;
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
       Direction stem_dir = CENTER;
       SCM stem_dir_scm = s->get_property_data (ly_symbol2scm ("direction"));
       if (is_direction (stem_dir_scm))
-       stem_dir = to_dir (stem_dir_scm);
+       {
+         stem_dir = to_dir (stem_dir_scm);
+         force_dir = true;
+       }
       else
-       stem_dir = Stem::get_default_dir (s);
+       stem_dir = to_dir (s->get_property ("default-direction"));
+
+      if (!stem_dir)
+       stem_dir = to_dir (s->get_property ("neutral-direction"));
 
       if (stem_dir)
        {
@@ -547,13 +576,24 @@ Beam::get_default_dir (Grob *me)
        }
     }
 
-  Direction dir = CENTER;
+
+  if (!force_dir)
+    {
+      if (abs (extremes[UP]) > -extremes[DOWN])
+       return DOWN;
+      else if (extremes[UP] < -extremes[DOWN])
+       return UP;
+    }
   
-  if (Direction d =  (Direction) sign (count[UP] - count[DOWN]))
+  Direction dir = CENTER;
+  Direction d = CENTER;
+  if ((d = (Direction) sign (count[UP] - count[DOWN])))
     dir = d;
-  else if (Direction d = (Direction)  sign (total[UP] / count[UP] - total[DOWN]/count[DOWN]))
+  else if (count[UP]
+          && count[DOWN]
+          && (d = (Direction)  sign (total[UP] / count[UP] - total[DOWN]/count[DOWN])))
     dir = d;
-  else if (Direction d = (Direction)  sign (total[UP] - total[DOWN]))
+  else if ((d = (Direction)  sign (total[UP] - total[DOWN])))
     dir = d;
   else
     dir = to_dir (me->get_property ("neutral-direction"));
@@ -569,7 +609,7 @@ Beam::set_stem_directions (Grob *me, Direction d)
 {
   extract_grob_set (me, "stems", stems);
 
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
 
@@ -607,8 +647,8 @@ Beam::consider_auto_knees (Grob *me)
   Grob *common = common_refpoint_of_array (stems, me, Y_AXIS);
   Real staff_space = Staff_symbol_referencer::staff_space (me);
 
-  Array<Interval> head_extents_array;
-  for (int i = 0; i < stems.size (); i++)
+  vector<Interval> head_extents_array;
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *stem = stems[i];
       if (Stem::is_invisible (stem))
@@ -634,7 +674,7 @@ Beam::consider_auto_knees (Grob *me)
              head_extents[-stemdir] = -stemdir * infinity_f;
            }
        }
-      head_extents_array.push (head_extents);
+      head_extents_array.push_back (head_extents);
 
       gaps.remove_interval (head_extents);
     }
@@ -642,7 +682,7 @@ Beam::consider_auto_knees (Grob *me)
   Interval max_gap;
   Real max_gap_len = 0.0;
 
-  for (int i = gaps.allowed_regions_.size () -1; i >= 0; i--)
+  for (vsize i = gaps.allowed_regions_.size () -1; i != VPOS ;i--)
     {
       Interval gap = gaps.allowed_regions_[i];
 
@@ -669,7 +709,7 @@ Beam::consider_auto_knees (Grob *me)
   if (max_gap_len > threshold)
     {
       int j = 0;
-      for (int i = 0; i < stems.size (); i++)
+      for (vsize i = 0; i < stems.size (); i++)
        {
          Grob *stem = stems[i];
          if (Stem::is_invisible (stem))
@@ -777,8 +817,8 @@ Beam::calc_least_squares_positions (SCM smob, SCM posns)
   Interval pos (0,0);
   if (count < 1)
     return ly_interval2scm (pos);
-  
-  Array<Real> x_posns;
+
+  vector<Real> x_posns;
   extract_grob_set (me, "stems", stems);
   Grob *commonx = common_refpoint_of_array (stems, me, X_AXIS);
   Grob *commony = common_refpoint_of_array (stems, me, Y_AXIS);
@@ -789,17 +829,17 @@ Beam::calc_least_squares_positions (SCM smob, SCM posns)
   Grob *lvs = last_visible_stem (me);
 
   Interval ideal (Stem::get_stem_info (fvs).ideal_y_
-                 + fvs->relative_coordinate (commony, Y_AXIS) -my_y,
+                 + fvs->relative_coordinate (commony, Y_AXIS) - my_y,
                  Stem::get_stem_info (lvs).ideal_y_
                  + lvs->relative_coordinate (commony, Y_AXIS) - my_y);
 
   Real x0 = first_visible_stem (me)->relative_coordinate (commonx, X_AXIS);
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
 
       Real x = s->relative_coordinate (commonx, X_AXIS) - x0;
-      x_posns.push (x);
+      x_posns.push_back (x);
     }
   Real dx = last_visible_stem (me)->relative_coordinate (commonx, X_AXIS) - x0;
 
@@ -838,13 +878,13 @@ Beam::calc_least_squares_positions (SCM smob, SCM posns)
     }
   else
     {
-      Array<Offset> ideals;
-      for (int i = 0; i < stems.size (); i++)
+      vector<Offset> ideals;
+      for (vsize i = 0; i < stems.size (); i++)
        {
          Grob *s = stems[i];
          if (Stem::is_invisible (s))
            continue;
-         ideals.push (Offset (x_posns[i],
+         ideals.push_back (Offset (x_posns[i],
                               Stem::get_stem_info (s).ideal_y_
                               + s->relative_coordinate (commony, Y_AXIS)
                               - my_y));
@@ -884,7 +924,7 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
   /*
     Code dup.
   */
-  Array<Real> x_posns;
+  vector<Real> x_posns;
   extract_grob_set (me, "stems", stems);
   Grob *commonx = common_refpoint_of_array (stems, me, X_AXIS);
   Grob *commony = common_refpoint_of_array (stems, me, Y_AXIS);
@@ -895,12 +935,12 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
     return posns;
 
   Real x0 = fvs->relative_coordinate (commonx, X_AXIS);
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
 
       Real x = s->relative_coordinate (commonx, X_AXIS) - x0;
-      x_posns.push (x);
+      x_posns.push_back (x);
     }
 
   Grob *lvs = last_visible_stem (me);
@@ -924,7 +964,7 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
   */
   Interval feasible_left_point;
   feasible_left_point.set_full ();
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
       if (Stem::is_invisible (s))
@@ -1047,13 +1087,13 @@ where_are_the_whole_beams (SCM beaming)
 /* Return the Y position of the stem-end, given the Y-left, Y-right
    in POS for stem S.  This Y position is relative to S. */
 Real
-Beam::calc_stem_y (Grob *me, Grob *s, Grob ** common,
+Beam::calc_stem_y (Grob *me, Grob *stem, Grob **common,
                   Real xl, Real xr,
                   Drul_array<Real> pos, bool french)
 {
   Real beam_translation = get_beam_translation (me);
 
-  Real r = s->relative_coordinate (common[X_AXIS], X_AXIS) - xl;
+  Real r = stem->relative_coordinate (common[X_AXIS], X_AXIS) - xl;
   Real dy = pos[RIGHT] - pos[LEFT];
   Real dx = xr - xl;
   Real stem_y_beam0 = (dy && dx
@@ -1061,8 +1101,8 @@ Beam::calc_stem_y (Grob *me, Grob *s, Grob ** common,
                       * dy
                       : 0) + pos[LEFT];
 
-  Direction my_dir = get_grob_direction (s);
-  SCM beaming = s->get_property ("beaming");
+  Direction my_dir = get_grob_direction (stem);
+  SCM beaming = stem->get_property ("beaming");
 
   Real stem_y = stem_y_beam0;
   if (french)
@@ -1073,13 +1113,13 @@ Beam::calc_stem_y (Grob *me, Grob *s, Grob ** common,
     }
   else
     {
-      Slice bm = Stem::beam_multiplicity (s);
+      Slice bm = Stem::beam_multiplicity (stem);
       if (!bm.is_empty ())
        stem_y += bm[my_dir] * beam_translation;
     }
 
   Real id = me->relative_coordinate (common[Y_AXIS], Y_AXIS)
-    - s->relative_coordinate (common[Y_AXIS], Y_AXIS);
+    - stem->relative_coordinate (common[Y_AXIS], Y_AXIS);
 
   return stem_y + id;
 }
@@ -1094,8 +1134,9 @@ Beam::set_stem_lengths (SCM smob)
 {
   Grob *me = unsmob_grob (smob);
 
-  /* trigger callback. */
+  /* trigger callbacks. */
   (void) me->get_property ("direction");
+  (void) me->get_property ("beaming");
 
   SCM posns = me->get_property ("positions");
   
@@ -1113,8 +1154,7 @@ Beam::set_stem_lengths (SCM smob)
 
   bool gap = false;
   Real thick = 0.0;
-  if (scm_is_number (me->get_property ("gap-count"))
-      && scm_to_int (me->get_property ("gap-count")))
+  if (robust_scm2int (me->get_property ("gap-count"), 0))
     {
       gap = true;
       thick = get_thickness (me);
@@ -1126,7 +1166,7 @@ Beam::set_stem_lengths (SCM smob)
   Real xl = fvs ? fvs->relative_coordinate (common[X_AXIS], X_AXIS) : 0.0;
   Real xr = lvs ? lvs->relative_coordinate (common[X_AXIS], X_AXIS) : 0.0;
 
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
       if (Stem::is_invisible (s))
@@ -1151,12 +1191,12 @@ Beam::set_stem_lengths (SCM smob)
 }
 
 void
-Beam::set_beaming (Grob *me, Beaming_info_list const *beaming)
+Beam::set_beaming (Grob *me, Beaming_pattern const *beaming)
 {
   extract_grob_set (me, "stems", stems);
 
   Direction d = LEFT;
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       /*
        Don't overwrite user settings.
@@ -1168,13 +1208,19 @@ Beam::set_beaming (Grob *me, Beaming_info_list const *beaming)
          if (beaming_prop == SCM_EOL
              || index_get_cell (beaming_prop, d) == SCM_EOL)
            {
-             int b = beaming->infos_.elem (i).beams_i_drul_[d];
+             int count = beaming->beamlet_count (i, d);
              if (i > 0
                  && i < stems.size () -1
                  && Stem::is_invisible (stem))
-               b = min (b, beaming->infos_.elem (i).beams_i_drul_[-d]);
+               count = min (count, beaming->beamlet_count (i,-d));
 
-             Stem::set_beaming (stem, b, d);
+             if ( ((i == 0 && d == LEFT)
+                   || (i == stems.size ()-1 && d == RIGHT))
+                  && stems.size () > 1
+                  && to_boolean (me->get_property ("clip-edges")))
+               count = 0;
+
+             Stem::set_beaming (stem, count, d);
            }
        }
       while (flip (&d) != LEFT);
@@ -1187,7 +1233,7 @@ Beam::forced_stem_count (Grob *me)
   extract_grob_set (me, "stems", stems);
 
   int f = 0;
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       Grob *s = stems[i];
 
@@ -1196,8 +1242,11 @@ Beam::forced_stem_count (Grob *me)
 
       /* I can imagine counting those boundaries as a half forced stem,
         but let's count them full for now. */
+      Direction defdir = to_dir (s->get_property ("default-direction"));
+      
       if (abs (Stem::chord_start_y (s)) > 0.1
-         && (get_grob_direction (s) != Stem::get_default_dir (s)))
+         && defdir
+         && get_grob_direction (s) != defdir)
        f++;
     }
   return f;
@@ -1208,7 +1257,7 @@ Beam::visible_stem_count (Grob *me)
 {
   extract_grob_set (me, "stems", stems);
   int c = 0;
-  for (int i = stems.size (); i--;)
+  for (vsize i = stems.size (); i--;)
     {
       if (!Stem::is_invisible (stems[i]))
        c++;
@@ -1221,7 +1270,7 @@ Beam::first_visible_stem (Grob *me)
 {
   extract_grob_set (me, "stems", stems);
 
-  for (int i = 0; i < stems.size (); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       if (!Stem::is_invisible (stems[i]))
        return stems[i];
@@ -1234,7 +1283,7 @@ Beam::last_visible_stem (Grob *me)
 {
   extract_grob_set (me, "stems", stems);
 
-  for (int i = stems.size (); i--;)
+  for (vsize i = stems.size (); i--;)
     {
       if (!Stem::is_invisible (stems[i]))
        return stems[i];
@@ -1349,7 +1398,7 @@ Beam::is_knee (Grob *me)
   bool knee = false;
   int d = 0;
   extract_grob_set (me, "stems", stems);
-  for (int i = stems.size (); i--;)
+  for (vsize i = stems.size (); i--;)
     {
       Direction dir = get_grob_direction (stems[i]);
       if (d && d != dir)
@@ -1371,7 +1420,7 @@ Beam::get_direction_beam_count (Grob *me, Direction d)
   extract_grob_set (me, "stems", stems);
   int bc = 0;
 
-  for (int i = stems.size (); i--;)
+  for (vsize i = stems.size (); i--;)
     {
       /*
        Should we take invisible stems into account?
@@ -1399,13 +1448,14 @@ ADD_INTERFACE (Beam,
               "beamed-stem-shorten "
               "beaming "
               "break-overshoot "
-              "chord-tremolo "
+              "clip-edges "
               "concaveness "
               "damping "
               "details "
               "direction " 
               "gap "
               "gap-count "
+              "grow-direction "
               "inspect-quants "
               "knee "
               "length-fraction "