LSST Applications g013ef56533+7c9321ec0f,g042eb84c57+c6cfa41bc3,g199a45376c+0ba108daf9,g1fd858c14a+fcad0d0313,g210f2d0738+c0f94c6586,g262e1987ae+a7e710680e,g29ae962dfc+fb55f2edb0,g2ac17093b6+61d6563b1e,g2b1d02342f+df6f932764,g2cef7863aa+aef1011c0b,g2f7ad74990+c0f94c6586,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+53cf87ae69,g47891489e3+4316d04fff,g511e8cfd20+baa56acf6c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+fd7ad03fde,g64539dfbff+c0f94c6586,g67b6fd64d1+4316d04fff,g67fd3c3899+c0f94c6586,g6985122a63+4316d04fff,g74acd417e5+ca833bee28,g786e29fd12+668abc6043,g81db2e9a8d+b2ec8e584f,g87389fa792+8856018cbb,g89139ef638+4316d04fff,g8d7436a09f+0a24083b20,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+11eb8fd5b8,gbf99507273+8c5ae1fdc5,gcdda8b9158+e4c84c9d5c,gce8aa8abaa+8c5ae1fdc5,gd7ef33dd92+4316d04fff,gdab6d2f7ff+ca833bee28,ge410e46f29+4316d04fff,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.40
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.