Comments (3)
Hi Gabriel,
I would not use a collision detection routine for this. I would just periodically check if the particles have swapped their position. You can do that while the simulation is running using a loop, or afterwards using the snapshots saved in the simulationarchive. Unless you check for this condition very often, it should not be a limiting factor for the runtime. You only get the time right approximately, but maybe that's good enough. You can always go back and check more often.
As a side note, when you calculate the semi-major axes of all particles, you probably want to use heliocentric coordinates rather than Jacobi coordinates.
I hope this helps.
Hanno
from rebound.
Hi, Hanno,
Thanks for the suggestions! I think that checking for around 50.000 output points (worst case scenario, since collisions can be detected earlier) won't be much a of a problem.
For larger integration times, the resolution at the lower end might be a little poor (on the order of 1e3 years, which is around 1e4 orbits of my innermost planet), but it's not an absolute problem for statistical purposes, I guess.
I will definitely pay more attention and make more use of the SimulationArchive, it's much better and easier than having to filter out which variables and analysis to make beforehand.
I will also keep in mind the change of coordinates and make some experiment with it! Thanks.
I will close it for now!
Kind regards,
Gabriel.
from rebound.
"Manually" check for collision condition after each timestep seems to be simpler for now. CLosing it!
from rebound.
Related Issues (20)
- `configure_ghostboxes` has no effect HOT 2
- Attribute error Question HOT 3
- Inconsistency in solution across inclination HOT 16
- Adding particles via Horizons system HOT 3
- There are quite significant differences in the orbital elements of asteroids and their clones HOT 9
- Soft sphere implementation HOT 3
- MEGNO time evolution HOT 17
- WHFast512 on macOS system HOT 1
- run solar system with and without planet to planet gravity forces HOT 1
- MOID calculation at a given timestep HOT 6
- How to read the binary file output? HOT 8
- how to read snapshots in a loop HOT 4
- viewing a SA from the command line ("sadump") HOT 6
- There don't appear to be any reference docs? HOT 3
- TRACE crashes after adding 2 particles after a merger event (on ARM-based Mac; [edit:] on Linux, too) HOT 3
- How to get the orbit elements of the first star? HOT 3
- 'Barycentre' primary convenience function HOT 3
- Create comets from particles in mean motion resonance with an eccentric planet HOT 13
- [TRACE+] Binary files size from simulation archive dependent on integrator HOT 11
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rebound.