]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/accidental-placement.cc
Run grand replace for 2015.
[lilypond.git] / lily / accidental-placement.cc
index c45947629e07484ac520ee7b6f21fb2c53b44c3f..c8c14d4bd980ce54cceba08d4c3f56eee003c937 100644 (file)
-/*   
-     accidental-placement.cc --  implement Accidental_placement
+/*
+  This file is part of LilyPond, the GNU music typesetter.
+
+  Copyright (C) 2002--2015 Han-Wen Nienhuys <hanwen@xs4all.nl>
 
-     source file of the GNU LilyPond music typesetter
+  LilyPond is free software: you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation, either version 3 of the License, or
+  (at your option) any later version.
 
-     (c) 2002--2004 Han-Wen Nienhuys <hanwen@cs.uu.nl>
+  LilyPond is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
 
+  You should have received a copy of the GNU General Public License
+  along with LilyPond.  If not, see <http://www.gnu.org/licenses/>.
 */
 
-#include <math.h>
+#include "accidental-placement.hh"
 
-#include "libc-extension.hh"   // isinf
+#include "accidental-interface.hh"
 #include "item.hh"
-#include "skyline.hh"
 #include "music.hh"
-#include "pitch.hh"
-#include "warn.hh"
-#include "accidental-placement.hh"
-#include "note-column.hh"
-#include "group-interface.hh"
 #include "note-collision.hh"
-#include "accidental-interface.hh"
+#include "note-column.hh"
+#include "pointer-group-interface.hh"
+#include "rhythmic-head.hh"
+#include "skyline.hh"
+#include "skyline-pair.hh"
+#include "stream-event.hh"
+#include "warn.hh"
 
-MAKE_SCHEME_CALLBACK (Accidental_placement,alignment_callback, 2);
-SCM
-Accidental_placement::alignment_callback (SCM s, SCM )
+static Pitch *
+accidental_pitch (Grob *acc)
 {
-  Grob * me =unsmob_grob (s);
+  SCM cause = acc->get_parent (Y_AXIS)->get_property ("cause");
 
-  Grob * par = me->get_parent (X_AXIS);
-  if (!to_boolean (par->get_property ("positioning-done")))
+  Stream_event *mcause = Stream_event::unsmob (cause);
+  if (!mcause)
     {
-      par->set_property ("positioning-done", SCM_BOOL_T);
-      position_accidentals (par);
+      programming_error ("note head has no event cause");
+      return 0;
     }
 
-  return scm_int2num (0);
+  return Pitch::unsmob (mcause->get_property ("pitch"));
 }
 
-
 void
-Accidental_placement::add_accidental (Grob* me, Grob* a)
+Accidental_placement::add_accidental (Grob *me, Grob *a, bool stagger, long context_hash)
 {
-  a->set_parent (me, X_AXIS);
-  a->add_offset_callback (alignment_callback_proc, X_AXIS);
-  SCM cause = a->get_parent (Y_AXIS)->get_property ("cause");
-
-  Music *mcause =unsmob_music (cause); 
-  if (!mcause)
-    {
-      programming_error ("Note head has no music cause!");
-      return; 
-    }
+  Pitch *p = accidental_pitch (a);
+  if (!p)
+    return;
 
-  Pitch *p= unsmob_pitch (mcause->get_property ("pitch"));
-
-  int n = p->get_notename ();
+  a->set_parent (me, X_AXIS);
+  a->set_property ("X-offset", Grob::x_parent_positioning_proc);
+  long n = p->get_notename ();
 
-  SCM accs = me->get_property ("accidental-grobs");
-  SCM key = scm_int2num (n);
-  SCM entry = scm_assq (key, accs);
+  SCM accs = me->get_object ("accidental-grobs");
+  SCM key = scm_cons (scm_from_int (n), scm_from_long  (stagger ? context_hash : 1));
+  // assoc because we're dealing with pairs
+  SCM entry = scm_assoc (key, accs);
   if (entry == SCM_BOOL_F)
-    {
-      entry = SCM_EOL;
-    }
+    entry = SCM_EOL;
   else
-    entry = ly_cdr (entry);
+    entry = scm_cdr (entry);
 
   entry = scm_cons (a->self_scm (), entry);
 
-  accs = scm_assq_set_x (accs,  key, entry);
+  accs = scm_assoc_set_x (accs, key, entry);
 
-  me->set_property ("accidental-grobs", accs);
+  me->set_object ("accidental-grobs", accs);
 }
 
 /*
   Split into break reminders.
- */
+*/
 void
