LSST Applications g013ef56533+b8d55c8942,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+7a3b874d60,g210f2d0738+7416ca6900,g262e1987ae+1d557ba9a3,g29ae962dfc+519d34895e,g2cef7863aa+aef1011c0b,g30d7c61c20+36d16ea71a,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+cb326ad149,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5a60e81ecd+6240c63dbc,g64539dfbff+7416ca6900,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g74acd417e5+0bae3c876a,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+dee7680868,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+14b8164b5b,g98a1a72a9c+8389601a76,g98df359435+fff771c62d,gb8cb2b794d+6728931916,gbf99507273+8c5ae1fdc5,gc2a301910b+7416ca6900,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+0bae3c876a,ge410e46f29+f459a6810c,ge41e95a9f2+7416ca6900,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
GaussHermiteEvaluator.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008, 2009, 2010, 2011 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24
25#ifndef LSST_AFW_MATH_SHAPELETS_HERMITEEVALUATOR_H
26#define LSST_AFW_MATH_SHAPELETS_HERMITEEVALUATOR_H
27
28#include "ndarray.h"
29#include "lsst/geom.h"
31#include "Eigen/Core"
32
33namespace lsst { namespace shapelet {
34
49public:
50
51 static int const computeOffset(int order) { return order*(order+1)/2; }
52 static int const computeIndex(int x, int y) { return computeOffset(x+y) + x; }
53
55 ++_i;
56 if (--_y < 0) {
57 _x = 0;
58 _y = ++_n;
59 } else {
60 ++_x;
61 }
62 return *this;
63 }
64
65 int const getOrder() const { return _n; }
66 int const getX() const { return _x; }
67 int const getY() const { return _y; }
68
69 int const getIndex() const { return _i; }
70
71 PackedIndex() : _n(0), _i(0), _x(0), _y(0) {}
72 PackedIndex(int const x, int const y) : _n(x+y), _i(computeOffset(_n) + x), _x(x), _y(y) {}
73
74private:
75 int _n;
76 int _i;
77 int _x;
78 int _y;
79};
80
85public:
86
95 static Eigen::MatrixXd computeInnerProductMatrix(
96 int rowOrder, int colOrder, double a, double b
97 );
98
99 int getOrder() const { return _xWorkspace.getSize<0>() - 1; }
100
106 Array1d const & target, double x, double y,
107 Array1d const & dx = Array1d(),
108 Array1d const & dy = Array1d()
109 ) const;
110
116 Array1d const & target, geom::Point2D const & point,
117 Array1d const & dx = Array1d(),
118 Array1d const & dy = Array1d()
119 ) const {
120 fillEvaluation(target, point.getX(), point.getY(), dx, dy);
121 }
122
128 Array1d const & target, geom::Extent2D const & point,
129 Array1d const & dx = Array1d(),
130 Array1d const & dy = Array1d()
131 ) const {
132 fillEvaluation(target, point.getX(), point.getY(), dx, dy);
133 }
134
139 void fillIntegration(Array1d const & target, int xMoment=0, int yMoment=0) const;
140
145 ndarray::Array<double const,1> const & coeff, double x, double y,
146 double * dx = 0, double * dy = 0
147 ) const;
148
153 ndarray::Array<double const,1> const & coeff, geom::Point2D const & point,
154 double * dx = 0, double * dy = 0
155 ) const {
156 return sumEvaluation(coeff, point.getX(), point.getY(), dx, dy);
157 }
158
163 ndarray::Array<double const,1> const & coeff, geom::Extent2D const & point,
164 double * dx = 0, double * dy = 0
165 ) const {
166 return sumEvaluation(coeff, point.getX(), point.getY(), dx, dy);
167 }
168
172 double sumIntegration(ndarray::Array<double const,1> const & coeff, int xMoment=0, int yMoment=0) const;
173
174 explicit GaussHermiteEvaluator(int order);
175
176private:
177
178 ndarray::Array<double,1,1> _xWorkspace;
179 ndarray::Array<double,1,1> _yWorkspace;
180 ndarray::Array<double,1,1> _dxWorkspace;
181 ndarray::Array<double,1,1> _dyWorkspace;
182};
183
184}} // namespace lsst::shapelet
185
186#endif // !defined(LSST_AFW_MATH_SHAPELETS_HERMITEEVALUATOR_H)
double sumEvaluation(ndarray::Array< double const, 1 > const &coeff, geom::Extent2D const &point, double *dx=0, double *dy=0) const
Evaluate a simple unscaled shapelet expansion at the given point.
static Eigen::MatrixXd computeInnerProductMatrix(int rowOrder, int colOrder, double a, double b)
Fill a matrix with the function inner products of two HERMITE shapelet basis functions with different...
void fillEvaluation(Array1d const &target, geom::Extent2D const &point, Array1d const &dx=Array1d(), Array1d const &dy=Array1d()) const
Fill a vector whose dot product with a HERMITE coefficient vector evaluates a simple unscaled shapele...
double sumIntegration(ndarray::Array< double const, 1 > const &coeff, int xMoment=0, int yMoment=0) const
Integrate a simple unscaled shapelet expansion.
void fillEvaluation(Array1d const &target, double x, double y, Array1d const &dx=Array1d(), Array1d const &dy=Array1d()) const
Fill a vector whose dot product with a HERMITE coefficient vector evaluates a simple unscaled shapele...
double sumEvaluation(ndarray::Array< double const, 1 > const &coeff, double x, double y, double *dx=0, double *dy=0) const
Evaluate a simple unscaled shapelet expansion at the given point.
void fillEvaluation(Array1d const &target, geom::Point2D const &point, Array1d const &dx=Array1d(), Array1d const &dy=Array1d()) const
Fill a vector whose dot product with a HERMITE coefficient vector evaluates a simple unscaled shapele...
double sumEvaluation(ndarray::Array< double const, 1 > const &coeff, geom::Point2D const &point, double *dx=0, double *dy=0) const
Evaluate a simple unscaled shapelet expansion at the given point.
void fillIntegration(Array1d const &target, int xMoment=0, int yMoment=0) const
Fill a vector whose dot product with a HERMITE coefficient vector integrates a simple unscaled shapel...
static int const computeIndex(int x, int y)
PackedIndex(int const x, int const y)
static int const computeOffset(int order)
Extent< double, 2 > Extent2D
Definition Extent.h:400
Point< double, 2 > Point2D
Definition Point.h:324
ndarray::Array< double, 1 > Array1d
Typedef for a commonly-used array type.
Definition constants.h:121
Constants, typedefs, and general-purpose functions for shapelets library.