lsst.astshim  14.0-8-g6380fcc
utils.h
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 #ifndef ASTSHIM_DETAIL_UTILS_H
23 #define ASTSHIM_DETAIL_UTILS_H
24 
25 #include <cctype>
26 #include <limits>
27 #include <sstream>
28 #include <stdexcept>
29 #include <vector>
30 
31 #include "astshim/base.h"
32 
33 namespace ast {
34 namespace detail {
35 
36 static const int FITSLEN = 80;
37 
39 inline void annulAstObject(AstObject *object) {
40  if (object != nullptr) {
41  astAnnul(object);
42  }
43 }
44 
45 template <typename T1, typename T2>
46 inline void assertEqual(T1 val1, std::string const &descr1, T2 val2, std::string const &descr2) {
47  if (val1 != val2) {
48  std::ostringstream os;
49  os << descr1 << " = " << val1 << " != " << descr2 << " = " << val2;
50  throw std::invalid_argument(os.str());
51  }
52 }
53 
57 inline void astBadToNan(std::vector<double> &p) {
58  for (auto &val : p) {
59  if (val == AST__BAD) {
60  val = std::numeric_limits<double>::quiet_NaN();
61  }
62  }
63 }
64 
68 void astBadToNan(ast::Array2D const &arr);
69 
77 inline std::string formatAxisAttr(std::string const &name, int axis) {
78  std::stringstream os;
79  os << name << "(" << axis << ")";
80  return os.str();
81 }
82 
86 bool isSeries(AstCmpMap const *cmpMap);
87 
93 std::string getClassName(AstObject const *rawObj);
94 
98 inline double safeDouble(double val) {
99  assertOK();
100  return val != AST__BAD ? val : std::numeric_limits<double>::quiet_NaN();
101 }
102 
106 inline std::string stringToUpper(std::string const &str) {
107  std::string upstr;
108  upstr.reserve(str.size());
109  for (char c : str) {
110  upstr += std::toupper(c);
111  }
112  return upstr;
113 }
114 
115 } // namespace detail
116 } // namespace ast
117 
118 #endif
ndarray::Array< double, 2, 2 > Array2D
Definition: base.h:42
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Definition: base.cc:49