]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/beam.cc
Merge branch 'lilypond/translation' of ssh://git.sv.gnu.org/srv/git/lilypond into...
[lilypond.git] / lily / beam.cc
index f39bf9f62c06e55f0914bf5e83004ad775b96bde..1ecd71a88ff5a87027c7aa987c67292a5954196d 100644 (file)
@@ -37,6 +37,7 @@
 
 #include "beam.hh"
 
+#include "align-interface.hh"
 #include "beam-scoring-problem.hh"
 #include "beaming-pattern.hh"
 #include "directional-element-interface.hh"
@@ -53,6 +54,7 @@
 #include "pointer-group-interface.hh"
 #include "rhythmic-head.hh"
 #include "spanner.hh"
+#include "staff-symbol.hh"
 #include "staff-symbol-referencer.hh"
 #include "stem.hh"
 #include "warn.hh"
 
 #include <map>
 
-
 Beam_stem_segment::Beam_stem_segment ()
 {
-  max_connect_ = 1000;         // infinity
+  max_connect_ = 1000;          // infinity
   stem_ = 0;
   width_ = 0.0;
   stem_x_ = 0.0;
@@ -77,7 +78,7 @@ Beam_stem_segment::Beam_stem_segment ()
 }
 
 bool
-beam_segment_less (Beam_segment const& a, Beam_segment const& b)
+beam_segment_less (Beam_segment const &a, Beam_segment const &b)
 {
   return a.horizontal_[LEFT] < b.horizontal_[LEFT];
 }
@@ -93,7 +94,7 @@ Beam::add_stem (Grob *me, Grob *s)
   if (Stem::get_beam (s))
     {
       programming_error ("Stem already has beam");
-      return ;
+      return;
     }
 
   Pointer_group_interface::add_grob (me, ly_symbol2scm ("stems"), s);
@@ -105,7 +106,7 @@ Real
 Beam::get_beam_thickness (Grob *me)
 {
   return robust_scm2double (me->get_property ("beam-thickness"), 0)
-    * Staff_symbol_referencer::staff_space (me);
+         * Staff_symbol_referencer::staff_space (me);
 }
 
 /* Return the translation between 2 adjoining beams. */
@@ -119,8 +120,8 @@ Beam::get_beam_translation (Grob *me)
   Real fract = robust_scm2double (me->get_property ("length-fraction"), 1.0);
 
   Real beam_translation = beam_count < 4
-    ? (2 * staff_space + line - beam_thickness) / 2.0
-    : (3 * staff_space + line - beam_thickness) / 3.0;
+                          ? (2 * staff_space + line - beam_thickness) / 2.0
+                          : (3 * staff_space + line - beam_thickness) / 3.0;
 
   return fract * beam_translation;
 }
@@ -149,7 +150,7 @@ Beam::calc_normal_stems (SCM smob)
   extract_grob_set (me, "stems", stems);
   SCM val = Grob_array::make_array ();
   Grob_array *ga = unsmob_grob_array (val);
-  for (vsize i = 0; i < stems.size ();  i++)
+  for (vsize i = 0; i < stems.size (); i++)
     if (Stem::is_normal_stem (stems[i]))
       ga->add (stems[i]);
 
@@ -176,33 +177,33 @@ Beam::calc_direction (SCM smob)
     {
       extract_grob_set (me, "stems", stems);
       if (stems.size () == 0)
-       {
-         me->warning (_ ("removing beam with no stems"));
-         me->suicide ();
+        {
+          me->warning (_ ("removing beam with no stems"));
+          me->suicide ();
 
-         return SCM_UNSPECIFIED;
-       }
+          return SCM_UNSPECIFIED;
+        }
       else
-       {
-         Grob *stem = first_normal_stem (me);
-
-         /*
-           This happens for chord tremolos.
-         */
-         if (!stem)
-           stem = stems[0];
-
-         if (is_direction (stem->get_property_data ("direction")))
-           dir = to_dir (stem->get_property_data ("direction"));
-         else
-           dir = to_dir (stem->get_property ("default-direction"));
-       }
+        {
+          Grob *stem = first_normal_stem (me);
+
+          /*
+            This happens for chord tremolos.
+          */
+          if (!stem)
+            stem = stems[0];
+
+          if (is_direction (stem->get_property_data ("direction")))
+            dir = to_dir (stem->get_property_data ("direction"));
+          else
+            dir = to_dir (stem->get_property ("default-direction"));
+        }
     }
 
   if (count >= 1)
     {
       if (!dir)
-       dir = get_default_dir (me);
+        dir = get_default_dir (me);
 
       consider_auto_knees (me);
     }
@@ -215,8 +216,6 @@ Beam::calc_direction (SCM smob)
   return scm_from_int (dir);
 }
 
-
-
 /* We want a maximal number of shared beams, but if there is choice, we
  * take the one that is closest to the end of the stem. This is for
  * situations like
@@ -231,8 +230,8 @@ Beam::calc_direction (SCM smob)
  */
 int
 position_with_maximal_common_beams (SCM left_beaming, SCM right_beaming,
-                                   Direction left_dir,
-                                   Direction right_dir)
+                                    Direction left_dir,
+                                    Direction right_dir)
 {
   Slice lslice = int_list_to_slice (scm_cdr (left_beaming));
 
@@ -243,17 +242,17 @@ position_with_maximal_common_beams (SCM left_beaming, SCM right_beaming,
     {
       int count = 0;
       for (SCM s = scm_car (right_beaming); scm_is_pair (s); s = scm_cdr (s))
-       {
-         int k = -right_dir * scm_to_int (scm_car (s)) + i;
-         if (scm_c_memq (scm_from_int (k), left_beaming) != SCM_BOOL_F)
-           count++;
-       }
+        {
+          int k = -right_dir * scm_to_int (scm_car (s)) + i;
+          if (scm_c_memq (scm_from_int (k), left_beaming) != SCM_BOOL_F)
+            count++;
+        }
 
       if (count >= best_count)
-       {
-         best_count = count;
-         best_start = i;
-       }
+        {
+          best_count = count;
+          best_start = i;
+        }
     }
 
   return best_start;
@@ -279,51 +278,51 @@ Beam::calc_beaming (SCM smob)
 
       Direction this_dir = get_grob_direction (this_stem);
       if (scm_is_pair (last_beaming) && scm_is_pair (this_beaming))
-       {
-         int start_point = position_with_maximal_common_beams
-           (last_beaming, this_beaming,
-            last_dir ? last_dir : this_dir,
-            this_dir);
-
-         Direction d = LEFT;
-         Slice new_slice;
-         do
-           {
-             new_slice.set_empty ();
-             SCM s = index_get_cell (this_beaming, d);
-             for (; scm_is_pair (s); s = scm_cdr (s))
-               {
-                 int new_beam_pos
-                   = start_point - this_dir * scm_to_int (scm_car (s));
-
-                 new_slice.add_point (new_beam_pos);
-                 scm_set_car_x (s, scm_from_int (new_beam_pos));
-               }
-           }
-         while (flip (&d) != LEFT);
-
-         if (!new_slice.is_empty ())
-           last_int = new_slice;
-       }
+        {
+          int start_point = position_with_maximal_common_beams
+                            (last_beaming, this_beaming,
+                             last_dir ? last_dir : this_dir,
+                             this_dir);
+
+          Direction d = LEFT;
+          Slice new_slice;
+          do
+            {
+              new_slice.set_empty ();
+              SCM s = index_get_cell (this_beaming, d);
+              for (; scm_is_pair (s); s = scm_cdr (s))
+                {
+                  int new_beam_pos
+                    = start_point - this_dir * scm_to_int (scm_car (s));
+
+                  new_slice.add_point (new_beam_pos);
+                  scm_set_car_x (s, scm_from_int (new_beam_pos));
+                }
+            }
+          while (flip (&d) != LEFT);
+
+          if (!new_slice.is_empty ())
+            last_int = new_slice;
+        }
       else
-       {
-         /*
-           FIXME: what's this for?
-          */
-         SCM s = scm_cdr (this_beaming);
-         for (; scm_is_pair (s); s = scm_cdr (s))
-           {
-             int np = -this_dir * scm_to_int (scm_car (s));
-             scm_set_car_x (s, scm_from_int (np));
-             last_int.add_point (np);
-           }
-       }
+        {
+          /*
+            FIXME: what's this for?
+           */
+          SCM s = scm_cdr (this_beaming);
+          for (; scm_is_pair (s); s = scm_cdr (s))
+            {
+              int np = -this_dir * scm_to_int (scm_car (s));
+              scm_set_car_x (s, scm_from_int (np));
+              last_int.add_point (np);
+            }
+        }
 
       if (scm_ilength (scm_cdr (this_beaming)) > 0)
-       {
-         last_beaming = this_beaming;
-         last_dir = this_dir;
-       }
+        {
+          last_beaming = this_beaming;
+          last_dir = this_dir;
+        }
     }
 
   return SCM_EOL;
@@ -331,12 +330,12 @@ Beam::calc_beaming (SCM smob)
 
 bool
 operator <(Beam_stem_segment const &a,
-          Beam_stem_segment const &b)
+           Beam_stem_segment const &b)
 {
   return a.rank_ < b.rank_;
 }
 
