]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/beam.cc
*** empty log message ***
[lilypond.git] / lily / beam.cc
index 655e760322dc8ccb4240e390318be11b01598f50..2f403f306722ad21b4996447596c04cfe99f8752 100644 (file)
@@ -3,17 +3,13 @@
   
   source file of the GNU LilyPond music typesetter
   
-  (c)  1997--2002 Han-Wen Nienhuys <hanwen@cs.uu.nl>
+  (c)  1997--2003 Han-Wen Nienhuys <hanwen@cs.uu.nl>
   Jan Nieuwenhuizen <janneke@gnu.org>
-  
 */
 
 /*
 TODO:
 
-
-  * Junk stem_info.
-
   * Use Number_pair i.s.o Interval to represent (yl, yr).
   
   - Determine auto knees based on positions if it's set by the user.
@@ -63,40 +59,49 @@ Beam::add_stem (Grob *me, Grob *s)
   
   s->add_dependency (me);
 
-  assert (!Stem::beam_l (s));
+  assert (!Stem::get_beam (s));
   s->set_grob_property ("beam", me->self_scm ());
 
   add_bound_item (dynamic_cast<Spanner*> (me), dynamic_cast<Item*> (s));
 }
 
 
-/*
-  this returns the translation between 2 adjoining beams.
- */
+Real
+Beam::get_thickness (Grob * me)
+{
+  SCM th = me->get_grob_property ("thickness");
+  if (gh_number_p (th))
+    return gh_scm2double (th)* Staff_symbol_referencer::staff_space (me);
+  else
+    return 0.0;
+}
+
+/* Return the translation between 2 adjoining beams. */
 Real
 Beam::get_beam_translation (Grob *me)
 {
   SCM func = me->get_grob_property ("space-function");
-  SCM s = gh_call2 (func, me->self_scm (), gh_int2scm (get_beam_count (me)));
+  SCM s = gh_call2 (func, me->self_scm (), scm_int2num (get_beam_count (me)));
   return gh_scm2double (s);
 }
 
-/*
-  Maximum beam_count.
- */
+/* Maximum beam_count. */
 int
 Beam::get_beam_count (Grob *me) 
 {
   int m = 0;
   for (SCM s = me->get_grob_property ("stems"); gh_pair_p (s); s = ly_cdr (s))
     {
-      Grob *sc = unsmob_grob (ly_car (s));
-      
-      m = m >? (Stem::beam_multiplicity (sc).length () + 1);
+      Grob *stem = unsmob_grob (ly_car (s));
+      m = m >? (Stem::beam_multiplicity (stem).length () + 1);
     }
   return m;
 }
 
+
+/*
+  Space return space between beams.
+ */
 MAKE_SCHEME_CALLBACK (Beam, space_function, 2);
 SCM
 Beam::space_function (SCM smob, SCM beam_count)
@@ -104,9 +109,8 @@ Beam::space_function (SCM smob, SCM beam_count)
   Grob *me = unsmob_grob (smob);
   
   Real staff_space = Staff_symbol_referencer::staff_space (me);
-  Real line = me->paper_l ()->get_var ("linethickness");
-  Real thickness = gh_scm2double (me->get_grob_property ("thickness"))
-    * staff_space;
+  Real line = me->get_paper ()->get_var ("linethickness");
+  Real thickness = get_thickness (me);
   
   Real beam_translation = gh_scm2int (beam_count) < 4
     ? (2*staff_space + line - thickness) / 2.0
