Loading [MathJax]/extensions/tex2jax.js
LSST Applications g04a91732dc+9666464c73,g0fba68d861+d6c6f70ffa,g1fd858c14a+94f68680cf,g208c678f98+1ca806343c,g271391ec13+ac98094cfc,g2c84ff76c0+120a924478,g2c9e612ef2+a92a2e6025,g35bb328faa+fcb1d3bbc8,g4d2262a081+7c332456db,g4e0f332c67+c58e4b632d,g53246c7159+fcb1d3bbc8,g60b5630c4e+a92a2e6025,g67b6fd64d1+9d1b2ab50a,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+506db7da85,g89139ef638+9d1b2ab50a,g8d6b6b353c+a92a2e6025,g9125e01d80+fcb1d3bbc8,g989de1cb63+9d1b2ab50a,g9f33ca652e+d1749da127,ga2b97cdc51+a92a2e6025,gabe3b4be73+1e0a283bba,gb1101e3267+6ecbd0580e,gb58c049af0+f03b321e39,gb89ab40317+9d1b2ab50a,gb90eeb9370+384e1fc23b,gcf25f946ba+506db7da85,gd315a588df+382ef11c06,gd6cbbdb0b4+75aa4b1db4,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+9984fb7d9e,ge278dab8ac+c61fbefdff,ge410e46f29+9d1b2ab50a,ge82c20c137+e12a08b75a,gf67bdafdda+9d1b2ab50a,gfd5510ef7b+864ead200a,v29.0.0.rc2
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
ScaledBasis2d.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2/*
3 * Developed for the LSST Data Management System.
4 * This product includes software developed by the LSST Project
5 * (https://www.lsst.org).
6 * See the COPYRIGHT file at the top-level directory of this distribution
7 * for details of code ownership.
8 *
9 * This program is free software: you can redistribute it and/or modify
10 * it under the terms of the GNU General Public License as published by
11 * the Free Software Foundation, either version 3 of the License, or
12 * (at your option) any later version.
13 *
14 * This program is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 * GNU General Public License for more details.
18 *
19 * You should have received a copy of the GNU General Public License
20 * along with this program. If not, see <https://www.gnu.org/licenses/>.
21 */
22#ifndef LSST_AFW_MATH_POLYNOMIALS_ScaledBasis2d_h_INCLUDED
23#define LSST_AFW_MATH_POLYNOMIALS_ScaledBasis2d_h_INCLUDED
24
26
27namespace lsst { namespace geom { namespace polynomials {
28
29template <typename Basis>
30class Function2d;
31
42template <typename Nested>
44public:
45
48
51
53 using Workspace = typename Nested::Workspace;
54
56 using IndexRange = typename Nested::IndexRange;
57
59 explicit ScaledBasis2d(Nested const & nested, Scaling2d const & scaling) :
60 _nested(nested),
61 _scaling(scaling)
62 {}
63
77 ScaledBasis2d(std::size_t order, Box2D const & box) :
78 _nested(order),
79 _scaling(makeUnitRangeScaling2d(box))
80 {}
81
83 ScaledBasis2d(ScaledBasis2d const &) = default;
84
87
90
93
95 Nested const & getNested() const noexcept { return _nested; }
96
98 Scaling2d const & getScaling() const noexcept { return _scaling; }
99
101 std::size_t getOrder() const { return getNested().getOrder(); }
102
104 std::size_t size() const { return getNested().size(); }
105
112 Scaled scaled(Scaling2d const & first) const {
113 return getNested().scaled(first.then(getScaling()));
114 }
115
117 int index(int x, int y) const { return getNested().index(x, y); }
118
139 IndexRange getIndices() const { return getNested().getIndices(); }
140
142 Workspace makeWorkspace() const { return getNested().makeWorkspace();}
143
162 template <typename Vector>
163 double sumWith(geom::Point2D const & point, Vector const & coefficients,
164 SumMode mode=SumMode::FAST) const {
165 return getNested().sumWith(getScaling().applyForward(point), coefficients, mode);
166 }
167
169 template <typename Vector>
170 double sumWith(geom::Point2D const & point, Vector const & coefficients,
171 Workspace & workspace, SumMode mode=SumMode::FAST) const {
172 return getNested().sumWith(getScaling().applyForward(point), coefficients, workspace, mode);
173 }
174
185 template <typename Vector>
186 void fill(geom::Point2D const & point, Vector && basis) const {
187 return getNested().fill(getScaling().applyForward(point),
188 std::forward<Vector>(basis));
189 }
190
192 template <typename Vector>
193 void fill(geom::Point2D const & point, Vector && basis, Workspace & workspace) const {
194 return getNested().fill(getScaling().applyForward(point),
196 workspace);
197 }
198
199private:
200 Nested _nested;
201 Scaling2d _scaling;
202};
203
204}}} // namespace lsst::geom::polynomials
205
206#endif // !LSST_AFW_MATH_POLYNOMIALS_ScaledBasis2d_h_INCLUDED
A floating-point coordinate rectangle geometry.
Definition Box.h:413
A 2-d function defined by a series expansion and its coefficients.
Definition Function2d.h:42
double sumWith(geom::Point2D const &point, Vector const &coefficients, Workspace &workspace, SumMode mode=SumMode::FAST) const
void fill(geom::Point2D const &point, Vector &&basis, Workspace &workspace) const
ScaledBasis2d(Nested const &nested, Scaling2d const &scaling)
Construct a scaled basis from a nested basis and a scaling transform.
Chebyshev1Basis2d< packing > const & getNested() const noexcept
void fill(geom::Point2D const &point, Vector &&basis) const
ScaledBasis2d & operator=(ScaledBasis2d const &)=default
double sumWith(geom::Point2D const &point, Vector const &coefficients, SumMode mode=SumMode::FAST) const
A 2-d separable affine transform that can be used to map one interval to another.
Definition Scaling2d.h:48
T forward(T... args)
Low-level polynomials (including special polynomials) in C++.
Scaling2d makeUnitRangeScaling2d(geom::Box2D const &box)
Return a Scaling1d that maps the given box to [-1, 1]x[-1, 1].
Definition Scaling2d.h:112
SumMode
Enum used to control how to sum polynomial terms.
Definition SafeSum.h:32
@ FAST
Summation using regular floating-point addition.
Definition SafeSum.h:34
Point< double, 2 > Point2D
Definition Point.h:324