-typedef map<int, vector<Beam_stem_segment> >  Position_stem_segments_map;
+typedef map<int, vector<Beam_stem_segment> > Position_stem_segments_map;
 
 // TODO - should store result in a property?
 vector<Beam_segment>
@@ -346,7 +345,7 @@ Beam::get_beam_segments (Grob *me_grob, Grob **common)
      Stem #'beaming is correct */
   (void) me_grob->get_property ("beaming");
 
-  Spanner *me = dynamic_cast<Spanner*> (me_grob);
+  Spanner *me = dynamic_cast<Spanner *> (me_grob);
 
   extract_grob_set (me, "stems", stems);
   Grob *commonx = common_refpoint_of_array (stems, me, X_AXIS);
@@ -377,50 +376,50 @@ Beam::get_beam_segments (Grob *me_grob, Grob **common)
       SCM beaming = stem->get_property ("beaming");
       Direction d = LEFT;
       do
-       {
-         // Find the maximum and minimum beam ranks.
-         // Given that RANKS is never reset to empty, the interval will always be
-         // smallest for the left beamlet of the first stem, and then it might grow.
-         // Do we really want this? (It only affects the tremolo gaps) --jneem
-         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);
-           }
-       }
+        {
+          // Find the maximum and minimum beam ranks.
+          // Given that RANKS is never reset to empty, the interval will always be
+          // smallest for the left beamlet of the first stem, and then it might grow.
+          // Do we really want this? (It only affects the tremolo gaps) --jneem
+          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));
+                       Drul_array<Real> (-0.5, 0.0));
 
   vector<Beam_segment> segments;
   for (Position_stem_segments_map::const_iterator i (stem_segments.begin ());
@@ -434,113 +433,113 @@ Beam::get_beam_segments (Grob *me_grob, Grob **common)
       // Iterate over all of the segments of the current beam rank,
       // merging the adjacent Beam_stem_segments into one Beam_segment
       // when appropriate.
-      int vertical_count =  (*i).first;
+      int vertical_count = (*i).first;
       for (vsize j = 0; j < segs.size (); j++)
-       {
-         // Keeping track of the different directions here is a little tricky.
-         // segs[j].dir_ is the direction of the beam segment relative to the stem
-         // (ie. segs[j].dir_ == LEFT if the beam segment sticks out to the left of
-         // its stem) whereas event_dir refers to the edge of the beam segment that
-         // we are currently looking at (ie. if segs[j].dir_ == event_dir then we
-         // are looking at that edge of the beam segment that is furthest from its
-         // stem).
-         Direction event_dir = LEFT;
-         Beam_stem_segment const& seg = segs[j];
-         do
-           {
-             Beam_stem_segment const& neighbor_seg = segs[j + event_dir];
-             // TODO: make names clearer? --jneem
-             // on_line_bound: whether the current segment is on the boundary of the WHOLE beam
-             // on_beam_bound: whether the current segment is on the boundary of just that part
-             //   of the beam with the current beam_rank
-             bool on_line_bound = (seg.dir_ == LEFT) ? seg.stem_index_ == 0
-               : seg.stem_index_ == stems.size() - 1;
-             bool on_beam_bound = (event_dir == LEFT) ? j == 0 :
-               j == segs.size () - 1;
-             bool inside_stem = (event_dir == LEFT)
-               ? seg.stem_index_ > 0
-               : seg.stem_index_ + 1 < stems.size () ;
-
-             bool event = on_beam_bound
-               || abs (seg.rank_ - neighbor_seg.rank_) > 1
-               || (abs (vertical_count) >= seg.max_connect_
-                   || abs (vertical_count) >= neighbor_seg.max_connect_);
-
-             if (!event)
-               // Then this edge of the current segment is irrelevent because it will
-               // be connected with the next segment in the event_dir direction.
-               continue;
-
-             current.vertical_count_ = vertical_count;
-             current.horizontal_[event_dir] = seg.stem_x_;
-             if (seg.dir_ == event_dir)
-               // then we are examining the edge of a beam segment that is furthest
-               // from its stem.
-               {
-                 if (on_line_bound
-                     && me->get_bound (event_dir)->break_status_dir ())
-                   {
-                     current.horizontal_[event_dir]
-                       = (robust_relative_extent (me->get_bound (event_dir),
-                                                  commonx, X_AXIS)[RIGHT]
-                          + event_dir * break_overshoot[event_dir]);
-                   }
-                 else
-                   {
-                     Grob *stem = stems[seg.stem_index_];
-                     Drul_array<Real> beamlet_length =
-                       robust_scm2interval (stem->get_property ("beamlet-default-length"), Interval (1.1, 1.1));
-                     Drul_array<Real> max_proportion =
-                       robust_scm2interval (stem->get_property ("beamlet-max-length-proportion"), Interval (0.75, 0.75));
-                     Real length = beamlet_length[seg.dir_];
-
-                     if (inside_stem)
-                       {
-                         Grob *neighbor_stem = stems[seg.stem_index_ + event_dir];
-                         Real neighbor_stem_x = neighbor_stem->relative_coordinate (commonx, X_AXIS);
-
-                         length = min (length,
-                                       fabs (neighbor_stem_x - seg.stem_x_) * max_proportion[seg.dir_]);
-                       }
-                     current.horizontal_[event_dir] += event_dir * length;
-                   }
-               }
-             else
-               // we are examining the edge of a beam segment that is closest
-               // (ie. touching, unless there is a gap) its stem.
-               {
-                 current.horizontal_[event_dir] += event_dir * seg.width_/2;
-                 if (seg.gapped_)
-                   {
-                     current.horizontal_[event_dir] -= event_dir * gap_length;
-
-                     if (Stem::is_invisible (seg.stem_))
-                       {
-                         /*
-                           Need to do this in case of whole notes. We don't want the
-                           heads to collide with the beams.
-                          */
-                         extract_grob_set (seg.stem_, "note-heads", heads);
-
-                         for (vsize k = 0; k < heads.size (); k ++)
-                           current.horizontal_[event_dir]
-                             = event_dir * min  (event_dir * current.horizontal_[event_dir],
-                                                 - gap_length/2
-                                                 + event_dir
-                                                   * heads[k]->extent (commonx,
-                                                                       X_AXIS)[-event_dir]);
-                       }
-                   }
-               }
-
-             if (event_dir == RIGHT)
-               {
-                 segments.push_back (current);
-                 current = Beam_segment ();
-               }
-           }
-         while (flip (&event_dir) != LEFT);
-       }
+        {
+          // Keeping track of the different directions here is a little tricky.
+          // segs[j].dir_ is the direction of the beam segment relative to the stem
+          // (ie. segs[j].dir_ == LEFT if the beam segment sticks out to the left of
+          // its stem) whereas event_dir refers to the edge of the beam segment that
+          // we are currently looking at (ie. if segs[j].dir_ == event_dir then we
+          // are looking at that edge of the beam segment that is furthest from its
+          // stem).
+          Direction event_dir = LEFT;
+          Beam_stem_segment const &seg = segs[j];
+          do
+            {
+              Beam_stem_segment const &neighbor_seg = segs[j + event_dir];
+              // TODO: make names clearer? --jneem
+              // on_line_bound: whether the current segment is on the boundary of the WHOLE beam
+              // on_beam_bound: whether the current segment is on the boundary of just that part
+              //   of the beam with the current beam_rank
+              bool on_line_bound = (seg.dir_ == LEFT) ? seg.stem_index_ == 0
+                                   : seg.stem_index_ == stems.size () - 1;
+              bool on_beam_bound = (event_dir == LEFT) ? j == 0
+                                   : j == segs.size () - 1;
+              bool inside_stem = (event_dir == LEFT)
+                                 ? seg.stem_index_ > 0
+                                 : seg.stem_index_ + 1 < stems.size ();
+
+              bool event = on_beam_bound
+                           || abs (seg.rank_ - neighbor_seg.rank_) > 1
+                           || (abs (vertical_count) >= seg.max_connect_
+                               || abs (vertical_count) >= neighbor_seg.max_connect_);
+
+              if (!event)
+                // Then this edge of the current segment is irrelevent because it will
+                // be connected with the next segment in the event_dir direction.
+                continue;
+
+              current.vertical_count_ = vertical_count;
+              current.horizontal_[event_dir] = seg.stem_x_;
+              if (seg.dir_ == event_dir)
+                // then we are examining the edge of a beam segment that is furthest
+                // from its stem.
+                {
+                  if (on_line_bound
+                      && me->get_bound (event_dir)->break_status_dir ())
+                    {
+                      current.horizontal_[event_dir]
+                        = (robust_relative_extent (me->get_bound (event_dir),
+                                                   commonx, X_AXIS)[RIGHT]
+                           + event_dir * break_overshoot[event_dir]);
+                    }
+                  else
+                    {
+                      Grob *stem = stems[seg.stem_index_];
+                      Drul_array<Real> beamlet_length
+                        = robust_scm2interval (stem->get_property ("beamlet-default-length"), Interval (1.1, 1.1));
+                      Drul_array<Real> max_proportion
+                        = robust_scm2interval (stem->get_property ("beamlet-max-length-proportion"), Interval (0.75, 0.75));
+                      Real length = beamlet_length[seg.dir_];
+
+                      if (inside_stem)
+                        {
+                          Grob *neighbor_stem = stems[seg.stem_index_ + event_dir];
+                          Real neighbor_stem_x = neighbor_stem->relative_coordinate (commonx, X_AXIS);
+
+                          length = min (length,
+                                        fabs (neighbor_stem_x - seg.stem_x_) * max_proportion[seg.dir_]);
+                        }
+                      current.horizontal_[event_dir] += event_dir * length;
+                    }
+                }
+              else
+                // we are examining the edge of a beam segment that is closest
+                // (ie. touching, unless there is a gap) its stem.
+                {
+                  current.horizontal_[event_dir] += event_dir * seg.width_ / 2;
+                  if (seg.gapped_)
+                    {
+                      current.horizontal_[event_dir] -= event_dir * gap_length;
+
+                      if (Stem::is_invisible (seg.stem_))
+                        {
+                          /*
+                            Need to do this in case of whole notes. We don't want the
+                            heads to collide with the beams.
+                           */
+                          extract_grob_set (seg.stem_, "note-heads", heads);
+
+                          for (vsize k = 0; k < heads.size (); k++)
+                            current.horizontal_[event_dir]
+                              = event_dir * min (event_dir * current.horizontal_[event_dir],
+                                                 - gap_length / 2
+                                                 + event_dir
+                                                 * heads[k]->extent (commonx,
+                                                                     X_AXIS)[-event_dir]);
+                        }
+                    }
+                }
+
+              if (event_dir == RIGHT)
+                {
+                  segments.push_back (current);
+                  current = Beam_segment ();
+                }
+            }
+          while (flip (&event_dir) != LEFT);
+        }
 
     }
 