-Accidental_placement::split_accidentals (Grob * accs,
-                                        Link_array<Grob> *break_reminder,
-                                        Link_array<Grob> *real_acc)
+Accidental_placement::split_accidentals (Grob *accs,
+                                         vector<Grob *> *break_reminder,
+                                         vector<Grob *> *real_acc)
 {
-  for (SCM acs =accs->get_property ("accidental-grobs"); is_pair (acs);
-       acs =ly_cdr (acs))
-    for (SCM s = ly_cdar (acs); is_pair (s); s = ly_cdr (s))
+  for (SCM acs = accs->get_object ("accidental-grobs"); scm_is_pair (acs);
+       acs = scm_cdr (acs))
+    for (SCM s = scm_cdar (acs); scm_is_pair (s); s = scm_cdr (s))
       {
-       Grob *a = unsmob_grob (ly_car (s));
+        Grob *a = Grob::unsmob (scm_car (s));
 
-       if (unsmob_grob (a->get_property ("tie")))
-         break_reminder->push (a);
-       else
-         real_acc->push (a);
+        if (Grob::is_smob (a->get_object ("tie"))
+            && !to_boolean (a->get_property ("forced")))
+          break_reminder->push_back (a);
+        else
+          real_acc->push_back (a);
       }
 }
 
-/*
-  Accidentals are special, because they appear and disappear after
-  ties at will.
-*/
-Interval
-Accidental_placement::get_relevant_accidental_extent (Grob *me,
-                                                     Item *item_col,
-                                                     Grob *left_object)
+vector<Grob *>
+Accidental_placement::get_relevant_accidentals (vector<Grob *> const &elts, Grob *left)
 {
-  Link_array<Grob> br, ra;
-  Link_array<Grob> *which = 0;
+  vector<Grob *> br;
+  vector<Grob *> ra;
+  vector<Grob *> ret;
+  bool right = dynamic_cast<Item *> (left)->break_status_dir () == RIGHT;
 
-  Accidental_placement::split_accidentals (me, &br, &ra);
-  br.concat (ra);
-  
-  if (dynamic_cast<Item*>(left_object)->break_status_dir () == RIGHT)
-    which = & br;
-  else
-    which = & ra;
-  
-  Interval extent;
-  for (int i = 0; i < which->size (); i++)
+  for (vsize i = 0; i < elts.size (); i++)
     {
-      extent.unite (which->elem (i)->extent (item_col, X_AXIS));
-    }
+      split_accidentals (elts[i], &br, &ra);
 
-  if (!extent.is_empty ())
-    {
-      Real p = robust_scm2double (me->get_property ("left-padding"), 0.2);
-      extent[LEFT] -= p;
+      ret.insert (ret.end (), ra.begin (), ra.end ());
+
+      if (right)
+        ret.insert (ret.end (), br.begin (), br.end ());
     }
-  
-  return extent;
+  return ret;
 }
 
-
-
 struct Accidental_placement_entry
 {
-  Array<Skyline_entry> left_skyline_;
-  Array<Skyline_entry> right_skyline_;
-  Interval vertical_extent_;
-  Array<Box> extents_;
-  Link_array<Grob> grobs_;
-  Real offset_; 
-  int notename_;
-  Accidental_placement_entry ()
-  {
-    offset_ =0.0;
-    notename_ = -1;
-  }
+  Skyline_pair horizontal_skylines_;
+  vector<Grob *> grobs_;
 };
 
-static Interval all_accidental_vertical_extent;
-Real ape_priority (Accidental_placement_entry const * a)
+Real ape_priority (Accidental_placement_entry const *a)
 {
-  return a->vertical_extent_[UP];
+  // right is up because we're horizontal
+  return a->horizontal_skylines_.right ();
 }
 
-  
-
-int ape_compare (Accidental_placement_entry *const &a,
-                Accidental_placement_entry *const &b)
+bool ape_less (Accidental_placement_entry *const &a,
+               Accidental_placement_entry *const &b)
 {
-  return sign (ape_priority (a) - ape_priority (b));
-}
+  vsize size_a = a->grobs_.size ();
+  vsize size_b = b->grobs_.size ();
+  if (size_a != size_b)
+    return size_b < size_a;
 
-int ape_rcompare (Accidental_placement_entry *const &a,
-                Accidental_placement_entry *const &b)
-{
-  return -sign (ape_priority (a) - ape_priority (b));
+  return ape_priority (a) < ape_priority (b);
 }
 
-
 /*
+  This function provides a method for sorting accidentals that belong to the
+  same note. The accidentals that this function considers to be "smallest"
+  will be placed to the left of the "larger" accidentals.
 
-TODO: should favor
-
-  b
- b
+  Naturals are the largest (so that they don't get confused with cancellation
+  naturals); apart from that, we order according to the alteration (so
+  double-flats are the smallest).
 
-placement
+  Precondition: the accidentals are attached to NoteHeads of the same note
+  name -- the octaves, however, may be different.
 */
-void
-stagger_apes (Link_array<Accidental_placement_entry> *apes)
+static bool
+acc_less (Grob *const &a, Grob *const &b)
 {
-  Link_array<Accidental_placement_entry> asc = *apes;
-
+  Pitch *p = accidental_pitch (a);
+  Pitch *q = accidental_pitch (b);
 
-  asc.sort (&ape_compare);
-
-  apes->clear ();
-
-  int i =0;
-  int parity = 1;
-  while (i < asc.size ())
+  if (!p || !q)
     {
-      Accidental_placement_entry * a = 0;      
-      if (parity)
-       a = asc.pop ();
-      else
-       a = asc[i++];
-
-      apes->push (a);
-      parity = !parity;
+      programming_error ("these accidentals do not have a pitch");
+      return false;
     }
 
-  apes->reverse ();
-}
+  if (p->get_octave () != q->get_octave ())
+    return p->get_octave () < q->get_octave ();
 
-  
+  if (p->get_alteration () == Rational (0))
+    return false;
+  if (q->get_alteration () == Rational (0))
+    return true;
 
-/*
-
-  This routine computes placements of accidentals. During
-  add_accidental (), accidentals are already grouped by note, so that
-  octaves are placed above each other; they form columns. Then the
-  columns are sorted: the biggest columns go closest to the note.
-  Then the columns are spaced as closely as possible (using skyline
-  spacing).
-  
-  
-  TODO: more advanced placement. Typically, the accs should be placed
-  to form a C shape, like this
-
-  
-           ##
-        b b
-       # #
-        b
-          b b
+  return p->get_alteration () < q->get_alteration ();
+}
 
-   The naturals should be left of the C as well; they should
-   be separate accs.
+/*
+  TODO: should favor
 
-   Note that this placement problem looks NP hard, so we just use a
-   simple strategy, not an optimal choice.
+  *  b
+  * b
 
+  placement
 */
