LSST Applications g04c3c9f7ca+2075667efa,g1e125bf412+5f448d5fcf,g2079a07aa2+3e9fd84d81,g2305ad1205+b635cf1488,g2bbee38e9b+6c6beb4891,g337abbeb29+6c6beb4891,g33d1c0ed96+6c6beb4891,g3a166c0a6a+6c6beb4891,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+42f171e1e6,g5c3423f6d4+d536b04327,g607f77f49a+d536b04327,g6f43f06aed+ca1339dc19,g858d7b2824+d536b04327,g8ee334c5b4+d7f9608c2f,g9963eaa53e+b3dc1655d3,g998f4353bf+d536b04327,g99cad8db69+8ef2408349,g9ddcbc5298+9a081db1e4,ga1e77700b3+2cbb763275,gadfd92a7e4+aec2f3b930,gae0086650b+585e252eca,gb0e22166c9+0e73c8378f,gb3b7280ab2+cb5fdb229e,gbb8dafda3b+a327199e22,gc120e1dc64+88074880ea,gc28159a63d+6c6beb4891,gcdd4ae20e8+bd241b2308,gcde1bda545+903e937d91,gcf0d15dbbd+bd241b2308,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+6c6beb4891,gfbcc870c63+b310236976,w.2024.23
LSST Data Management Base Package
Loading...
Searching...
No Matches
MultiShapeletBasis.h
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2/*
3 * LSST Data Management System
4 * Copyright 2008-2013 LSST Corporation.
5 *
6 * This product includes software developed by the
7 * LSST Project (http://www.lsst.org/).
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 LSST License Statement and
20 * the GNU General Public License along with this program. If not,
21 * see <http://www.lsstcorp.org/LegalNotices/>.
22 */
23#ifndef LSST_SHAPELET_MultiShapeletBasis_h_INCLUDED
24#define LSST_SHAPELET_MultiShapeletBasis_h_INCLUDED
25
27
28namespace lsst { namespace shapelet {
29
41public:
42
58 MultiShapeletBasisComponent(double radius, int order, ndarray::Array<double const,2,2> const & matrix);
59
61 double getRadius() const { return _radius; }
62
64 int getOrder() const { return _order; }
65
67 ndarray::Array<double const,2,2> getMatrix() const { return _matrix; }
68
69private:
70 friend class MultiShapeletBasis;
71
72 double _radius;
73 int _order;
74 ndarray::Array<double const,2,2> _matrix;
75};
76
94public:
97 typedef ComponentVector::const_iterator Iterator;
98
100 explicit MultiShapeletBasis(int size);
101
103 int getSize() const { return _size; }
104
106 int getComponentCount() const { return _components.size(); }
107
109
110 Iterator begin() const { return _components.begin(); }
111 Iterator end() const { return _components.end(); }
113
119 void addComponent(double radius, int order, ndarray::Array<double const,2,2> const & matrix);
120
122 void scale(double factor);
123
125 void normalize();
126
128 void merge(MultiShapeletBasis const & other);
129
138 afw::geom::ellipses::Ellipse const & ellipse,
139 ndarray::Array<double const,1,1> const & coefficients
140 ) const;
141
142private:
143 int _size;
144 ComponentVector _components;
145};
146
147}} // namespace lsst::shapelet
148
149#endif // !LSST_SHAPELET_MultiShapeletBasis_h_INCLUDED
ndarray::Array< double const, 2, 2 > coefficients
T begin(T... args)
An ellipse defined by an arbitrary BaseCore and a center point.
Definition Ellipse.h:51
Simple struct that represents one shapelet expansion in a MultiShapeletBasis.
MultiShapeletBasisComponent(double radius, int order, ndarray::Array< double const, 2, 2 > const &matrix)
Main constructor for MultiShapeletBasisComponent.
int getOrder() const
Order of this shapelet expansion.
double getRadius() const
Return the radius of this shapelet expansion.
ndarray::Array< double const, 2, 2 > getMatrix() const
Matrix whose elements [i,j] map MultiShapeletBasis elements j to shapelet terms i.
A basis formed from a linear combination of shapelet bases that differ only in radius.
void merge(MultiShapeletBasis const &other)
Combine the given basis with this (in place), by appending its elements.
MultiShapeletBasis(int size)
Construct a MultiShapeletBasis with the given number of elements (i.e. free amplitudes).
void normalize()
Rescale all matrices so each element has unit flux.
Iterator begin() const
Iterator over the components (distinct shapelet bases) of the MultiShapeletBasis.
ComponentVector::const_iterator Iterator
int getComponentCount() const
Return the number of components (distinct shapelet bases) in the MultiShapeletBasis.
int getSize() const
Return the number of elements (i.e. free amplitudes) in the MultiShapeletBasis.
MultiShapeletBasisComponent Component
void scale(double factor)
Multiply the radius of all basis elements by the given factor.
std::vector< Component > ComponentVector
MultiShapeletFunction makeFunction(afw::geom::ellipses::Ellipse const &ellipse, ndarray::Array< double const, 1, 1 > const &coefficients) const
Create a MultiShapeletFunction from the basis.
void addComponent(double radius, int order, ndarray::Array< double const, 2, 2 > const &matrix)
Add a new component (shapelet basis) to the MultiShapeletBasis.
A multi-scale shapelet function.
T end(T... args)
T size(T... args)
table::Key< int > order