lsst.astshim g45278abf61+052faf71bd
PolyMap.h
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2017 AURA/LSST.
4 *
5 * This product includes software developed by the
6 * LSST Project (http://www.lsst.org/).
7 *
8 * This program is free software: you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation, either version 3 of the License, or
11 * (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the LSST License Statement and
19 * the GNU General Public License along with this program. If not,
20 * see <https://www.lsstcorp.org/LegalNotices/>.
21 */
22#ifndef ASTSHIM_POLYMAP_H
23#define ASTSHIM_POLYMAP_H
24
25#include <algorithm>
26#include <memory>
27#include <vector>
28
29#include "astshim/base.h"
30#include "astshim/Mapping.h"
31
32namespace ast {
33
49class PolyMap : public Mapping {
50 friend class Object;
51
52public:
99 explicit PolyMap(ConstArray2D const &coeff_f, ConstArray2D const &coeff_i,
100 std::string const &options = "IterInverse=0")
101 : Mapping(reinterpret_cast<AstMapping *>(_makeRawPolyMap(coeff_f, coeff_i, options))) {}
102
126 explicit PolyMap(ConstArray2D const &coeff_f, int nout, std::string const &options = "IterInverse=0")
127 : Mapping(reinterpret_cast<AstMapping *>(_makeRawPolyMap(coeff_f, nout, options))) {}
128
129 virtual ~PolyMap() {}
130
132 PolyMap(PolyMap const &) = default;
133 PolyMap(PolyMap &&) = default;
134 PolyMap &operator=(PolyMap const &) = delete;
135 PolyMap &operator=(PolyMap &&) = default;
136
138 std::shared_ptr<PolyMap> copy() const { return std::static_pointer_cast<PolyMap>(copyPolymorphic()); }
139
141 bool getIterInverse() const { return getB("IterInverse"); }
142
144 int getNIterInverse() const { return getI("NIterInverse"); }
145
147 double getTolInverse() const { return getD("TolInverse"); }
148
232 PolyMap polyTran(bool forward, double acc, double maxacc, int maxorder, std::vector<double> const &lbnd,
233 std::vector<double> const &ubnd) const;
234
235protected:
236 virtual std::shared_ptr<Object> copyPolymorphic() const override {
237 return copyImpl<PolyMap, AstPolyMap>();
238 }
239
241 PolyMap(AstPolyMap *map);
242
243private:
245 AstPolyMap *_makeRawPolyMap(ConstArray2D const &coeff_f, ConstArray2D const &coeff_i,
246 std::string const &options = "") const;
247
249 AstPolyMap *_makeRawPolyMap(ConstArray2D const &coeff_f, int nout, std::string const &options = "") const;
250};
251
252} // namespace ast
253
254#endif
Definition: Mapping.h:59
Definition: Object.h:51
double getD(std::string const &attrib) const
Definition: Object.h:374
bool getB(std::string const &attrib) const
Definition: Object.h:348
int getI(std::string const &attrib) const
Definition: Object.h:400
Definition: PolyMap.h:49
PolyMap(PolyMap &&)=default
PolyMap(ConstArray2D const &coeff_f, ConstArray2D const &coeff_i, std::string const &options="IterInverse=0")
Definition: PolyMap.h:99
int getNIterInverse() const
Get NIterInverse: maximum number of iterations for iterative inverse.
Definition: PolyMap.h:144
virtual ~PolyMap()
Definition: PolyMap.h:129
PolyMap(PolyMap const &)=default
Copy constructor: make a deep copy.
bool getIterInverse() const
Get IterInverse: does this provide an iterative inverse transformation?
Definition: PolyMap.h:141
virtual std::shared_ptr< Object > copyPolymorphic() const override
Definition: PolyMap.h:236
double getTolInverse() const
Get TolInverse: target relative error for iterative inverse.
Definition: PolyMap.h:147
PolyMap(ConstArray2D const &coeff_f, int nout, std::string const &options="IterInverse=0")
Definition: PolyMap.h:126
std::shared_ptr< PolyMap > copy() const
Return a deep copy of this object.
Definition: PolyMap.h:138
PolyMap & operator=(PolyMap &&)=default
PolyMap & operator=(PolyMap const &)=delete
PolyMap polyTran(bool forward, double acc, double maxacc, int maxorder, std::vector< double > const &lbnd, std::vector< double > const &ubnd) const
Definition: PolyMap.cc:30
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
ndarray::Array< const double, 2, 2 > ConstArray2D
Definition: base.h:46