@@ -146,7 +150,7 @@ Beam::before_line_breaking (SCM smob)
        {
          me->warning (_ ("Beam has less than two stems. Removing beam."));
 
-         unsmob_grob (gh_car (stems))->remove_grob_property ("beam");
+         unsmob_grob (gh_car (stems))->set_grob_property ("beam", SCM_EOL);
          me->suicide ();
 
          return SCM_UNSPECIFIED;
@@ -161,7 +165,7 @@ Beam::before_line_breaking (SCM smob)
     {
       Direction d = get_default_dir (me);
 
-      consider_auto_knees (me, d);
+      consider_auto_knees (me);
       set_stem_directions (me, d);
 
       connect_beams (me);
@@ -173,6 +177,49 @@ Beam::before_line_breaking (SCM smob)
 }
 
 
+/*
+  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
+
+       x
+      |
+      |
+  |===|
+  |=
+  |
+  x
+
+  
+ */
+int
+position_with_maximal_common_beams (SCM left_beaming, SCM right_beaming,
+                                   Direction left_dir,
+                                   Direction right_dir)
+{
+  Slice lslice = int_list_to_slice (gh_cdr (left_beaming));
+
+  int best_count = 0;
+  int best_start = 0;
+  for (int i = lslice[-left_dir];
+       (i - lslice[left_dir])* left_dir <= 0 ; i+= left_dir) 
+    {
+      int count =0;
+      for ( SCM s = gh_car (right_beaming); gh_pair_p (s); s = gh_cdr (s))
+       {
+         int k = - right_dir * gh_scm2int (gh_car (s)) + i;
+         if (scm_memq (scm_int2num (k), left_beaming) != SCM_BOOL_F)
+           count ++;
+       }
+
+      if (count >= best_count)
+       {
+         best_count = count; 
+         best_start = i;
+       }
+    }
+
+  return best_start;
+}
 
 void
 Beam::connect_beams (Grob *me)
@@ -182,15 +229,19 @@ Beam::connect_beams (Grob *me)
 
   Slice last_int;
   last_int.set_empty();
+  SCM last_beaming = SCM_EOL;
+  Direction last_dir = CENTER;
   for (int i = 0; i< stems.size(); i++)
     {
       Grob *this_stem = stems[i];
       SCM this_beaming = this_stem->get_grob_property ("beaming");
 
       Direction this_dir = Directional_element_interface::get(this_stem);
-      if (i > 0)
+      if (gh_pair_p (last_beaming) && gh_pair_p (this_beaming))
        {
-         int start_point = last_int [this_dir];
+         int start_point = position_with_maximal_common_beams
+           (last_beaming, this_beaming,
+            last_dir, this_dir);
          
          Direction d = LEFT;
          Slice new_slice ; 
@@ -207,8 +258,10 @@ Beam::connect_beams (Grob *me)
                    start_point - this_dir * gh_scm2int (gh_car (s));
 
                  new_slice.add_point (new_beam_pos);
-                 gh_set_car_x (s, gh_int2scm (new_beam_pos));
+                 gh_set_car_x (s, scm_int2num (new_beam_pos));
                }
+
+
            }
          while (flip (&d) != LEFT);
 
@@ -222,14 +275,20 @@ Beam::connect_beams (Grob *me)
          for (; gh_pair_p (s); s = gh_cdr (s))
            {
              int np = - this_dir * gh_scm2int (gh_car(s));
-             gh_set_car_x (s, gh_int2scm (np));
+             gh_set_car_x (s, scm_int2num (np));
              last_int.add_point (np);
            }
        }
 
       if (i == stems.size () -1)
        {
-         gh_set_cdr_x ( this_beaming, SCM_EOL);
+         gh_set_cdr_x (this_beaming, SCM_EOL);
+       }
+
+      if (scm_ilength (gh_cdr (this_beaming)) > 0)
+       {
+         last_beaming = this_beaming;
+         last_dir = this_dir;
        }
     }
  }
@@ -269,7 +328,7 @@ Beam::brew_molecule (SCM grob)
   Real dy = pos.delta ();
   Real dydx = dy && dx ? dy/dx : 0;
   
-  Real thick = gh_scm2double (me->get_grob_property ("thickness"));
+  Real thick = get_thickness (me);
   Real bdy = get_beam_translation (me);
 
   SCM last_beaming = SCM_EOL;;
@@ -279,7 +338,8 @@ Beam::brew_molecule (SCM grob)
 
   SCM gap = me->get_grob_property ("gap");
   Molecule the_beam;
