transforms.h

Go to the documentation of this file.

/*
 * Authors:
 *   ? <?@?.?>
 *   Krzysztof Kosiński <tweenk.pl@gmail.com>
 *   Johan Engelen
 * 
 * Copyright ?-2012 Authors
 *
 * This library is free software; you can redistribute it and/or
 * modify it either under the terms of the GNU Lesser General Public
 * License version 2.1 as published by the Free Software Foundation
 * (the "LGPL") or, at your option, under the terms of the Mozilla
 * Public License Version 1.1 (the "MPL"). If you do not alter this
 * notice, a recipient may use your version of this file under either
 * the MPL or the LGPL.
 *
 * You should have received a copy of the LGPL along with this library
 * in the file COPYING-LGPL-2.1; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
 * You should have received a copy of the MPL along with this library
 * in the file COPYING-MPL-1.1
 *
 * The contents of this file are subject to the Mozilla Public License
 * Version 1.1 (the "License"); you may not use this file except in
 * compliance with the License. You may obtain a copy of the License at
 * http://www.mozilla.org/MPL/
 *
 * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY
 * OF ANY KIND, either express or implied. See the LGPL or the MPL for
 * the specific language governing rights and limitations.
 */
 
#ifndef LIB2GEOM_SEEN_TRANSFORMS_H
#define LIB2GEOM_SEEN_TRANSFORMS_H
 
#include <cmath>
#include <2geom/forward.h>
#include <2geom/affine.h>
#include <2geom/angle.h>
#include <boost/concept/assert.hpp>
 
