/*
- offset.cc -- implement Offset
+ 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>
- (c) 1997--1999 Han-Wen Nienhuys <hanwen@cs.uu.nl>
+ 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>
-#ifndef STANDALONE
-#include "string.hh"
-#endif
#include "offset.hh"
-
#ifndef STANDALONE
-String
-Offset::str () const
+string
+Offset::to_string () const
{
- String s;
- s = String("(") + to_str (coordinate_a_[X_AXIS]) + ", "
- + to_str (coordinate_a_[Y_AXIS]) + ")";
+ string s;
+ s = string (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", "
+ + ::to_string (coordinate_a_[Y_AXIS]) + ")";
return s;
}
#endif
+/*
+ free bsd fix by John Galbraith
+*/
+
Offset
complex_multiply (Offset z1, Offset z2)
{
Offset z;
- z[X_AXIS] = z1[X_AXIS] * z2[X_AXIS] - z1[Y_AXIS]*z2[Y_AXIS];
- z[Y_AXIS] = z1[X_AXIS] * z2[Y_AXIS] + z1[Y_AXIS] * z2[X_AXIS];
+ if (!isinf (z2[Y_AXIS]))
+ {
+ z[X_AXIS] = z1[X_AXIS] * z2[X_AXIS] - z1[Y_AXIS] * z2[Y_AXIS];
+ z[Y_AXIS] = z1[X_AXIS] * z2[Y_AXIS] + z1[Y_AXIS] * z2[X_AXIS];
+ }
+ return z;
+}
+
+Offset
+complex_conjugate (Offset o)
+{
+ o[Y_AXIS] = -o[Y_AXIS];
+ return o;
+}
+
+Offset
+complex_divide (Offset z1, Offset z2)
+{
+ z2 = complex_conjugate (z2);
+ Offset z = complex_multiply (z1, z2);
+ z *= 1 / z2.length ();
return z;
}
{
Real s = sin (o[Y_AXIS]);
Real c = cos (o[Y_AXIS]);
-
+
Real r = exp (o[X_AXIS]);
- return Offset(r*c, r*s);
+ return Offset (r * c, r * s);
}
Real
Offset::arg () const
{
- return atan2 (y (), x());
+ return atan2 (coordinate_a_[Y_AXIS], coordinate_a_[X_AXIS]);
}
+Real
+Offset::angle_degrees () const
+{
+ return arg () * 180 / M_PI;
+}
/**
euclidian vector length / complex modulus
- */
+*/
Real
Offset::length () const
{
- return sqrt (sqr (x()) + sqr (y()));
+ return hypot (coordinate_a_[X_AXIS], coordinate_a_[Y_AXIS]);
+}
+
+bool
+Offset::is_sane () const
+{
+ return !isnan (coordinate_a_[X_AXIS])
+ && !isnan (coordinate_a_ [Y_AXIS])
+ && !isinf (coordinate_a_[X_AXIS])
+ && !isinf (coordinate_a_[Y_AXIS]);
+}
+
+Offset
+Offset::direction () const
+{
+ Offset d = *this;
+ if (isinf (d[X_AXIS]))
+ {
+ if (!isinf (d[Y_AXIS]))
+ return Offset ((d[X_AXIS] > 0.0 ? 1.0 : -1.0), 0.0);
+ }
+ else if (isinf (d[Y_AXIS]))
+ return Offset (0.0, (d[Y_AXIS] > 0.0 ? 1.0 : -1.0));
+ else if (d[X_AXIS] == 0.0 && d[Y_AXIS] == 0.0)
+ return d;
+ // The other cases propagate or produce NaN as appropriate.
+
+ d /= length ();
+ return d;
}
-void
-Offset::mirror (Axis a)
+
+Offset
+Offset::swapped () const
{
- coordinate_a_[a] = - coordinate_a_[a];
+ return Offset (coordinate_a_[Y_AXIS], coordinate_a_[X_AXIS]);
}