This is the group project of Emilie Grégoire and Mitja Desmet for the course of computational physics in first master of physics at the VUB.
The structure is as follows:
- project.py contains all definitions necessary to numerically integrate the geodesic equation and write the results to a file.
- animation.py is used to plot the results from the integration using VPython. Best is to open this file in VIdle, the IDE which is installed alongside VPython. The functions are written specifically for the metrics we analysed to more accurately represent the situation, but the one for Minkowski should be good for every metric.
- testing.py was used to actually do all the numerical analysis and testing.
- einstein.py uses the EinsteinPy package do also calculate geodesics. We used this to compare our results to.
All .txt files are results from the different computations. The ones denote with solveEinstein come from the EinsteinPy package. The ones which start with solveGE come from calculations done with the solveGE function (see project.py). Similarly, the ones with RK4 are the ones which come from RK4 (also see project.py for implementation).
- If there appears Minskowski in the name, the considered spacetime was Minkowski spactime. Otherwise it was Schwarzschild spacetime. For this last one we always considered a mass of 5.972e24 kg.
- The files with phi or theta are results from calculations with as initial condition only had a velocity in the phi or theta direction respectively.
- mass refers to the only calculation with a different mass: the mass considerd was 5.972e30 kg.
- Finally randphi had as initial condition a velocity in both the r and phi direction.
To end, interesting-links contains some links to resources we used for our project.