@@ -554,6 +553,8 @@ Beam::print (SCM grob)
   Spanner *me = unsmob_spanner (grob);
   Grob *commonx = 0;
   vector<Beam_segment> segments = get_beam_segments (me, &commonx);
+  if (!segments.size ())
+    return SCM_EOL;
 
   Interval span;
   if (normal_stem_count (me))
@@ -583,7 +584,7 @@ Beam::print (SCM grob)
   scale_drul (&pos, Staff_symbol_referencer::staff_space (me));
 
   Real dy = pos[RIGHT] - pos[LEFT];
-  Real slope = (dy && span.length ()) ? dy / span.length ()  : 0;
+  Real slope = (dy && span.length ()) ? dy / span.length () : 0;
 
   Real beam_thickness = get_beam_thickness (me);
   Real beam_dy = get_beam_translation (me);
@@ -598,7 +599,7 @@ Beam::print (SCM grob)
                  ? segments[0].vertical_count_
                  : segments.back ().vertical_count_);
 
-  for (vsize i = 0; i < segments.size (); i ++)
+  for (vsize i = 0; i < segments.size (); i++)
     {
       Real local_slope = slope;
       /*
@@ -607,8 +608,8 @@ Beam::print (SCM grob)
       */
       if (feather_dir)
         local_slope += (feather_dir * segments[i].vertical_count_
-                                    * beam_dy
-                                    * placements.length ()
+                        * beam_dy
+                        * placements.length ()
                         / span.length ());
 
       Stencil b = Lookup::beam (local_slope, segments[i].horizontal_.length (), beam_thickness, blot);
@@ -665,15 +666,15 @@ Beam::print (SCM grob)
       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.
+        This code prints the demerits for each beam. Perhaps this
+        should be switchable for those who want to twiddle with the
+        parameters.
       */
       string str;
       SCM properties = Font_interface::text_font_alist_chain (me);
 
-      properties = scm_cons(scm_acons (ly_symbol2scm ("font-size"), scm_from_int (-5), SCM_EOL),
-                            properties);
+      properties = scm_cons (scm_acons (ly_symbol2scm ("font-size"), scm_from_int (-5), SCM_EOL),
+                             properties);
 
       Direction stem_dir = stems.size () ? to_dir (stems[0]->get_property ("direction")) : UP;
 
@@ -681,10 +682,10 @@ Beam::print (SCM grob)
                                        (me->layout ()->self_scm (), properties, annotation));
 
       if (!score.is_empty ())
