Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Increase precision of Duration to cover 10^-29 seconds to 10^28 seconds #186

Open
ChristopherRabotin opened this issue Dec 3, 2022 · 11 comments · May be fixed by #326
Open

Increase precision of Duration to cover 10^-29 seconds to 10^28 seconds #186

ChristopherRabotin opened this issue Dec 3, 2022 · 11 comments · May be fixed by #326

Comments

@ChristopherRabotin
Copy link
Member

ChristopherRabotin commented Dec 3, 2022

This is really just for fun and I don't imagine anyone ever needing this, but behind an ultra feature crate, the Duration struct could change the centuries from i16 to i64 and add a subnanos field of type u64. This would make the structure relatively big (192 bits / 24 octets) but it would allow representations of duration from 5.4e-29 seconds (or 54 quectoseconds) to 2.9e+28 seconds (2.9e+4 yottaseconds or 29 qettaseconds).

According to https://en.wikipedia.org/wiki/Orders_of_magnitude_(time), this means that the smallest duration would be 162 thousand times longer than the mean lifetime of W and Z bosons and the largest time would be just short of the radioactive half-life of tellurium-128, the longest known half-life of any elemental isotope.

All in one library! \o/

Note: Currently, the maximum duration is 1.03e+14 seconds, or just short of the typical lifetime of a biological species on Earth (130 teraseconds -- Duration::MAX is 103 teraseconds).

Math:

In [7]: f"{Duration.init_from_max().to_seconds():e}"
Out[7]: '1.034079e+14'
In [8]: 1.0/2**64 * 1e-9 # they are sub-nanoseconds hence the 1e-9
Out[8]: 5.421010862427523e-29
(...)
In [14]: (Unit.Century * 1).to_seconds()
Out[14]: 3155760000.0

In [15]: 2**(64-1) * 3155760000.0 # we can't do the math with Durations since they max out
Out[15]: 2.9106748539024827e+28
@ChristopherRabotin ChristopherRabotin changed the title Increase precisions from 10^-29 seconds to 10^28 seconds Increase precision of Duration to cover 10^-29 seconds to 10^28 seconds Dec 3, 2022
@ChristopherRabotin
Copy link
Member Author

More seriously, it might make sense to have a crate feature called attosecond where simply a u32 would be added as subnanos.

Adding a u32 as subnanos would support 0.2328... attoseconds. Note that to support all of the specs of the Precise Time Protocol, femtoseconds may be needed since the maximum allowable Jitter is 250 fm: https://github.com/opencomputeproject/Time-Appliance-Project/tree/master/Open-Time-Server/#pps-out .

In [3]: 1.0/2**32 * 1e-9
Out[3]: 2.3283064365386964e-19

@ChristopherRabotin
Copy link
Member Author

I should run some benchmarks to see if there is any impact in changing the centuries to a 32 bit integer and adding that subnanos field. The fields of Duration are not exposed publicly, so changing those won't be a breaking change.

@gwbres
Copy link
Collaborator

gwbres commented Jun 18, 2024

Hello,

I'm curious as to where this topic lies on your roadmap ?

You have stated in several places that the current accuracy is 1E-9s.
My current application is RADAR positioning, which is then theoretically limited to 30cm accuracy.
When you say 1E-9, I presume it is not the actual accuracy but more like the smallest error the testbench tolerates, therefore the "garanteed" value, is that correct ? Don't get me wrong, it's great to have one, very few people go this far because it's usually the most complex thing to express, correctly define and target. And I say this because I work in the industry, use and make instruments 🤣

So far my experience is that hifitime has never been the limiting factor, but I'm also thinking about the future and what we could do better and hope for. IMO it would be the minimum to support what PTP requires. I'm pretty sure PTP will once again evolve very soon, most likely when PTP3/White rabbit becomes the standard, which already achieves 1E-11 accuracy (somewhere around that). Lab clocks from the 90s have an accuracy of 1E-15 at 1second. The best clock these days are located around 1E-18/1E-19 (optical clocks). Anyways, I would argue that it is not something you want to change regularly and should be well thought ahead and scaled properly (and indeed, should remain focused on physics). And also, a 10x magnitude margin is always minimum :)

Side note (a)

I have been trying to keep track of this project for a while, but just like you I am far too involved in my own stuff, it is hard to take part in anything else. I'm curious as to how they express time. If they did not choose hifitime, it might be a good idea to discuss and see if they could join the adventure

Side note (b)

Me laughing at this proposal, is not me laughing at the actual proposal but more at the idea of making this toolkit better, while it is already outstanding. And yes, emojis are not and will never be proper means of communications 😆

@ChristopherRabotin
Copy link
Member Author

So I've yet to start any work on this, but if needed, I don't think it'll be a huge amount of work (especially once #296 is implemented since lots of the testing can be formally verified). As long as there is no backward incompatibility, we could also plan this for version 4.1.