-  Real lt = me->paper_l ()->get_var ("linethickness");
+  Real lt = me->get_paper ()->get_var ("linethickness");
+  
   for (int i = 0; i< stems.size(); i++)
     {
       Grob * st =stems[i];
@@ -288,6 +348,11 @@ Beam::brew_molecule (SCM grob)
       Real xposn = st->relative_coordinate (xcommon, X_AXIS);
       Real stem_width = gh_scm2double (st->get_grob_property ("thickness")) *lt;
 
+      /*
+       We do the space left of ST, with lfliebertjes pointing to the
+       right from the left stem, and rfliebertjes pointing left from
+       right stem.
+       */
       if (i > 0)
        {
          SCM left = gh_cdr (last_beaming);
@@ -353,18 +418,14 @@ Beam::brew_molecule (SCM grob)
 
          if (lfliebertjes.size() || rfliebertjes.size())
            {
-
              Real nw_f;
-             if (!Stem::first_head (st))
-               nw_f = 0;
-             else
-               {
-                 int t = Stem::duration_log (st); 
 
-                 SCM proc = me->get_grob_property ("flag-width-function");
-                 SCM result = gh_call1 (proc, gh_int2scm (t));
-                 nw_f = gh_scm2double (result);
-               }
+             int t = Stem::duration_log (st); 
+
+             SCM proc = me->get_grob_property ("flag-width-function");
+             SCM result = gh_call1 (proc, scm_int2num (t));
+             nw_f = gh_scm2double (result);
+               
              
              /* Half beam should be one note-width,
                 but let's make sure two half-beams never touch */
@@ -408,17 +469,17 @@ Beam::brew_molecule (SCM grob)
       String str;
       if (1)
        {
-         str += to_str (gh_scm2int (me->get_grob_property ("best-idx")));
+         str += to_string (gh_scm2int (me->get_grob_property ("best-idx")));
          str += ":";
        }
-      str += to_str (gh_scm2double (me->get_grob_property ("quant-score")),
+      str += to_string (gh_scm2double (me->get_grob_property ("quant-score")),
                     "%.2f");
 
       SCM properties = Font_interface::font_alist_chain (me);
 
       
-      Molecule tm = Text_item::text2molecule (me, ly_str02scm (str.ch_C ()), properties);
-      the_beam.add_at_edge (Y_AXIS, UP, tm, 5.0);
+      Molecule tm = Text_item::text2molecule (me, scm_makfrom0str (str.to_str0 ()), properties);
+      the_beam.add_at_edge (Y_AXIS, UP, tm, 5.0, 0);
     }
 #endif
     
@@ -439,8 +500,8 @@ Beam::get_default_dir (Grob *me)
   count[UP]  = count[DOWN] = 0;
   Direction d = DOWN;
 
-  Link_array<Item> stems=
-       Pointer_group_interface__extract_grobs (me, (Item*)0, "stems");
+  Link_array<Grob> stems=
+       Pointer_group_interface__extract_grobs (me, (Grob*)0, "stems");
 
   for (int i=0; i <stems.size (); i++)
     do {
@@ -459,10 +520,10 @@ Beam::get_default_dir (Grob *me)
   
   SCM func = me->get_grob_property ("dir-function");
   SCM s = gh_call2 (func,
-                   gh_cons (gh_int2scm (count[UP]),
-                            gh_int2scm (count[DOWN])),
-                   gh_cons (gh_int2scm (total[UP]),
-                            gh_int2scm (total[DOWN])));
+                   gh_cons (scm_int2num (count[UP]),
+                            scm_int2num (count[DOWN])),
+                   gh_cons (scm_int2num (total[UP]),
+                            scm_int2num (total[DOWN])));
 
   if (gh_number_p (s) && gh_scm2int (s))
     return to_dir (s);
@@ -478,97 +539,181 @@ Beam::get_default_dir (Grob *me)
 void
 Beam::set_stem_directions (Grob *me, Direction d)
 {
-  Link_array<Item> stems
-    =Pointer_group_interface__extract_grobs (me, (Item*) 0, "stems");
+  Link_array<Grob> stems
+    =Pointer_group_interface__extract_grobs (me, (Grob*) 0, "stems");
   
   for (int i=0; i <stems.size (); i++)
     {
       Grob *s = stems[i];
-      /* For knees, non-forced stems should probably have their
-        natural direction. In any case, when knee, beam direction is
-        foe. */
-      if (knee_b(me))
-       Stem::get_direction (s); // this actually sets it, if necessary
-      else
-       {
-         SCM force = s->remove_grob_property ("dir-forced");
-         if (!gh_boolean_p (force) || !gh_scm2bool (force))
-           Directional_element_interface::set (s, d);
-       }
+  
+      SCM forcedir = s->get_grob_property ("direction");
+      if (!to_dir (forcedir))
+       Directional_element_interface::set (s, d);
     }
-} 
+}
+
+/*
+  A union of intervals in the real line.
 
-/* Simplistic auto-knees; only consider vertical gap between two
-   adjacent chords.
+  Abysmal performance (quadratic) for large N, hopefully we don't have
+  that large N. In any case, this should probably be rewritten to use
+  a balanced tree.
+ */
+struct Int_set
+{
+  Array<Interval> allowed_regions_;
+
+  Int_set()
+  {
+    set_full();
+  }
 
-   This may decide for a knee that's impossible to fit sane scoring
-   criteria (eg, stem lengths).  We may need something smarter. */
+  void set_full()
+  {
+    allowed_regions_.clear();
+    Interval s;
+    s.set_full ();
+    allowed_regions_.push (s);
+  }
+
+  void remove_interval (Interval rm)
+  {
+    for (int i = 0; i < allowed_regions_.size(); )
+      {
+       Interval s = rm;
+
+       s.intersect (allowed_regions_[i]);
+
+       if (!s.empty_b ())
+         {
+           Interval before = allowed_regions_[i];
+           Interval after = allowed_regions_[i];
+
+           before[RIGHT] = s[LEFT];
+           after[LEFT] = s[RIGHT];
+
+           if (!before.empty_b() && before.length () > 0.0)
+             {
+               allowed_regions_.insert (before, i);
+               i++;
+             }
+           allowed_regions_.del (i);
+           if (!after.empty_b () && after.length () > 0.0)
+             {
+               allowed_regions_.insert (after, i);
+               i++;
+             }
+         }
+       else
+         i++;
+      }
+  }
+};
+
+
+/*
+  Only try horizontal beams for knees.  No reliable detection of
+  anything else is possible here, since we don't know funky-beaming
+  settings, or X-distances (slopes!)  People that want sloped
+  knee-beams, should set the directions manually.
+ */
 void
-Beam::consider_auto_knees (Grob *me, Direction d)
+Beam::consider_auto_knees (Grob* me)
 {
   SCM scm = me->get_grob_property ("auto-knee-gap");
-
   if (!gh_number_p (scm))
-    return;
-  
-  bool knee_b = false;
+    return ;
+
+  Real threshold = gh_scm2double (scm);
   
-  Real staff_space = Staff_symbol_referencer::staff_space (me);
-  Real gap = gh_scm2double (scm) / staff_space;
+  Int_set gaps;
+
+  gaps.set_full ();
 
   Link_array<Grob> stems=
     Pointer_group_interface__extract_grobs (me, (Grob*)0, "stems");
       
   Grob *common = common_refpoint_of_array (stems, me,  Y_AXIS);
-
-  int l = 0;
-  for (int r=1; r < stems.size (); r++)
+  Real staff_space = Staff_symbol_referencer::staff_space (me);
+  
+  Array<Interval> hps_array;  
+  for (int i=0; i < stems.size (); i++)
     {
-      if (!Stem::invisible_b (stems[r-1]))
-       l = r - 1;
-      Grob *right = stems[r];
-      Grob *left = stems[l];
-      if (Stem::invisible_b (left))
+      Grob* stem = stems[i];
+      if (Stem::invisible_b (stem))
        continue;
-      if (Stem::invisible_b (right))
-       continue;
-         
-      Real left_y = Stem::extremal_heads (left)[d]
-       ->relative_coordinate (common, Y_AXIS);
-      Real right_y = Stem::extremal_heads (right)[-d]
-       ->relative_coordinate (common, Y_AXIS);
 
-      Real dy = right_y - left_y;
-
-      if (abs (dy) >= gap)
+      Interval hps = Stem::head_positions (stem);
+      if(!hps.empty_b())
        {
-         knee_b = true;
-         Direction knee_dir = (right_y > left_y ? UP : DOWN);
-         if (!Stem::invisible_b (left)
-             && left->get_grob_property ("dir-forced") != SCM_BOOL_T)
-           {
-             Directional_element_interface::set (left, knee_dir);
-             left->set_grob_property ("dir-forced", SCM_BOOL_T);
+         hps[LEFT] += -1;
+         hps[RIGHT] += 1; 
+         hps *= staff_space * 0.5 ;
 
-           }
-         if (!Stem::invisible_b (right)
-             && stems[r]->get_grob_property ("dir-forced") != SCM_BOOL_T)
+         /*
+           We could subtract beam Y position, but this routine only
+           sets stem directions, a constant shift does not have an
+           influence.
+           
+          */
+         hps += stem->relative_coordinate (common, Y_AXIS);
+
+         if (to_dir (stem->get_grob_property ("direction")))
            {
-             Directional_element_interface::set (right, -knee_dir);
-             right->set_grob_property ("dir-forced", SCM_BOOL_T);
+             Direction stemdir = to_dir (stem->get_grob_property ("direction"));
+             hps[-stemdir] = - stemdir * infinity_f;
            }
        }
+      hps_array.push (hps);
+
+      gaps.remove_interval (hps);
     }
 
-  if (knee_b)
+  Interval max_gap;
+  Real max_gap_len =0.0;
+
+  for (int i  = gaps.allowed_regions_.size() -1;  i >=  0 ; i--)
     {
-      me->set_grob_property ("knee", SCM_BOOL_T);
-       
-      for (int i=0; i < stems.size (); i++)
-       stems[i]->set_grob_property ("stem-info", SCM_EOL);
+      Interval gap = gaps.allowed_regions_[i];
+
+      /*
+       the outer gaps are not knees.
+       */
+      if (isinf (gap[LEFT]) || isinf(gap[RIGHT]))
+       continue;
+      
+      if (gap.length () >= max_gap_len)
+       {
+         max_gap_len = gap.length();
+         max_gap = gap;
+       }
+    }
+
+  if (max_gap_len > threshold)
+    {
+      int j = 0;
+      for (int i = 0; i < stems.size(); i++)
+       {
+         Grob* stem = stems[i];
+         if (Stem::invisible_b (stem))
+           continue;
+
+         Interval hps = hps_array[j++];
+
+
+         Direction d =  (hps.center () < max_gap.center()) ?
+           UP : DOWN ;
+         
+         stem->set_grob_property ("direction", scm_int2num (d));
+         
+         hps.intersect (max_gap);
+         assert (hps.empty_b () || hps.length () < 1e-6 );
+       }
     }
 }
 
+
+
 /* Set stem's shorten property if unset.
 
  TODO:
@@ -588,19 +733,19 @@ Beam::set_stem_shorten (Grob *me)
   if (knee_b(me))
     return ;
   
-  Real forced_fraction = forced_stem_count (me) / visible_stem_count (me);
+  Real forced_fraction = 1.0 * forced_stem_count (me)
+    / visible_stem_count (me);
 
   int beam_count = get_beam_count (me);
 
-  SCM shorten = me->get_grob_property ("beamed-stem-shorten");
-  if (shorten == SCM_EOL)
+  SCM shorten_list = me->get_grob_property ("beamed-stem-shorten");
+  if (shorten_list == SCM_EOL)
     return;
 
-  int sz = scm_ilength (shorten);
-  
   Real staff_space = Staff_symbol_referencer::staff_space (me);
-  SCM shorten_elt = scm_list_ref (shorten,
-                                 gh_int2scm (beam_count <? (sz - 1)));
+  
+  SCM shorten_elt =
+    robust_list_ref (beam_count -1, shorten_list);
   Real shorten_f = gh_scm2double (shorten_elt) * staff_space;
 
   /* your similar cute comment here */
@@ -637,6 +782,10 @@ Beam::after_line_breaking (SCM smob)
   return SCM_UNSPECIFIED;
 }
 
+
+/*
+  Compute  a first approximation to the beam slope.
+ */
 MAKE_SCHEME_CALLBACK (Beam, least_squares, 1);
 SCM
 Beam::least_squares (SCM smob)
@@ -664,9 +813,9 @@ Beam::least_squares (SCM smob)
   Grob *fvs  = first_visible_stem (me);
   Grob *lvs  = last_visible_stem (me);
   
-  Interval ideal (Stem::calc_stem_info (fvs).ideal_y_
+  Interval ideal (Stem::get_stem_info (fvs).ideal_y_
                  + fvs->relative_coordinate (commony, Y_AXIS) -my_y,
-                 Stem::calc_stem_info (lvs).ideal_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);
@@ -688,23 +837,18 @@ Beam::least_squares (SCM smob)
       Interval chord (Stem::chord_start_y (first_visible_stem (me)),
                      Stem::chord_start_y (last_visible_stem (me)));
 
+      /* Simple beams (2 stems) on middle line should be allowed to be
+        slightly sloped.
+        
+        However, if both stems reach middle line,
+        ideal[LEFT] == ideal[RIGHT] and ideal.delta () == 0.
 
-      /*
-       TODO -- use scoring for this.
-
-       complicated, because we take stem-info.ideal for determining
-       beam slopes.
-       */
-      /* Make simple beam on middle line have small tilt */
+        For that case, we apply artificial slope */
       if (!ideal[LEFT] && chord.delta () && count == 2)
        {
-
-         /*
-           FIXME. -> UP
-         */
+         /* FIXME. -> UP */
          Direction d = (Direction) (sign (chord.delta ()) * UP);
          pos[d] = gh_scm2double (me->get_grob_property ("thickness")) / 2;
-         //                * dir;
          pos[-d] = - pos[d];
        }
       else
@@ -725,7 +869,7 @@ Beam::least_squares (SCM smob)
          if (Stem::invisible_b (s))
            continue;
          ideals.push (Offset (x_posns[i],
-                              Stem::calc_stem_info (s).ideal_y_
+                              Stem::get_stem_info (s).ideal_y_
                               + s->relative_coordinate (commony, Y_AXIS)
                               - my_y));
        }
@@ -801,7 +945,7 @@ Beam::shift_region_to_valid (SCM grob)
       Direction d = Stem::get_direction (s);
 
       Real left_y =
-       Stem::calc_stem_info (s).shortest_y_
+       Stem::get_stem_info (s).shortest_y_
        - dydx * x_posns [i];
 
       /*
@@ -932,9 +1076,13 @@ Beam::check_concave (SCM smob)
       concave *= dir;
       concaveness2 = concave / (stems.size () - 2);
       
-      /* ugh: this is the a kludge to get
-        input/regression/beam-concave.ly to behave as
-        baerenreiter. */
+      /*
+
+      ugh: this is the a kludge to get
+      input/regression/beam-concave.ly to behave as
+      baerenreiter.
+
+      */
 
       /*
        huh? we're dividing twice (which is not scalable) meaning that
@@ -999,6 +1147,10 @@ Beam::slope_damping (SCM smob)
   return SCM_UNSPECIFIED;
 }
 
+/*
+  Report slice containing the numbers that are both in (car BEAMING)
+  and (cdr BEAMING)
+ */
 Slice
 where_are_the_whole_beams(SCM beaming)
 {
@@ -1014,10 +1166,8 @@ where_are_the_whole_beams(SCM beaming)
   return l;
 }
 
-/*
-  Calculate 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.
- */
+/* 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,
                   Real xl, Real xr,
@@ -1086,8 +1236,7 @@ Beam::set_stem_lengths (Grob *me)
       &&gh_scm2double (me->get_grob_property ("gap")))
   {
     gap = true;
-    thick = gh_scm2double (me->get_grob_property ("thickness"))
-      * Staff_symbol_referencer::staff_space(me);
+    thick = get_thickness(me);
   }
       
   // ugh -> use commonx
@@ -1138,13 +1287,18 @@ Beam::set_beaming (Grob *me, Beaming_info_list *beaming)
              ||(d == RIGHT && i == stems.size () -1))
            continue;
 
-
-         SCM beaming_prop = stems[i]->get_grob_property ("beaming");
+         Grob *st =  stems[i];
+         SCM beaming_prop = st->get_grob_property ("beaming");
          if (beaming_prop == SCM_EOL ||
              index_get_cell (beaming_prop, d) == SCM_EOL)
            {
              int b = beaming->infos_.elem (i).beams_i_drul_[d];
-             Stem::set_beaming (stems[i], b, d);
+             if (i>0
+                 && i < stems.size() -1
+                 && Stem::invisible_b (st))
+               b = b <? beaming->infos_.elem(i).beams_i_drul_[-d];
+             
+             Stem::set_beaming (st, b, d);
            }
        }
       while (flip (&d) != LEFT);
@@ -1164,7 +1318,9 @@ Beam::forced_stem_count (Grob *me)
       if (Stem::invisible_b (s))
        continue;
 
-      if (((int)Stem::chord_start_y (s)) 
+      /* I can imagine counting those boundaries as a half forced stem,
+        but let's count them full for now. */
+      if (abs (Stem::chord_start_y (s)) > 0.1
         && (Stem::get_direction (s) != Stem::get_default_dir (s)))
         f++;
     }
@@ -1245,10 +1401,6 @@ Beam::rest_collision_callback (SCM element_smob, SCM axis)
       || !Beam::visible_stem_count (beam))
     return gh_double2scm (0.0);
 
-  // make callback for rest from this.
-  // todo: make sure this calced already.
-
-  //  Interval pos = ly_scm2interval (beam->get_grob_property ("positions"));
   Interval pos (0, 0);
   SCM s = beam->get_grob_property ("positions");
   if (gh_pair_p (s) && gh_number_p (ly_car (s)))
@@ -1261,32 +1413,51 @@ Beam::rest_collision_callback (SCM element_smob, SCM axis)
   Real dydx = dy && dx ? dy/dx : 0;
   
   Direction d = Stem::get_direction (stem);
-  Real beamy = (stem->relative_coordinate (0, X_AXIS) - x0) * dydx + pos[LEFT];
+  Real stem_y = (pos[LEFT]
+                + (stem->relative_coordinate (0, X_AXIS) - x0) * dydx)
+    * d;
+  
+  Real beam_translation = get_beam_translation (beam);
+  Real beam_thickness = gh_scm2double (beam->get_grob_property ("thickness"));
+  int beam_count = get_direction_beam_count (beam, d);
+  Real height_of_my_beams = beam_thickness
+    + (beam_count - 1) * beam_translation;
+  Real beam_y = stem_y - height_of_my_beams + beam_thickness / 2.0;
 
   Real staff_space = Staff_symbol_referencer::staff_space (rest);
 
