sbasis-curve.h revision d37634d73670180f99a3e0ea583621373d90ec4f
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm/**
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * \file
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * \brief Symmetric power basis curve
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *//*
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * Authors:
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * MenTaLguY <mental@rydia.net>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * Marco Cecchetti <mrcekets at gmail.com>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * Krzysztof KosiƄski <tweenk.pl@gmail.com>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * Copyright 2007-2009 Authors
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * This library is free software; you can redistribute it and/or
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * modify it either under the terms of the GNU Lesser General Public
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * License version 2.1 as published by the Free Software Foundation
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * (the "LGPL") or, at your option, under the terms of the Mozilla
64caa91f2899a6648503a75dc7310841955b74fdJucaBlues * Public License Version 1.1 (the "MPL"). If you do not alter this
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * notice, a recipient may use your version of this file under either
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * the MPL or the LGPL.
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * You should have received a copy of the LGPL along with this library
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * in the file COPYING-LGPL-2.1; if not, write to the Free Software
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * You should have received a copy of the MPL along with this library
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * in the file COPYING-MPL-1.1
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen *
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * The contents of this file are subject to the Mozilla Public License
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * Version 1.1 (the "License"); you may not use this file except in
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * compliance with the License. You may obtain a copy of the License at
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen * http://www.mozilla.org/MPL/
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * OF ANY KIND, either express or implied. See the LGPL or the MPL for
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * the specific language governing rights and limitations.
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm */
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#ifndef _2GEOM_SBASIS_CURVE_H_
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#define _2GEOM_SBASIS_CURVE_H_
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#include <2geom/curve.h>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#include <2geom/nearest-point.h>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#include <2geom/sbasis-geometric.h>
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrmnamespace Geom
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen{
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen
d431763a9ec8059aa4962688de8144319969fb0fjohanengelen/** @brief Symmetric power basis curve.
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * Symmetric power basis (S-basis for short) polynomials are a versatile numeric representation
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * of arbitrary continuous curves. They combine the properties of Bezier curves
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * (geometric interpretation of parameters, numerical stability near ends of the curve)
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * and the monomial basis (fast evaluation). They are the main representation of curves
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * in 2Geom.
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm *
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm * @ingroup Curves
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm */
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrmclass SBasisCurve : public Curve {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrmprivate:
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm SBasisCurve();
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm D2<SBasis> inner;
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrmpublic:
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm explicit SBasisCurve(D2<SBasis> const &sb) : inner(sb) {}
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm explicit SBasisCurve(Curve const &other) : inner(other.toSBasis()) {}
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Curve *duplicate() const { return new SBasisCurve(*this); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Point initialPoint() const { return inner.at0(); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Point finalPoint() const { return inner.at1(); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual bool isDegenerate() const { return inner.isConstant(); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Point pointAt(Coord t) const { return inner.valueAt(t); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual std::vector<Point> pointAndDerivatives(Coord t, unsigned n) const {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm return inner.valueAndDerivatives(t, n);
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Coord valueAt(Coord t, Dim2 d) const { return inner[d].valueAt(t); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual void setInitial(Point const &v) {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm for (unsigned d = 0; d < 2; d++) { inner[d][0][0] = v[d]; }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual void setFinal(Point const &v) {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm for (unsigned d = 0; d < 2; d++) { inner[d][0][1] = v[d]; }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Rect boundsFast() const { return *bounds_fast(inner); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Rect boundsExact() const { return *bounds_exact(inner); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual OptRect boundsLocal(OptInterval const &i, unsigned deg) const {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm return bounds_local(inner, i, deg);
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm }
c0cd5511d3b975ebe07d019c1f5528108725e438johanengelen virtual std::vector<Coord> roots(Coord v, Dim2 d) const { return Geom::roots(inner[d] - v); }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm virtual Coord nearestPoint( Point const& p, Coord from = 0, Coord to = 1 ) const {
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen return nearest_point(p, inner, from, to);
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual std::vector<Coord> allNearestPoints( Point const& p, Coord from = 0,
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen Coord to = 1 ) const
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen {
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen return all_nearest_points(p, inner, from, to);
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual Coord length(Coord tolerance) const { return ::Geom::length(inner, tolerance); }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual Curve *portion(Coord f, Coord t) const {
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen return new SBasisCurve(Geom::portion(inner, f, t));
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual Curve *transformed(Affine const &m) const { return new SBasisCurve(inner * m); }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual Curve *derivative() const {
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen return new SBasisCurve(Geom::derivative(inner));
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual D2<SBasis> toSBasis() const { return inner; }
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen virtual int degreesOfFreedom() const {
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm return inner[0].degreesOfFreedom() + inner[1].degreesOfFreedom();
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen }
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm};
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen} // end namespace Geom
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm#endif // _2GEOM_SBASIS_CURVE_H_
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
7073d105e612f7dc898c292742bee9655d2a51b2johanengelen
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm/*
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm Local Variables:
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm mode:c++
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm c-file-style:"stroustrup"
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm indent-tabs-mode:nil
f07bfd5a05d43a6d11f7cd442f085149092dea88pjrm fill-column:99
0903335a0099bd7ee779925f43a15a2216a0e863johanengelen End:
0903335a0099bd7ee779925f43a15a2216a0e863johanengelen*/
0903335a0099bd7ee779925f43a15a2216a0e863johanengelen// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99 :
0903335a0099bd7ee779925f43a15a2216a0e863johanengelen