LSSTApplications  19.0.0-10-g920eed2,19.0.0-11-g48a0200+2,19.0.0-18-gfc4e62b+13,19.0.0-2-g3b2f90d+2,19.0.0-2-gd671419+5,19.0.0-20-g5a5a17ab+11,19.0.0-21-g2644856+13,19.0.0-23-g84eeccb+1,19.0.0-24-g878c510+1,19.0.0-25-g6c8df7140,19.0.0-25-gb330496+1,19.0.0-3-g2b32d65+5,19.0.0-3-g8227491+12,19.0.0-3-g9c54d0d+12,19.0.0-3-gca68e65+8,19.0.0-3-gcfc5f51+5,19.0.0-3-ge110943+11,19.0.0-3-ge74d124,19.0.0-3-gfe04aa6+13,19.0.0-30-g9c3fd16+1,19.0.0-4-g06f5963+5,19.0.0-4-g3d16501+13,19.0.0-4-g4a9c019+5,19.0.0-4-g5a8b323,19.0.0-4-g66397f0+1,19.0.0-4-g8278b9b+1,19.0.0-4-g8557e14,19.0.0-4-g8964aba+13,19.0.0-4-ge404a01+12,19.0.0-5-g40f3a5a,19.0.0-5-g4db63b3,19.0.0-5-gfb03ce7+13,19.0.0-6-gbaebbfb+12,19.0.0-61-gec4c6e08+1,19.0.0-7-g039c0b5+11,19.0.0-7-gbea9075+4,19.0.0-7-gc567de5+13,19.0.0-71-g41c0270,19.0.0-9-g2f02add+1,19.0.0-9-g463f923+12,w.2020.22
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:78
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:76
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