]> 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 49acbbd541189fbc8b05b55a1f9b9d955ce8274e..9dd66cef9a4124658bf447b14d81dc1f9ecac1d1 100644 (file)
@@ -26,8 +26,9 @@
 
 #include "beam.hh"
 
-#include "beaming.hh"
+#include "beaming-pattern.hh"
 #include "directional-element-interface.hh"
+#include "main.hh"
 #include "international.hh"
 #include "interval-set.hh"
 #include "item.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)
 {
@@ -100,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)
@@ -115,29 +129,21 @@ 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)
-       {
-         me->suicide ();
          return SCM_UNSPECIFIED;
        }
       else 
@@ -147,24 +153,24 @@ Beam::calc_direction (SCM smob)
          /*
            ugh: stems[0] case happens for chord tremolo.
          */
-         d = to_dir ((stem ? stem : stems[0])->get_property ("default-direction"));
+         dir = to_dir ((stem ? stem : stems[0])->get_property ("default-direction"));
        }
     }
 
   if (count >= 1)
     {
-      if (!d)
-       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);
 }
 
 
@@ -278,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;
 
-  xcommon = me->get_bound (LEFT)->common_refpoint (xcommon, X_AXIS);
-  xcommon = me->get_bound (RIGHT)->common_refpoint (xcommon, X_AXIS);
+             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);
+       }
+      
+    }
+
+  return segments;
+}
+
+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);
 
-  Real x0, dx;
+  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.back ()->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?");
@@ -322,184 +475,39 @@ 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 ("line-thickness"));
-
-  for (vsize 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;
-
-      vector<int> full_beams;
-      vector<int> lfliebertjes;
-      vector<int> rfliebertjes;
-
-      for (SCM s = left;
-          scm_is_pair (s); s = scm_cdr (s))
-       {
-         int b = scm_to_int (scm_car (s));
-         if (scm_c_memq (scm_car (s), right) != SCM_BOOL_F)
-           full_beams.push_back (b);
-         else
-           lfliebertjes.push_back (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_back (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 ("blot-diameter"));
-      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);
-
-         vector_sort (full_beams, default_compare);
-         if (stem_dir == UP)
-           reverse (full_beams);
-       }
-
-      int k = 0;
-      for (vsize 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 local_slope = slope;
+      if (feather_dir)
        {
-         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 (vsize 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 (vsize 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
@@ -519,26 +527,42 @@ 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);
 
+  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 = to_dir (s->get_property ("default-direction"));
 
@@ -552,6 +576,15 @@ Beam::get_default_dir (Grob *me)
        }
     }
 
+
+  if (!force_dir)
+    {
+      if (abs (extremes[UP]) > -extremes[DOWN])
+       return DOWN;
+      else if (extremes[UP] < -extremes[DOWN])
+       return UP;
+    }
+  
   Direction dir = CENTER;
   Direction d = CENTER;
   if ((d = (Direction) sign (count[UP] - count[DOWN])))
@@ -784,7 +817,7 @@ Beam::calc_least_squares_positions (SCM smob, SCM posns)
   Interval pos (0,0);
   if (count < 1)
     return ly_interval2scm (pos);
-  
+
   vector<Real> x_posns;
   extract_grob_set (me, "stems", stems);
   Grob *commonx = common_refpoint_of_array (stems, me, X_AXIS);
@@ -796,7 +829,7 @@ 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);
 
@@ -1054,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
@@ -1068,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)
@@ -1080,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;
 }
@@ -1101,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");
   
@@ -1120,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);
@@ -1158,7 +1191,7 @@ 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);
 
@@ -1175,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_.at (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_.at (i).beams_i_drul_[-d]);
+               count = min (count, beaming->beamlet_count (i,-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, b, d);
+             Stem::set_beaming (stem, count, d);
            }
        }
       while (flip (&d) != LEFT);
@@ -1409,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 "