-       {
-         score.translate_axis (me->relative_coordinate(commonx, X_AXIS), X_AXIS);
-         the_beam.add_at_edge (Y_AXIS, stem_dir, score, 1.0);
-       }
+        {
+          score.translate_axis (me->relative_coordinate (commonx, X_AXIS), X_AXIS);
+          the_beam.add_at_edge (Y_AXIS, stem_dir, score, 1.0);
+        }
     }
 #endif
 
@@ -703,10 +704,10 @@ Beam::get_default_dir (Grob *me)
       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]);
-       }
+        {
+          if (sign (positions[d]) == d)
+            extremes[d] = d * max (d * positions[d], d * extremes[d]);
+        }
       while (flip (&d) != DOWN);
     }
 
@@ -720,30 +721,29 @@ Beam::get_default_dir (Grob *me)
       Direction stem_dir = CENTER;
       SCM stem_dir_scm = s->get_property_data ("direction");
       if (is_direction (stem_dir_scm))
-       {
-         stem_dir = to_dir (stem_dir_scm);
-         force_dir = true;
-       }
+        {
+          stem_dir = to_dir (stem_dir_scm);
+          force_dir = true;
+        }
       else
-       stem_dir = to_dir (s->get_property ("default-direction"));
+        stem_dir = to_dir (s->get_property ("default-direction"));
 
       if (!stem_dir)
-       stem_dir = to_dir (s->get_property ("neutral-direction"));
+        stem_dir = to_dir (s->get_property ("neutral-direction"));
 
       if (stem_dir)
-       {
-         count[stem_dir] ++;
-         total[stem_dir] += max (int (- stem_dir * Stem::head_positions (s) [-stem_dir]), 0);
-       }
+        {
+          count[stem_dir]++;
+          total[stem_dir] += max (int (- stem_dir * Stem::head_positions (s) [-stem_dir]), 0);
+        }
     }
 
-
   if (!force_dir)
     {
       if (abs (extremes[UP]) > -extremes[DOWN])
-       return DOWN;
+        return DOWN;
       else if (extremes[UP] < -extremes[DOWN])
-       return UP;
+        return UP;
     }
 
   Direction dir = CENTER;
@@ -751,10 +751,10 @@ Beam::get_default_dir (Grob *me)
   if ((d = (Direction) sign (count[UP] - count[DOWN])))
     dir = d;
   else if (count[UP]
-          && count[DOWN]
-          && (d = (Direction)  sign (total[UP] / count[UP] - total[DOWN]/count[DOWN])))
+           && count[DOWN]
+           && (d = (Direction) sign (total[UP] / count[UP] - total[DOWN] / count[DOWN])))
     dir = d;
-  else if ((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"));
@@ -776,7 +776,7 @@ Beam::set_stem_directions (Grob *me, Direction d)
 
       SCM forcedir = s->get_property_data ("direction");
       if (!to_dir (forcedir))
-       set_grob_direction (s, d);
+        set_grob_direction (s, d);
     }
 }
 
@@ -815,24 +815,24 @@ Beam::consider_auto_knees (Grob *me)
 
       Interval head_extents = Stem::head_positions (stem);
       if (!head_extents.is_empty ())
-       {
-         head_extents[LEFT] += -1;
-         head_extents[RIGHT] += 1;
-         head_extents *= staff_space * 0.5;
-
-         /*
-           We could subtract beam Y position, but this routine only
-           sets stem directions, a constant shift does not have an
-           influence.
-         */
-         head_extents += stem->pure_relative_y_coordinate (common, 0, INT_MAX);
-
-         if (to_dir (stem->get_property_data ("direction")))
-           {
-             Direction stemdir = to_dir (stem->get_property ("direction"));
-             head_extents[-stemdir] = -stemdir * infinity_f;
-           }
-       }
+        {
+          head_extents[LEFT] += -1;
+          head_extents[RIGHT] += 1;
+          head_extents *= staff_space * 0.5;
+
+          /*
+            We could subtract beam Y position, but this routine only
+            sets stem directions, a constant shift does not have an
+            influence.
+          */
+          head_extents += stem->pure_relative_y_coordinate (common, 0, INT_MAX);
+
+          if (to_dir (stem->get_property_data ("direction")))
+            {
+              Direction stemdir = to_dir (stem->get_property ("direction"));
+              head_extents[-stemdir] = -stemdir * infinity_f;
+            }
+        }
       head_extents_array.push_back (head_extents);
 
       gaps.remove_interval (head_extents);
@@ -841,46 +841,46 @@ Beam::consider_auto_knees (Grob *me)
   Interval max_gap;
   Real max_gap_len = 0.0;
 
-  for (vsize i = gaps.allowed_regions_.size () -1; i != VPOS ;i--)
+  for (vsize i = gaps.allowed_regions_.size () - 1; i != VPOS; i--)
     {
       Interval gap = gaps.allowed_regions_[i];
 
       /*
-       the outer gaps are not knees.
+        the outer gaps are not knees.
       */
       if (isinf (gap[LEFT]) || isinf (gap[RIGHT]))
-       continue;
+        continue;
 
       if (gap.length () >= max_gap_len)
-       {
-         max_gap_len = gap.length ();
-         max_gap = gap;
-       }
+        {
+          max_gap_len = gap.length ();
+          max_gap = gap;
+        }
     }
 
   Real beam_translation = get_beam_translation (me);
   Real beam_thickness = Beam::get_beam_thickness (me);
   int beam_count = Beam::get_beam_count (me);
   Real height_of_beams = beam_thickness / 2
-    + (beam_count - 1) * beam_translation;
+                         + (beam_count - 1) * beam_translation;
   Real threshold = scm_to_double (scm) + height_of_beams;
 
   if (max_gap_len > threshold)
     {
       int j = 0;
       for (vsize i = 0; i < stems.size (); i++)
-       {
-         Grob *stem = stems[i];
-         Interval head_extents = head_extents_array[j++];
+        {
+          Grob *stem = stems[i];
+          Interval head_extents = head_extents_array[j++];
 
-         Direction d = (head_extents.center () < max_gap.center ())
-           ? UP : DOWN;
+          Direction d = (head_extents.center () < max_gap.center ())
+                        ? UP : DOWN;
 
-         stem->set_property ("direction", scm_from_int (d));
+          stem->set_property ("direction", scm_from_int (d));
 
-         head_extents.intersect (max_gap);
-         assert (head_extents.is_empty () || head_extents.length () < 1e-6);
-       }
+          head_extents.intersect (max_gap);
+          assert (head_extents.is_empty () || head_extents.length () < 1e-6);
+        }
     }
 }
 
@@ -894,16 +894,14 @@ This is done in beam because the shorten has to be uniform over the
 entire beam.
 */
 
-
-
 void
 set_minimum_dy (Grob *me, Real *dy)
 {
   if (*dy)
     {
       /*
-       If dy is smaller than the smallest quant, we
-       get absurd direction-sign penalties.
+        If dy is smaller than the smallest quant, we
+        get absurd direction-sign penalties.
       */
 
       Real ss = Staff_symbol_referencer::staff_space (me);
@@ -914,12 +912,10 @@ set_minimum_dy (Grob *me, Real *dy)
       Real hang = 1.0 - (beam_thickness - slt) / 2;
 
       *dy = sign (*dy) * max (fabs (*dy),
-                             min (min (sit, inter), hang));
+                              min (min (sit, inter), hang));
     }
 }
 
