Memory issue with satellite.at()?

I have an application where I am trying to calculate the position of a satellite on a very fine time grid (say millions of points).

The simple call to geocentric = satellite.at(t) seems to run into performance issue, where for large n scales much worse than linearly, and by the time n is O(1 M) the thing is essentially requiring >10 GB or RAM and bringing my terminal down.

I can make a more detailed report, including all the related information, but I want to make a quick check first of whether I am doing something wrong and/or you have any immediate insight into this.

Thanks in advance!

[lucabaldini]

Wow! The t._nutation_angles = iau2000b(t.tt) really saved my day.

I haven’t made extensive tests, yet, but at the first glance the time/memory footprint is well below the level at which I would notice it.

I don’t know if mine is really a peculiar use case, but it might be worth noting this in the docs.

Thanks!

[brandon-rhodes]

The new code has now been released in Skyfield 1.21. I’m going to close the issue for now, but please re-open if you run into further problems!