-  
-  Real rest_dim = rest->extent (rest, Y_AXIS)[d]*2.0 / staff_space; // refp??
+  /* Better calculate relative-distance directly, rather than using
+     rest_dim? */
+  Grob *common_x = rest->common_refpoint (beam, Y_AXIS);
+  Real rest_dim = rest->extent (common_x, Y_AXIS)[d] / staff_space * d;
 
-  Real minimum_dist
-    = gh_scm2double (rest->get_grob_property ("minimum-beam-collision-distance"));
-  Real dist =
-    minimum_dist +  -d  * (beamy - rest_dim) >? 0;
+  Real minimum_distance = gh_scm2double
+    (rest->get_grob_property ("minimum-beam-collision-distance"));
 
+  Real distance = beam_y - rest_dim;
+  Real shift = 0;
+  if (distance < 0)
+    shift = minimum_distance - distance;
+  else if (minimum_distance > distance)
+    shift = minimum_distance - distance;
+      
   int stafflines = Staff_symbol_referencer::line_count (rest);
 
-  // move discretely by half spaces.
-  int discrete_dist = int (ceil (dist));
+  /* Always move discretely by half spaces */
+  Real discrete_shift = ceil (shift * 2.0) / 2.0;
 
-  // move by whole spaces inside the staff.
-  if (discrete_dist < stafflines+1)
-    discrete_dist = int (ceil (discrete_dist / 2.0)* 2.0);
+  /* Inside staff, move by whole spaces*/
+  if ((rest->extent (common_x, Y_AXIS)[d] + discrete_shift) * d
+      < stafflines / 2.0
+      ||(rest->extent (common_x, Y_AXIS)[-d] + discrete_shift) * -d
+      < stafflines / 2.0)
+    discrete_shift = ceil (discrete_shift);
 