-
-
 MAKE_SCHEME_CALLBACK (Beam, calc_stem_shorten, 1)
 SCM
 Beam::calc_stem_shorten (SCM smob)
@@ -933,7 +929,7 @@ Beam::calc_stem_shorten (SCM smob)
     return scm_from_int (0);
 
   Real forced_fraction = 1.0 * forced_stem_count (me)
-    / normal_stem_count (me);
+                         / normal_stem_count (me);
 
   int beam_count = get_beam_count (me);
 
@@ -944,19 +940,17 @@ Beam::calc_stem_shorten (SCM smob)
   Real staff_space = Staff_symbol_referencer::staff_space (me);
 
   SCM shorten_elt
-    = robust_list_ref (beam_count -1, shorten_list);
+    = robust_list_ref (beam_count - 1, shorten_list);
   Real shorten = scm_to_double (shorten_elt) * staff_space;
 
   shorten *= forced_fraction;
 
-
   if (shorten)
     return scm_from_double (shorten);
 
   return scm_from_double (0.0);
 }
 
-
 Interval
 Beam::no_visible_stem_positions (Grob *me, Interval default_value)
 {
@@ -966,22 +960,25 @@ Beam::no_visible_stem_positions (Grob *me, Interval default_value)
 
   Interval head_positions;
   Slice multiplicity;
-  for (vsize i = 0; i < stems.size(); i++)
+  for (vsize i = 0; i < stems.size (); i++)
     {
       head_positions.unite (Stem::head_positions (stems[i]));
       multiplicity.unite (Stem::beam_multiplicity (stems[i]));
     }
 
   Direction dir = get_grob_direction (me);
-  Real y = head_positions[dir]
-    * 0.5 * Staff_symbol_referencer::staff_space (me)
-    + dir * get_beam_translation (me) * (multiplicity.length () + 1);
+
+  if (!dir)
+    programming_error ("The beam should have a direction by now.");
+
+  Real y = head_positions.linear_combination (dir)
+           * 0.5 * Staff_symbol_referencer::staff_space (me)
+           + dir * get_beam_translation (me) * (multiplicity.length () + 1);
 
   y /= Staff_symbol_referencer::staff_space (me);
-  return Interval (y,y);
+  return Interval (y, y);
 }
 
-
 /*
   Compute a first approximation to the beam slope.
 */
@@ -992,7 +989,7 @@ Beam::calc_least_squares_positions (SCM smob, SCM /* posns */)
   Grob *me = unsmob_grob (smob);
 
   int count = normal_stem_count (me);
-  Interval pos (0,0);
+  Interval pos (0, 0);
   if (count < 1)
     return ly_interval2scm (no_visible_stem_positions (me, pos));
 
@@ -1007,9 +1004,9 @@ Beam::calc_least_squares_positions (SCM smob, SCM /* posns */)
   Grob *lvs = last_normal_stem (me);
 
   Interval ideal (Stem::get_stem_info (fvs).ideal_y_
-                 + fvs->relative_coordinate (commony, Y_AXIS) - my_y,
-                 Stem::get_stem_info (lvs).ideal_y_
-                 + lvs->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_normal_stem (me)->relative_coordinate (commonx, X_AXIS);
   for (vsize i = 0; i < stems.size (); i++)
@@ -1028,29 +1025,29 @@ Beam::calc_least_squares_positions (SCM smob, SCM /* posns */)
   if (!ideal.delta ())
     {
       Interval chord (Stem::chord_start_y (stems[0]),
-                     Stem::chord_start_y (stems.back ()));
+                      Stem::chord_start_y (stems.back ()));
 
       /* Simple beams (2 stems) on middle line should be allowed to be
-        slightly sloped.
+         slightly sloped.
 
-        However, if both stems reach middle line,
-        ideal[LEFT] == ideal[RIGHT] and ideal.delta () == 0.
+         However, if both stems reach middle line,
+         ideal[LEFT] == ideal[RIGHT] and ideal.delta () == 0.
 
-        For that case, we apply artificial slope */
+         For that case, we apply artificial slope */
       if (!ideal[LEFT] && chord.delta () && count == 2)
-       {
-         /* FIXME. -> UP */
-         Direction d = (Direction) (sign (chord.delta ()) * UP);
-         pos[d] = get_beam_thickness (me) / 2;
-         pos[-d] = -pos[d];
-       }
+        {
+          /* FIXME. -> UP */
+          Direction d = (Direction) (sign (chord.delta ()) * UP);
+          pos[d] = get_beam_thickness (me) / 2;
+          pos[-d] = -pos[d];
+        }
       else
-       pos = ideal;
+        pos = ideal;
 
       /*
-       For broken beams this doesn't work well. In this case, the
-       slope esp. of the first part of a broken beam should predict
-       where the second part goes.
+        For broken beams this doesn't work well. In this case, the
+        slope esp. of the first part of a broken beam should predict
+        where the second part goes.
       */
       ldy = pos[RIGHT] - pos[LEFT];
     }
@@ -1058,13 +1055,13 @@ Beam::calc_least_squares_positions (SCM smob, SCM /* posns */)
     {
       vector<Offset> ideals;
       for (vsize i = 0; i < stems.size (); i++)
-       {
-         Grob *s = stems[i];
-         ideals.push_back (Offset (x_posns[i],
-                              Stem::get_stem_info (s).ideal_y_
-                              + s->relative_coordinate (commony, Y_AXIS)
-                              - my_y));
-       }
+        {
+          Grob *s = stems[i];
+          ideals.push_back (Offset (x_posns[i],
+                                    Stem::get_stem_info (s).ideal_y_
+                                    + s->relative_coordinate (commony, Y_AXIS)
+                                    - my_y));
+        }
 
       minimise_least_squares (&slope, &y, ideals);
 
@@ -1081,11 +1078,10 @@ Beam::calc_least_squares_positions (SCM smob, SCM /* posns */)
   */
   scale_drul (&pos, 1 / Staff_symbol_referencer::staff_space (me));
 
-  me->set_property ("least-squares-dy",  scm_from_double (ldy));
+  me->set_property ("least-squares-dy", scm_from_double (ldy));
   return ly_interval2scm (pos);
 }
 
-
 // Assuming V is not empty, pick a 'reasonable' point inside V.
 static Real
 point_in_interval (Interval v, Real dist)
@@ -1119,10 +1115,11 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
   extract_grob_set (me, "covered-grobs", covered);
 
   Grob *common[NO_AXES] = { me, me };
-  for (Axis a = X_AXIS; a < NO_AXES; incr (a)) {
-    common[a] = common_refpoint_of_array (stems, me, a);
-    common[a] = common_refpoint_of_array (covered, common[a], a);
-  }
+  for (Axis a = X_AXIS; a < NO_AXES; incr (a))
+    {
+      common[a] = common_refpoint_of_array (stems, me, a);
+      common[a] = common_refpoint_of_array (covered, common[a], a);
+    }
   Grob *fvs = first_normal_stem (me);
 
   if (!fvs)
