LSSTApplications  19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
LSSTDataManagementBasePackage
functional.cc
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 #include <stdexcept>
23 #include <string>
24 
25 #include "astshim/functional.h"
26 #include "astshim/UnitMap.h"
27 #include "astshim/UnitNormMap.h"
28 #include "astshim/ParallelMap.h"
29 #include "astshim/SeriesMap.h"
30 
31 namespace ast {
32 
34  std::shared_ptr<FrameSet> const merged = first.copy();
35  std::shared_ptr<FrameSet> const newFrames = second.copy();
36 
37  newFrames->setCurrent(FrameSet::BASE);
38  int const joinNAxes = first.getFrame(FrameSet::CURRENT)->getNAxes();
39  merged->addFrame(FrameSet::CURRENT, UnitMap(joinNAxes), *newFrames);
40 
41  // All frame numbers from `second` have been offset in `merged` by number of frames in `first`
42  int const mergedCurrent = first.getNFrame() + second.getCurrent();
43  merged->setCurrent(mergedCurrent);
44 
45  return merged;
46 }
47 
49  auto naxes = center.size();
50  if (mapping1d.getNIn() != 1) {
51  throw std::invalid_argument("mapping1d has " + std::to_string(mapping1d.getNIn()) +
52  " inputs, instead of 1");
53  }
54  if (mapping1d.getNOut() != 1) {
55  throw std::invalid_argument("mapping1d has " + std::to_string(mapping1d.getNOut()) +
56  " outputs, instead of 1");
57  }
58  auto unitNormMap = UnitNormMap(center);
59  return std::make_shared<Mapping>(
60  unitNormMap.then(UnitMap(naxes).under(mapping1d)).then(*unitNormMap.inverted()));
61 }
62 
63 } // namespace ast
ast::FrameSet
A FrameSet consists of a set of one or more Frames (which describe coordinate systems),...
Definition: FrameSet.h:99
std::shared_ptr
STL class.
lsst::afw::table._match.second
second
Definition: _match.py:76
std::vector< double >
std::vector::size
T size(T... args)
UnitNormMap.h
ast::Mapping
An abstract base class for objects which transform one set of coordinates to another.
Definition: Mapping.h:59
SeriesMap.h
ParallelMap.h
lsst::afw::table._match.first
first
Definition: _match.py:74
ast::append
std::shared_ptr< FrameSet > append(FrameSet const &first, FrameSet const &second)
Construct a FrameSet that performs two transformations in series.
Definition: functional.cc:33
UnitMap.h
ast::Mapping::getNOut
int getNOut() const
Get NOut: the number of output axes.
Definition: Mapping.h:82
ast::FrameSet::CURRENT
static constexpr int CURRENT
index of current frame
Definition: FrameSet.h:105
ast::makeRadialMapping
std::shared_ptr< Mapping > makeRadialMapping(std::vector< double > const &center, Mapping const &mapping1d)
Construct a radially symmetric mapping from a 1-dimensional mapping.
Definition: functional.cc:48
std::to_string
T to_string(T... args)
ast::UnitNormMap
The forward transformation of a UnitNormMap subtracts the specified centre and then transforms the re...
Definition: UnitNormMap.h:53
std::invalid_argument
STL class.
functional.h
ast::FrameSet::setCurrent
void setCurrent(int ind)
Set Current: index of current Frame, starting from 1.
Definition: FrameSet.h:476
ast
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
ast::Mapping::getNIn
int getNIn() const
Get NIn: the number of input axes.
Definition: Mapping.h:77
ast::FrameSet::addFrame
virtual void addFrame(int iframe, Mapping const &map, Frame const &frame)
Add a new Frame and an associated Mapping to this FrameSet so as to define a new coordinate system,...
Definition: FrameSet.h:210
ast::FrameSet::BASE
static constexpr int BASE
index of base frame
Definition: FrameSet.h:104
lsst::afw::geom.python.transform.then
then
Definition: transform.py:113
ast::UnitMap
A UnitMap is a unit (null) Mapping that has no effect on the coordinates supplied to it.
Definition: UnitMap.h:44