]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/stem.cc
release: 1.3.83
[lilypond.git] / lily / stem.cc
index 60f7b1bad336edfdce8a3d64213fa2667e925665..b64f89bafacbf81b40311789b1855ab428700ef2 100644 (file)
 
   source file of the GNU LilyPond music typesetter
 
-  (c) 1996,1997 Han-Wen Nienhuys <hanwen@stack.nl>
+  (c) 1996--2000 Han-Wen Nienhuys <hanwen@cs.uu.nl>
+    Jan Nieuwenhuizen <janneke@gnu.org>
 
   TODO: This is way too hairy
 */
+#include <math.h>              // m_pi
 
+#include "directional-element-interface.hh"
+#include "note-head.hh"
 #include "stem.hh"
-#include "dimen.hh" 
 #include "debug.hh"
 #include "paper-def.hh"
-#include "note-head.hh"
+#include "rhythmic-head.hh"
 #include "lookup.hh"
 #include "molecule.hh"
-#include "p-col.hh"
+#include "paper-column.hh"
 #include "misc.hh"
 #include "beam.hh"
 #include "rest.hh"
+#include "group-interface.hh"
+#include "cross-staff.hh"
+#include "staff-symbol-referencer.hh"
+#include "spanner.hh"
 
-const int STEMLEN=7;
-
-IMPLEMENT_IS_TYPE_B1 (Stem,Item);
 
-Stem::Stem () 
+void
+Stem::set_beaming (Score_element*me ,int i,  Direction d )
 {
-  /*
-    TODO: staff-size
-   */
-  abbrev_flag_i_ = 0;
-  beam_l_ = 0;
-  beams_left_i_ = 0;
-  beams_right_i_ = 0;
-
-  stem_bottom_f_ = stem_top_f_ = 0;
-  flag_i_ = 2;
-  dir_ = CENTER;
-  staff_size_i_ = 8;
-
-  stem_xoffset_f_ =0;
-  beam_gap_i_ = 0;
+  SCM pair = me->get_elt_property ("beaming");
+  
+  if (!gh_pair_p (pair))
+    {
+      pair = gh_cons (gh_int2scm (0),gh_int2scm (0));
+      me->      set_elt_property ("beaming", pair);
+    }
+  index_set_cell (pair, d, gh_int2scm (i));
 }
 
 int
-Stem::min_head_i () const
+Stem::beam_count (Score_element*me,Direction d)
 {
-  int m = 1000;
-  for (int i =0; i < head_l_arr_.size (); i++)
-    m = m <? head_l_arr_[i]->position_i_;
-  return m;
+  SCM p=me->get_elt_property ("beaming");
+  if (gh_pair_p (p))
+    return gh_scm2int (index_cell (p,d));
+  else
+    return 0;
 }
 
-int
-Stem::max_head_i () const
+Interval
+Stem::head_positions (Score_element*me) 
 {
-  int m = -1000;
-  for (int i =0; i < head_l_arr_.size (); i++)
-    m = m >? head_l_arr_[i]->position_i_;
-  return m;
-  
-}
+  if (!heads_i (me))
+    {
+      Interval iv;
+      return iv;
+    }
 
-void
-Stem::do_print () const
-{
-#ifndef NPRINT
-  DOUT << "flag "<< flag_i_ << "abbrev_flag_i_" << abbrev_flag_i_;
-  if (beam_l_)
-    DOUT << "beamed";
-#endif
+  Drul_array<Score_element*> e (extremal_heads (me));
+
+  return Interval (Staff_symbol_referencer::position_f (e[DOWN]),
+                  Staff_symbol_referencer::position_f ( e[UP]));
 }
 
-Real 
-Stem::stem_length_f () const
+
+Real
+Stem::chord_start_f (Score_element*me) 
 {
-  return stem_top_f_-stem_bottom_f_ ;
+  return head_positions(me)[get_direction (me)]
+    * Staff_symbol_referencer::staff_space (me)/2.0;
 }
 
 Real