-  return gh_double2scm (-d *  discrete_dist);
+  return gh_double2scm (-d * discrete_shift);
 }
 
 bool
-Beam::knee_b (Grob*me)
+Beam::knee_b (Grob* me)
 {
   SCM k = me->get_grob_property ("knee");
   if (gh_boolean_p (k))
@@ -1295,36 +1466,59 @@ Beam::knee_b (Grob*me)
   bool knee = false;
   int d = 0;
   for (SCM s = me->get_grob_property ("stems"); gh_pair_p (s); s = ly_cdr (s))
-    if (d != Directional_element_interface::get (unsmob_grob (ly_car (s))))
-      {
-       knee = true;
-       break;
-      }
+    {
+      Direction dir = Directional_element_interface::get
+       (unsmob_grob (ly_car (s)));
+      if (d && d != dir)
+       {
+         knee = true;
+         break;
+       }
+      d = dir;
+    }
   
   me->set_grob_property ("knee", gh_bool2scm (knee));
 
   return knee;
 }
 
+int
+Beam::get_direction_beam_count (Grob *me, Direction d )
+{
+  Link_array<Grob>stems = 
+    Pointer_group_interface__extract_grobs (me, (Grob*) 0, "stems");
+  int bc = 0;
+  
+  for (int i = stems.size (); i--;)
+    {
+      /*
+       Should we take invisible stems into account?
+       */
+      if (Stem::get_direction (stems[i]) == d)
+        bc = bc >? (Stem::beam_multiplicity (stems[i]).length () + 1);
+    }
 
-ADD_INTERFACE (Beam, "beam-interface",
-  "A beam.
-
-#'thickness= weight of beams, in staffspace
-
-
-We take the least squares line through the ideal-length stems, and
-then damp that using
-
-       damped = tanh (slope)
+  return bc;
+}
 
-this gives an unquantized left and right position for the beam end.
-Then we take all combinations of quantings near these left and right
-positions, and give them a score (according to how close they are to
-the ideal slope, how close the result is to the ideal stems, etc.). We
-take the best scoring combination.
 
-",
-  "french-beaming position-callbacks concaveness-gap concaveness-threshold dir-function quant-score auto-knee-gap gap chord-tremolo beamed-stem-shorten shorten least-squares-dy damping flag-width-function neutral-direction positions space-function thickness");
+ADD_INTERFACE (Beam, "beam-interface",
+  "A beam. \n\n"
+" "
+"#'thickness= weight of beams, in staffspace "
+" "
+" "
+"We take the least squares line through the ideal-length stems, and "
+"then damp that using "
+" \n"
+"      damped = tanh (slope) \n"
+" \n"
+"this gives an unquantized left and right position for the beam end. "
+"Then we take all combinations of quantings near these left and right "
+"positions, and give them a score (according to how close they are to "
+"the ideal slope, how close the result is to the ideal stems, etc.). We "
+"take the best scoring combination. "
+,
+  "knee french-beaming position-callbacks concaveness-gap concaveness-threshold dir-function quant-score auto-knee-gap gap chord-tremolo beamed-stem-shorten shorten least-squares-dy damping flag-width-function neutral-direction positions space-function thickness");