oem.hpp

Go to the documentation of this file.

#pragma once
 
#include "ffi_core.hpp"
 
#include <array>
#include <string>
#include <string_view>
#include <vector>
 
namespace siderust {
 
namespace oem {
 
namespace detail {
 
struct OemStatesGuard {
  SiderustOemState *ptr = nullptr;
  unsigned long count = 0;
 
  ~OemStatesGuard() { siderust_oem_states_free(ptr, count); }
 
  OemStatesGuard() = default;
  OemStatesGuard(const OemStatesGuard &) = delete;
  OemStatesGuard &operator=(const OemStatesGuard &) = delete;
  OemStatesGuard(OemStatesGuard &&) = delete;
  OemStatesGuard &operator=(OemStatesGuard &&) = delete;
};
 
} // namespace detail
 
struct StateVector {
  double epoch_jd;               
  std::array<double, 3> pos_km;  
  std::array<double, 3> vel_kms; 
};
 
inline std::vector<StateVector> parse(std::string_view text) {
  const std::string buf{text};
  SiderustOemState *raw_ptr = nullptr;
  unsigned long count = 0;
 
  check_status(siderust_oem_parse_str(buf.c_str(), &raw_ptr, &count), "oem::parse");
 
  detail::OemStatesGuard guard{};
  guard.ptr = raw_ptr;
  guard.count = count;
 
  std::vector<StateVector> result;
  result.reserve(static_cast<std::size_t>(count));
  for (unsigned long i = 0; i < count; ++i) {
    const auto &s = raw_ptr[i];
    result.push_back({s.epoch_jd,
                      {s.pos_km[0], s.pos_km[1], s.pos_km[2]},
                      {s.vel_kms[0], s.vel_kms[1], s.vel_kms[2]}});
  }
 
  return result;
}
 
} // namespace oem
 
} // namespace siderust