]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/beam.cc
test
[lilypond.git] / lily / beam.cc
index 7c40e0758f3b6a21343ab3197acffea0a5b79758..c2d024df6adba40cc5b8df8ef5126f07a7c4d469 100644 (file)
@@ -24,14 +24,13 @@ Notes:
 
  - Stems run to the Y-center of the beam.
   
- - beam_space is the offset between Y centers of the beam.
+ - beam_translation is the offset between Y centers of the beam.
 
 */
 
 
 #include <math.h> // tanh.
 
-#include "align-interface.hh"
 #include "molecule.hh" 
 #include "directional-element-interface.hh"
 #include "beaming.hh"
@@ -57,23 +56,6 @@ Notes:
 #endif
 
 
-const int INTER_QUANT_PENALTY = 1000; 
-const int SECONDARY_BEAM_DEMERIT  = 15;
-const int STEM_LENGTH_DEMERIT_FACTOR = 5;
-// possibly ridiculous, but too short stems just won't do
-const int STEM_LENGTH_LIMIT_PENALTY = 5000;
-const int DAMPING_DIRECTIION_PENALTY = 800;
-const int MUSICAL_DIRECTION_FACTOR = 400;
-const int IDEAL_SLOPE_FACTOR = 10;
-const int REGION_SIZE = 2;
-
-
-static Real
-shrink_extra_weight (Real x)
-{
-  return fabs (x) * ((x < 0) ? 1.5 : 1.0);
-}
-
 void
 Beam::add_stem (Grob *me, Grob *s)
 {
@@ -81,14 +63,18 @@ 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_beam_space (Grob *me)
+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)));
@@ -118,15 +104,15 @@ 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 line = me->get_paper ()->get_var ("linethickness");
   Real thickness = gh_scm2double (me->get_grob_property ("thickness"))
     * staff_space;
   
-  Real beam_space = gh_scm2int (beam_count) < 4
+  Real beam_translation = gh_scm2int (beam_count) < 4
     ? (2*staff_space + line - thickness) / 2.0
     : (3*staff_space + line - thickness) / 3.0;
   
-  return gh_double2scm (beam_space);
+  return gh_double2scm (beam_translation);
 }
 
 
@@ -231,6 +217,7 @@ Beam::connect_beams (Grob *me)
        }
       else
        {
+         gh_set_car_x ( this_beaming, SCM_EOL);
          SCM s = gh_cdr (this_beaming);
          for (; gh_pair_p (s); s = gh_cdr (s))
            {
@@ -239,8 +226,13 @@ Beam::connect_beams (Grob *me)
              last_int.add_point (np);
            }
        }
+
+      if (i == stems.size () -1)
+       {
+         gh_set_cdr_x ( this_beaming, SCM_EOL);
+       }
     }
-}
+ }
 
 MAKE_SCHEME_CALLBACK (Beam, brew_molecule, 1);
 SCM
@@ -278,7 +270,7 @@ Beam::brew_molecule (SCM grob)
   Real dydx = dy && dx ? dy/dx : 0;
   
   Real thick = gh_scm2double (me->get_grob_property ("thickness"));
-  Real bdy = get_beam_space (me);
+  Real bdy = get_beam_translation (me);
 
   SCM last_beaming = SCM_EOL;;
   Real last_xposn = -1;
@@ -287,7 +279,7 @@ 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];
@@ -416,16 +408,16 @@ 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);
+      Molecule tm = Text_item::text2molecule (me, ly_str02scm (str.to_str0 ()), properties);
       the_beam.add_at_edge (Y_AXIS, UP, tm, 5.0);
     }
 #endif
