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
MapSplit.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 <cmath>
23 #include <limits>
24 #include <sstream>
25 #include <stdexcept>
26 
27 #include "astshim/base.h"
28 #include "astshim/detail/utils.h"
29 #include "astshim/Object.h"
30 #include "astshim/Mapping.h"
31 #include "astshim/MapSplit.h"
32 
33 namespace ast {
34 
36  std::vector<int> locOut;
37  locOut.reserve(map.getNOut()); // the max # of elements astMapSplit may set
38  AstMapping *rawSplitMap;
39  astMapSplit(map.getRawPtr(), in.size(), in.data(), locOut.data(), &rawSplitMap);
40  assertOK();
41  if (!rawSplitMap) {
42  throw std::runtime_error("Could not split the map");
43  }
44  splitMap = Object::fromAstObject<Mapping>(reinterpret_cast<AstObject *>(rawSplitMap), false);
45  origIn = in;
46  // copy the splitMap->getNOut() elements of `locOut` that contain data to `origOut`
47  const int newNOut = splitMap->getNOut();
48  for (int i = 0; i < newNOut; ++i) {
49  origOut.push_back(locOut[i]);
50  }
51 }
52 
53 } // namespace ast
ast::assertOK
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition: base.cc:49
std::vector::reserve
T reserve(T... args)
std::vector< int >
std::vector::size
T size(T... args)
ast::Mapping
An abstract base class for objects which transform one set of coordinates to another.
Definition: Mapping.h:59
ast::Mapping::getNOut
int getNOut() const
Get NOut: the number of output axes.
Definition: Mapping.h:82
std::vector::push_back
T push_back(T... args)
ast::MapSplit::origIn
std::vector< int > origIn
Indices of the inputs of the original mapping were picked for the split mapping.
Definition: MapSplit.h:73
MapSplit.h
ast::MapSplit::splitMap
std::shared_ptr< Mapping > splitMap
The Mapping that was split off.
Definition: MapSplit.h:67
std::runtime_error
STL class.
utils.h
base.h
ast
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
Mapping.h
ast::Object::getRawPtr
AstObject const * getRawPtr() const
Get the raw AST pointer.
Definition: Object.h:292
ast::MapSplit::origOut
std::vector< int > origOut
Indices of the outputs of the original mapping which are fed by the picked inputs.
Definition: MapSplit.h:83
ast::MapSplit::MapSplit
MapSplit(Mapping const &map, std::vector< int > const &in)
Construct a MapSplit.
Definition: MapSplit.cc:35
Object.h
std::vector::data
T data(T... args)