22 #ifndef ASTSHIM_FRAME_H
23 #define ASTSHIM_FRAME_H
30 #include "astshim/Mapping.h"
31 #include "astshim/Object.h"
32 #include "astshim/base.h"
33 #include "astshim/detail/utils.h"
100 : frame(frame), mapping(mapping) {}
168 explicit Frame(
int naxes, std::string
const &options =
"")
169 :
Mapping(reinterpret_cast<AstMapping *>(astFrame(naxes,
"%s", options.c_str()))) {}
199 assertPointLength(a,
"a");
200 assertPointLength(b,
"b");
201 assertPointLength(c,
"c");
202 return detail::safeDouble(astAngle(
getRawPtr(), a.data(), b.data(), c.data()));
227 assertPointLength(a,
"a");
228 assertPointLength(b,
"b");
229 return detail::safeDouble(astAxAngle(
getRawPtr(), a.data(), b.data(), axis));
245 return detail::safeDouble(astAxDistance(
getRawPtr(), axis, v1, v2));
260 double axOffset(
int axis,
double v1,
double dist)
const {
261 return detail::safeDouble(astAxOffset(
getRawPtr(), axis, v1, dist));
458 std::shared_ptr<FrameSet>
convert(
Frame const &to, std::string
const &domainlist =
"");
474 assertPointLength(point1,
"point1");
475 assertPointLength(point2,
"point2");
476 return detail::safeDouble(astDistance(
getRawPtr(), point1.data(), point2.data()));
769 std::shared_ptr<FrameSet>
findFrame(
Frame const &tmplt, std::string
const &domainlist =
"");
782 std::string
format(
int axis,
double value)
const {
783 char const *rawstr = astFormat(
getRawPtr(), axis, value);
785 return std::string(rawstr);
793 bool ret = astGetActiveUnit(
getRawPtr());
807 double getBottom(
int axis)
const {
return getD(detail::formatAxisAttr(
"Bottom", axis)); }
817 int getDigits(
int axis)
const {
return getI(detail::formatAxisAttr(
"Digits", axis)); }
822 bool getDirection(
int axis)
const {
return getB(detail::formatAxisAttr(
"Direction", axis)); }
842 std::string
getFormat(
int axis)
const {
return getC(detail::formatAxisAttr(
"Format", axis)); }
848 std::string
getInternalUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"InternalUnit", axis)); }
853 std::string
getLabel(
int axis)
const {
return getC(detail::formatAxisAttr(
"Label", axis)); }
882 std::string
getNormUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"NormUnit", axis)); }
912 std::string
getSymbol(
int axis)
const {
return getC(detail::formatAxisAttr(
"Symbol", axis)); }
928 double getTop(
int axis)
const {
return getD(detail::formatAxisAttr(
"Top", axis)); }
933 std::string
getUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"Unit", axis)); }
966 std::vector<double>
intersect(std::vector<double>
const &a1, std::vector<double>
const &a2,
967 std::vector<double>
const &b1, std::vector<double>
const &b2)
const;
986 std::vector<int> ret(other.
getNIn());
1045 detail::astBadToNan(value);
1073 assertPointLength(point1,
"point1");
1074 assertPointLength(point2,
"point2");
1076 astOffset(
getRawPtr(), point1.data(), point2.data(),
offset, ret.data());
1078 detail::astBadToNan(ret);
1115 detail::assertEqual(
getNIn(),
"naxes", 2,
" cannot call offset2");
1116 assertPointLength(point1,
"point1");
1120 detail::astBadToNan(point2);
1136 detail::assertEqual(perm.size(),
"perm.size()",
static_cast<std::size_t
>(
getNAxes()),
"naxes");
1194 std::vector<double>
const &point3)
const;
1205 void setBottom(
int axis,
double bottom) {
setD(detail::formatAxisAttr(
"Bottom", axis), bottom); }
1215 void setDigits(
int axis,
int digits) {
setD(detail::formatAxisAttr(
"Digits", axis), digits); }
1221 setB(detail::formatAxisAttr(
"Direction", axis), direction);
1248 setC(detail::formatAxisAttr(
"Format", axis), format);
1254 void setLabel(
int axis, std::string
const &label) {
setC(detail::formatAxisAttr(
"Label", axis), label); }
1311 setC(detail::formatAxisAttr(
"Symbol", axis), symbol);
1328 void setTop(
int axis,
double top) {
setD(detail::formatAxisAttr(
"Top", axis), top); }
1333 void setUnit(
int axis, std::string
const &unit) {
setC(detail::formatAxisAttr(
"Unit", axis), unit); }
1501 int nread = astUnformat(
getRawPtr(), axis, str.c_str(), &value);
1503 return NReadValue(nread, detail::safeDouble(value));
1516 explicit Frame(AstFrame *rawPtr) :
Mapping(reinterpret_cast<AstMapping *>(rawPtr)) {
1518 std::ostringstream os;
1519 os <<
"This is a " <<
getClassName() <<
", which is not a Frame";
1520 throw std::invalid_argument(os.str());
1524 virtual std::shared_ptr<Object>
copyPolymorphic()
const override {
return copyImpl<Frame, AstFrame>(); }
1533 template <
typename T>
1534 void assertPointLength(T
const &p,
char const *name)
const {
1535 if (static_cast<int>(p.size()) !=
getNIn()) {
1536 std::ostringstream os;
1537 os <<
"point " << name <<
" has " << p.size() <<
" axes, but " <<
getNIn() <<
" required";
1538 throw std::invalid_argument(os.str());
std::string getInternalUnit(int axis) const
Definition: Frame.h:848
std::shared_ptr< Mapping > mapping
Mapping.
Definition: Frame.h:102
double angle(PointD const &a, PointD const &b, PointD const &c) const
Definition: Frame.h:198
PointD norm(PointD value) const
Definition: Frame.h:1042
FrameMapping(std::shared_ptr< Frame > frame, std::shared_ptr< Mapping > mapping)
Definition: Frame.h:99
std::string getNormUnit(int axis) const
Definition: Frame.h:882
DirectionPoint offset2(PointD const &point1, double angle, double offset) const
Definition: Frame.h:1114
std::string getTitle() const
Definition: Frame.h:923
void setI(std::string const &attrib, int value)
Definition: Object.h:480
int getDigits(int axis) const
Definition: Frame.h:817
int getNAxes() const
Definition: Frame.h:876
double distance(PointD const &point1, PointD const &point2) const
Definition: Frame.h:473
void setSystem(std::string const &system)
Definition: Frame.h:1318
std::string getSymbol(int axis) const
Definition: Frame.h:912
void setTop(int axis, double top)
Definition: Frame.h:1328
std::shared_ptr< FrameSet > convert(Frame const &to, std::string const &domainlist="")
Definition: Frame.cc:31
std::shared_ptr< Frame > copy() const
Return a deep copy of this object.
Definition: Frame.h:179
CmpFrame under(Frame const &next) const
Definition: Frame.cc:66
double axDistance(int axis, double v1, double v2) const
Definition: Frame.h:244
double getObsAlt() const
Definition: Frame.h:887
void setBottom(int axis, double bottom)
Definition: Frame.h:1205
void setMatchEnd(bool match)
Definition: Frame.h:1259
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Definition: base.cc:49
void setObsAlt(double alt)
Definition: Frame.h:1276
std::string format(int axis, double value) const
Definition: Frame.h:782
void setMaxAxes(int maxAxes)
Definition: Frame.h:1265
PointD offset(PointD point1, PointD point2, double offset) const
Definition: Frame.h:1072
Definition: FrameSet.h:99
std::string getObsLat() const
Definition: Frame.h:892
FrameMapping pickAxes(std::vector< int > const &axes) const
Definition: Frame.cc:68
std::string const getC(std::string const &attrib) const
Definition: Object.h:344
void setPermute(bool permute)
Definition: Frame.h:1300
ResolvedPoint resolve(std::vector< double > const &point1, std::vector< double > const &point2, std::vector< double > const &point3) const
Definition: Frame.cc:84
void setB(std::string const &attrib, bool value)
Definition: Object.h:432
void setDigits(int axis, int digits)
Definition: Frame.h:1215
bool getActiveUnit() const
Definition: Frame.h:792
Frame(AstFrame *rawPtr)
Definition: Frame.h:1516
std::vector< int > matchAxes(Frame const &other) const
Definition: Frame.h:985
std::string getLabel(int axis) const
Definition: Frame.h:853
std::shared_ptr< Frame > frame
Frame.
Definition: Frame.h:101
double value
Value that was read.
Definition: Frame.h:69
int getNIn() const
Definition: Mapping.h:76
std::string getAlignSystem() const
Definition: Frame.h:802
void setEpoch(std::string const &epoch)
Definition: Frame.h:1242
void setObsLon(std::string const &lon)
Definition: Frame.h:1286
virtual std::shared_ptr< Object > copyPolymorphic() const override
Definition: Frame.h:1524
void setTitle(std::string const &title)
Definition: Frame.h:1323
bool getPermute() const
Definition: Frame.h:902
void setSymbol(int axis, std::string const &symbol)
Definition: Frame.h:1310
Frame(int naxes, std::string const &options="")
Definition: Frame.h:168
std::string getUnit(int axis) const
Definition: Frame.h:933
void permAxes(std::vector< int > perm)
Definition: Frame.h:1135
void setC(std::string const &attrib, std::string const &value)
Definition: Object.h:444
double d1
Resolved vector component 1.
Definition: Frame.h:84
std::string getObsLon() const
Definition: Frame.h:897
double getBottom(int axis) const
Definition: Frame.h:807
void setActiveUnit(bool enable)
Definition: Frame.h:1292
bool getDirection(int axis) const
Definition: Frame.h:822
std::string getDomain() const
Definition: Frame.h:827
void setFormat(int axis, std::string const &format)
Definition: Frame.h:1247
double axAngle(PointD const &a, PointD const &b, int axis) const
Definition: Frame.h:226
PointD point
Point.
Definition: Frame.h:53
double getD(std::string const &attrib) const
Definition: Object.h:357
std::vector< double > PointD
Definition: base.h:57
double getDut1() const
Definition: Frame.h:832
NReadValue unformat(int axis, std::string const &str) const
Definition: Frame.h:1499
int getMinAxes() const
Definition: Frame.h:870
bool getPreserveAxes() const
Definition: Frame.h:907
void setLabel(int axis, std::string const &label)
Definition: Frame.h:1254
NReadValue(int nread, double value)
Definition: Frame.h:67
double d2
Resolved vector component 2.
Definition: Frame.h:85
void setPreserveAxes(bool preserve)
Definition: Frame.h:1305
void setD(std::string const &attrib, double value)
Definition: Object.h:456
Definition: CmpFrame.h:60
DirectionPoint(double direction, PointD const &point)
Definition: Frame.h:51
void setDigits(int digits)
Definition: Frame.h:1210
void setDomain(std::string const &domain)
Definition: Frame.h:1227
void setEpoch(double epoch)
Definition: Frame.h:1237
std::shared_ptr< FrameSet > findFrame(Frame const &tmplt, std::string const &domainlist="")
Definition: Frame.cc:41
double direction
Direction, an angle in radians.
Definition: Frame.h:51
double getEpoch() const
Definition: Frame.h:837
double axOffset(int axis, double v1, double dist) const
Definition: Frame.h:260
std::string getSystem() const
Definition: Frame.h:918
bool getMatchEnd() const
Definition: Frame.h:858
void setDirection(bool direction, int axis)
Definition: Frame.h:1220
int getI(std::string const &attrib) const
Definition: Object.h:383
void setAlignSystem(std::string const &system)
Definition: Frame.h:1200
std::vector< double > point
Point.
Definition: Frame.h:83
int nread
Number of characters that was read.
Definition: Frame.h:67
bool getB(std::string const &attrib) const
Definition: Object.h:331
void setUnit(int axis, std::string const &unit)
Definition: Frame.h:1333
std::string getFormat(int axis) const
Definition: Frame.h:842
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
Definition: Frame.cc:51
void setMinAxes(int minAxes)
Definition: Frame.h:1271
int getMaxAxes() const
Definition: Frame.h:864
double getTop(int axis) const
Definition: Frame.h:928
int getDigits() const
Definition: Frame.h:812
void setObsLat(std::string const &lat)
Definition: Frame.h:1281
void setDut1(double dut1)
Definition: Frame.h:1232
std::string getClassName() const
Definition: Object.h:118
ResolvedPoint(int naxes)
Definition: Frame.h:82
AstObject const * getRawPtr() const
Definition: Object.h:270