LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
base.cc
Go to the documentation of this file.
1 #include "astshim/base.h"
2 
3 #include <sstream>
4 #include <stdexcept>
5 #include <string>
6 #include <vector>
7 
8 namespace ast {
9 namespace {
10 
11 static std::ostringstream errorMsgStream;
12 
13 /*
14 Write an error message to `errorMsgStream`
15 
16 Intended to be registered as an error handler to AST by calling `astSetPutErr(reportError)`.
17 */
18 void reportError(int errNum, const char *errMsg) { errorMsgStream << errMsg; }
19 
20 /*
21 Instantiate this class to register `reportError` as an AST error handler.
22 */
23 class ErrorHandler {
24 public:
25  ErrorHandler() { astSetPutErr(reportError); }
26 
27  ErrorHandler(ErrorHandler const &) = delete;
28  ErrorHandler(ErrorHandler &&) = delete;
29  ErrorHandler &operator=(ErrorHandler const &) = delete;
30  ErrorHandler &operator=(ErrorHandler &&) = delete;
31 
32  static std::string getErrMsg() {
33  auto errMsg = errorMsgStream.str();
34  // clear status bits
35  errorMsgStream.clear();
36  if (errMsg.empty()) {
37  errMsg = "Failed with AST status = " + std::to_string(astStatus);
38  } else {
39  // empty the stream
40  errorMsgStream.str("");
41  }
42  astClearStatus;
43  return errMsg;
44  }
45 };
46 
47 } // namespace
48 
49 void assertOK(AstObject *rawPtr1, AstObject *rawPtr2) {
50  // Construct ErrorHandler once, the first time this function is called.
51  // This is done to initialize `errorMsgStream` and register `reportError` as the AST error handler.
52  // See https://isocpp.org/wiki/faq/ctors#static-init-order-on-first-use
53  static ErrorHandler *errHandler = new ErrorHandler();
54  if (!astOK) {
55  if (rawPtr1) {
56  astAnnul(rawPtr1);
57  }
58  if (rawPtr2) {
59  astAnnul(rawPtr2);
60  }
61  throw std::runtime_error(errHandler->getErrMsg());
62  }
63 }
64 
66  return static_cast<ConstArray2D>(arrayFromVector(const_cast<std::vector<double> &>(vec), nAxes));
67 }
68 
70  int nPoints = vec.size() / nAxes;
71  if (nPoints * nAxes != vec.size()) {
73  os << "vec length = " << vec.size() << " not a multiple of nAxes = " << nAxes;
74  throw std::runtime_error(os.str());
75  }
76  Array2D::Index shape = ndarray::makeVector(nAxes, nPoints);
77  Array2D::Index strides = ndarray::makeVector(nPoints, 1);
78  return external(vec.data(), shape, strides);
79 }
80 
81 } // namespace ast
std::ostream * os
Definition: Schema.cc:557
T clear(T... args)
T data(T... args)
AST wrapper classes and functions.
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition: base.cc:49
ConstArray2D arrayFromVector(std::vector< double > const &vec, int nAxes)
Reshape a vector as a 2-dimensional array that shares the same memory.
Definition: base.cc:65
ndarray::Array< double, 2, 2 > Array2D
2D array of const double; typically used for lists of const points
Definition: base.h:42
ndarray::Array< const double, 2, 2 > ConstArray2D
2D array of const double; typically used for lists of const points
Definition: base.h:46
FilterProperty & operator=(FilterProperty const &)=default
T size(T... args)
T str(T... args)
T to_string(T... args)