@@ -1159,20 +1156,20 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
     {
       Grob *s = stems[i];
       if (Stem::is_invisible (s))
-       continue;
+        continue;
 
       Direction d = get_grob_direction (s);
       Real left_y
-       = Stem::get_stem_info (s).shortest_y_
-       - slope * x_posns [i];
+        = Stem::get_stem_info (s).shortest_y_
+          - slope * x_posns [i];
 
       /*
-       left_y is now relative to the stem S. We want relative to
-       ourselves, so translate:
+        left_y is now relative to the stem S. We want relative to
+        ourselves, so translate:
       */
       left_y
-       += + s->relative_coordinate (common[Y_AXIS], Y_AXIS)
-       - me->relative_coordinate (common[Y_AXIS], Y_AXIS);
+      += + s->relative_coordinate (common[Y_AXIS], Y_AXIS)
+         - me->relative_coordinate (common[Y_AXIS], Y_AXIS);
 
       Interval flp;
       flp.set_full ();
@@ -1181,14 +1178,7 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
       feasible_left_point.intersect (flp);
     }
 
-  /*
-    We have two intervals here, one for the up variant (beams goes
-    over the collision) one for the down.
-  */
-  Drul_array<Interval> collision_free (feasible_left_point,
-                                       feasible_left_point);
-
-  vector<Grob*> filtered;
+  vector<Grob *> filtered;
   /*
     We only update these for objects that are too large for quanting
     to find a workaround.  Typically, these are notes with
@@ -1200,11 +1190,18 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
     take care of computing the impact those exactly.
   */
   Real min_y_size = 2.0;
-  for (vsize i = 0; i < covered.size(); i++)
+
+  // A list of intervals into which beams may not fall
+  vector<Interval> forbidden_intervals;
+
+  for (vsize i = 0; i < covered.size (); i++)
     {
-      if (!covered[i]->is_live())
+      if (!covered[i]->is_live ())
         continue;
-      
+
+      if (Beam::has_interface (covered[i]) && is_cross_staff (covered[i]))
+        continue;
+
       Box b;
       for (Axis a = X_AXIS; a < NO_AXES; incr (a))
         b[a] = covered[i]->extent (common[a], a);
@@ -1233,19 +1230,19 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
                 would resolve the problem, eg.
 
                  x    x
-                |    | 
+                |    |
                 =====
 
                 =====
-                |   |  
+                |   |
                 x   x
-                
+
                 Such beams would need a coordinating grob to resolve
                 the collision, since both will likely want to occupy
                 the centerline.
               */
               Direction stemdir = get_grob_direction (head_stem);
-              b[Y_AXIS][stemdir] = stemdir * infinity_f; 
+              b[Y_AXIS][stemdir] = stemdir * infinity_f;
             }
           else
             {
@@ -1263,65 +1260,97 @@ Beam::shift_region_to_valid (SCM grob, SCM posns)
           Real dy = slope * x;
 
           Direction yd = DOWN;
+          Interval disallowed;
           do
             {
               Real left_y = b[Y_AXIS][yd];
 
-              if (left_y == yd * infinity_f)
-                {
-                  collision_free[yd].set_empty ();
-                  continue;
-                }
-
               left_y -= dy;
 
               // Translate back to beam as ref point.
               left_y -= me->relative_coordinate (common[Y_AXIS], Y_AXIS);
-            
-              Interval allowed;
-              allowed.set_full ();
 
-              allowed[-yd] = left_y;
-              collision_free[yd].intersect (allowed);
+              disallowed[yd] = left_y;
             }
           while (flip (&yd) != DOWN);
+
+          forbidden_intervals.push_back (disallowed);
         }
       while (flip (&d) != LEFT);
     }
 
-  Grob_array *arr = 
-    Pointer_group_interface::get_grob_array (me,
-                                             ly_symbol2scm ("covered-grobs"));
+  Grob_array *arr
+    Pointer_group_interface::get_grob_array (me,
+                                               ly_symbol2scm ("covered-grobs"));
   arr->set_array (filtered);
 
-  if (collision_free[DOWN].contains (beam_left_y)
-      || collision_free[UP].contains (beam_left_y))
+  vector_sort (forbidden_intervals, Interval::left_less);
+  Real epsilon = 1.0e-10;
+  Interval feasible_beam_placements (beam_left_y, beam_left_y);
+
+  /*
+    forbidden_intervals contains a vector of intervals in which
+    the beam cannot start.  it iterates through these intervals,
+    pushing feasible_beam_placements epsilon over or epsilon under a
+    collision.  when this type of change happens, the loop is marked
+    as "dirty" and re-iterated.
+
+    TODO: figure out a faster ways that this loop can happen via
+    a better search algorithm and/or OOP.
+  */
+
+  bool dirty = false;
+  do
     {
-      // We're good to go. Do nothing.
+      dirty = false;
+      for (vsize i = 0; i < forbidden_intervals.size (); i++)
+        {
+          Direction d = DOWN;
+          do
+            {
+              if (forbidden_intervals[i][d] == d * infinity_f)
+                feasible_beam_placements[d] = d * infinity_f;
+              else if (forbidden_intervals[i].contains (feasible_beam_placements[d]))
+                {
+                  feasible_beam_placements[d] = d * epsilon + forbidden_intervals[i][d];
+                  dirty = true;
+                }
+            }
+          while (flip (&d) != DOWN);
+        }
     }
-  else if (!collision_free[DOWN].is_empty ()
-           || !collision_free[UP].is_empty ())
-    {
-      // We have space above or below collisions (or, no collisions at
-      // all).
-      Interval best =  
-        (collision_free[DOWN].length () > collision_free[UP].length ()) ?
-        collision_free[DOWN] : collision_free[UP];
+  while (dirty);
 
-      beam_left_y = point_in_interval (best, 2.0);
+  // if the beam placement falls out of the feasible region, we push it
+  // to infinity so that it can never be a feasible candidate below
+  Direction d = DOWN;
+  do
+    {
+      if (!feasible_left_point.contains (feasible_beam_placements[d]))
+        feasible_beam_placements[d] = d * infinity_f;
     }
-  else if (!feasible_left_point.is_empty ())
+  while (flip (&d) != DOWN);
+
+  if ((feasible_beam_placements[UP] == infinity_f && feasible_beam_placements[DOWN] == -infinity_f) && !feasible_left_point.is_empty ())
     {
       // We are somewhat screwed: we have a collision, but at least
       // there is a way to satisfy stem length constraints.
       beam_left_y = point_in_interval (feasible_left_point, 2.0);
     }
+  else if (!feasible_left_point.is_empty ())
+    {
+      // Only one of them offers is feasible solution. Pick that one.
+      if (abs (beam_left_y - feasible_beam_placements[DOWN]) > abs (beam_left_y - feasible_beam_placements[UP]))
+        beam_left_y = feasible_beam_placements[UP];
+      else
+        beam_left_y = feasible_beam_placements[DOWN];
+    }
   else
     {
       // We are completely screwed.
-      warning (_ ("no viable initial configuration found: may not find good beam slope"));
+      me->warning (_ ("no viable initial configuration found: may not find good beam slope"));
     }
-  
+
   pos = Drul_array<Real> (beam_left_y, (beam_left_y + beam_dy));
   scale_drul (&pos, 1 / Staff_symbol_referencer::staff_space (me));
 
@@ -1363,7 +1392,7 @@ Beam::slope_damping (SCM smob, SCM posns)
       Grob *commonx = fvs->common_refpoint (lvs, X_AXIS);
 
       Real dx = last_normal_stem (me)->relative_coordinate (commonx, X_AXIS)
-       - first_normal_stem (me)->relative_coordinate (commonx, X_AXIS);
+                - first_normal_stem (me)->relative_coordinate (commonx, X_AXIS);
 
       Real slope = dy && dx ? dy / dx : 0;
 
@@ -1382,13 +1411,12 @@ Beam::slope_damping (SCM smob, SCM posns)
   return ly_interval2scm (pos);
 }
 
