Конвертировать boost :: posix_time :: ptime в NTP datetamp

Мне нужно конвертировать boost::posix_time::ptime в NTP Датапечать в соответствии с
RFC 5905 представлен следующей структурой:

struct NtpDatestamp {
std::int32_t era_number;
std::uint32_t seconds_since_era_epoch;
std::uint64_t fraction_of_second;
};

RFC 5905 гласит следующее:

Для преобразования системного времени в любом формате в формат даты и времени NTP
требует, чтобы количество секунд s от эпохи расцвета до системы
время будет определено. Определить целое число era а также timestamp дано s,

era = s / 2^(32) and timestamp = s - era * 2^(32),

который работает для положительных и отрицательных дат. Определить s учитывая эпоху
и отметка времени,

s = era * 2^(32) + timestamp.

Поэтому я попробовал следующее:

const auto system_time = boost::posix_time::time_from_string("1899-12-31 00:00:00.000");

const boost::posix_time::ptime prime_epoch{boost::gregorian::date{1900, 1, 1}};
// Calculate the number of seconds from the prime epoch to the system time.
const boost::posix_time::time_duration time_duration{system_time - prime_epoch};
const std::int64_t s{time_duration.total_seconds()};
const std::int32_t era_number{static_cast<std::int32_t>(s / std::pow(2, 32))};
const std::uint64_t seconds_since_era_epoch{static_cast<std::uint64_t>(s - s / std::pow(2, 32) * std::pow(2, 32))};
// The fraction of a NTP Datestamp is measured in Attoseconds.
const std::uint64_t fraction_of_second{static_cast<std::uint64_t>(time_duration.total_microseconds() * 1e12)};

Но это дает неверные результаты.

Я полностью озадачен этой (на самом деле простой) проблемой в данный момент.

Может ли кто-нибудь направить меня в правильном направлении? Как я могу получить номер эпохи, смещение эпохи а также доля даты NTP от boost::posix_time::ptime?

Редактировать: Либо расчеты в RFC 5905 недостаточно точны, либо я их неверно истолковал. Благодаря комментариям я изменил расчет на следующий (на этот раз полный пример):

#include <cmath>
#include <cstdint>
#include <iostream>

#include <boost/date_time.hpp>

int main() {
const auto system_time =
boost::posix_time::time_from_string("1899-12-31 00:00:00.000");

const boost::posix_time::ptime prime_epoch{
boost::gregorian::date{1900, 1, 1}};
// Calculate the number of seconds from the prime epoch to the system time.
const boost::posix_time::time_duration time_duration{prime_epoch -
system_time};

// s is correctly determined now.
std::int64_t s{time_duration.total_seconds()};
if (prime_epoch > system_time) {
// boost::posix_time::time_duration does not take the sign into account.
s *= -1;
}

// TODO(wolters): The following calculations do not return the correct
// results, but the RFC 5905 states them
const std::int32_t era{static_cast<std::int32_t>(s / std::pow(2, 32))};
const std::uint64_t timestamp{
static_cast<std::uint64_t>(s - era * std::pow(2, 32))};
// The fraction of a NTP Datestamp is measured in Attoseconds.
// TODO(wolters): `boost::posix_time::ptime` does NOT resolve to attoseconds,
// but doesn't the target format expect the value to be specified as
// attoseconds? Doesn't the following depend on Boost compile options?
const std::uint64_t fraction{
static_cast<std::uint64_t>(time_duration.fractional_seconds())};

std::cout << "s = " << std::dec << s << '\n';
// TODO(wolters): This does still not match the expected results; taken from
// Figure 4 of https://www.ietf.org/rfc/rfc5905.txt
std::cout << "Era (expected: -1) = " << std::dec << era << '\n';
std::cout << "Timestamp (expected: 4294880896) = " << std::dec << timestamp
<< '\n';
std::cout << "Fraction (expected: 0) = " << std::dec << fraction << '\n';
}

s сейчас рассчитывается правильно, но остальные расчеты неверны. Я думаю, что я пропускаю что-то важное полностью …