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 | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | List of all members
lsst::afw::math::PolynomialFunction2< ReturnT > Class Template Reference

2-dimensional polynomial function with cross terms More...

#include <FunctionLibrary.h>

Inheritance diagram for lsst::afw::math::PolynomialFunction2< ReturnT >:
lsst::afw::math::BasePolynomialFunction2< ReturnT > lsst::afw::math::Function2< ReturnT > lsst::afw::table::io::PersistableFacade< Function2< ReturnT > > lsst::afw::math::Function< ReturnT > lsst::afw::table::io::PersistableFacade< Function< ReturnT > > lsst::afw::table::io::Persistable

Public Member Functions

 PolynomialFunction2 (unsigned int order)
 Construct a polynomial function of specified order. More...
 
 PolynomialFunction2 (std::vector< double > params)
 Construct a polynomial function with specified parameters. More...
 
 PolynomialFunction2 (PolynomialFunction2 const &)=default
 
 PolynomialFunction2 (PolynomialFunction2 &&)=default
 
PolynomialFunction2operator= (PolynomialFunction2 const &)=default
 
PolynomialFunction2operator= (PolynomialFunction2 &&)=default
 
 ~PolynomialFunction2 () noexcept override=default
 
std::shared_ptr< Function2< ReturnT > > clone () const override
 Return a pointer to a deep copy of this function. More...
 
ReturnT operator() (double x, double y) const noexcept(IS_NOTHROW_INIT< ReturnT >) override
 
std::vector< double > getDFuncDParameters (double x, double y) const override
 Return the coefficients of the Function's parameters, evaluated at (x, y) I.e. More...
 
std::string toString (std::string const &prefix) const override
 Return a string representation of the function. More...
 
bool isPersistable () const noexcept override
 Return true if this particular object can be persisted using afw::table::io. More...
 
int getOrder () const noexcept
 Get the polynomial order. More...
 
bool isLinearCombination () const noexcept override
 Is the function a linear combination of its parameters? More...
 
unsigned int getNParameters () const noexcept
 Return the number of function parameters. More...
 
virtual double getParameter (unsigned int ind) const
 Get one function parameter without range checking. More...
 
std::vector< double > const & getParameters () const noexcept
 Return all function parameters. More...
 
void setParameter (unsigned int ind, double newValue)
 Set one function parameter without range checking. More...
 
void setParameters (std::vector< double > const &params)
 Set all function parameters. More...
 
void writeFits (std::string const &fileName, std::string const &mode="w") const
 Write the object to a regular FITS file. More...
 
void writeFits (fits::MemFileManager &manager, std::string const &mode="w") const
 Write the object to a FITS image in memory. More...
 
void writeFits (fits::Fits &fitsfile) const
 Write the object to an already-open FITS object. More...
 

Static Public Member Functions

static int nParametersFromOrder (int order)
 Compute number of parameters from polynomial order. More...
 
static int orderFromNParameters (int nParameters)
 Compute polynomial order from the number of parameters. More...
 
static std::shared_ptr< Function2< ReturnT > > readFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< Function2< ReturnT > > readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< Function2< ReturnT > > readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< Function< ReturnT > > readFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< Function< ReturnT > > readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< Function< ReturnT > > readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< Function2< ReturnT > > dynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 
static std::shared_ptr< Function< ReturnT > > dynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 

Protected Types

using OutputArchiveHandle = io::OutputArchiveHandle
 

Protected Member Functions

std::string getPersistenceName () const override
 Return the unique name used to persist this object and look up its factory. More...
 
void write (afw::table::io::OutputArchiveHandle &handle) const override
 Write the object to one or more catalogs. More...
 
 PolynomialFunction2 ()
 
std::string getPythonModule () const override
 Return the fully-qualified Python module that should be imported to guarantee that its factory is registered. More...
 

Protected Attributes

int _order
 order of polynomial More...
 
std::vector< double > _params
 
bool _isCacheValid
 

Detailed Description

template<typename ReturnT>
class lsst::afw::math::PolynomialFunction2< ReturnT >

2-dimensional polynomial function with cross terms

f(x,y) = c0 (0th order)