-
 MAKE_SCHEME_CALLBACK (Beam, quanting, 2);
 SCM
 Beam::quanting (SCM smob, SCM posns)
 {
   Grob *me = unsmob_grob (smob);
-  Drul_array<Real> ys(0, 0);
+  Drul_array<Real> ys (0, 0);
   ys = robust_scm2drul (posns, ys);
   Beam_scoring_problem problem (me, ys);
 
@@ -1396,7 +1424,6 @@ Beam::quanting (SCM smob, SCM posns)
   return ly_interval2scm (ys);
 }
 
-
 /*
   Report slice containing the numbers that are both in (car BEAMING)
   and (cdr BEAMING)
@@ -1410,7 +1437,7 @@ where_are_the_whole_beams (SCM beaming)
     {
       if (scm_c_memq (scm_car (s), scm_cdr (beaming)) != SCM_BOOL_F)
 
-       l.add_point (scm_to_int (scm_car (s)));
+        l.add_point (scm_to_int (scm_car (s)));
     }
 
   return l;
@@ -1420,27 +1447,27 @@ where_are_the_whole_beams (SCM beaming)
    in POS for stem S.  This Y position is relative to S. */
 Real
 Beam::calc_stem_y (Grob *me, Grob *stem, Grob **common,
-                  Real xl, Real xr, Direction feather_dir,
-                  Drul_array<Real> pos, bool french)
+                   Real xl, Real xr, Direction feather_dir,
+                   Drul_array<Real> pos, bool french)
 {
   Real beam_translation = get_beam_translation (me);
   Direction stem_dir = get_grob_direction (stem);
 
   Real dx = xr - xl;
-  Real relx = dx ? (stem->relative_coordinate (common[X_AXIS], X_AXIS) - xl)/dx : 0;
-  Real xdir = 2*relx-1;
+  Real relx = dx ? (stem->relative_coordinate (common[X_AXIS], X_AXIS) - xl) / dx : 0;
+  Real xdir = 2 * relx - 1;
 
-  Real stem_y = linear_combination(pos, xdir);
+  Real stem_y = linear_combination (pos, xdir);
 
   SCM beaming = stem->get_property ("beaming");
 
   Slice beam_slice (french
-                   ? where_are_the_whole_beams (beaming)
-                   : Stem::beam_multiplicity (stem));
+                    ? where_are_the_whole_beams (beaming)
+                    : Stem::beam_multiplicity (stem));
   if (beam_slice.is_empty ())
-    beam_slice = Slice (0,0);
-  Interval beam_multiplicity(beam_slice[LEFT],
-                            beam_slice[RIGHT]);
+    beam_slice = Slice (0, 0);
+  Interval beam_multiplicity (beam_slice[LEFT],
+                              beam_slice[RIGHT]);
 
   /*
     feather dir = 1 , relx 0->1 : factor 0 -> 1
@@ -1454,9 +1481,9 @@ Beam::calc_stem_y (Grob *me, Grob *stem, Grob **common,
     feather_factor = 1 - relx;
 
   stem_y += feather_factor * beam_translation
-    * beam_multiplicity[Direction(((french) ? DOWN : UP)*stem_dir)];
+            * beam_multiplicity[Direction (((french) ? DOWN : UP) * stem_dir)];
   Real id = me->relative_coordinate (common[Y_AXIS], Y_AXIS)
-    - stem->relative_coordinate (common[Y_AXIS], Y_AXIS);
+            - stem->relative_coordinate (common[Y_AXIS], Y_AXIS);
 
   return stem_y + id;
 }
@@ -1510,22 +1537,22 @@ Beam::set_stem_lengths (SCM smob)
 
       bool french = to_boolean (s->get_property ("french-beaming"));
       Real stem_y = calc_stem_y (me, s, common,
-                                xl, xr, feather_dir,
-                                pos, french && s != lvs && s!= fvs);
+                                 xl, xr, feather_dir,
+                                 pos, french && s != lvs && s != fvs);
 
       /*
-       Make the stems go up to the end of the beam. This doesn't matter
-       for normal beams, but for tremolo beams it looks silly otherwise.
+        Make the stems go up to the end of the beam. This doesn't matter
+        for normal beams, but for tremolo beams it looks silly otherwise.
       */
       if (gap
-         && !Stem::is_invisible (s))
-       stem_y += thick * 0.5 * get_grob_direction (s);
+          && !Stem::is_invisible (s))
+        stem_y += thick * 0.5 * get_grob_direction (s);
 
       /*
-       Do set_stemend for invisible stems too, so tuplet brackets
-       have a reference point for sloping
+        Do set_stem_positions for invisible stems too, so tuplet brackets
+        have a reference point for sloping
        */
-      Stem::set_stemend (s, 2 * stem_y / staff_space);
+      Stem::set_stem_positions (s, 2 * stem_y / staff_space);
     }
 
   return posns;
@@ -1540,30 +1567,30 @@ Beam::set_beaming (Grob *me, Beaming_pattern const *beaming)
   for (vsize i = 0; i < stems.size (); i++)
     {
       /*
-       Don't overwrite user settings.
+        Don't overwrite user settings.
       */
       do
-       {
-         Grob *stem = stems[i];
-         SCM beaming_prop = stem->get_property ("beaming");
-         if (beaming_prop == SCM_EOL
-             || index_get_cell (beaming_prop, d) == SCM_EOL)
-           {
-             int count = beaming->beamlet_count (i, d);
-             if (i > 0
-                 && i + 1 < stems.size ()
-                 && Stem::is_invisible (stem))
-               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, count, d);
-           }
-       }
+        {
+          Grob *stem = stems[i];
+          SCM beaming_prop = stem->get_property ("beaming");
+          if (beaming_prop == SCM_EOL
+              || index_get_cell (beaming_prop, d) == SCM_EOL)
+            {
+              int count = beaming->beamlet_count (i, d);
+              if (i > 0
+                  && i + 1 < stems.size ()
+                  && Stem::is_invisible (stem))
+                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, count, d);
+            }
+        }
       while (flip (&d) != LEFT);
     }
 }
@@ -1579,13 +1606,13 @@ Beam::forced_stem_count (Grob *me)
       Grob *s = stems[i];
 
       /* I can imagine counting those boundaries as a half forced stem,
-        but let's count them full for now. */
+         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
-         && defdir
-         && get_grob_direction (s) != defdir)
-       f++;
+          && defdir
+          && get_grob_direction (s) != defdir)
+        f++;
     }
   return f;
 }
@@ -1642,7 +1669,7 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
     return scm_from_double (0.0);
 
   Drul_array<Real> pos (robust_scm2drul (beam->get_property ("positions"),
-                                        Drul_array<Real> (0,0)));
+                                         Drul_array<Real> (0, 0)));
 
   Real staff_space = Staff_symbol_referencer::staff_space (rest);
 
@@ -1650,8 +1677,8 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
 
   Real dy = pos[RIGHT] - pos[LEFT];
 
