LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
Public Member Functions | Public Attributes | List of all members
ast::QuadApprox Class Reference

A quadratic approximation to a 2D Mapping. More...

#include <QuadApprox.h>

Public Member Functions

 QuadApprox (Mapping const &map, std::vector< double > const &lbnd, std::vector< double > const &ubnd, int nx=3, int ny=3)
 Obtain a quadratic approximation to a 2D Mapping. More...
 
 QuadApprox (QuadApprox const &)=default
 
 QuadApprox (QuadApprox &&)=default
 
QuadApproxoperator= (QuadApprox const &)=default
 
QuadApproxoperator= (QuadApprox &&)=default
 

Public Attributes

std::vector< double > fit
 A vector of coefficients of the quadratic approximation to the specified transformation. More...
 
double rms
 The RMS residual between the fit and the Mapping, summed over all Mapping outputs. More...
 

Detailed Description

A quadratic approximation to a 2D Mapping.

Construct the class to compute the contained fields:

Definition at line 36 of file QuadApprox.h.

Constructor & Destructor Documentation

◆ QuadApprox() [1/3]

ast::QuadApprox::QuadApprox ( Mapping const &  map,
std::vector< double > const &  lbnd,
std::vector< double > const &  ubnd,
int  nx = 3,
int  ny = 3 
)
explicit

Obtain a quadratic approximation to a 2D Mapping.

Compute the coefficients of a quadratic fit to the supplied Mapping over the input area specified by lbnd and ubnd. The Mapping must have 2 inputs, but may have any number of outputs. The i'th Mapping output is modeled as a quadratic function of the 2 inputs (x,y):

output_i = a_i_0 + a_i_1*x + a_i_2*y + a_i_3*x*y + a_i_4*x*x + a_i_5*y*y

The fit vector is set tothe values of the co-efficients a_0_0, a_0_1, etc.

Parameters
[in]mapMapping to fit.
[in]lbndThe lower bounds of a box defined within the input coordinate system of the Mapping. The number of elements in this vector should equal map.getNIn(). This box should specify the region over which the fit is to be performed.
[in]ubndThe upper bounds of the box specifying the region over which the fit is to be performed.
[in]nxThe number of points to place along the first Mapping input. The first point is at lbnd[0] and the last is at ubnd[0]". If a value less than three is supplied a value of three will be used.
[in]nyThe number of points to place along the second Mapping input. The first point is at lbnd[1] and the last is at ubnd[1]. If a value less than three is supplied a value of three will be used.
Exceptions
std::invalid_argumentif the mapping does not have 2 inputs, or if lbnd or ubnd do not each contain 2 elements.
std::runtime_errorif the fit cannot be computed.

Definition at line 31 of file QuadApprox.cc.

33  : fit(6 * map.getNOut()), rms(0) {
34  int const nIn = map.getNIn();
35  detail::assertEqual(nIn, "map.getNIn()", 2, "required nIn");
36  detail::assertEqual(lbnd.size(), "lbnd.size", static_cast<std::size_t>(nIn), "nIn");
37  detail::assertEqual(ubnd.size(), "ubnd.size", static_cast<std::size_t>(nIn), "nIn");
38  fit.reserve(6 * map.getNOut());
39  bool isok = astQuadApprox(map.getRawPtr(), lbnd.data(), ubnd.data(), nx, ny, fit.data(), &rms);
40  assertOK();
41  if (!isok) {
42  throw std::runtime_error("Failed to fit a quadratic approximation");
43  }
44 }
double rms
The RMS residual between the fit and the Mapping, summed over all Mapping outputs.
Definition: QuadApprox.h:102
std::vector< double > fit
A vector of coefficients of the quadratic approximation to the specified transformation.
Definition: QuadApprox.h:98
T data(T... args)
void assertEqual(T1 val1, std::string const &descr1, T2 val2, std::string const &descr2)
Definition: utils.h:48
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition: base.cc:49
T reserve(T... args)
T size(T... args)

◆ QuadApprox() [2/3]

ast::QuadApprox::QuadApprox ( QuadApprox const &  )
default

◆ QuadApprox() [3/3]

ast::QuadApprox::QuadApprox ( QuadApprox &&  )
default

Member Function Documentation

◆ operator=() [1/2]

QuadApprox& ast::QuadApprox::operator= ( QuadApprox &&  )
default

◆ operator=() [2/2]

QuadApprox& ast::QuadApprox::operator= ( QuadApprox const &  )
default

Member Data Documentation

◆ fit

std::vector<double> ast::QuadApprox::fit

A vector of coefficients of the quadratic approximation to the specified transformation.

This vector will contain "6*NOut", elements: the first 6 elements hold the fit to the first Mapping output, the next 6 elements hold the fit to the second Mapping output, etc. So if the Mapping has 2 inputs and 2 outputs the quadratic approximation to the forward transformation is:

X_out = fit[0] + fit[1]*X_in + fit[2]*Y_in + fit[3]*X_in*Y_in +
        fit[4]*X_in*X_in + fit[5]*Y_in*Y_in
Y_out = fit[6] + fit[7]*X_in + fit[8]*Y_in + fit[9]*X_in*Y_in +
        fit[10]*X_in*X_in + fit[11]*Y_in*Y_in
X_out = fit(1) + fit(2)*X_in + fit(3)*Y_in + fit(4)*X_in*Y_in +
        fit(5)*X_in*X_in + fit(6)*Y_in*Y_in
Y_out = fit(7) + fit(8)*X_in + fit(9)*Y_in + fit(10)*X_in*Y_in +
        fit(11)*X_in*X_in + fit(12)*Y_in*Y_in

Definition at line 98 of file QuadApprox.h.

◆ rms

double ast::QuadApprox::rms

The RMS residual between the fit and the Mapping, summed over all Mapping outputs.

Definition at line 102 of file QuadApprox.h.


The documentation for this class was generated from the following files: