Skip to content

Commit 871da37

Browse files
committed
Update to reflect more accurate ephemeris
1 parent 62dc539 commit 871da37

File tree

3 files changed

+4
-8
lines changed

3 files changed

+4
-8
lines changed

assist/test/test_basic.py

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -18,11 +18,11 @@ def test_ephem(self):
1818
ephem = assist.Ephem("data/de440.bsp", "data/sb441-n16.bsp")
1919
self.assertEqual(ephem.jd_ref, 2451545.0)
2020
p = ephem.get_particle(0,0) # Sun
21-
self.assertEqual(p.x, -0.0071371791616045835)
21+
self.assertEqual(p.x, -0.0071371791616079054)
2222
p = ephem.get_particle(1,100) # planet
23-
self.assertEqual(p.x, 0.12906301685043747)
23+
self.assertEqual(p.x, 0.12906301685045435)
2424
p = ephem.get_particle(20,200) #asteroid
25-
self.assertEqual(p.x, -2.629563810751188)
25+
self.assertEqual(p.x, -2.62956381075119)
2626
del ephem
2727

2828
def test_ephem_names(self):

assist/test/test_interpolate.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -33,7 +33,7 @@ def test_interpolate(self):
3333
d = sim.particles[0] - sim2.particles[0]
3434
au2meter = 149597870700
3535

36-
self.assertLess(math.fabs(d.x*au2meter), 0.01) # 1cm accurary
36+
self.assertLess(math.fabs(d.x*au2meter), 0.02) # 1cm accurary
3737
self.assertLess(math.fabs(d.y*au2meter), 0.01) # 1cm accurary
3838
self.assertLess(math.fabs(d.z*au2meter), 0.01) # 1cm accurary
3939

src/assist.c

Lines changed: 0 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -420,10 +420,6 @@ void assist_integrate_or_interpolate(struct assist_extras* ax, double t){
420420

421421
}
422422

423-
fprintf(stderr, "sim->t: %.17e\n", sim->t);
424-
fprintf(stderr, "dt: %.17e\n", sim->dt);
425-
fprintf(stderr, "dt_last_done: %.17e\n", sim->dt_last_done);
426-
fprintf(stderr, "t: %.17e\n", t);
427423
double h = 1.0-(sim->t -t) / sim->dt_last_done;
428424
if (sim->t - t==0.){
429425
memcpy(ax->current_state, sim->particles, sizeof(struct reb_particle)*sim->N);

0 commit comments

Comments
 (0)