-  Drul_array<Grob*> visible_stems (first_normal_stem (beam),
-                                  last_normal_stem (beam));
+  Drul_array<Grob *> visible_stems (first_normal_stem (beam),
+                                    last_normal_stem (beam));
   extract_grob_set (beam, "stems", stems);
 
   Grob *common = common_refpoint_of_array (stems, beam, X_AXIS);
@@ -1662,7 +1689,7 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
 
   Direction d = get_grob_direction (stem);
   Real stem_y = pos[LEFT]
-    + (stem->relative_coordinate (common, X_AXIS) - x0) * slope;
+                + (stem->relative_coordinate (common, X_AXIS) - x0) * slope;
 
   Real beam_translation = get_beam_translation (beam);
   Real beam_thickness = Beam::get_beam_thickness (beam);
@@ -1674,7 +1701,7 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
     = Stem::beam_multiplicity (stem).length () + 1;
 
   Real height_of_my_beams = beam_thickness / 2
-    + (beam_count - 1) * beam_translation;
+                            + (beam_count - 1) * beam_translation;
   Real beam_y = stem_y - d * height_of_my_beams;
 
   Grob *common_y = rest->common_refpoint (beam, Y_AXIS);
@@ -1685,7 +1712,7 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
   Real rest_dim = rest_extent[d];
   Real minimum_distance
     = staff_space * (robust_scm2double (stem->get_property ("stemlet-length"), 0.0)
-                    + robust_scm2double (rest->get_property ("minimum-distance"), 0.0));
+                     + robust_scm2double (rest->get_property ("minimum-distance"), 0.0));
 
   Real shift = d * min (d * (beam_y - d * minimum_distance - rest_dim), 0.0);
 
@@ -1705,6 +1732,77 @@ Beam::rest_collision_callback (SCM smob, SCM prev_offset)
   return scm_from_double (offset + staff_space * shift);
 }
 
+MAKE_SCHEME_CALLBACK_WITH_OPTARGS (Beam, pure_rest_collision_callback, 4, 1, "");
+SCM
+Beam::pure_rest_collision_callback (SCM smob,
+                                    SCM, /* prev_offset */
+                                    SCM, /* start */
+                                    SCM /* end */)
+{
+  Real amount = 0.0;
+
+  Grob *me = unsmob_grob (smob);
+  Grob *stem = unsmob_grob (me->get_object ("stem"));
+  if (!stem)
+    return scm_from_double (amount);
+  Grob *beam = unsmob_grob (stem->get_object ("beam"));
+  if (!beam
+      || !Beam::normal_stem_count (beam))
+    return scm_from_double (amount);
+
+  Real ss = Staff_symbol_referencer::staff_space (me);
+
+  /*
+    This gives the extrema of rest positions.
+    In general, beams are never typeset more than one staff space away
+    from the staff in either direction.
+  */
+  Grob *staff = Staff_symbol_referencer::get_staff_symbol (me);
+  Interval rest_max_pos = staff ? Staff_symbol::line_span (staff) : Interval (0.0, 0.0);
+  rest_max_pos.widen (1);
+  rest_max_pos *= ss / 2;
+
+  extract_grob_set (beam, "stems", stems);
+  vector<Grob *> my_stems;
+
+  for (vsize i = 0; i < stems.size (); i++)
+    if (Stem::head_count (stems[i]) || stems[i] == stem)
+      my_stems.push_back (stems[i]);
+
+  vsize idx = -1;
+
+  for (vsize i = 0; i < my_stems.size (); i++)
+    if (my_stems[i] == stem)
+      {
+        idx = i;
+        break;
+      }
+  Grob *left;
+  Grob *right;
+
+  if (idx == (vsize)-1 || my_stems.size () == 1)
+    return scm_from_double (amount);
+  else if (idx == 0)
+    left = right = my_stems[1];
+  else if (idx == my_stems.size () - 1)
+    left = right = my_stems[idx - 1];
+  else
+    {
+      left = my_stems[idx - 1];
+      right = my_stems[idx + 1];
+    }
+  Direction beamdir = get_grob_direction (beam);
+  /*
+    Take the position between the two bounding head_positions,
+    then bound it by the minimum and maximum positions outside the staff.
+    4.0 = 2.0 to get out of staff space * 2.0 for the average
+  */
+  amount = min (max ((Stem::head_positions (left)[beamdir] + Stem::head_positions (right)[beamdir]) / 4.0, rest_max_pos[DOWN]), rest_max_pos[UP]);
+
+  return scm_from_double (amount);
+}
+
+
 bool
 Beam::is_knee (Grob *me)
 {
@@ -1719,10 +1817,10 @@ Beam::is_knee (Grob *me)
     {
       Direction dir = get_grob_direction (stems[i]);
       if (d && d != dir)
-       {
-         knee = true;
-         break;
-       }
+        {
+          knee = true;
+          break;
+        }
       d = dir;
     }
 
@@ -1758,22 +1856,22 @@ Beam::get_direction_beam_count (Grob *me, Direction d)
   for (vsize i = stems.size (); i--;)
     {
       /*
-       Should we take invisible stems into account?
+        Should we take invisible stems into account?
       */
       if (get_grob_direction (stems[i]) == d)
-       bc = max (bc, (Stem::beam_multiplicity (stems[i]).length () + 1));
+        bc = max (bc, (Stem::beam_multiplicity (stems[i]).length () + 1));
     }
 
   return bc;
 }
 
 ADD_INTERFACE (Beam,
-              "A beam.\n"
-              "\n"
-              "The @code{beam-thickness} property is the weight of beams,"
-              " measured in staffspace.  The @code{direction} property is"
-              " not user-serviceable.  Use the @code{direction} property"
-              " of @code{Stem} instead.\n"
+               "A beam.\n"
+               "\n"
+               "The @code{beam-thickness} property is the weight of beams,"
+               " measured in staffspace.  The @code{direction} property is"
+               " not user-serviceable.  Use the @code{direction} property"
+               " of @code{Stem} instead.\n"
                "\n"
                "The following properties may be set in the @code{details}"
                " list.\n"
@@ -1809,32 +1907,32 @@ ADD_INTERFACE (Beam,
                " calculating direction penalties.\n"
                "@end table\n",
 
-              /* properties */
-              "annotation "
-              "auto-knee-gap "
-              "beamed-stem-shorten "
-              "beaming "
-              "beam-thickness "
-              "break-overshoot "
-              "clip-edges "
-              "concaveness "
-              "collision-interfaces "
-              "collision-voice-only "
-              "covered-grobs "
-              "damping "
-              "details "
-              "direction "
-              "gap "
-              "gap-count "
-              "grow-direction "
-              "inspect-quants "
-              "knee "
-              "length-fraction "
-              "least-squares-dy "
-              "neutral-direction "
-              "normal-stems "
-              "positions "
-              "quantized-positions "
-              "shorten "
-              "stems "
-              );
+               /* properties */
+               "annotation "
+               "auto-knee-gap "
+               "beamed-stem-shorten "
+               "beaming "
+               "beam-thickness "
+               "break-overshoot "
+               "clip-edges "
+               "concaveness "
+               "collision-interfaces "
+               "collision-voice-only "
+               "covered-grobs "
+               "damping "
+               "details "
+               "direction "
+               "gap "
+               "gap-count "
+               "grow-direction "
+               "inspect-quants "
+               "knee "
+               "length-fraction "
+               "least-squares-dy "
+               "neutral-direction "
+               "normal-stems "
+               "positions "
+               "quantized-positions "
+               "shorten "
+               "stems "
+              );