X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=flower%2Foffset.cc;h=bb97b9407e4eebbd62dd963f3b635a4049ec9325;hb=638ad29c79717fc1ba895b14aee256dd3a03f25b;hp=1fae9a3b3cedc99660e67cae5460dfd7b024970a;hpb=c659cb200486c2f908703696a1b2873e78c8160a;p=lilypond.git diff --git a/flower/offset.cc b/flower/offset.cc index 1fae9a3b3c..bb97b9407e 100644 --- a/flower/offset.cc +++ b/flower/offset.cc @@ -1,54 +1,55 @@ /* - 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--2011 Han-Wen Nienhuys - (c) 1997--2004 Han-Wen Nienhuys + 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 . */ #include "offset.hh" -#include - #ifndef STANDALONE -String +string Offset::to_string () const { - String s; - s = String (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", " - + ::to_string (coordinate_a_[Y_AXIS]) + ")"; + string s; + s = string (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", " + + ::to_string (coordinate_a_[Y_AXIS]) + ")"; return s; } #endif - -bool -isinf_b (Real r) -{ - return (fabs (r) > 1e20); -} - /* free bsd fix by John Galbraith - */ - +*/ + Offset complex_multiply (Offset z1, Offset z2) { Offset z; - if (!isinf_b (z2[Y_AXIS])) - { - z[X_AXIS] = z1[X_AXIS] * z2[X_AXIS] - z1[Y_AXIS]*z2[Y_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]; + o[Y_AXIS] = -o[Y_AXIS]; return o; } @@ -57,21 +58,19 @@ complex_divide (Offset z1, Offset z2) { z2 = complex_conjugate (z2); Offset z = complex_multiply (z1, z2); - z *= 1/z2.length (); + z *= 1 / z2.length (); return z; } - - Offset complex_exp (Offset o) { 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 @@ -80,11 +79,40 @@ Offset::arg () const 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 (coordinate_a_[X_AXIS]) + sqr (coordinate_a_[Y_AXIS])); + return sqrt (sqr (coordinate_a_[X_AXIS]) + + sqr (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; + d /= length (); + return d; +} + +Offset +Offset::swapped () const +{ + return Offset (coordinate_a_[Y_AXIS], coordinate_a_[X_AXIS]); }