-Stem::stem_start_f () const
+Stem::stem_end_position (Score_element*me) 
 {
-  return (dir_ < 0)? stem_top_f_ : stem_bottom_f_;
+  SCM p =me->get_elt_property ("stem-end-position");
+  Real pos;
+  if (!gh_number_p (p))
+    {
+
+      pos = get_default_stem_end_position (me);
+      me->set_elt_property ("stem-end-position", gh_double2scm (pos));
+    }
+  else
+    pos = gh_scm2double (p);
+
+  return pos;
 }
 
-Real
-Stem::stem_end_f () const
+Direction
+Stem::get_direction (Score_element*me)
 {
-  return (dir_ < 0)? stem_bottom_f_ : stem_top_f_;
+  Direction d = Directional_element_interface::get (me);
+
+  if (!d)
+    {
+       d = get_default_dir (me);
+       // urg, AAARGH!
+       Directional_element_interface::set (me, d);
+    }
+  return d ;
 }
 
 
 void
-Stem::set_stemend (Real se)
+Stem::set_stemend (Score_element*me, Real se)
 {
   // todo: margins
-  if (!  ((dir_ > 0 && se >= max_head_i ()) || 
-         (se <= min_head_i () && dir_ <0)))    
-    warning ("Weird stem size; check for narrow beams");
+  Direction d= get_direction (me);
+  
+  if (d && d * head_positions(me)[get_direction (me)] >= se*d)
+    warning (_ ("Weird stem size; check for narrow beams"));
 
-  stem_top_f_  = (dir_ < 0) ? max_head_i () : se;
-  stem_bottom_f_  = (dir_ < 0) ? se  : min_head_i ();
+  me->set_elt_property ("stem-end-position", gh_double2scm (se));
 }
 
 int
-Stem::type_i () const
+Stem::type_i (Score_element*me) 
 {
-  return head_l_arr_[0]->balltype_i_;
+  return first_head (me) ?  Rhythmic_head::balltype_i (first_head (me)) : 2;
 }
 
-void
-Stem::add (Note_head *n)
+/*
+  Note head that determines hshift for upstems
+ */ 
+Score_element*
+Stem::support_head (Score_element*me)
+{
+  SCM h = me->get_elt_property ("support-head");
+  Score_element * nh = unsmob_element (h);
+  if (nh)
+    return nh;
+  else if (heads_i (me) == 1)
+    {
+      /*
+       UGH.
+       */
+      
+      return unsmob_element (gh_car (me->get_elt_property ("heads")));
+    }
+  else
+    return first_head (me);
+}
+
+
+int
+Stem::heads_i (Score_element*me)
 {
-  n->add_dependency (this);
-  head_l_arr_.push (n);
+  Pointer_group_interface gi (me, "heads");
+  return gi.count ();
 }
-void
-Stem::add (Rest *r)
+
+/*
+  The note head which forms one end of the stem.  
+ */
+Score_element*
+Stem::first_head (Score_element*me)
 {
-  rest_l_arr_.push (r);
-  r->add_dependency (this);    // ?
+  return extremal_heads (me)[-get_direction (me)];
 }
 
-bool
-Stem::invisible_b () const
+/*
+  START is part where stem reaches `last' head. 
+ */
+Drul_array<Score_element*>
+Stem::extremal_heads (Score_element*me) 
 {
+  const int inf = 1000000;
+  Drul_array<int> extpos;
+  extpos[DOWN] = inf;
+  extpos[UP] = -inf;  
   
-  return (!head_l_arr_.size () || 
-    head_l_arr_[0]->balltype_i_ <= 0);
+  Drul_array<Score_element *> exthead;
+  exthead[LEFT] = exthead[RIGHT] =0;
+  
+  for (SCM s = me->get_elt_property ("heads"); gh_pair_p (s); s = gh_cdr (s))
+    {
+      Score_element * n = unsmob_element (gh_car (s));
+
+      
+      int p = int(Staff_symbol_referencer::position_f (n));
+
+      Direction d = LEFT;
+      do {
+      if (d* p > d* extpos[d])
+       {
+         exthead[d] = n;
+         extpos[d] = p;
+       }
+      } while (flip (&d) != DOWN);
+    }
 
+  return exthead;
 }
 
