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" 67 NReadValue(
int nread,
double value) : nread(nread), value(value){};
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()))) {}
179 std::shared_ptr<Frame>
copy()
const {
return std::static_pointer_cast<
Frame>(copyPolymorphic()); }
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));
457 std::shared_ptr<FrameSet> convert(
Frame const &to, std::string
const &domainlist =
"");
473 assertPointLength(point1,
"point1");
474 assertPointLength(point2,
"point2");
475 return detail::safeDouble(astDistance(getRawPtr(), point1.data(), point2.data()));
767 std::shared_ptr<FrameSet> findFrame(
Frame const &tmplt, std::string
const &domainlist =
"");
780 std::string
format(
int axis,
double value)
const {
781 char const *rawstr = astFormat(getRawPtr(), axis, value);
783 return std::string(rawstr);
791 bool ret = astGetActiveUnit(getRawPtr());
805 double getBottom(
int axis)
const {
return getD(detail::formatAxisAttr(
"Bottom", axis)); }
815 int getDigits(
int axis)
const {
return getI(detail::formatAxisAttr(
"Digits", axis)); }
820 bool getDirection(
int axis)
const {
return getB(detail::formatAxisAttr(
"Direction", axis)); }
825 std::string
getDomain()
const {
return getC(
"Domain"); }
830 double getDut1()
const {
return getD(
"Dut1"); }
840 std::string
getFormat(
int axis)
const {
return getC(detail::formatAxisAttr(
"Format", axis)); }
846 std::string
getInternalUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"InternalUnit", axis)); }
851 std::string
getLabel(
int axis)
const {
return getC(detail::formatAxisAttr(
"Label", axis)); }
880 std::string
getNormUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"NormUnit", axis)); }
890 std::string
getObsLat()
const {
return getC(
"ObsLat"); }
895 std::string
getObsLon()
const {
return getC(
"ObsLon"); }
910 std::string
getSymbol(
int axis)
const {
return getC(detail::formatAxisAttr(
"Symbol", axis)); }
916 std::string
getSystem()
const {
return getC(
"System"); }
921 std::string
getTitle()
const {
return getC(
"Title"); }
926 double getTop(
int axis)
const {
return getD(detail::formatAxisAttr(
"Top", axis)); }
931 std::string
getUnit(
int axis)
const {
return getC(detail::formatAxisAttr(
"Unit", axis)); }
964 std::vector<double> intersect(std::vector<double>
const &a1, std::vector<double>
const &a2,
965 std::vector<double>
const &b1, std::vector<double>
const &b2)
const;
984 std::vector<int> ret(other.
getNIn());
985 astMatchAxes(getRawPtr(), other.
getRawPtr(), ret.data());
1044 astNorm(getRawPtr(), value.data());
1046 detail::astBadToNan(value);
1074 assertPointLength(point1,
"point1");
1075 assertPointLength(point2,
"point2");
1077 astOffset(getRawPtr(), point1.data(), point2.data(), offset, ret.data());
1079 detail::astBadToNan(ret);
1116 detail::assertEqual(getNIn(),
"naxes", 2,
" cannot call offset2");
1117 assertPointLength(point1,
"point1");
1119 double offsetAngle = astOffset2(getRawPtr(), point1.data(), angle, offset, point2.data());
1121 detail::astBadToNan(point2);
1137 detail::assertEqual(perm.size(),
"perm.size()",
static_cast<std::size_t
>(getNAxes()),
"naxes");
1138 astPermAxes(getRawPtr(), perm.data());
1161 FrameMapping pickAxes(std::vector<int>
const &axes)
const;
1194 ResolvedPoint resolve(std::vector<double>
const &point1, std::vector<double>
const &point2,
1195 std::vector<double>
const &point3)
const;
1206 void setBottom(
int axis,
double bottom) { setD(detail::formatAxisAttr(
"Bottom", axis), bottom); }
1216 void setDigits(
int axis,
int digits) { setD(detail::formatAxisAttr(
"Digits", axis), digits); }
1222 setB(detail::formatAxisAttr(
"Direction", axis), direction);
1228 void setDomain(std::string
const &domain) { setC(
"Domain", domain); }
1243 void setEpoch(std::string
const &epoch) { setC(
"Epoch", epoch); }
1249 setC(detail::formatAxisAttr(
"Format", axis), format);
1255 void setLabel(
int axis, std::string
const &label) { setC(detail::formatAxisAttr(
"Label", axis), label); }
1282 void setObsLat(std::string
const &lat) { setC(
"ObsLat", lat); }
1287 void setObsLon(std::string
const &lon) { setC(
"ObsLon", lon); }
1294 astSetActiveUnit(getRawPtr(), enable);
1312 setC(detail::formatAxisAttr(
"Symbol", axis), symbol);
1319 void setSystem(std::string
const &system) { setC(
"System", system); }
1324 void setTitle(std::string
const &title) { setC(
"Title", title); }
1329 void setTop(
int axis,
double top) { setD(detail::formatAxisAttr(
"Top", axis), top); }
1334 void setUnit(
int axis, std::string
const &unit) { setC(detail::formatAxisAttr(
"Unit", axis), unit); }
1502 int nread = astUnformat(getRawPtr(), axis, str.c_str(), &value);
1504 return NReadValue(nread, detail::safeDouble(value));
1517 explicit Frame(AstFrame *rawPtr) :
Mapping(reinterpret_cast<AstMapping *>(rawPtr)) {
1518 if (!astIsAFrame(getRawPtr())) {
1519 std::ostringstream os;
1520 os <<
"This is a " << getClassName() <<
", which is not a Frame";
1521 throw std::invalid_argument(os.str());
1525 virtual std::shared_ptr<Object>
copyPolymorphic()
const override {
return copyImpl<Frame, AstFrame>(); }
1534 template <
typename T>
1535 void assertPointLength(T
const &p,
char const *name)
const {
1536 if (static_cast<int>(p.size()) != getNIn()) {
1537 std::ostringstream os;
1538 os <<
"point " << name <<
" has " << p.size() <<
" axes, but " << getNIn() <<
" required";
1539 throw std::invalid_argument(os.str());
double axOffset(int axis, double v1, double dist) const
Definition: Frame.h:260
Definition: CmpFrame.h:60
Frame(int naxes, std::string const &options="")
Definition: Frame.h:168
std::string getDomain() const
Definition: Frame.h:825
void setDigits(int digits)
Definition: Frame.h:1211
void setObsLat(std::string const &lat)
Definition: Frame.h:1282
double getTop(int axis) const
Definition: Frame.h:926
AstObject const * getRawPtr() const
Definition: Object.h:286
double value
Value that was read.
Definition: Frame.h:69
void setEpoch(std::string const &epoch)
Definition: Frame.h:1243
void setPreserveAxes(bool preserve)
Definition: Frame.h:1306
double distance(PointD const &point1, PointD const &point2) const
Definition: Frame.h:472
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
double axDistance(int axis, double v1, double v2) const
Definition: Frame.h:244
void setMaxAxes(int maxAxes)
Definition: Frame.h:1266
void setDut1(double dut1)
Definition: Frame.h:1233
std::shared_ptr< Frame > copy() const
Return a deep copy of this object.
Definition: Frame.h:179
double d2
Resolved vector component 2.
Definition: Frame.h:85
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Definition: base.cc:49
void setSymbol(int axis, std::string const &symbol)
Definition: Frame.h:1311
DirectionPoint(double direction, PointD const &point)
Definition: Frame.h:51
double getBottom(int axis) const
Definition: Frame.h:805
std::vector< double > point
Point.
Definition: Frame.h:83
double axAngle(PointD const &a, PointD const &b, int axis) const
Definition: Frame.h:226
void setFormat(int axis, std::string const &format)
Definition: Frame.h:1248
void permAxes(std::vector< int > perm)
Definition: Frame.h:1136
double direction
Direction, an angle in radians.
Definition: Frame.h:51
void setTop(int axis, double top)
Definition: Frame.h:1329
bool getMatchEnd() const
Definition: Frame.h:856
std::string getFormat(int axis) const
Definition: Frame.h:840
PointD point
Point.
Definition: Frame.h:53
void setEpoch(double epoch)
Definition: Frame.h:1238
void setDomain(std::string const &domain)
Definition: Frame.h:1228
std::string getSystem() const
Definition: Frame.h:916
int getNAxes() const
Definition: Frame.h:874
std::shared_ptr< Mapping > mapping
Mapping.
Definition: Frame.h:102
std::string getUnit(int axis) const
Definition: Frame.h:931
void setBottom(int axis, double bottom)
Definition: Frame.h:1206
Frame(AstFrame *rawPtr)
Definition: Frame.h:1517
double d1
Resolved vector component 1.
Definition: Frame.h:84
void setLabel(int axis, std::string const &label)
Definition: Frame.h:1255
std::shared_ptr< Frame > frame
Frame.
Definition: Frame.h:101
void setObsAlt(double alt)
Definition: Frame.h:1277
PointD offset(PointD point1, PointD point2, double offset) const
Definition: Frame.h:1073
int getDigits() const
Definition: Frame.h:810
bool getDirection(int axis) const
Definition: Frame.h:820
std::string getObsLon() const
Definition: Frame.h:895
int getDigits(int axis) const
Definition: Frame.h:815
double getObsAlt() const
Definition: Frame.h:885
std::vector< double > PointD
Definition: base.h:57
double angle(PointD const &a, PointD const &b, PointD const &c) const
Definition: Frame.h:198
int getMaxAxes() const
Definition: Frame.h:862
std::string getInternalUnit(int axis) const
Definition: Frame.h:846
double getEpoch() const
Definition: Frame.h:835
std::string getTitle() const
Definition: Frame.h:921
int getNIn() const
Definition: Mapping.h:76
bool getPreserveAxes() const
Definition: Frame.h:905
std::string getLabel(int axis) const
Definition: Frame.h:851
ResolvedPoint(int naxes)
Definition: Frame.h:82
double getDut1() const
Definition: Frame.h:830
void setObsLon(std::string const &lon)
Definition: Frame.h:1287
void setMatchEnd(bool match)
Definition: Frame.h:1260
bool getPermute() const
Definition: Frame.h:900
std::string getAlignSystem() const
Definition: Frame.h:800
void setUnit(int axis, std::string const &unit)
Definition: Frame.h:1334
void setSystem(std::string const &system)
Definition: Frame.h:1319
std::string getObsLat() const
Definition: Frame.h:890
FrameMapping(std::shared_ptr< Frame > frame, std::shared_ptr< Mapping > mapping)
Definition: Frame.h:99
void setDirection(bool direction, int axis)
Definition: Frame.h:1221
DirectionPoint offset2(PointD const &point1, double angle, double offset) const
Definition: Frame.h:1115
NReadValue unformat(int axis, std::string const &str) const
Definition: Frame.h:1500
std::string format(int axis, double value) const
Definition: Frame.h:780
std::string getSymbol(int axis) const
Definition: Frame.h:910
void setMinAxes(int minAxes)
Definition: Frame.h:1272
Definition: FrameSet.h:99
void setActiveUnit(bool enable)
Definition: Frame.h:1293
void setAlignSystem(std::string const &system)
Definition: Frame.h:1201
NReadValue(int nread, double value)
Definition: Frame.h:67
bool getActiveUnit() const
Definition: Frame.h:790
virtual std::shared_ptr< Object > copyPolymorphic() const override
Definition: Frame.h:1525
void setPermute(bool permute)
Definition: Frame.h:1301
void setDigits(int axis, int digits)
Definition: Frame.h:1216
void setTitle(std::string const &title)
Definition: Frame.h:1324
std::vector< int > matchAxes(Frame const &other) const
Definition: Frame.h:983
PointD norm(PointD value) const
Definition: Frame.h:1043
std::string getNormUnit(int axis) const
Definition: Frame.h:880
int getMinAxes() const
Definition: Frame.h:868