Intermediate products for the most recent y are cached, so when computing for a set of x, y it is more efficient to change x before you change y.

Definition at line 449 of file FunctionLibrary.h.

Member Typedef Documentation

◆ OutputArchiveHandle

using lsst::afw::table::io::Persistable::OutputArchiveHandle = io::OutputArchiveHandle
protectedinherited

Definition at line 108 of file Persistable.h.

Constructor & Destructor Documentation

◆ PolynomialFunction2() [1/5]

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::PolynomialFunction2 ( unsigned int  order)
inlineexplicit

Construct a polynomial function of specified order.

The polynomial will have (order + 1) * (order + 2) / 2 coefficients

The parameters are initialized to zero.

Parameters
orderorder of polynomial (0 for constant)

Definition at line 458 of file FunctionLibrary.h.

459  : BasePolynomialFunction2<ReturnT>(order), _oldY(0), _xCoeffs(this->_order + 1) {}
table::Key< int > order

◆ PolynomialFunction2() [2/5]

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::PolynomialFunction2 ( std::vector< double >  params)
inlineexplicit

Construct a polynomial function with specified parameters.

The order of the polynomial is determined from the length of the params vector: order = (sqrt(1 + 8 * length) - 3) / 2 and if this is not an integer then the length is unsuitable

Exceptions
lsst::pex::exceptions::InvalidParameterErrorif params length is unsuitable
lsst::pex::exceptions::Exceptionif an internal sanity check fails
Parameters
paramspolynomial coefficients (const, x, y, x^2, xy, y^2...); length must be one of 1, 3, 6, 10, 15...

Definition at line 471 of file FunctionLibrary.h.

474  : BasePolynomialFunction2<ReturnT>(params), _oldY(0), _xCoeffs(this->_order + 1) {}

◆ PolynomialFunction2() [3/5]

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::PolynomialFunction2 ( PolynomialFunction2< ReturnT > const &  )
default

◆ PolynomialFunction2() [4/5]

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::PolynomialFunction2 ( PolynomialFunction2< ReturnT > &&  )
default

◆ ~PolynomialFunction2()

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::~PolynomialFunction2 ( )
overridedefaultnoexcept

◆ PolynomialFunction2() [5/5]

template<typename ReturnT >
lsst::afw::math::PolynomialFunction2< ReturnT >::PolynomialFunction2 ( )
inlineexplicitprotected

Definition at line 562 of file FunctionLibrary.h.

562 : BasePolynomialFunction2<ReturnT>(), _oldY(0), _xCoeffs(0) {}

Member Function Documentation

◆ clone()

template<typename ReturnT >
std::shared_ptr<Function2<ReturnT> > lsst::afw::math::PolynomialFunction2< ReturnT >::clone ( ) const
inlineoverridevirtual

Return a pointer to a deep copy of this function.

This function exists instead of a copy constructor so one can obtain a copy of an actual function instead of a useless copy of the base class.

Every non-virtual function must override this method.

Returns
a pointer to a deep copy of the function

Implements lsst::afw::math::Function2< ReturnT >.

Definition at line 482 of file FunctionLibrary.h.

◆ dynamicCast() [1/2]

template std::shared_ptr< meas::extensions::psfex::PsfexPsf > lsst::afw::table::io::PersistableFacade< Function2< ReturnT > >::dynamicCast ( std::shared_ptr< Persistable > const &  ptr)
staticinherited

Dynamically cast a shared_ptr.

Dynamically cast a shared pointer and raise on failure.

You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;

Exceptions
lsst::pex::exceptions::LogicErrorif the cast fails

param[in] ptr The pointer to be cast.

Returns
The cast pointer.
Exceptions
lsst::pex::exceptions::TypeErrorIf the dynamic cast fails.

Definition at line 218 of file Persistable.cc.

18  {
19  auto result = std::dynamic_pointer_cast<T>(ptr);
20  if (!result) {
21  throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
22  }
23  return result;
24 }
py::object result
Definition: _schema.cc:429
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
uint64_t * ptr
Definition: RangeSet.cc:88

◆ dynamicCast() [2/2]