-// if dir_ is set we return fake values.
-int
-Stem::get_center_distance_from_top ()
+void
+Stem::add_head (Score_element*me, Score_element *n)
 {
-  if (dir_)
-    return (dir_ > 0) ? 0 : 1;
+  n->set_elt_property ("stem", me->self_scm ());
+  n->add_dependency (me);
+
+  if (Note_head::has_interface (n))
+    {
+      Pointer_group_interface (me, "heads").add_element (n);
+    }
+  else
+    {
+      n->set_elt_property ("rest", n->self_scm ());
+    }
+}
 
-  int staff_center = staff_size_i_ / 2;
-  int max = max_head_i () - staff_center;
-  return max >? 0;
+bool
+Stem::invisible_b (Score_element*me)
+{
+  return !(heads_i (me) && Rhythmic_head::balltype_i (support_head (me)) >= 1);
 }
 
-// if dir_ is set we return fake values.
 int
-Stem::get_center_distance_from_bottom ()
+Stem::get_center_distance (Score_element*me, Direction d)
 {
-  if (dir_)
-    return (dir_ > 0) ? 1 : 0;
-
-  int staff_center = staff_size_i_ / 2;
-  int min = staff_center - min_head_i ();
-  return min >? 0;
+  int staff_center = 0;
+  int distance = (int) (d*(head_positions(me)[d] - staff_center));
+  return distance >? 0;
 }
 
 Direction
-Stem::get_default_dir ()
+Stem::get_default_dir (Score_element*me) 
 {
-  if (dir_)
-    return dir_;
-  return (get_center_distance_from_top () >=
-         get_center_distance_from_bottom ()) ? 
-    (Direction)-1 : (Direction)1;
-}
+  int du = get_center_distance (me,UP);
+  int dd = get_center_distance (me,DOWN);
 
+  if (sign (dd - du))
+    return Direction (sign (dd -du));
 
-void
-Stem::set_default_dir ()
-{
-  dir_ = get_default_dir ();
+  return to_dir (me->get_elt_property ("default-neutral-direction"));
 }
 
-void
-Stem::set_default_stemlen ()
+/*
+  ugh. A is used for different purposes. This functionality should be
+  moved into scheme at some point to get rid of the silly
+  conversions. (but lets wait till we have namespaces in SCM)
+ */
+Real
+Stem::get_default_stem_end_position (Score_element*me) 
 {
-  if (!dir_)
-    set_default_dir ();
-
-  // ugh... how about non 5-line staffs?
-  bool on_ledger_line_b = ((max_head_i () < -2 && dir_ == 1)
-                          //    || (min_head_i () > staff_size_i_ && dir_ == -1));
-                          || (min_head_i () > staff_size_i_ + 3 && dir_ == -1));
-  if (on_ledger_line_b)
+  bool grace_b = to_boolean (me->get_elt_property ("grace"));
+  String type_str = grace_b ? "grace-" : "";
+  SCM s;
+  Array<Real> a;
+
+  Real length_f = 0.;
+  SCM scm_len = me->get_elt_property("length");
+  if (gh_number_p (scm_len))
     {
-      set_stemend (staff_size_i_ / 2 - 1);
+      length_f = gh_scm2double (scm_len);
     }
-  else 
+  else
     {
-      Real dy = paper ()->interbeam_f ();
-      Real len = STEMLEN;
-      // ugh, should get nice *rule* for this
-      if (abbrev_flag_i_ > 1)
-       len += (abbrev_flag_i_ - 1)* dy / 2;
-      set_stemend ((dir_ > 0) ? max_head_i () + len :
-                  min_head_i () - len);
+      s = scm_eval2 (ly_symbol2scm ((type_str + "stem-length").ch_C()),
+                    SCM_EOL);
+      for (SCM q = s; q != SCM_EOL; q = gh_cdr (q))
+       a.push (gh_scm2double (gh_car (q)));
+               
+      // stem uses half-spaces
+      length_f = a[((flag_i (me) - 2) >? 0) <? (a.size () - 1)] * 2;
     }
-}
 