+void
+stagger_apes (vector<Accidental_placement_entry *> *apes)
+{
+  vector<Accidental_placement_entry *> asc = *apes;
+
+  vector_sort (asc, &ape_less);
+  // we do the staggering below based on size
+  // this ensures that if a placement has 4 entries, it will
+  // always be closer to the NoteColumn than a placement with 1
+  // this allows accidentals to be on-average closer to notes
+  // while still preserving octave alignment
+  vector<vector<Accidental_placement_entry *> > ascs;
+
+  vsize sz = INT_MAX;
+  for (vsize i = 0; i < asc.size (); i++)
+    {
+      vsize my_sz = asc[i]->grobs_.size ();
+      if (sz != my_sz)
+        ascs.push_back (vector<Accidental_placement_entry *> ());
+      ascs.back ().push_back (asc[i]);
+      sz = my_sz;
+    }
 
-/*
-  TODO: there should be more space in the following situation
+  apes->clear ();
 
+  for (vsize i = 0; i < ascs.size (); i++)
+    {
+      int parity = 1;
+      for (vsize j = 0; j < ascs[i].size ();)
+        {
+          Accidental_placement_entry *a = 0;
+          if (parity)
+            {
+              a = ascs[i].back ();
+              ascs[i].pop_back ();
+            }
+          else
+            a = ascs[i][j++];
+
+          apes->push_back (a);
+          parity = !parity;
+        }
+    }
 
-    Natural + downstem
+  reverse (*apes);
+}
 
-    |_       
-    | |    X
-    |_|   | 
-      |   |
-      
- */
-SCM
-Accidental_placement::position_accidentals (Grob * me)
+static vector<Accidental_placement_entry *>
+build_apes (SCM accs)
 {
-  if (!me->live ())
-    return SCM_UNSPECIFIED;
-  
-  SCM accs = me->get_property ("accidental-grobs");
-
-  /*
-    TODO: there is a bug in this code. If two accs are on the same
-    Y-position, they share an Ape, and will be printed in overstrike.
-   */
-  Link_array<Accidental_placement_entry> apes;
-  for (SCM s = accs; is_pair (s); s =ly_cdr (s))
+  vector<Accidental_placement_entry *> apes;
+  for (SCM s = accs; scm_is_pair (s); s = scm_cdr (s))
     {
       Accidental_placement_entry *ape = new Accidental_placement_entry;
-      ape->notename_ = ly_scm2int (ly_caar (s));
-      
-      for (SCM t = ly_cdar (s); is_pair (t); t =ly_cdr (t))
-       ape->grobs_.push (unsmob_grob (ly_car (t)));
 
-      apes.push (ape);
+      for (SCM t = scm_cdar (s); scm_is_pair (t); t = scm_cdr (t))
+        ape->grobs_.push_back (Grob::unsmob (scm_car (t)));
+
+      apes.push_back (ape);
     }
 
+  return apes;
+}
 
-  Grob *common[] = {me, 0};
+static void
+set_ape_skylines (Accidental_placement_entry *ape,
+                  Grob **common, Real padding)
+{
+  vector<Grob *> accs (ape->grobs_);
+  vector_sort (accs, &acc_less);
+
+  /* We know that each accidental has the same note name and we assume that
+     accidentals in different octaves won't collide. If two or more
+     accidentals are in the same octave:
+     1) if they are the same accidental, print them in overstrike
+     2) otherwise, shift one to the left so they don't overlap. */
+  int last_octave = 0;
+  Real offset = 0;
+  Real last_offset = 0;
+  Rational last_alteration (0);
+  for (vsize i = accs.size (); i--;)
+    {
+      Grob *a = accs[i];
+      Pitch *p = accidental_pitch (a);
+
+      if (!p)
+        continue;
+
+      if (i == accs.size () - 1 || p->get_octave () != last_octave)
+        {
+          last_offset = 0;
+          offset = a->extent (a, X_AXIS)[LEFT] - padding;
+        }
+      else if (p->get_alteration () == last_alteration)
+        a->translate_axis (last_offset, X_AXIS);
+      else /* Our alteration is different from the last one */
+        {
+          Real this_offset = offset - a->extent (a, X_AXIS)[RIGHT];
+          a->translate_axis (this_offset, X_AXIS);
+
+          last_offset = this_offset;
+          offset -= a->extent (a, X_AXIS).length () + padding;
+        }
+
+      if (Skyline_pair *sky = Skyline_pair::unsmob (a->get_property ("horizontal-skylines")))
+        {
+          Skyline_pair copy (*sky);
+          copy.raise (a->relative_coordinate (common[X_AXIS], X_AXIS));
+          copy.shift (a->relative_coordinate (common[Y_AXIS], Y_AXIS));
+          ape->horizontal_skylines_.merge (copy);
+        }
+
+      last_octave = p->get_octave ();
+      last_alteration = p->get_alteration ();
+    }
+}
 
-  /*
-    First we must extract *all* pointers. We can only determine
-    extents if we're sure that we've found the right common refpoint
-   */
-  Link_array<Grob> note_cols, heads;
-  for (int i= apes.size (); i--;)
-    { 
-      Accidental_placement_entry * ape = apes[i];
-      for (int j = ape->grobs_.size (); j--;)
-       {
-         Grob * a = ape->grobs_[j];
-
-         if (common[Y_AXIS])
-           common[Y_AXIS] = common[Y_AXIS]->common_refpoint (a, Y_AXIS);
-         else
-           common[Y_AXIS] = a;
-         
-         Grob *head = a->get_parent (Y_AXIS);
-
-         Grob * col = head->get_parent (X_AXIS);
-         if (Note_column::has_interface (col))
-           note_cols.push (col);
-         else
-           heads.push (head);
-       }
+static vector<Grob *>
+extract_heads_and_stems (vector<Accidental_placement_entry *> const &apes)
+{
+  vector<Grob *> note_cols;
+  vector<Grob *> ret;
+
+  for (vsize i = apes.size (); i--;)
+    {
+      Accidental_placement_entry *ape = apes[i];
+      for (vsize j = ape->grobs_.size (); j--;)
+        {
+          Grob *acc = ape->grobs_[j];
+          Grob *head = acc->get_parent (Y_AXIS);
+          Grob *col = head->get_parent (X_AXIS);
+
+          if (Note_column::has_interface (col))
+            note_cols.push_back (col);
+          else
+            ret.push_back (head);
+        }
     }
 
   /*
-    This is a little kludgy: to get all notes, we look if there are
-    collisions as well.
-   */
-  for (int i = note_cols.size () ; i--;)
+    This is a little kludgy: in case there are note columns without
+    accidentals, we get them from the Note_collision objects.
+  */
+  for (vsize i = note_cols.size (); i--;)
     {
       Grob *c = note_cols[i]->get_parent (X_AXIS);
       if (Note_collision_interface::has_interface (c))
-       {
-         Link_array<Grob> gs =
-           Pointer_group_interface__extract_grobs (c, (Grob*)0, "elements");
-      
-         note_cols.concat (gs);
-       }
-    }
-  
-  for (int i = note_cols.size () ; i--;)
-    {
-      heads.concat (Pointer_group_interface__extract_grobs (note_cols[i],
-                                                           (Grob*)0,
-                                                           "note-heads"));
-      
+        {
+          extract_grob_set (c, "elements", columns);
+          concat (note_cols, columns);
+        }
     }
-  heads.default_sort ();
-  heads.uniq ();
-  common[Y_AXIS] = common_refpoint_of_array (heads, common[Y_AXIS], Y_AXIS);
 
+  /* Now that we have all of the columns, grab all of the note-heads */
+  for (vsize i = note_cols.size (); i--;)
+    concat (ret, extract_grob_array (note_cols[i], "note-heads"));
+
+  /* Now that we have all of the heads, grab all of the stems */
+  for (vsize i = ret.size (); i--;)
+    if (Grob *s = Rhythmic_head::get_stem (ret[i]))
+      ret.push_back (s);
   
-  for (int i= apes.size (); i--;)
-    {
-      Accidental_placement_entry * ape = apes[i];
-      ape->left_skyline_ = empty_skyline (LEFT);
-      ape->right_skyline_ = empty_skyline (RIGHT);
-   
-      for (int j = apes[i]->grobs_.size (); j--;)
-       {
-         Grob * a = apes[i]->grobs_[j];
-
-         Array<Box> boxes = Accidental_interface::accurate_boxes (a, common);
-         
-         ape->extents_.concat (boxes);
-         for (int j  = boxes.size (); j--;)
-           {
-             insert_extent_into_skyline (&ape->left_skyline_, boxes[j], Y_AXIS, LEFT);
-             insert_extent_into_skyline (&ape->right_skyline_ , boxes[j], Y_AXIS, RIGHT);
-           }
-       }
-    }
+  uniquify (ret);
+  return ret;
+}
 
+static Grob *
+common_refpoint_of_accidentals (vector<Accidental_placement_entry *> const &apes, Axis a)
+{
+  Grob *ret = 0;
 
-  Interval total;
-  for (int i = apes.size (); i--;)
-    {
-      Interval y ;
-      
-      for (int j = apes[i]->extents_.size (); j--;)
-       {
-         y.unite (apes[i]->extents_[j][Y_AXIS]);
-       }
-      apes[i]->vertical_extent_ = y;
-      total.unite (y);
-    }
-  all_accidental_vertical_extent = total;
-  stagger_apes (&apes);
+  for (vsize i = apes.size (); i--;)
+    for (vsize j = apes[i]->grobs_.size (); j--;)
+      {
+        if (!ret)
+          ret = apes[i]->grobs_[j];
+        else
+          ret = ret->common_refpoint (apes[i]->grobs_[j], a);
+      }
 
-  Accidental_placement_entry * head_ape = new Accidental_placement_entry;
-  common[X_AXIS] = common_refpoint_of_array (heads, common[X_AXIS], X_AXIS);  
-  Array<Skyline_entry> head_skyline (empty_skyline (LEFT));
-  Array<Box> head_extents;
-  for (int i = heads.size (); i--;)
-    {
-      Box b (heads[i]->extent (common[X_AXIS] , X_AXIS),
-           heads[i]->extent (common[Y_AXIS], Y_AXIS));
+  return ret;
+}
 
-      insert_extent_into_skyline (&head_skyline, b , Y_AXIS, LEFT);
-    }
+static Skyline
+build_heads_skyline (vector<Grob *> const &heads_and_stems,
+                     Grob **common)
+{
+  vector<Box> head_extents;
+  for (vsize i = heads_and_stems.size (); i--;)
+    head_extents.push_back (Box (heads_and_stems[i]->extent (common[X_AXIS], X_AXIS),
+                                 heads_and_stems[i]->pure_height (common[Y_AXIS], 0, INT_MAX)));
 
-  head_ape-> left_skyline_ = head_skyline;
-  head_ape->offset_ = 0.0;
+  return Skyline (head_extents, Y_AXIS, LEFT);
+}
 
-  Real padding = robust_scm2double (me->get_property ("padding"),0.2);
+/*
+  Position the apes, starting from the right, so that they don't collide.
+  Return the total width.
+*/
+static Interval
+position_apes (Grob *me,
+               vector<Accidental_placement_entry *> const &apes,
+               Skyline const &heads_skyline)
+{
+  Real padding = robust_scm2double (me->get_property ("padding"), 0.2);
+  Skyline left_skyline = heads_skyline;
+  left_skyline.raise (-robust_scm2double (me->get_property ("right-padding"), 0));
 
-  Array<Skyline_entry> left_skyline = head_ape->left_skyline_;
-  heighten_skyline (&left_skyline,
-                   -robust_scm2double (me->get_property ("right-padding"), 0));
   /*
     Add accs entries right-to-left.
-   */
-  for (int i= apes.size (); i-- > 0;)
+  */
+  Interval width;
+  Real last_offset = 0.0;
+  for (vsize i = apes.size (); i-- > 0;)
     {
-      Real offset =
-       -skyline_meshing_distance (apes[i]->right_skyline_, left_skyline);
+      Accidental_placement_entry *ape = apes[i];
+
+      Real offset = -ape->horizontal_skylines_[RIGHT]
+                    .distance (left_skyline, 0.1);
       if (isinf (offset))
-       offset = (i < apes.size () - 1) ? apes[i+1]->offset_ : 0.0;
+        offset = last_offset;
       else
-       offset -= padding;
-
-      apes[i]->offset_ = offset;
+        offset -= padding;
 
-      Array<Skyline_entry> new_left_skyline = apes[i]->left_skyline_;
-      heighten_skyline (&new_left_skyline, apes[i]->offset_);
-      merge_skyline (&new_left_skyline, left_skyline, LEFT);
+      Skyline new_left_skyline = ape->horizontal_skylines_[LEFT];
+      new_left_skyline.raise (offset);
+      new_left_skyline.merge (left_skyline);
       left_skyline = new_left_skyline;
-    }      
 
-  for (int i = apes.size (); i--;)
-    {
-      Accidental_placement_entry* ape = apes[i];
-      for (int j  = ape->grobs_.size (); j--;)
-       {
-         ape->grobs_[j]->translate_axis (ape->offset_, X_AXIS);
-       }
+      /* Shift all of the accidentals in this ape */
+      for (vsize j = ape->grobs_.size (); j--;)
+        ape->grobs_[j]->translate_axis (offset, X_AXIS);
+
+      for (LEFT_and_RIGHT (d))
+        {
+          Real mh = ape->horizontal_skylines_[d].max_height ();
+          if (!isinf (mh))
+            width.add_point (mh + offset);
+        }
+
+      last_offset = offset;
     }
 
-  
-  Interval left_extent, right_extent;
-  Accidental_placement_entry *ape = apes[0];
+  return width;
+}
+
+/*
+  This routine computes placements of accidentals. During
+  add_accidental (), accidentals are already grouped by note, so that
+  octaves are placed above each other; they form columns. Then the
+  columns are sorted: the biggest columns go closest to the note.
+  Then the columns are spaced as closely as possible (using skyline
+  spacing).
 
-  for (int i = ape->extents_.size (); i--;)
-    left_extent.unite (ape->offset_ +  ape->extents_[i][X_AXIS]);
 
-  ape = apes.top ();
-  for (int i = ape->extents_.size (); i--;)
-    right_extent.unite (ape->offset_  + ape->extents_[i][X_AXIS]);
+  TODO: more advanced placement. Typically, the accs should be placed
+  to form a C shape, like this
 
-  
-  left_extent[LEFT] -= robust_scm2double (me->get_property ("left-padding"), 0);
+  *     ##
+  *  b b
+  * # #
+  *  b
+  *    b b
 
-  
-  Interval width (left_extent[LEFT], right_extent[RIGHT]);
+  The naturals should be left of the C as well; they should
+  be separate accs.
 
-  SCM scm_width = ly_interval2scm (width);
-  me->set_extent (scm_width, X_AXIS);
-  
-  for (int i = apes.size (); i--;)
-    delete apes[i];
+  Note that this placement problem looks NP hard, so we just use a
+  simple strategy, not an optimal choice.
+*/
+
+/*
+  TODO: there should be more space in the following situation
+
+
+  Natural + downstem
+
+  *
+  *  |_
+  *  | |    X
+  *  |_|   |
+  *    |   |
+  *
+*/
+
+MAKE_SCHEME_CALLBACK (Accidental_placement, calc_positioning_done, 1);
+SCM
+Accidental_placement::calc_positioning_done (SCM smob)
+{
+  Grob *me = Grob::unsmob (smob);
+  if (!me->is_live ())
+    return SCM_BOOL_T;
+
+  me->set_property ("positioning-done", SCM_BOOL_T);
+
+  SCM accs = me->get_object ("accidental-grobs");
+  if (!scm_is_pair (accs))
+    return SCM_BOOL_T;
+
+  vector<Accidental_placement_entry *> apes = build_apes (accs);
+
+  Grob *common[] = {me, 0};
+
+  vector<Grob *> heads_and_stems = extract_heads_and_stems (apes);
+
+  common[Y_AXIS] = common_refpoint_of_accidentals (apes, Y_AXIS);
+  common[Y_AXIS] = common_refpoint_of_array (heads_and_stems, common[Y_AXIS], Y_AXIS);
+  common[X_AXIS] = common_refpoint_of_array (heads_and_stems, me, X_AXIS);
+  Real padding = robust_scm2double (me->get_property ("padding"), 0.2);
+
+  for (vsize i = apes.size (); i--;)
+    set_ape_skylines (apes[i], common, padding);
+  Skyline heads_skyline = build_heads_skyline (heads_and_stems, common);
+
+  stagger_apes (&apes);
+  Interval width = position_apes (me, apes, heads_skyline);
+
+  me->flush_extent_cache (X_AXIS);
+  me->set_property ("X-extent", ly_interval2scm (width));
+
+  junk_pointers (apes);
 
-  return SCM_UNSPECIFIED;
+  return SCM_BOOL_T;
 }
 
 ADD_INTERFACE (Accidental_placement,
-             "accidental-placement-interface",
-             "Resolve accidental collisions.",
-             "left-padding padding right-padding accidental-grobs positioning-done")
+               "Resolve accidental collisions.",
+
+               /* properties */
+               "accidental-grobs "
+               "direction "
+               "padding "
+               "positioning-done "
+               "right-padding "
+               "script-priority "
+              );