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

Integrat Error Message #26

Open
Mira-Ao opened this issue Feb 21, 2024 · 5 comments
Open

Integrat Error Message #26

Mira-Ao opened this issue Feb 21, 2024 · 5 comments

Comments

@Mira-Ao
Copy link

Mira-Ao commented Feb 21, 2024

When I do the following as described on the your website

./integrat MPCORB.DAT new_mpcorb.dat 2019 11 01 -f ~/.find_orb/*430t

This kind of output is returned and I cannot proceed.

Using an update
Got 0 hashes
Integrat version Feb 20 2024 11:55:54
Integrating to 1 Jan 2019 0:00:00 = JD 2458484.50000
Got '0'
integrat: integrat.cpp:714: double try_to_integrate(const char*, char*, double, double, double): Assertion `got_it' failed.
Aborted

@Bill-Gray
Copy link
Owner

Ah, I see the problem... basically, the documentation is lagging behind the program.

I ran into numerous issues with the MPCORB format. In particular, we had cases where an orbit would start out with eccentricity just barely below 1.0, and after integration would have eccentricity just barely above 1.0. MPCORB can't handle parabolic or hyperbolic orbits.

So the program now uses the "standard" orbit format (.sof), which can handle all orbits (and has numerous other advantages). I'll modify the documentation to reflect this.

@Mira-Ao
Copy link
Author

Mira-Ao commented Feb 22, 2024

Thank you, please let me know when the relevant part of the document has been reflected.
One more thing, regarding the comet orbit file, the JPL link in the documentation is already invalid. Which file should I use in this case?
https://www.minorplanetcenter.net/iau/MPCORB/CometEls.txt
Is this correct?

@Mira-Ao
Copy link
Author

Mira-Ao commented Feb 22, 2024

What I did before I ran Integrat was download MPC's AllCometEls.txt and rename the file to ELEMENTS.COMET,
https://projectpluto.com/astcheck.htm
I ran mpc2sof by referring to the web at So the sof format should be used.

So it is not only https://projectpluto.com/find_sou.htm#run_integrat that I saw.

@Mira-Ao
Copy link
Author

Mira-Ao commented Mar 21, 2024

Perhaps you updated some web pages, I don't think the Integrat exe file (or the link to the dedicated Integrat page) was there before, am I wrong?

@Mira-Ao
Copy link
Author

Mira-Ao commented Mar 23, 2024

OK, I finally understand what you are saying. So now Integrat is for integrating .sof files from the beginning.

But even using jpl's 430t is terribly time consuming.
ephem start 2460398.500000
Perturber (1) Ceres calculated
Perturber (2) Pallas calculatedconds remain; 1 done 0
Perturber (4) Vesta calculatedeconds remain; 3 done 0
Hi! I've got process number 0, PID 5158, parent's is 27
Hi! I've got process number 1, PID 5315, parent's is 51
Hi! I've got process number 2, PID 5316, parent's is 53
Hi! I've got process number 3, PID 5317, parent's is 53
114 seconds elapsed; 746567 seconds remain; 204 done 0

Does this indication mean that it is frightening that it takes so much time just to calculate Vesta, which is only the fourth of a million asteroids?

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

No branches or pull requests

2 participants