-void
-Stem::set_default_extents ()
-{
-  if (!stem_length_f ())
-    set_default_stemlen ();
-
-  set_stemend ((dir_< 0) ? 
-              max_head_i ()-stem_length_f (): min_head_i () + stem_length_f ());
-  // ugh, a whole ball is wider
-  if (head_l_arr_[0]->balltype_i_ <= 0)
-    stem_xoffset_f_ = paper ()->note_width () / 2;
-  else if (dir_ > 0)   
-    stem_xoffset_f_ = paper ()->note_width () - paper ()->rule_thickness ();
-  else
-    stem_xoffset_f_ = 0;
+
+  a.clear ();
+  s = scm_eval2 (ly_symbol2scm ((type_str + "stem-shorten").ch_C()),
+                SCM_EOL);
+  for (SCM q = s; q != SCM_EOL; q = gh_cdr (q))
+    a.push (gh_scm2double (gh_car (q)));
+
+
+  // stem uses half-spaces
+
+  // fixme: use gh_list_ref () iso. array[]
+  Real shorten_f = a[((flag_i (me) - 2) >? 0) <? (a.size () - 1)] * 2;
+
+  /* URGURGURG
+     'set-default-stemlen' sets direction too
+   */
+  Direction dir = get_direction (me);
+  if (!dir)
+    {
+      dir = get_default_dir (me);
+      Directional_element_interface::set (me, dir);
+    }
+  
+  /* 
+    stems in unnatural (forced) direction should be shortened, 
+    according to [Roush & Gourlay]
+   */
+  if (((int)chord_start_f (me))
+      && (get_direction (me) != get_default_dir (me)))
+    length_f -= shorten_f;
+
+
+   Real st = head_positions(me)[dir] + dir * length_f;
+  
+   bool no_extend_b = to_boolean (me->get_elt_property ("no-stem-extend"));
+    if (!grace_b && !no_extend_b && dir * st < 0)
+      st = 0.0;
+
+  return st;
 }
 
 /*
-  TODO
-  
-  move into note_column.cc
+  FIXME: wrong name
+ */
+int
+Stem::flag_i (Score_element*me) 
+{
+  SCM s = me->get_elt_property ("duration-log");
+  return  (gh_number_p (s)) ? gh_scm2int (s) : 2;
+}
 
-  */
 void
-Stem::set_noteheads ()
+Stem::position_noteheads (Score_element*me)
 {
-  if (!head_l_arr_.size ())
+  if (!heads_i (me))
     return;
-  head_l_arr_.sort (Note_head::compare);
-  if (dir_ < 0) 
-    head_l_arr_.reverse ();
   
-  head_l_arr_[0]->extremal_i_ = -1;
-  head_l_arr_.top ()->extremal_i_ = 1;
-  int parity=1;
-  int lastpos = head_l_arr_[0]->position_i_;
-  for (int i=1; i < head_l_arr_.size (); i ++) 
+  Link_array<Score_element> heads =
+    Pointer_group_interface__extract_elements (me, (Score_element*)0, "heads");
+
+  heads.sort (compare_position);
+  Direction dir =get_direction (me);
+  
+  if (dir < 0)
+    heads.reverse ();
+
+
+  Real w = support_head (me)->extent (X_AXIS)[dir];
+  for (int i=0; i < heads.size (); i++)
     {
-      int dy =abs (lastpos- head_l_arr_[i]->position_i_);
-       
-      if (dy <= 1) 
+      heads[i]->translate_axis (w - heads[i]->extent (X_AXIS)[dir], X_AXIS);
+    }
+  
+  bool parity= true;           // todo: make me settable.
+  int lastpos = int (Staff_symbol_referencer::position_f (heads[0]));
+  for (int i=1; i < heads.size (); i ++)
+    {
+      Real p = Staff_symbol_referencer::position_f (heads[i]);
+      int dy =abs (lastpos- (int)p);
+
+      if (dy <= 1)
        {
          if (parity)
-           head_l_arr_[i]->x_dir_ = (stem_xoffset_f_>0) ? UP:DOWN;
+           {
+             Real l  = heads[i]->extent (X_AXIS).length ();
+             heads[i]->translate_axis (l * get_direction (me), X_AXIS);
+           }
          parity = !parity;
        }
       else
-       parity = 0;
-      lastpos = head_l_arr_[i]->position_i_;
+       parity = true;
+      
+      lastpos = int (p);
+    }
+}
+
+MAKE_SCHEME_CALLBACK(Stem,before_line_breaking);
+SCM
+Stem::before_line_breaking (SCM smob)
+{
+  Score_element*me = unsmob_element (smob);
+  stem_end_position (me);      // ugh. Trigger direction calc.
+  position_noteheads (me);
+
+  if (invisible_b (me))
+    {
+      me->remove_elt_property ("molecule-callback");
+      // suicide();
     }
+  
+  set_spacing_hints (me);
+  return SCM_UNSPECIFIED;
 }
 