Concerning stattime, I reached out to them a while ago and never heard back: pendulum-project/statime#333 . They only closed the issue without comment. So I assume that they are more interested in rewriting the code and opening the door for mistakes that have been solved with hifitime -- but it isn't my problem ;-) I know of at least one other project that didn't want to use hifitime although they have exactly the same use case I had originally, and specifically because they said that they want to try to make something similar.

@gwbres
Copy link
Collaborator

gwbres commented Jun 19, 2024

So I've yet to start any work on this, but if needed, I don't think it'll be a huge amount of work (especially once #296 is implemented since lots of the testing can be formally verified). As long as there is no backward incompatibility, we could also plan this for version 4.1.

Sounds great

Concerning stattime, I reached out to them a while ago and never heard back: pendulum-project/statime#333 . They only closed the issue without comment. So I assume that they are more interested in rewriting the code and opening the door for mistakes that have been solved with hifitime -- but it isn't my problem ;-) I know of at least one other project that didn't want to use hifitime although they have exactly the same use case I had originally, and specifically because they said that they want to try to make something similar.

it looks like a private company (although I don't really quite grasp these private companies working on github). That's their decision, we'll see what they achieve. I see that they use std::time shipped by tokio. As a side note, all my GNSS related libraries (which I consider my "main" libraries these days) use hifitime and not std anymore

@ChristopherRabotin
Copy link
Member Author

Should Duration be a single i128?

Currently Unit supports multiplications with one float type (f64) and one signed integer type (i64). The advantage of supporting only one type for each is that when the user requests something like 0.159 * Unit::Second, the compiler knows to convert that 0.159 to a f64. Similarly, when the user specifies -159 * Unit::Centuries, the compiler knows to convert that to a i64.

Using a single i128 integer would allow the Duration to cover 5e19 centuries if that was an integer counter on the nanoseconds. More usefully though, we could change this integer to count the number of attoseconds (1e-18 seconds) (or smaller!) and that would still cover 53,914,487,622 centuries.

I'd probably just recommend that the Duration be a counter of zepto seconds (1e-21 seconds) which would still allow the representation of dates for up to 53,914,487 centuries, or roughly 1.7e17 seconds, or a bit more than the age of the Earth (and about 25% of the age of the universe).

This kind of seems like a massive overkill, but it would also greatly simplify the code for Duration, getting rid of the normalize() function all together, and enabling lots of const fn functions. It would extend the list of possible units and frequencies (a functionality I actually never use). I also expect that it would massively speed up the computations of durations because we wouldn't need to normalize any of the durations.

@gwbres , what do you think?

@pwnorbitals, you have a lot of experience in embedded systems: would using a signed 128-bit integer be an issue for supporting some of the platforms you've demonstrated hifitime on? As of March 2024, the Rust compiler has fixed the alignment difference for i128/u128 between Rust and C, if I understand correctly.

@ChristopherRabotin
Copy link
Member Author

The i128 is part of core rust. Soon, I'll try to implement this change and see how it affects performance. Switching to a supported type would simplify a lot of the testing and complexity of the Duration type.

@gwbres
Copy link
Collaborator

gwbres commented Jul 25, 2024

The i128 is part of core rust. Soon, I'll try to implement this change and see how it affects performance. Switching to a supported type would simplify a lot of the testing and complexity of the Duration type.

looking forward to it, keep us updated

@gwbres
Copy link
Collaborator

gwbres commented Aug 5, 2024

Hello,

to illustrate this topic, Hifitime-4 reaches its limit when we survey laboratories that have stable clocks. I say stable to differentiate with Ultra stable clocks. The following performances are probably in the region of MASER-H clocks (about 1E-13@1sec). They are beyond GNSS performances but also more than millions of times "poorer" than the best technology, yet already very far from general public technology.

As an example, here we survey a station for 48 hours, you can see that it is static

alt

This is its absolute offset to GST that we obtained with GNSS-RTK + Hifitime-4.
The 1ns granularity (1m on the RADAR) is very clear. The clock is stable enough we barely see any drift

clk
drift

Here I activate +cggtts (extra smoothing on timing)

srsys

@ChristopherRabotin
Copy link
Member Author

The change to a zeptosecond counter allows the Duration to max out at 53914487 centuries 22747 days 2 h 20 min 41 s 43 ms 369 μs 864 ns (plus some more sub nanoseconds which aren't shown here).

@ChristopherRabotin
Copy link
Member Author

It seems to me that Hifitime v4 has been in development for too long. I think it's best if we release version 4.0 now, which has been in use in several projects on the main branch for quite some time, and move the remaining version 4 issues to a version 5. We would then switch to the development of version 5 soon. What do you think, @gwbres ?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants