LSSTApplications  19.0.0-14-gb0260a2+55ebc804df,20.0.0+1b2fec4223,20.0.0+291834023b,20.0.0+3aed93c9a0,20.0.0+445d2e5176,20.0.0+44e6c265c0,20.0.0+5f0d44a82c,20.0.0+dbd2281351,20.0.0+e2e26847c2,20.0.0+f45b7d88f4,20.0.0+f48bb53858,20.0.0-1-g10df615+65b272dc7c,20.0.0-1-g253301a+1b2fec4223,20.0.0-1-g2b7511a+f48bb53858,20.0.0-1-g4d801e7+990eb492d1,20.0.0-1-g5b95a8c+535fc78ebf,20.0.0-1-gc96f8cb+92cff2790a,20.0.0-1-gd1c87d7+add1f556b4,20.0.0-1-gedffbd8+f862251818,20.0.0-10-gd6f3b0e+42e798760b,20.0.0-14-g199c894+1c6164e6a1,20.0.0-15-gf5872d0+92cff2790a,20.0.0-16-g6fe20a1+367edb3a36,20.0.0-2-g4dae9ad+92cff2790a,20.0.0-2-g5ad0983+42e798760b,20.0.0-2-g7818986+add1f556b4,20.0.0-2-gb095acb+63636aeed8,20.0.0-2-gdaeb0e8+e2e179d1a0,20.0.0-2-gec03fae+de92faa26f,20.0.0-22-g629db56c+016cf65f00,20.0.0-4-gfea843c+f45b7d88f4,20.0.0-46-g00fa051+95f8c0b24e,20.0.0-5-g357b56b+f45b7d88f4,20.0.0-5-gfcebe35+5d9781fc8a,20.0.0-7-gc818c6e6+f7ac80044f,20.0.0-7-gcda7bf1+85dc274a4b,20.0.0-9-g61a2a9a3d+63f5489896,w.2020.38
LSSTDataManagementBasePackage
Frame.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 <vector>
24 
25 #include "astshim/CmpFrame.h"
26 #include "astshim/Frame.h"
27 #include "astshim/FrameSet.h"
28 
29 namespace ast {
30 
32  auto *rawFrameSet =
33  reinterpret_cast<AstFrameSet *>(astConvert(getRawPtr(), to.getRawPtr(), domainlist.c_str()));
34  assertOK(reinterpret_cast<AstObject *>(rawFrameSet));
35  if (!rawFrameSet) {
37  }
38  return Object::fromAstObject<FrameSet>(reinterpret_cast<AstObject *>(rawFrameSet), false);
39 }
40 
41 std::shared_ptr<FrameSet> Frame::findFrame(Frame const &tmplt, std::string const &domainlist) {
42  auto *rawFrameSet =
43  reinterpret_cast<AstFrameSet *>(astFindFrame(getRawPtr(), tmplt.getRawPtr(), domainlist.c_str()));
44  assertOK(reinterpret_cast<AstObject *>(rawFrameSet));
45  if (!rawFrameSet) {
47  }
48  return Object::fromAstObject<FrameSet>(reinterpret_cast<AstObject *>(rawFrameSet), false);
49 }
50 
52  std::vector<double> const &b1, std::vector<double> const &b2) const {
53  int const naxes = 2;
54  detail::assertEqual(getNAxes(), "# axes", naxes, "");
55  detail::assertEqual(a1.size(), "a1.size()", static_cast<std::size_t>(naxes), "");
56  detail::assertEqual(a2.size(), "a2.size()", static_cast<std::size_t>(naxes), "");
57  detail::assertEqual(b1.size(), "b1.size()", static_cast<std::size_t>(naxes), "");
58  detail::assertEqual(b2.size(), "b2.size()", static_cast<std::size_t>(naxes), "");
59  std::vector<double> ret(naxes);
60  astIntersect(getRawPtr(), a1.data(), a2.data(), b1.data(), b2.data(), ret.data());
61  assertOK();
63  return ret;
64 }
65 
66 CmpFrame Frame::under(Frame const &next) const { return CmpFrame(*this, next); }
67 
69  AstMapping *rawMap;
70  auto *rawFrame =
71  reinterpret_cast<AstFrame *>(astPickAxes(getRawPtr(), axes.size(), axes.data(), &rawMap));
72  assertOK(reinterpret_cast<AstObject *>(rawFrame), reinterpret_cast<AstObject *>(rawMap));
74  try {
75  frame = Object::fromAstObject<Frame>(reinterpret_cast<AstObject *>(rawFrame), false);
76  } catch (...) {
77  astAnnul(rawMap);
78  throw;
79  }
80  auto map = Object::fromAstObject<Mapping>(reinterpret_cast<AstObject *>(rawMap), false);
81  return FrameMapping(frame, map);
82 }
83 
85  std::vector<double> const &point3) const {
86  int const naxes = getNAxes();
87  detail::assertEqual(point1.size(), "a1.size()", static_cast<std::size_t>(naxes), "");
88  detail::assertEqual(point2.size(), "a2.size()", static_cast<std::size_t>(naxes), "");
89  detail::assertEqual(point3.size(), "b1.size()", static_cast<std::size_t>(naxes), "");
90  ResolvedPoint ret(naxes);
91  astResolve(getRawPtr(), point1.data(), point2.data(), point3.data(), ret.point.data(), &ret.d1, &ret.d2);
92  assertOK();
94  return ret;
95 }
96 
97 } // namespace ast
astshim.fitsChanContinued.next
def next(self)
Definition: fitsChanContinued.py:105
FrameSet.h
ast::FrameMapping
Struct returned by Frame::pickAxes containing a frame and a mapping.
Definition: Frame.h:91
std::string
STL class.
std::shared_ptr
STL class.
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< double >
std::vector::size
T size(T... args)
ast::Frame::convert
std::shared_ptr< FrameSet > convert(Frame const &to, std::string const &domainlist="")
Compute a frameset that describes the conversion between this frame and another frame.
Definition: Frame.cc:31
ast::detail::astBadToNan
void astBadToNan(std::vector< double > &p)
Replace AST__BAD with a quiet NaN in a vector.
Definition: utils.h:59
ast::Frame::getNAxes
int getNAxes() const
Get NAxes: the number of axes in the frame (i.e.
Definition: Frame.h:876
ast::Frame
Frame is used to represent a coordinate system.
Definition: Frame.h:157
ast::ResolvedPoint
Struct returned by Frame::resolve containing a point and the resolved vector components.
Definition: Frame.h:75
ast::Frame::pickAxes
FrameMapping pickAxes(std::vector< int > const &axes) const
Create a new Frame whose axes are copied from an existing Frame along with other Frame attributes,...
Definition: Frame.cc:68
ast::Frame::resolve
ResolvedPoint resolve(std::vector< double > const &point1, std::vector< double > const &point2, std::vector< double > const &point3) const
Resolve a vector into two orthogonal components.
Definition: Frame.cc:84
std::string::c_str
T c_str(T... args)
CmpFrame.h
ast::ResolvedPoint::d2
double d2
Resolved vector component 2.
Definition: Frame.h:85
ast::Frame::under
CmpFrame under(Frame const &next) const
Combine this frame with another to form a compound frame (CmpFrame), with the axes of this frame foll...
Definition: Frame.cc:66
ast::Frame::intersect
std::vector< double > intersect(std::vector< double > const &a1, std::vector< double > const &a2, std::vector< double > const &b1, std::vector< double > const &b2) const
Find the point of intersection between two geodesic curves.
Definition: Frame.cc:51
to
table::Key< int > to
Definition: TransformMap.cc:349
ast
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
ast::ResolvedPoint::point
std::vector< double > point
Point.
Definition: Frame.h:83
ast::Object::getRawPtr
AstObject const * getRawPtr() const
Get the raw AST pointer.
Definition: Object.h:292
std::size_t
Frame.h
ast::ResolvedPoint::d1
double d1
Resolved vector component 1.
Definition: Frame.h:84
ast::CmpFrame
A CmpFrame is a compound Frame which allows two component Frames (of any class) to be merged together...
Definition: CmpFrame.h:60
ast::detail::assertEqual
void assertEqual(T1 val1, std::string const &descr1, T2 val2, std::string const &descr2)
Definition: utils.h:48
std::vector::data
T data(T... args)
ast::Frame::findFrame
std::shared_ptr< FrameSet > findFrame(Frame const &tmplt, std::string const &domainlist="")
Find a coordinate system with specified characteristics.
Definition: Frame.cc:41