LSST Applications g0aad566f14+ff8335139e,g2079a07aa2+86d27d4dc4,g2305ad1205+f692bcb212,g2bbee38e9b+047b288a59,g337abbeb29+047b288a59,g33d1c0ed96+047b288a59,g3a166c0a6a+047b288a59,g3d1719c13e+ff8335139e,g487adcacf7+0ce5160e95,g4be5004598+ff8335139e,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+d1f8abcd0b,g63cd9335cc+585e252eca,g858d7b2824+ff8335139e,g88963caddf+5804a2ccce,g88c60e1ba4+f263127a16,g99cad8db69+5fb0aeae0c,g9ddcbc5298+9a081db1e4,ga1e77700b3+a912195c07,ga57fefb910+66aa5c8fea,gab192c5e6e+cead8cd1b9,gae0086650b+585e252eca,gb0e22166c9+60f28cb32d,gb1fda7adca+d2c2560e0a,gb3a676b8dc+b4feba26a1,gb4b16eec92+fe5f3eb638,gba4ed39666+c2a2e4ac27,gbb8dafda3b+d91932a454,gc120e1dc64+a1259a995c,gc28159a63d+047b288a59,gc3e9b769f7+c20394caa9,gcf0d15dbbd+66aa5c8fea,gdaeeff99f8+f9a426f77a,ge79ae78c31+047b288a59,w.2024.20
LSST Data Management Base Package
Loading...
Searching...
No Matches
constants.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_CONSTANTS_H
26#define LSST_AFW_MATH_SHAPELETS_CONSTANTS_H
27
37#include "ndarray.h"
38#include "lsst/pex/exceptions.h"
39
40namespace lsst { namespace shapelet {
41
42extern double const BASIS_NORMALIZATION;
43
45
53 HERMITE,
71 LAGUERRE,
91};
92
94inline int computeOffset(int order) { return order * (order + 1) / 2; }
95
97inline int computeSize(int order) { return computeOffset(order + 1); }
98
104inline int computeOrder(int size) {
105 int order = (std::sqrt(8*size + 1) - 3)/2;
106 if (computeSize(order) != size) {
107 throw LSST_EXCEPT(
109 "Invalid size for shapelet coefficient matrix"
110 );
111 }
112 return order;
113}
114
121typedef ndarray::Array<double,1> Array1d;
122
124inline double intSqrt(int n) {
125 return std::sqrt(double(n));
126}
127
129inline double rationalSqrt(int n, int d) {
130 return std::sqrt(double(n) / double(d));
131}
132
133}} // namespace lsst::shapelet
134
135#endif // !defined(LSST_AFW_MATH_SHAPELETS_CONSTANTS_H)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
An ellipse core with quadrupole moments as parameters.
Definition Quadrupole.h:47
Reports invalid arguments.
Definition Runtime.h:66
int computeOrder(int size)
Infer the order of a shapelet expansion from the number of coefficients.
Definition constants.h:104
int computeOffset(int order)
Return the offset of the given order in a coefficient vector.
Definition constants.h:94
double intSqrt(int n)
Compute the square root of an integer number.
Definition constants.h:124
int computeSize(int order)
Return the size of the coefficient vector for the given order.
Definition constants.h:97
double const BASIS_NORMALIZATION
Normalization factor for 1-d orthonormal shapelets: pi^(-1/4)
double rationalSqrt(int n, int d)
Compute the square root of a rational number i.e. sqrt(n/d)
Definition constants.h:129
ndarray::Array< double, 1 > Array1d
Typedef for a commonly-used array type.
Definition constants.h:121
@ HERMITE
Cartesian shapelets or Gauss-Hermite functions, as defined in Refregier, 2003.
Definition constants.h:53
@ LAGUERRE
Polar shapelets or Gauss-Laguerre functions, as defined in Bernstein and Jarvis, 2002.
Definition constants.h:71
afw::geom::ellipses::Quadrupole EllipseCore
Definition constants.h:44
T sqrt(T... args)
table::Key< int > order