+
+
+/**
+   set stem directions for hinting the optical spacing correction.
+
+   Modifies DIR_LIST property of the Stem's Paper_column
+
+   TODO: more advanced: supply height of noteheads as well, for more advanced spacing possibilities
+ */
 void
-Stem::do_pre_processing ()
+Stem::set_spacing_hints (Score_element*me) 
 {
-  if (stem_bottom_f_== stem_top_f_)
-    set_default_extents ();
-  set_noteheads ();
-  flag_i_ = flag_i_;
-  transparent_b_ = invisible_b ();
-  set_empty (invisible_b ());
+  if (!invisible_b (me))
+    {
+      SCM scmdir  = gh_int2scm (get_direction (me));
+
+      Item* item = dynamic_cast<Item*> (me);
+      Item * col =  item->column_l ();
+      SCM dirlist =col->get_elt_property ("dir-list");
+      if (dirlist == SCM_UNDEFINED)
+       dirlist = SCM_EOL;
+
+      if (scm_sloppy_memq (scmdir, dirlist) == SCM_EOL)
+       {
+         dirlist = gh_cons (scmdir, dirlist);
+         col->set_elt_property ("dir-list", dirlist);
+       }
+    }
 }
 
+Molecule
+Stem::flag (Score_element*me)
+{
+  String style;
+  SCM st = me->get_elt_property ("flag-style");
+  if ( gh_string_p (st))
+    {
+      style = ly_scm2string (st);
+    }
+
+  char c = (get_direction (me) == UP) ? 'u' : 'd';
+  Molecule m = me->lookup_l ()->afm_find (String ("flags-") + to_str (c) + 
+                                     to_str (flag_i (me)));
+  if (!style.empty_b ())
+    m.add_molecule(me->lookup_l ()->afm_find (String ("flags-") + to_str (c) + style));
+  return m;
+}
 
 Interval
-Stem::do_width () const
+Stem::dim_callback (Score_element *se, Axis ) 
 {
   Interval r (0, 0);
-  if (abbrev_flag_i_)
-    {
-      r = abbrev_mol ().extent ().x ();
-    }
-  else if (beam_l_ || abs (flag_i_) <= 2)
+  if (unsmob_element (se->get_elt_property ("beam")) || abs (flag_i (se)) <= 2)
     ;  // TODO!
   else
     {
-      Paper_def*p= paper ();
-      r = p->lookup_l ()->flag (flag_i_, dir_).dim.x ();
-      r+= stem_xoffset_f_;
+      r = flag (se).extent (X_AXIS);
     }
   return r;
 }
 
 
