2 This file is part of LilyPond, the GNU music typesetter.
4 Copyright (C) 1997--2015 Han-Wen Nienhuys <hanwen@xs4all.nl>
6 LilyPond is free software: you can redistribute it and/or modify
7 it under the terms of the GNU General Public License as published by
8 the Free Software Foundation, either version 3 of the License, or
9 (at your option) any later version.
11 LilyPond is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
16 You should have received a copy of the GNU General Public License
17 along with LilyPond. If not, see <http://www.gnu.org/licenses/>.
26 Offset::to_string () const
29 s = string (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", "
30 + ::to_string (coordinate_a_[Y_AXIS]) + ")";
36 free bsd fix by John Galbraith
40 complex_multiply (Offset z1, Offset z2)
43 if (!isinf (z2[Y_AXIS]))
45 z[X_AXIS] = z1[X_AXIS] * z2[X_AXIS] - z1[Y_AXIS] * z2[Y_AXIS];
46 z[Y_AXIS] = z1[X_AXIS] * z2[Y_AXIS] + z1[Y_AXIS] * z2[X_AXIS];
52 complex_conjugate (Offset o)
54 o[Y_AXIS] = -o[Y_AXIS];
59 complex_divide (Offset z1, Offset z2)
61 z2 = complex_conjugate (z2);
62 Offset z = complex_multiply (z1, z2);
63 z *= 1 / z2.length ();
68 complex_exp (Offset o)
70 Real s = sin (o[Y_AXIS]);
71 Real c = cos (o[Y_AXIS]);
73 Real r = exp (o[X_AXIS]);
75 return Offset (r * c, r * s);
81 return atan2 (coordinate_a_[Y_AXIS], coordinate_a_[X_AXIS]);
85 Offset::angle_degrees () const
87 return arg () * 180 / M_PI;
90 euclidian vector length / complex modulus
93 Offset::length () const
95 return hypot (coordinate_a_[X_AXIS], coordinate_a_[Y_AXIS]);
99 Offset::is_sane () const
101 return !isnan (coordinate_a_[X_AXIS])
102 && !isnan (coordinate_a_ [Y_AXIS])
103 && !isinf (coordinate_a_[X_AXIS])
104 && !isinf (coordinate_a_[Y_AXIS]);
108 Offset::direction () const
111 if (isinf (d[X_AXIS]))
113 if (!isinf (d[Y_AXIS]))
114 return Offset ((d[X_AXIS] > 0.0 ? 1.0 : -1.0), 0.0);
116 else if (isinf (d[Y_AXIS]))
117 return Offset (0.0, (d[Y_AXIS] > 0.0 ? 1.0 : -1.0));
118 else if (d[X_AXIS] == 0.0 && d[Y_AXIS] == 0.0)
120 // The other cases propagate or produce NaN as appropriate.
127 Offset::swapped () const
129 return Offset (coordinate_a_[Y_AXIS], coordinate_a_[X_AXIS]);