template std::shared_ptr< meas::extensions::psfex::PsfexPsf > lsst::afw::table::io::PersistableFacade< Function< ReturnT > >::dynamicCast ( std::shared_ptr< Persistable > const &  ptr)
staticinherited

Dynamically cast a shared_ptr.

Dynamically cast a shared pointer and raise on failure.

You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;

Exceptions
lsst::pex::exceptions::LogicErrorif the cast fails

param[in] ptr The pointer to be cast.

Returns
The cast pointer.
Exceptions
lsst::pex::exceptions::TypeErrorIf the dynamic cast fails.

Definition at line 218 of file Persistable.cc.

18  {
19  auto result = std::dynamic_pointer_cast<T>(ptr);
20  if (!result) {
21  throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
22  }
23  return result;
24 }

◆ getDFuncDParameters()

template<typename ReturnT >
std::vector< double > lsst::afw::math::PolynomialFunction2< ReturnT >::getDFuncDParameters ( double  x,
double  y 
) const
overridevirtual

Return the coefficients of the Function's parameters, evaluated at (x, y) I.e.

given c0, c1, c2, c3 ... return 1, x, y, x^2 ...

Reimplemented from lsst::afw::math::BasePolynomialFunction2< ReturnT >.

Definition at line 34 of file Function.cc.

34  {
35  std::vector<double> coeffs(this->getNParameters());
36 
37  //
38  // Go through params order by order, evaluating x^r y^s; we do this by first evaluating
39  // y^s for a complete order, then going through again multiplying by x^r
40  //
41  int i0 = 0; // starting index for this order's coefficients
42  for (int order = 0; order <= this->_order; ++order) {
43  coeffs[i0] = 1;
44  double zn = y; // y^s
45  for (int i = 1; i <= order; ++i) {
46  coeffs[i0 + i] = zn;
47  zn *= y;
48  }
49 
50  zn = x; // x^r
51  for (int i = order - 1; i >= 0; --i) {
52  coeffs[i0 + i] *= zn;
53  zn *= x;
54  }
55 
56  i0 += order + 1;
57  }
58 
59  assert(i0 == static_cast<int>(coeffs.size()));
60 
61  return coeffs;
62 }
double x
int y
Definition: SpanSet.cc:48
unsigned int getNParameters() const noexcept
Return the number of function parameters.
Definition: Function.h:112

◆ getNParameters()

template<typename ReturnT >
unsigned int lsst::afw::math::Function< ReturnT >::getNParameters ( ) const
inlinenoexceptinherited

Return the number of function parameters.

Returns
the number of function parameters

Definition at line 112 of file Function.h.

112 { return _params.size(); }
T size(T... args)

◆ getOrder()

template<typename ReturnT >
int lsst::afw::math::BasePolynomialFunction2< ReturnT >::getOrder ( ) const
inlinenoexceptinherited

Get the polynomial order.

Definition at line 360 of file Function.h.

360 { return _order; }

◆ getParameter()

template<typename ReturnT >
virtual double lsst::afw::math::Function< ReturnT >::getParameter ( unsigned int  ind) const
inlinevirtualinherited

Get one function parameter without range checking.

Returns
the specified function parameter
Parameters
indindex of parameter

Definition at line 119 of file Function.h.

120  {
121  return _params[ind];
122  }

◆ getParameters()

template<typename ReturnT >
std::vector<double> const& lsst::afw::math::Function< ReturnT >::getParameters ( ) const
inlinenoexceptinherited

Return all function parameters.

Returns
the function parameters as a vector

Definition at line 129 of file Function.h.

129 { return _params; }

◆ getPersistenceName()

template<typename ReturnT >
std::string lsst::afw::math::PolynomialFunction2< ReturnT >::getPersistenceName
overrideprotectedvirtual

Return the unique name used to persist this object and look up its factory.

Must be less than ArchiveIndexSchema::MAX_NAME_LENGTH characters.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 207 of file FunctionLibrary.cc.

207  {
208  return "PolynomialFunction2" + Suffix<ReturnT>::get();
209 }

◆ getPythonModule()

template<typename ReturnT >
std::string lsst::afw::math::Function< ReturnT >::getPythonModule ( ) const
inlineoverrideprotectedvirtualinherited

Return the fully-qualified Python module that should be imported to guarantee that its factory is registered.

Must be less than ArchiveIndexSchema::MAX_MODULE_LENGTH characters.

Will be ignored if empty.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 188 of file Function.h.

188 { return "lsst.afw.math"; }

◆ isLinearCombination()

template<typename ReturnT >
bool lsst::afw::math::BasePolynomialFunction2< ReturnT >::isLinearCombination ( ) const
inlineoverridevirtualnoexceptinherited

Is the function a linear combination of its parameters?

Returns
true if the function can be expressed as: sum over i of parameter_i * function_i(args)
Warning
: subclasses must override if true.

Reimplemented from lsst::afw::math::Function< ReturnT >.

Definition at line 362 of file Function.h.

362 { return true; }

◆ isPersistable()

template<typename ReturnT >
bool lsst::afw::math::PolynomialFunction2< ReturnT >::isPersistable ( ) const
inlineoverridevirtualnoexcept

Return true if this particular object can be persisted using afw::table::io.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 549 of file FunctionLibrary.h.

549 { return true; }

◆ nParametersFromOrder()

template<typename ReturnT >
static int lsst::afw::math::BasePolynomialFunction2< ReturnT >::nParametersFromOrder ( int  order)
inlinestaticinherited

Compute number of parameters from polynomial order.

Exceptions
lsst::pex::exceptions::InvalidParameterErrorif order < 0

Definition at line 369 of file Function.h.

369  {
370  if (order < 0) {
372  os << "order=" << order << " invalid: must be >= 0";
374  }
375  return (order + 1) * (order + 2) / 2;
376  }
std::ostream * os
Definition: Schema.cc:557
Reports invalid arguments.
Definition: Runtime.h:66

◆ operator()()

template<typename ReturnT >
ReturnT lsst::afw::math::PolynomialFunction2< ReturnT >::operator() ( double  x,
double  y 
) const
inlineoverridevirtualnoexcept

Implements lsst::afw::math::Function2< ReturnT >.

Definition at line 486 of file FunctionLibrary.h.

486  {
487  /* Solve as follows:
488  - f(x,y) = Cx0 + Cx1 x + Cx2 x^2 + Cx3 x^3 + ...
489  where:
490  Cx0 = P0 + P2 y + P5 y^2 + P9 y^3 + ...
491  Cx1 = P1 + P4 y + P8 y2 + ...
492  Cx2 = P3 + P7 y + ...
493  Cx3 = P6 + ...
494  ...
495 
496  Compute Cx0, Cx1...Cxn by solving 1-d polynomials in y in the usual way.
497  These values are cached and only recomputed for new values of Y or if the parameters change.
498 
499  Then compute f(x,y) by solving the 1-d polynomial in x in the usual way.
500  */
501  const int maxXCoeffInd = this->_order;
502 
503  if ((y != _oldY) || !this->_isCacheValid) {
504  // update _xCoeffs cache
505  // note: paramInd is decremented in both of the following loops
506  int paramInd = static_cast<int>(this->_params.size()) - 1;
507 
508  // initialize _xCoeffs to coeffs for pure y^n; e.g. for 3rd order:
509  // _xCoeffs[0] = _params[9], _xCoeffs[1] = _params[8], ... _xCoeffs[3] = _params[6]
510  for (int xCoeffInd = 0; xCoeffInd <= maxXCoeffInd; ++xCoeffInd, --paramInd) {
511  _xCoeffs[xCoeffInd] = this->_params[paramInd];
512  }
513 
514  // finish computing _xCoeffs
515  for (int xCoeffInd = 0, endXCoeffInd = maxXCoeffInd; paramInd >= 0; --paramInd) {
516  _xCoeffs[xCoeffInd] = (_xCoeffs[xCoeffInd] * y) + this->_params[paramInd];
517  ++xCoeffInd;
518  if (xCoeffInd >= endXCoeffInd) {
519  xCoeffInd = 0;
520  --endXCoeffInd;
521  }
522  }
523 
524  _oldY = y;
525  this->_isCacheValid = true;
526  }
527 
528  // use _xCoeffs to compute result
529  double retVal = _xCoeffs[maxXCoeffInd];
530  for (int xCoeffInd = maxXCoeffInd - 1; xCoeffInd >= 0; --xCoeffInd) {
531  retVal = (retVal * x) + _xCoeffs[xCoeffInd];
532  }
533  return static_cast<ReturnT>(retVal);
534  }