-  
-Molecule
- Stem::abbrev_mol () const
+const Real ANGLE = 20* (2.0*M_PI/360.0); // ugh! Should be settable.
+
+
+MAKE_SCHEME_CALLBACK(Stem,brew_molecule);
+
+SCM
+Stem::brew_molecule (SCM smob) 
 {
-  Real dy = paper ()->interbeam_f ();
-  Real w = 1.5 * paper ()->lookup_l ()->ball (2).dim.x ().length ();
-  Real beamdy = paper ()->interline_f () / 2;
-  int beams_i = 0;
-  Real slope = paper ()->internote_f () / 4;
-  if (beam_l_) {
-    // huh?
-      slope = 2 * beam_l_->slope;
-    // ugh, rather calc from Abbreviation_req
-      beams_i = beams_right_i_ >? beams_left_i_; 
-  }
-  paper ()->lookup_l ()->beam (slope, 20 PT);
+  Score_element*me = unsmob_element (smob);
+  Molecule mol;
+  Direction d = get_direction (me);
+  
+  
+  Real y1 = Staff_symbol_referencer::position_f (first_head (me));
+  Real y2 = stem_end_position (me);
   
-  Molecule beams;
-  Atom a (paper ()->lookup_l ()->beam (slope, w));
-  a.translate (Offset(- w / 2, stem_end_f () - (w / 2 * slope)));
-  // ugh
-    if (!beams_i)
-      a.translate (dy + beamdy - dir_ * dy, Y_AXIS);
-    else
-      a.translate (2 * beamdy - dir_ * (beamdy - dy), Y_AXIS);
+  Interval stem_y(y1,y2);
+  stem_y.unite (Interval (y2,y1));
+
+  Real dy = Staff_symbol_referencer::staff_space (me)/2.0;
+  Real head_wid = 0;
+  if (support_head (me))
+    head_wid = support_head (me)->extent (X_AXIS).length ();
+  stem_y[Direction(-d)] += d * head_wid * tan(ANGLE)/(2*dy);
   
-  for (int i = 0; i < abbrev_flag_i_; i++) 
+  if (!invisible_b (me))
     {
-      Atom b (a);
-      b.translate (-dir_ * dy * (beams_i + i), Y_AXIS);
-      beams.add (b);
+      Real stem_width = me->paper_l ()->get_var ("stemthickness");
+      Molecule ss =me->lookup_l ()->filledbox (Box (Interval (-stem_width/2, stem_width/2),
+                                                Interval (stem_y[DOWN]*dy, stem_y[UP]*dy)));
+      mol.add_molecule (ss);
     }
-  
-  return beams;
+
+  if (!beam_l (me) && abs (flag_i (me)) > 2)
+    {
+      Molecule fl = flag (me);
+      fl.translate_axis(stem_y[d]*dy, Y_AXIS);
+      mol.add_molecule (fl);
+    }
+
+  return mol.create_scheme();
 }
 
