File GaussHermiteEvaluator.h

namespace lsst

Class for a simple mapping implementing a generic AstrometryTransform.

Remove all non-astronomical counts from the Chunk Exposure’s pixels.

Forward declarations for lsst::utils::Cache

For details on the Cache class, see the Cache.h file.

It uses a template rather than a pointer so that the derived classes can use the specifics of the transform. The class simplePolyMapping overloads a few routines.

A base class for image defects

Numeric constants used by the Integrate.h integrator routines.

Compute Image Statistics

Note

Gauss-Kronrod-Patterson quadrature coefficients for use in quadpack routine qng. These coefficients were calculated with 101 decimal digit arithmetic by L. W. Fullerton, Bell Labs, Nov 1981.

Note

The Statistics class itself can only handle lsst::afw::image::MaskedImage() types. The philosophy has been to handle other types by making them look like lsst::afw::image::MaskedImage() and reusing that code. Users should have no need to instantiate a Statistics object directly, but should use the overloaded makeStatistics() factory functions.

namespace shapelet

class GaussHermiteEvaluator

#include <GaussHermiteEvaluator.h>

A class to evaluate HERMITE shapelet-related quantities.

Public Functions

int getOrder() const

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 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 shapelet expansion at the given point.

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 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 shapelet expansion.

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.

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.

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.

double sumIntegration(ndarray::Array<double const, 1> const &coeff, int xMoment = 0, int yMoment = 0) const

Integrate a simple unscaled shapelet expansion.

GaussHermiteEvaluator(int order)

Public Static Functions

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 scales. \( M_{\mathbf{i},\mathbf{j}} = \int d^2 \mathbf{x} \psi_\mathbf{i}(a\mathbf{x})\phi_\mathbf{j}(b\mathbf{x}) \).

Private Members

ndarray::Array<double, 1, 1> _xWorkspace

ndarray::Array<double, 1, 1> _yWorkspace

ndarray::Array<double, 1, 1> _dxWorkspace

ndarray::Array<double, 1, 1> _dyWorkspace

class PackedIndex

#include <GaussHermiteEvaluator.h>

An iterator-like object to help in traversing “packed” shapelet or Hermite polynomial matrix or vector dimensions.

A pair of indices (x,y) is mapped to the packed position i = (x+y)(x+y+1)/2 + x.

Typical usage is in a nested loop of the form:

for (PackedIndex i; i.getOrder() <= order; ++i) {
    // utilize i
}

Public Functions

PackedIndex &operator++()

int const getOrder() const

int const getX() const

int const getY() const

int const getIndex() const

PackedIndex()

PackedIndex(int const x, int const y)

Public Static Functions

static int const computeOffset(int order)

static int const computeIndex(int x, int y)

Private Members

int _n

int _i

int _x

int _y