@@ -492,76 +484,95 @@ Beam::set_stem_directions (Grob *me, Direction d)
   for (int i=0; i <stems.size (); i++)
     {
       Grob *s = stems[i];
-      SCM force = s->remove_grob_property ("dir-forced");
-      if (!gh_boolean_p (force) || !gh_scm2bool (force))
-       Directional_element_interface::set (s, d);
+      /* For knees, non-forced stems should probably have their
+        natural direction. In any case, when knee, beam direction is
+        foe.
+        
+        TODO: for x staff knees, set direction pointing to 'the
+        other' staff, rather than natural.
+      */
+      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);
+       }
     }
 } 
 
 /* Simplistic auto-knees; only consider vertical gap between two
    adjacent chords.
 
-  `Forced' stem directions are ignored.  If you don't want auto-knees,
-  don't set, or unset auto-knee-gap. */
+   This may decide for a knee that's impossible to fit sane scoring
+   criteria (eg, stem lengths).  We may need something smarter. */
 void
 Beam::consider_auto_knees (Grob *me, Direction d)
 {
   SCM scm = me->get_grob_property ("auto-knee-gap");
 
-  if (gh_number_p (scm))
-    {
-      bool knee_b = false;
-      Real knee_y = 0;
-      Real staff_space = Staff_symbol_referencer::staff_space (me);
-      Real gap = gh_scm2double (scm) / staff_space;
-
+  if (!gh_number_p (scm))
+    return;
+  
+  bool knee_b = false;
+  
+  Real staff_space = Staff_symbol_referencer::staff_space (me);
+  Real gap = gh_scm2double (scm) / staff_space;
 
-      Link_array<Grob> stems=
-       Pointer_group_interface__extract_grobs (me, (Grob*)0, "stems");
+  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 i=1; i < stems.size (); i++)
-        {
-         if (!Stem::invisible_b (stems[i-1]))
-           l = i - 1;
-         if (Stem::invisible_b (stems[l]))
-           continue;
-         if (Stem::invisible_b (stems[i]))
-           continue;
+  Grob *common = common_refpoint_of_array (stems, me,  Y_AXIS);
+
+  int l = 0;
+  for (int r=1; r < stems.size (); r++)
+    {
+      if (!Stem::invisible_b (stems[r-1]))
+       l = r - 1;
+      Grob *right = stems[r];
+      Grob *left = stems[l];
+      if (Stem::invisible_b (left))
+       continue;
+      if (Stem::invisible_b (right))
+       continue;
          
-         Real left = Stem::extremal_heads (stems[l])[d]
-           ->relative_coordinate (common, Y_AXIS);
-         Real right = Stem::extremal_heads (stems[i])[-d]
-           ->relative_coordinate (common, Y_AXIS);
+      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 - left;
+      Real dy = right_y - left_y;
 
-         if (abs (dy) >= gap)
-           {
-             knee_y = (right + left) / 2;
-             knee_b = true;
-             break;
-           }
-       }
-      
-      if (knee_b)
+      if (abs (dy) >= gap)
        {
-         for (int i=0; i < stems.size (); i++)
+         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)
            {
-             Grob *s = stems[i];         
-             if (Stem::invisible_b (s) || 
-                 s->get_grob_property ("dir-forced") == SCM_BOOL_T)
-               continue;
-             Real y = Stem::extremal_heads (stems[i])[d]
-               ->relative_coordinate (common, Y_AXIS);
+             Directional_element_interface::set (left, knee_dir);
+             left->set_grob_property ("dir-forced", SCM_BOOL_T);
 
-             Directional_element_interface::set (s, y < knee_y ? UP : DOWN);
-             s->set_grob_property ("dir-forced", SCM_BOOL_T);
+           }
+         if (!Stem::invisible_b (right)
+             && stems[r]->get_grob_property ("dir-forced") != SCM_BOOL_T)
+           {
+             Directional_element_interface::set (right, -knee_dir);
+             right->set_grob_property ("dir-forced", SCM_BOOL_T);
            }
        }
     }
+
+  if (knee_b)
+    {
+      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);
+    }
 }
 
 /* Set stem's shorten property if unset.
@@ -569,17 +580,20 @@ Beam::consider_auto_knees (Grob *me, Direction d)
  TODO:
    take some y-position (chord/beam/nearest?) into account
    scmify forced-fraction
-
-   TODO:
-   
-   why is shorten stored in beam, and not directly in stem?
+  This is done in beam because the shorten has to be uniform over the
+  entire beam.
 
 */
 void
-Beam::set_stem_shorten (Grob *m)
+Beam::set_stem_shorten (Grob *me)
 {
-  Spanner*me = dynamic_cast<Spanner*> (m);
-
+  /*
+    shortening looks silly for x staff beams
+   */
+  if (knee_b(me))
+    return ;
+  
   Real forced_fraction = forced_stem_count (me) / visible_stem_count (me);
 
   int beam_count = get_beam_count (me);
@@ -629,378 +643,6 @@ Beam::after_line_breaking (SCM smob)
   return SCM_UNSPECIFIED;
 }
 
-struct Quant_score
-{
-  Real yl;
-  Real yr;
-  Real demerits;
-};
-
-
-/*
-  TODO:
-  
-   - Make all demerits customisable
-
-   - One sensible check per demerit (what's this --hwn)
-
-   - Add demerits for quants per se, as to forbid a specific quant
-     entirely
-
-*/
-MAKE_SCHEME_CALLBACK (Beam, quanting, 1);
-SCM
-Beam::quanting (SCM smob)
-{
-  Grob *me = unsmob_grob (smob);
-
-  SCM s = me->get_grob_property ("positions");
-  Real yl = gh_scm2double (gh_car (s));
-  Real yr = gh_scm2double (gh_cdr (s));
-
-  Real ss = Staff_symbol_referencer::staff_space (me);
-  Real thickness = gh_scm2double (me->get_grob_property ("thickness")) / ss;
-  Real slt = me->paper_l ()->get_var ("linethickness") / ss;
-
-
-  SCM sdy = me->get_grob_property ("least-squares-dy");
-  Real dy_mus = gh_number_p (sdy) ? gh_scm2double (sdy) : 0.0;
-  
-  Real straddle = 0.0;
-  Real sit = (thickness - slt) / 2;
-  Real inter = 0.5;
-  Real hang = 1.0 - (thickness - slt) / 2;
-  Real quants [] = {straddle, sit, inter, hang };
-  
-  int num_quants = int (sizeof (quants)/sizeof (Real));
-  Array<Real> quantsl;
-  Array<Real> quantsr;
-
-  /*
-    going to REGION_SIZE == 2, yields another 0.6 second with
-    wtk1-fugue2.
-
-
-    (result indexes between 70 and 575)  ? --hwn. 
-
-  */
-
-
-  
-  /*
-    Do stem computations.  These depend on YL and YR linearly, so we can
-    precompute for every stem 2 factors.
-   */
-  Link_array<Grob> stems=
-    Pointer_group_interface__extract_grobs (me, (Grob*)0, "stems");
-  Array<Stem_info> stem_infos;
-  Array<Real> lbase_lengths;
-  Array<Real> rbase_lengths;  
-
-  Drul_array<bool> dirs_found(0,0);
-  Grob *common_y = common_refpoint_of_array (stems, me, Y_AXIS);
-  
-  bool french = to_boolean (me->get_grob_property ("french-beaming"));
-  for (int i= 0; i < stems.size(); i++)
-    {
-      Grob*s = stems[i];
-      stem_infos.push (Stem::calc_stem_info (s));
-      dirs_found[stem_infos.top ().dir_] = true;
-
-      Real b = calc_stem_y (me, s, common_y , Interval (1,0), french && i > 0&& (i < stems.size  () -1));
-      lbase_lengths.push (b);
-
-      Real a = calc_stem_y (me, s, common_y , Interval (0,1),  french && i > 0&& (i < stems.size  () -1));
-      rbase_lengths.push (a);
-    }
-
-  Direction ldir = Direction (stem_infos[0].dir_);
-  Direction rdir = Direction (stem_infos.top ().dir_);
-  bool knee_b = dirs_found[LEFT] && dirs_found[RIGHT];
-
-
-  int region_size = REGION_SIZE;
-  /*
-    Knees are harder, lets try some more possibilities for knees. 
-   */
-  if (knee_b)
-    region_size += 2;
-  
-  for (int i = -region_size ; i < region_size; i++)
-    for (int j = 0; j < num_quants; j++)
-      {
-       quantsl.push (i + quants[j] + int (yl));
-       quantsr.push (i + quants[j] + int (yr));
-      }
-
-  Array<Quant_score> qscores;
-  
-  for (int l =0; l < quantsl.size (); l++)  
-    for (int r =0; r < quantsr.size (); r++)
-      {
-       Quant_score qs;
-       qs.yl = quantsl[l];
-       qs.yr = quantsr[r];
-       qs.demerits = 0.0;
-       
-       qscores.push (qs);
-      }
-
-
-  /*
-    This is a longish function, but we don't separate this out into
-    neat modular separate subfunctions, as the subfunctions would be
-    called for many values of YL, YR. By precomputing various
-    parameters outside of the loop, we can save a lot of time.
-
-  */
-
-
-  Grob *fvs  = first_visible_stem (me);
-  Grob *lvs  = last_visible_stem (me);
-
-  bool xstaff= false;
-  if (lvs && fvs)
-    {
-      Grob *commony = fvs->common_refpoint (lvs, Y_AXIS);
-      xstaff = Align_interface::has_interface (commony);
-    }
-
-  for (int i = qscores.size (); i--;)
-    if (qscores[i].demerits < 100)
-      {
-       qscores[i].demerits
-         += score_slopes_dy (me, qscores[i].yl, qscores[i].yr,
-                             dy_mus, yr- yl, xstaff); 
-      }
-
-  Real rad = Staff_symbol_referencer::staff_radius (me);
-  int beam_count = get_beam_count (me);
-  Real beam_space = beam_count < 4
-    ? (2*ss + slt - thickness) / 2.0
-     : (3*ss + slt - thickness) / 3.0;
-
-  for (int i = qscores.size (); i--;)
-    if (qscores[i].demerits < 100)
-      {
-       qscores[i].demerits
-         += score_forbidden_quants (me, qscores[i].yl, qscores[i].yr,
-                                    rad, slt, thickness, beam_space,
-                                    beam_count, ldir, rdir); 
-      }
-
-
-  for (int i = qscores.size (); i--;)
-    if (qscores[i].demerits < 100)
-      {
-       qscores[i].demerits
-         += score_stem_lengths (stems, stem_infos,
-                                lbase_lengths, rbase_lengths,
-                                knee_b,
-                                me, qscores[i].yl, qscores[i].yr);
-      }
-
-
-  Real best = 1e6;
-  int best_idx = -1;
-  for (int i = qscores.size (); i--;)
-    {
-      if (qscores[i].demerits < best)
-       {
-         best = qscores [i].demerits ;
-         best_idx = i;
-       }
-    }
-
-  
-  me->set_grob_property ("positions",
-                        gh_cons (gh_double2scm (qscores[best_idx].yl),
-                                 gh_double2scm (qscores[best_idx].yr))
-                        );
-
-#if DEBUG_QUANTING
-
-  // debug quanting
-  me->set_grob_property ("quant-score",
-                        gh_double2scm (qscores[best_idx].demerits));
-  me->set_grob_property ("best-idx", gh_int2scm (best_idx));
-#endif
-
-  return SCM_UNSPECIFIED;
-}
-
-Real
-Beam::score_stem_lengths (Link_array<Grob>stems,
-                         Array<Stem_info> stem_infos,
-                         Array<Real> left_factor,
-                         Array<Real> right_factor,
-                         bool knee, 
-                         Grob*me,
-                         Real yl, Real yr)
-{
-  Real demerit_score = 0.0 ;
-  Real pen = STEM_LENGTH_LIMIT_PENALTY;
-  
-  for (int i=0; i < stems.size (); i++)
-    {
-      Grob* s = stems[i];
-      if (Stem::invisible_b (s))
-       continue;
-
-      Real current_y =
-       yl * left_factor[i] + right_factor[i]* yr;
-
-      Stem_info info = stem_infos[i];
-      Direction d = info.dir_;
-
-      demerit_score += pen
-       * ( 0 >? (info.dir_ * (info.shortest_y_ - current_y)));
-      
-      demerit_score += STEM_LENGTH_DEMERIT_FACTOR
-       * shrink_extra_weight (d * current_y  - info.dir_ * info.ideal_y_);
-    }
-
-  demerit_score *= 2.0 / stems.size (); 
-
-  return demerit_score;
-}
-
-Real
-Beam::score_slopes_dy (Grob *me,
-                      Real yl, Real yr,
-                      Real dy_mus, Real dy_damp,
-                      bool xstaff)
-{
-  Real dy = yr - yl;
-
-  Real dem = 0.0;
-  if (sign (dy_damp) != sign (dy))
-    {
-      dem += DAMPING_DIRECTIION_PENALTY;
-    }
-
-   dem += MUSICAL_DIRECTION_FACTOR * (0 >? (fabs (dy) - fabs (dy_mus)));
-
-
-   Real slope_penalty = IDEAL_SLOPE_FACTOR;
-
-   /*
-     Xstaff beams tend to use extreme slopes to get short stems. We
-     put in a penalty here.
-   */
-   if (xstaff)
-     slope_penalty *= 10;
-
-   dem += shrink_extra_weight (fabs (dy_damp) - fabs (dy))* slope_penalty;
-   return dem;
-}
-
-static Real
-my_modf (Real x)
-{
-  return x - floor (x);
-}
-
-Real
-Beam::score_forbidden_quants (Grob*me,
-                             Real yl, Real yr,
-                             Real rad,
-                             Real slt,
-                             Real thickness, Real beam_space,
-                             int beam_count,
-                             Direction ldir, Direction rdir)
-{
-  Real dy = yr - yl;
-
-  Real dem = 0.0;
-  if (fabs (yl) < rad && fabs ( my_modf (yl) - 0.5) < 1e-3)
-    dem += INTER_QUANT_PENALTY;
-  if (fabs (yr) < rad && fabs ( my_modf (yr) - 0.5) < 1e-3)
-    dem += INTER_QUANT_PENALTY;
-
-  // todo: use beam_count of outer stems.
-  if (beam_count >= 2)
-    {
-     
-      Real straddle = 0.0;
-      Real sit = (thickness - slt) / 2;
-      Real inter = 0.5;
-      Real hang = 1.0 - (thickness - slt) / 2;
-      
-
-      if (fabs (yl - ldir * beam_space) < rad
-         && fabs (my_modf (yl) - inter) < 1e-3)
-       dem += SECONDARY_BEAM_DEMERIT;
-      if (fabs (yr - rdir * beam_space) < rad
-         && fabs (my_modf (yr) - inter) < 1e-3)
-       dem += SECONDARY_BEAM_DEMERIT;
-
-      Real eps = 1e-3;
-
-      /*
-       Can't we simply compute the distance between the nearest
-       staffline and the secondary beam? That would get rid of the
-       silly case analysis here (which is probably not when we have
-       different beam-thicknesses.)
-
-       --hwn
-       */
-
-
-      // hmm, without Interval/Drul_array, you get ~ 4x same code...
-      if (fabs (yl - ldir * beam_space) < rad + inter)
-       {
-         if (ldir == UP && dy <= eps
-             && fabs (my_modf (yl) - sit) < eps)
-           dem += SECONDARY_BEAM_DEMERIT;
-         
-         if (ldir == DOWN && dy >= eps
-             && fabs (my_modf (yl) - hang) < eps)
-           dem += SECONDARY_BEAM_DEMERIT;
-       }
-
-      if (fabs (yr - rdir * beam_space) < rad + inter)
-       {
-         if (rdir == UP && dy >= eps
-             && fabs (my_modf (yr) - sit) < eps)
-           dem += SECONDARY_BEAM_DEMERIT;
-         
-         if (rdir == DOWN && dy <= eps
-             && fabs (my_modf (yr) - hang) < eps)
-           dem += SECONDARY_BEAM_DEMERIT;
-       }
-      
-      if (beam_count >= 3)
-       {
-         if (fabs (yl - 2 * ldir * beam_space) < rad + inter)
-           {
-             if (ldir == UP && dy <= eps
-                 && fabs (my_modf (yl) - straddle) < eps)
-               dem += SECONDARY_BEAM_DEMERIT;
-             
-             if (ldir == DOWN && dy >= eps
-                 && fabs (my_modf (yl) - straddle) < eps)
-               dem += SECONDARY_BEAM_DEMERIT;
-       }
-         
-         if (fabs (yr - 2 * rdir * beam_space) < rad + inter)
-           {
-             if (rdir == UP && dy >= eps
-                 && fabs (my_modf (yr) - straddle) < eps)
-               dem += SECONDARY_BEAM_DEMERIT;
-             
-             if (rdir == DOWN && dy <= eps
-                 && fabs (my_modf (yr) - straddle) < eps)
-               dem += SECONDARY_BEAM_DEMERIT;
-           }
-       }
-    }
-  
-  return dem;
-}
-
-  
-
 MAKE_SCHEME_CALLBACK (Beam, least_squares, 1);
 SCM
 Beam::least_squares (SCM smob)
@@ -1380,27 +1022,23 @@ where_are_the_whole_beams(SCM beaming)
 
 /*
   Calculate the Y position of the stem-end, given the Y-left, Y-right
-  in POS, and for stem S.
+  in POS for stem S. This Y position is relative to S.
  */
 Real
-Beam::calc_stem_y (Grob *me, Grob* s, Grob * common_y, Interval pos, bool french) 
+Beam::calc_stem_y (Grob *me, Grob* s, Grob ** common,
+                  Real xl, Real xr,
+                  Interval pos, bool french) 
 {
-  Real beam_space = get_beam_space (me);
+  Real beam_translation = get_beam_translation (me);
 
-  // ugh -> use commonx
-  Grob * fvs = first_visible_stem (me);
-  Grob *lvs = last_visible_stem (me);
     
-  Real x0 = fvs ? fvs->relative_coordinate (0, X_AXIS) : 0.0;
-  Real dx = fvs ? lvs->relative_coordinate (0, X_AXIS) - x0 : 0.0;
-  Real r = s->relative_coordinate (0, X_AXIS) - x0;
+  Real r = s->relative_coordinate (common[X_AXIS], X_AXIS) - xl;
   Real dy = pos.delta ();
+  Real dx = xr - xl;
   Real stem_y_beam0 = (dy && dx
                       ? r / dx
                       * dy
                       : 0) + pos[LEFT];
-
-
   
   Direction my_dir = Directional_element_interface::get (s);
   SCM beaming = s->get_grob_property ("beaming");
@@ -1410,17 +1048,17 @@ Beam::calc_stem_y (Grob *me, Grob* s, Grob * common_y, Interval pos, bool french
     {
       Slice bm = where_are_the_whole_beams (beaming);
       if (!bm.empty_b())
-       stem_y += beam_space * bm[-my_dir];
+       stem_y += beam_translation * bm[-my_dir];
     }
   else
     {
       Slice bm = Stem::beam_multiplicity(s);
       if (!bm.empty_b())
-       stem_y +=bm[my_dir] * beam_space;
+       stem_y +=bm[my_dir] * beam_translation;
     }
   
-  Real id = me->relative_coordinate (common_y, Y_AXIS)
-    - s->relative_coordinate (common_y, Y_AXIS);
+  Real id = me->relative_coordinate (common[Y_AXIS], Y_AXIS)
+    - s->relative_coordinate (common[Y_AXIS], Y_AXIS);
   
   return stem_y + id;
 }
@@ -1438,21 +1076,51 @@ Beam::set_stem_lengths (Grob *me)
   if (stems.size () <= 1)
     return;
   
-  Grob *common = common_refpoint_of_array (stems, me, Y_AXIS);
+  Grob *common[2];
+  for (int a = 2; a--;)
+    common[a] = common_refpoint_of_array (stems, me, Axis(a));
+  
   Interval pos = ly_scm2interval (me->get_grob_property ("positions"));
   Real staff_space = Staff_symbol_referencer::staff_space (me);
 
   bool french = to_boolean (me->get_grob_property ("french-beaming"));
+
+  
+  bool gap = false;
+  Real thick =0.0;
+  if (gh_number_p (me->get_grob_property ("gap"))
+      &&gh_scm2double (me->get_grob_property ("gap")))
+  {
+    gap = true;
+    thick = gh_scm2double (me->get_grob_property ("thickness"))
+      * Staff_symbol_referencer::staff_space(me);
+  }
+      
+  // ugh -> use commonx
+  Grob * fvs = first_visible_stem (me);
+  Grob *lvs = last_visible_stem (me);
+    
+  Real xl = fvs ? fvs->relative_coordinate (common[X_AXIS], X_AXIS) : 0.0;
+  Real xr = lvs ? lvs->relative_coordinate (common[X_AXIS], X_AXIS) : 0.0;
+  
   for (int i=0; i < stems.size (); i++)
     {
       Grob* s = stems[i];
       if (Stem::invisible_b (s))
        continue;
 
-      Real stem_y = calc_stem_y (me, s, common , pos, french && i > 0&& (i < stems.size  () -1));
+      Real stem_y = calc_stem_y (me, s, common,
+                                xl, xr,
+                                pos, french && i > 0&& (i < stems.size  () -1));
 
-        Stem::set_stemend (s, 2* stem_y / staff_space);
+      /*
+       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_y += thick * 0.5 * Directional_element_interface::get(s);
+      
+      Stem::set_stemend (s, 2* stem_y / staff_space);
     }
 }
 
@@ -1623,8 +1291,31 @@ Beam::rest_collision_callback (SCM element_smob, SCM axis)
   return gh_double2scm (-d *  discrete_dist);
 }
 
+bool
+Beam::knee_b (Grob*me)
+{
+  SCM k = me->get_grob_property ("knee");
+  if (gh_boolean_p (k))
+    return gh_scm2bool (k);
 
+  bool knee = false;
+  int d = 0;
+  for (SCM s = me->get_grob_property ("stems"); gh_pair_p (s); s = ly_cdr (s))
+    {
+      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;
+}
 
 ADD_INTERFACE (Beam, "beam-interface",
   "A beam.