-Molecule*
- Stem::brew_molecule_p () const 
+Real
+Stem::off_callback (Score_element * me, Axis)
 {
-  Molecule *mol_p =new Molecule;
-  
-  Real bot  = stem_bottom_f_;
-  Real top = stem_top_f_;
-  
-  assert (bot!=top);
+  Real r=0;
+  if (Score_element * f = first_head (me))
+    {
+      Interval head_wid(0, f->extent (X_AXIS).length ());
+
+      if (to_boolean (me->get_elt_property ("stem-centered")))
+       return head_wid.center ();
+      
+      Real rule_thick = me->paper_l ()->get_var ("stemthickness");
+      Direction d = get_direction (me);
+      r = head_wid[d] - d * rule_thick ;
+    }
+  return r;
+}
+
+
+
+Score_element*
+Stem::beam_l (Score_element*me)
+{
+  SCM b=  me->get_elt_property ("beam");
+  return unsmob_element (b);
+}
+
+
+// ugh still very long.
+Stem_info
+Stem::calc_stem_info (Score_element*me) 
+{
+  Score_element * beam = beam_l (me);
+
+  Direction beam_dir = Directional_element_interface::get (beam);
+  if (!beam_dir)
+    {
+      programming_error ("Beam dir not set.");
+      beam_dir = UP;
+    }
+    
+
+  Real staff_space = Staff_symbol_referencer::staff_space (me);
+  Real half_space = staff_space / 2;
+  Real interbeam_f = me->paper_l ()->interbeam_f (Beam::get_multiplicity (beam));
+  Real thick = gh_scm2double (beam->get_elt_property ("beam-thickness"));
+  int multiplicity = Beam::get_multiplicity (beam);
+
+  Stem_info info; 
+  info.idealy_f_ = chord_start_f (me);
+
+  // for simplicity, we calculate as if dir == UP
+  info.idealy_f_ *= beam_dir;
+  SCM grace_prop = me->get_elt_property ("grace");
+
+  bool grace_b = to_boolean (grace_prop);
   
-  Paper_def *p =paper ();
+  Array<Real> a;
+  SCM s;
+  String type_str = grace_b ? "grace-" : "";
   
-  Real dy = p->internote_f ();
-  if (!invisible_b ())
+  s = scm_eval2 (ly_symbol2scm ((type_str + "beamed-stem-minimum-length").ch_C()), SCM_EOL);
+  a.clear ();
+  for (SCM q = s; q != SCM_EOL; q = gh_cdr (q))
+    a.push (gh_scm2double (gh_car (q)));
+
+
+  Real minimum_length = a[multiplicity <? (a.size () - 1)] * staff_space;
+  s = scm_eval2 (ly_symbol2scm ((type_str + "beamed-stem-length").ch_C()),
+                SCM_EOL);
+
+  a.clear();
+  for (SCM q = s; q != SCM_EOL; q = gh_cdr (q))
+    a.push (gh_scm2double (gh_car (q)));
+
+  Real stem_length =  a[multiplicity <? (a.size () - 1)] * staff_space;
+
+  if (!beam_dir || (beam_dir == Directional_element_interface::get (me)))
+    /* normal beamed stem */
     {
-      Symbol ss =p->lookup_l ()->stem (bot*dy,top*dy);
-      mol_p->add (Atom (ss));
+      if (multiplicity)
+       {
+         info.idealy_f_ += thick + (multiplicity - 1) * interbeam_f;
+       }
+      info.miny_f_ = info.idealy_f_;
+      info.maxy_f_ = INT_MAX;
+
+      info.idealy_f_ += stem_length;
+      info.miny_f_ += minimum_length;
+
+      /*
+       lowest beam of (UP) beam must never be lower than second staffline
+
+       Hmm, reference (Wanske?)
+
+       Although this (additional) rule is probably correct,
+       I expect that highest beam (UP) should also never be lower
+       than middle staffline, just as normal stems.
+       
+      */
+      bool no_extend_b = to_boolean (me->get_elt_property ("no-stem-extend"));
+      if (!grace_b && !no_extend_b)
+       {
+         /* highest beam of (UP) beam must never be lower than middle
+            staffline
+            lowest beam of (UP) beam must never be lower than second staffline
+          */
+         info.miny_f_ =
+           info.miny_f_ >? 0
+           >? (- 2 * half_space - thick
+               + (multiplicity > 0) * thick
+               + interbeam_f * (multiplicity - 1));
+       }
     }
-  
-  if (!beam_l_ &&abs (flag_i_) > 2)
+  else
+    /* knee */
     {
-      Symbol fl = p->lookup_l ()->flag (flag_i_, dir_);
-      mol_p->add_at_edge (Y_AXIS, dir_, Molecule (Atom (fl)));
-      assert (!abbrev_flag_i_);
+      info.idealy_f_ -= thick;
+      info.maxy_f_ = info.idealy_f_;
+      info.miny_f_ = -INT_MAX;
+
+      info.idealy_f_ -= stem_length;
+      info.maxy_f_ -= minimum_length;
     }
   
-  if (abbrev_flag_i_)
-    mol_p->add (abbrev_mol ());
-  
-  mol_p->translate (stem_xoffset_f_, X_AXIS);
-  return mol_p;
+  info.idealy_f_ = (info.maxy_f_ <? info.idealy_f_) >? info.miny_f_;
+
+  s = beam->get_elt_property ("shorten");
+  if (gh_number_p (s))
+    info.idealy_f_ -= gh_scm2double (s);
+
+  Real interstaff_f = -beam_dir* calc_interstaff_dist (dynamic_cast<Item*> (me), dynamic_cast<Spanner*> (beam));
+
+  info.idealy_f_ += interstaff_f;
+  info.miny_f_ += interstaff_f;
+  info.maxy_f_ += interstaff_f ;
+
+  return info;
 }
 
-Real
- Stem::hpos_f () const
+bool
+Stem::has_interface (Score_element*m)
 {
-  return Item::hpos_f () + stem_xoffset_f_;
+  return m && m->has_interface (ly_symbol2scm ("stem-interface"));
 }
 
-
 void
- Stem::do_substitute_dependency (Score_elem*o,Score_elem*n)
-{
-  Item * o_l = o->item ();
-  Item * n_l = n? n->item () : 0;
-  head_l_arr_.substitute ((Note_head*)o_l, (Note_head*)n_l);
-  rest_l_arr_.substitute ((Rest*)o_l, (Rest*)n_l);
+Stem::set_interface (Score_element*me)
+{    
+  me->set_elt_property ("heads", SCM_EOL);
+  me->add_offset_callback ( &Stem::off_callback, X_AXIS);
+  me->set_interface (ly_symbol2scm ("stem-interface"));
 }