◆ operator=() [1/2]

template<typename ReturnT >
PolynomialFunction2& lsst::afw::math::PolynomialFunction2< ReturnT >::operator= ( PolynomialFunction2< ReturnT > &&  )
default

◆ operator=() [2/2]

template<typename ReturnT >
PolynomialFunction2& lsst::afw::math::PolynomialFunction2< ReturnT >::operator= ( PolynomialFunction2< ReturnT > const &  )
default

◆ orderFromNParameters()

template<typename ReturnT >
static int lsst::afw::math::BasePolynomialFunction2< ReturnT >::orderFromNParameters ( int  nParameters)
inlinestaticinherited

Compute polynomial order from the number of parameters.

Only certain values of nParameters are acceptable, including: nParameters order 1 0 3 1 6 2 10 3 15 4 ...

Exceptions
lsst::pex::exceptions::InvalidParameterErrorif nParameters is invalid

Definition at line 392 of file Function.h.

392  {
393  int order = static_cast<int>(
394  0.5 + ((-3.0 + (std::sqrt(1.0 + (8.0 * static_cast<double>(nParameters))))) / 2.0));
397  os << "nParameters=" << nParameters << " invalid: order is not an integer";
399  }
400  return order;
401  }
static int nParametersFromOrder(int order)
Compute number of parameters from polynomial order.
Definition: Function.h:369
T sqrt(T... args)

◆ readFits() [1/6]

static std::shared_ptr<Function< ReturnT > > lsst::afw::table::io::PersistableFacade< Function< ReturnT > >::readFits ( fits::Fits fitsfile)
inlinestaticinherited

Read an object from an already open FITS object.

Parameters
[in]fitsfileFITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }
static std::shared_ptr< Function< ReturnT > > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [2/6]

static std::shared_ptr<Function2< ReturnT > > lsst::afw::table::io::PersistableFacade< Function2< ReturnT > >::readFits ( fits::Fits fitsfile)
inlinestaticinherited

Read an object from an already open FITS object.

Parameters
[in]fitsfileFITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }

◆ readFits() [3/6]

static std::shared_ptr<Function2< ReturnT > > lsst::afw::table::io::PersistableFacade< Function2< ReturnT > >::readFits ( fits::MemFileManager manager,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a FITS file in memory.

Parameters
[in]managerManager for the memory to read from.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 205 of file Persistable.h.

205  {
206  return dynamicCast(Persistable::_readFits(manager, hdu));
207  }

◆ readFits() [4/6]

static std::shared_ptr<Function< ReturnT > > lsst::afw::table::io::PersistableFacade< Function< ReturnT > >::readFits ( fits::MemFileManager manager,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a FITS file in memory.

Parameters
[in]managerManager for the memory to read from.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 205 of file Persistable.h.

205  {
206  return dynamicCast(Persistable::_readFits(manager, hdu));
207  }

◆ readFits() [5/6]

static std::shared_ptr<Function2< ReturnT > > lsst::afw::table::io::PersistableFacade< Function2< ReturnT > >::readFits ( std::string const &  fileName,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a regular FITS file.

Parameters
[in]fileNameName of the file to read.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 194 of file Persistable.h.

194  {
195  return dynamicCast(Persistable::_readFits(fileName, hdu));
196  }

◆ readFits() [6/6]

static std::shared_ptr<Function< ReturnT > > lsst::afw::table::io::PersistableFacade< Function< ReturnT > >::readFits ( std::string const &  fileName,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a regular FITS file.

Parameters
[in]fileNameName of the file to read.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 194 of file Persistable.h.

194  {
195  return dynamicCast(Persistable::_readFits(fileName, hdu));
196  }

◆ setParameter()

template<typename ReturnT >
void lsst::afw::math::Function< ReturnT >::setParameter ( unsigned int  ind,
double  newValue 
)
inlineinherited

Set one function parameter without range checking.

Parameters
indindex of parameter
newValuenew value for parameter

Definition at line 143 of file Function.h.

145  {
146  _isCacheValid = false;
147  _params[ind] = newValue;
148  }

◆ setParameters()

template<typename ReturnT >
void lsst::afw::math::Function< ReturnT >::setParameters ( std::vector< double > const &  params)
inlineinherited

Set all function parameters.

Exceptions
lsst::pex::exceptions::InvalidParameterErrorif the wrong number of parameters is supplied.
Parameters
paramsvector of function parameters

Definition at line 156 of file Function.h.

157  {
158  if (_params.size() != params.size()) {
159  throw LSST_EXCEPT(
161  (boost::format("params has %d entries instead of %d") % params.size() % _params.size())
162  .str());
163  }
164  _isCacheValid = false;
165  _params = params;
166  }
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174

◆ toString()

template<typename ReturnT >
std::string lsst::afw::math::PolynomialFunction2< ReturnT >::toString ( std::string const &  ) const
inlineoverridevirtual

Return a string representation of the function.

Returns
a string representation of the function

Reimplemented from lsst::afw::math::Function2< ReturnT >.

Definition at line 542 of file FunctionLibrary.h.

542  {
544  os << "PolynomialFunction2 [" << this->_order << "]: ";
545  os << Function2<ReturnT>::toString(prefix);
546  return os.str();
547  }
std::string prefix
Definition: SchemaMapper.cc:72

◆ write()

template<typename ReturnT >
void lsst::afw::math::PolynomialFunction2< ReturnT >::write ( afw::table::io::OutputArchiveHandle handle) const
overrideprotectedvirtual

Write the object to one or more catalogs.

The handle object passed to this function provides an interface for adding new catalogs and adding nested objects to the same archive (while checking for duplicates). See OutputArchiveHandle for more information.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 239 of file FunctionLibrary.cc.

239  {
240  PolynomialFunction2PersistenceHelper const keys(this->getNParameters());
241  table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
242  keys.coefficients.assignVector(*catalog.addNew(), this->getParameters());
243  handle.saveCatalog(catalog);
244 }
CatalogT< BaseRecord > BaseCatalog
Definition: fwd.h:72

◆ writeFits() [1/3]

void lsst::afw::table::io::Persistable::writeFits ( fits::Fits fitsfile) const
inherited

Write the object to an already-open FITS object.

Parameters
[in]fitsfileOpen FITS object to write to.

Definition at line 18 of file Persistable.cc.

18  {
19  OutputArchive archive;
20  archive.put(this);
21  archive.writeFits(fitsfile);
22 }

◆ writeFits() [2/3]

void lsst::afw::table::io::Persistable::writeFits ( fits::MemFileManager manager,
std::string const &  mode = "w" 
) const
inherited

Write the object to a FITS image in memory.

Parameters
[in]managerName of the file to write to.
[in]modeIf "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file.

Definition at line 29 of file Persistable.cc.

29  {
30  fits::Fits fitsfile(manager, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
31  writeFits(fitsfile);
32 }
void writeFits(std::string const &fileName, std::string const &mode="w") const
Write the object to a regular FITS file.
Definition: Persistable.cc:24

◆ writeFits() [3/3]

void lsst::afw::table::io::Persistable::writeFits ( std::string const &  fileName,
std::string const &  mode = "w" 
) const
inherited

Write the object to a regular FITS file.

Parameters
[in]fileNameName of the file to write to.
[in]modeIf "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file.

Definition at line 24 of file Persistable.cc.

24  {
25  fits::Fits fitsfile(fileName, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
26  writeFits(fitsfile);
27 }

Member Data Documentation

◆ _isCacheValid

template<typename ReturnT >
bool lsst::afw::math::Function< ReturnT >::_isCacheValid
mutableprotectedinherited

Definition at line 186 of file Function.h.

◆ _order

template<typename ReturnT >
int lsst::afw::math::BasePolynomialFunction2< ReturnT >::_order
protectedinherited

order of polynomial

Definition at line 431 of file Function.h.

◆ _params

template<typename ReturnT >
std::vector<double> lsst::afw::math::Function< ReturnT >::_params
protectedinherited

Definition at line 185 of file Function.h.


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