namespace Geom {
 
template <typename T>
struct TransformConcept {
    T t, t2;
    Affine m;
    Point p;
    bool bool_;
    Coord epsilon;
    void constraints() {
        m = t;  //implicit conversion
        m *= t;
        m = m * t;
        m = t * m;
        p *= t;
        p = p * t;
        t *= t;
        t = t * t;
        t = pow(t, 3);
        bool_ = (t == t);
        bool_ = (t != t);
        t = T::identity();
        t = t.inverse();
        bool_ = are_near(t, t2);
        bool_ = are_near(t, t2, epsilon);
    }
};
 
template <typename T>
class TransformOperations
    : boost::equality_comparable< T
    , boost::multipliable< T
      > >
{
public:
    template <typename T2>
    Affine operator*(T2 const &t) const {
        Affine ret(*static_cast<T const*>(this)); ret *= t; return ret;
    }
};
 
template <typename T>
T pow(T const &t, int n) {
    BOOST_CONCEPT_ASSERT((TransformConcept<T>));
    if (n == 0) return T::identity();
    T result(T::identity());
    T x(n < 0 ? t.inverse() : t);
    if (n < 0) n = -n;
    while ( n ) { // binary exponentiation - fast
        if ( n & 1 ) { result *= x; --n; }
        x *= x; n /= 2;
    }
    return result;
}
 
class Translate
    : public TransformOperations< Translate >
{
    Point vec;
public:
    Translate() = default;
    explicit Translate(Point const &p) : vec(p) {}
    Translate(Coord x, Coord y) : vec(x, y) {}
 
    operator Affine() const { return Affine(1, 0, 0, 1, vec[X], vec[Y]); }
    Coord operator[](Dim2 dim) const { return vec[dim]; }
    Coord operator[](unsigned dim) const { return vec[dim]; }
    Translate &operator*=(Translate const &o) { vec += o.vec; return *this; }
    bool operator==(Translate const &o) const { return vec == o.vec; }
 
    Point vector() const { return vec; }
    Translate inverse() const { return Translate(-vec); }
    static Translate identity() { return {}; }
 
    friend class Point;
};
 
inline bool are_near(Translate const &a, Translate const &b, Coord eps = EPSILON) {
    return are_near(a[X], b[X], eps) && are_near(a[Y], b[Y], eps);
}
 
class Scale
    : public TransformOperations< Scale >
{
    Point vec = { 1, 1 };
public:
    Scale() = default;
    explicit Scale(Point const &p) : vec(p) {}
    Scale(Coord x, Coord y) : vec(x, y) {}
    explicit Scale(Coord s) : vec(s, s) {}
    inline operator Affine() const { return Affine(vec[X], 0, 0, vec[Y], 0, 0); }
 
    Coord operator[](Dim2 d) const { return vec[d]; }
    Coord operator[](unsigned d) const { return vec[d]; }
    //TODO: should we keep these mutators? add them to the other transforms?
    Coord &operator[](Dim2 d) { return vec[d]; }
    Coord &operator[](unsigned d) { return vec[d]; }
    Scale &operator*=(Scale const &b) { vec[X] *= b[X]; vec[Y] *= b[Y]; return *this; }
    bool operator==(Scale const &o) const { return vec == o.vec; }
 
    Point vector() const { return vec; }
    Scale inverse() const { return Scale(1./vec[0], 1./vec[1]); }
    static Scale identity() { return {}; }
 
    friend class Point;
};
 
inline bool are_near(Scale const &a, Scale const &b, Coord eps=EPSILON) {
    return are_near(a[X], b[X], eps) && are_near(a[Y], b[Y], eps);
}
 
class Rotate
    : public TransformOperations< Rotate >
{
    Point vec = { 1, 0 };
public:
    Rotate() = default;
    explicit Rotate(Coord theta) : vec(Point::polar(theta)) {}
    explicit Rotate(Point const &p) : vec(p.normalized()) {}
    explicit Rotate(Coord x, Coord y) : Rotate(Point(x, y)) {}
    operator Affine() const { return Affine(vec[X], vec[Y], -vec[Y], vec[X], 0, 0); }
 
    Point const &vector() const { return vec; }
    Coord angle() const { return atan2(vec); }
    Coord operator[](Dim2 dim) const { return vec[dim]; }
    Coord operator[](unsigned dim) const { return vec[dim]; }
    Rotate &operator*=(Rotate const &o) { vec *= o; return *this; }
    bool operator==(Rotate const &o) const { return vec == o.vec; }
    Rotate inverse() const {
        Rotate r;
        r.vec = Point(vec[X], -vec[Y]);
        return r;
    }
    static Rotate identity() { return {}; }
    static Rotate from_degrees(Coord deg) { return Rotate(rad_from_deg(deg)); }
    static Affine around(Point const &p, Coord angle);
 
    friend class Point;
};
 
inline bool are_near(Rotate const &a, Rotate const &b, Coord eps = EPSILON) {
    return are_near(a[X], b[X], eps) && are_near(a[Y], b[Y], eps);
}
 
template <typename S>
class ShearBase
    : public TransformOperations<S>
{
protected:
    Coord f = 0;
    ShearBase() = default;
    explicit ShearBase(Coord _f) : f(_f) {}
public:
    Coord factor() const { return f; }
    void setFactor(Coord nf) { f = nf; }
    S &operator*=(S const &s) { f += s.f; return static_cast<S &>(*this); }
    bool operator==(ShearBase<S> const &s) const { return f == s.f; }
    S inverse() const { return S(-f); }
    static S identity() { return {}; }
 
    friend class Point;
    friend class Affine;
};
 
class HShear
    : public ShearBase<HShear>
{
public:
    HShear() = default;
    explicit HShear(Coord h) : ShearBase<HShear>(h) {}
    operator Affine() const { return Affine(1, 0, f, 1, 0, 0); }
};
 
inline bool are_near(HShear const &a, HShear const &b, Coord eps=EPSILON) {
    return are_near(a.factor(), b.factor(), eps);
}
 
class VShear
    : public ShearBase<VShear>
{
public:
    VShear() = default;
    explicit VShear(Coord h) : ShearBase<VShear>(h) {}
    operator Affine() const { return Affine(1, f, 0, 1, 0, 0); }
};
 
inline bool are_near(VShear const &a, VShear const &b, Coord eps = EPSILON) {
    return are_near(a.factor(), b.factor(), eps);
}
 
class Zoom
    : public TransformOperations< Zoom >
{
    Coord _scale = 1;
    Point _trans;
public:
    Zoom() = default;
    explicit Zoom(Coord s) : _scale(s) {}
    explicit Zoom(Point const &t) : _trans(t) {}
    explicit Zoom(Translate const &t) : Zoom(t.vector()) {}
    Zoom(Coord s, Point const &t) : _scale(s), _trans(t) {}
    Zoom(Coord s, Translate const &t) : Zoom(s, t.vector()) {}
 
    operator Affine() const {
        return Affine(_scale, 0, 0, _scale, _trans[X] * _scale, _trans[Y] * _scale);
    }
    Zoom &operator*=(Zoom const &z) {
        _trans += z._trans / _scale;
        _scale *= z._scale;
        return *this;
    }
    bool operator==(Zoom const &z) const { return _scale == z._scale && _trans == z._trans; }
 
    Coord scale() const { return _scale; }
    void setScale(Coord s) { _scale = s; }
    Point translation() const { return _trans; }
    void setTranslation(Point const &p) { _trans = p; }
    Zoom inverse() const { return Zoom(1 / _scale, Translate(-_trans * _scale)); }
    static Zoom identity() { return {}; }
    static Zoom map_rect(Rect const &old_r, Rect const &new_r);
    
    friend class Point;
    friend class Affine;
};
 
inline bool are_near(Zoom const &a, Zoom const &b, Coord eps = EPSILON) {
    return are_near(a.scale(), b.scale(), eps) &&
           are_near(a.translation(), b.translation(), eps);
}
 
template<>
inline Scale pow(Scale const &s, int n) {
    return Scale(::pow(s[X], n), ::pow(s[Y], n));
}
template<>
inline Translate pow(Translate const &t, int n) {
    return Translate(t[X] * n, t[Y] * n);
}
 
Affine reflection(Point const & vector, Point const & origin);
 
//TODO: decomposition of Affine into some finite combination of the above classes
 
} // namespace Geom
 
#endif // LIB2GEOM_SEEN_TRANSFORMS_H
 
/*
  Local Variables:
  mode:c++
  c-file-style:"stroustrup"
  c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
  indent-tabs-mode:nil
  fill-column:99
  End:
*/
// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99 :