/*
- misc.cc -- implement various stuff
+ This file is part of LilyPond, the GNU music typesetter.
- source file of the GNU LilyPond music typesetter
+ Copyright (C) 1997--2015 Han-Wen Nienhuys <hanwen@xs4all.nl>
+ Jan Nieuwenhuizen <janneke@gnu.org>
- (c) 1997--1998, 1998 Han-Wen Nienhuys <hanwen@cs.uu.nl>
- Jan Nieuwenhuizen <janneke@gnu.org>
+ 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.
+
+ 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 <complex>
#include "misc.hh"
+#include "offset.hh"
+#include "warn.hh"
-#ifndef STANDALONE
-#include "item.hh"
-#endif
-
+/*
+ Return the 2-log, rounded down
+*/
int
-intlog2(int d) {
- int i=0;
- while (!(d&1))
+intlog2 (int d)
+{
+ if (d <= 0)
+ error ("intlog2 with negative argument: " + ::to_string (d));
+ int i = 0;
+ while ((d != 1))
{
- d/= 2;
- i++;
+ d /= 2;
+ i++;
}
- assert (!(d/2));
+
+ assert (! (d / 2));
return i;
}
double
-log_2(double x) {
- return log (x) /log (2.0);
+log_2 (double x)
+{
+ return log (x) / log (2.0);
}
-#ifndef STANDALONE
-Interval
-itemlist_width (const Array<Item*> &its)
+Real
+directed_round (Real f, Direction d)
{
- Interval iv ;
- iv.set_empty();
-
- for (int j =0; j < its.size(); j++)
- {
- iv.unite (its[j]->width());
-
- }
- return iv;
+ if (d < 0)
+ return floor (f);
+ else
+ return ceil (f);
}
-#endif
-
-
/*
- TODO
- group in some Array_*
- make more generic / templatise
+ 0 at threshold, 1 at 0, with 1/x falloff.
*/
-int
-get_lower_bound (Array<Real> const& positions, Real x)
+Real
+peak_around (Real epsilon, Real threshold, Real x)
{
- if (x < positions[0])
- return 0;
- for (int i = 1; i < positions.size (); i++)
- if (x < positions[i])
- return i - 1;
- return positions.size () - 1;
+ if (x < 0)
+ return 1.0;
+ return max (- epsilon * (x - threshold) / ((x + epsilon) * threshold), 0.0);
}
-Slice
-get_bounds_slice (Array<Real> const& positions, Real x)
+/*
+ 0 at 0, 1 at standard_x, and increasing thereafter.
+ */
+Real
+convex_amplifier (Real standard_x, Real increase_factor, Real x)
{
- int l = get_lower_bound (positions, x);
- int u = positions.size () - 1 <? l + 1;
- if (x < positions[l])
- u = l;
- return Slice (l, u);
+ return (exp (increase_factor * x / standard_x) - 1.0) / (exp (increase_factor) - 1.0);
}
-Interval
-get_bounds_iv (Array<Real> const& positions, Real x)
+string
+camel_case_to_lisp_identifier (const string &in)
{
- Slice slice = get_bounds_slice (positions, x);
- return Interval (positions[slice.min ()], positions[slice.max ()]);
-}
+ vector<char> out;
-// silly name
-Interval
-quantise_iv (Array<Real> const& positions, Real period, Real x)
-{
- /*
- ugh
- assume that
- * positions are sorted,
- * positions are nonnegative
- * period starts at zero
- */
-
- int n = (int)(x / period);
- Real frac = (x / period - n) * period;
- if (frac < 0)
+ /* don't add '-' before first character */
+ out.push_back (char (tolower (in[0])));
+
+ for (size_t inpos = 1; inpos < in.size (); inpos++)
{
- frac += period;
- n--;
+ if (isupper (in[inpos]))
+ out.push_back ('-');
+ out.push_back ( char (tolower (in[inpos])));
}
- Slice slice = get_bounds_slice (positions, frac);
- Interval iv(positions[slice.min ()], positions[slice.max ()]);
+ string result (&out[0], out.size ());
+ replace_all (&result, '_', '-');
- if (slice.min () == slice.max ())
- {
- if (slice.min () == 0)
- iv.min () = - period + positions.top ();
- else
- iv.max () = period + positions[0];
- }
+ return result;
+}
+
+Offset
+get_point_in_y_direction (Offset orig, Real slope, Real dist, Direction dir)
+{
+ if (slope == infinity_f)
+ return orig + Offset (dir * dist, 0.0);
+
+ Real x = slope == 0.0 ? 1.0 * dir : 1.0 * sign (slope) * dir;
+ Real y = slope * x;
+ Real angle = atan2 (y, x);
- iv += period * n;
+ complex<Real> orig_c (orig[X_AXIS], orig[Y_AXIS]);
+ complex<Real> to_move = polar (dist, angle);
+ complex<Real> res = orig_c + to_move;
- return iv;
+ return Offset (real (res), imag (res));
}