Remove accuracy argument from grav_undo_planets()

Author: attipaci

include/novas.h

@@ -1155,7 +1155,7 @@ int grav_undef(double jd_tdb, enum novas_accuracy accuracy, const double *pos_ap
 
 int grav_planets(const double *pos_src, const double *pos_obs, const novas_planet_bundle *planets, double *out);
 
-int grav_undo_planets(const double *pos_app, const double *pos_obs, enum novas_accuracy accuracy, const novas_planet_bundle *planets, double *out);
+int grav_undo_planets(const double *pos_app, const double *pos_obs, const novas_planet_bundle *planets, double *out);
 
 int make_airborne_observer(const on_surface *location, const double *vel, observer *obs);
 

src/frames.c

@@ -923,7 +923,7 @@ int novas_app_to_geom(const novas_frame *frame, enum novas_reference_system sys,
   frame_aberration(frame, APP_TO_GEOM, app_pos);
 
   // Undo gravitational deflection and aberration.
-  prop_error(fn, grav_undo_planets(app_pos, frame->obs_pos, frame->accuracy, &frame->planets, geom_icrs), 0);
+  prop_error(fn, grav_undo_planets(app_pos, frame->obs_pos, &frame->planets, geom_icrs), 0);
 
   return 0;
 }

src/novas.c

@@ -4388,8 +4388,6 @@ int grav_planets(const double *pos_src, const double *pos_obs, const novas_plane
  * @param pos_obs     [AU] Position 3-vector of observer (or the geocenter), with respect to
  *                    origin at solar system barycenter, referred to ICRS axes,
  *                    components in AU.
- * @param accuracy    NOVAS_FULL_ACCURACY (0) for sub-μas accuracy or NOVAS_REDUCED_ACCURACY (1)
- *                    for sub-mas accuracy.
  * @param planets     Apparent planet data containing positions and velocities for the major
  *                    gravitating bodies in the solar-system.
  * @param[out] out    [AU] Nominal position vector of observed object, with respect to origin at
@@ -4405,12 +4403,12 @@ int grav_planets(const double *pos_src, const double *pos_obs, const novas_plane
  * @since 1.1
  * @author Attila Kovacs
  */
-int grav_undo_planets(const double *pos_app, const double *pos_obs, enum novas_accuracy accuracy, const novas_planet_bundle *planets, double *out) {
+int grav_undo_planets(const double *pos_app, const double *pos_obs, const novas_planet_bundle *planets, double *out) {
   static const char *fn = "grav_undo_planets";
 
+  const double tol = 1e-13;
   double pos_def[3] = { }, pos0[3] = { };
   double l;
-  double tol = accuracy == NOVAS_FULL_ACCURACY ? 1e-13 : 1e-10;
   int i;
 
   if(!pos_app || !pos_obs)
@@ -4499,7 +4497,6 @@ int obs_planets(double jd_tdb, enum novas_accuracy accuracy, const double *pos_o
   if(!pos_obs)
     return novas_error(-1, EINVAL, fn, "NULL observer position parameter");
 
-
   // Set up the structures of type 'object' containing the body information.
   if(!initialized) {
     for(i = 0; i < NOVAS_PLANETS; i++)
@@ -4666,7 +4663,7 @@ int grav_undef(double jd_tdb, enum novas_accuracy accuracy, const double *pos_ap
     return novas_error(-1, EINVAL, fn, "NULL source position 3-vector: pos_app=%p, out=%p", pos_app, out);
 
   prop_error(fn, obs_planets(jd_tdb, accuracy, pos_obs, pl_mask, &planets), 0);
-  prop_error(fn, grav_undo_planets(pos_app, pos_obs, accuracy, &planets, out), 0);
+  prop_error(fn, grav_undo_planets(pos_app, pos_obs, &planets, out), 0);
   return 0;
 }
 

test/src/test-errors.c

@@ -917,13 +917,13 @@ static int test_grav_undo_planets() {
   double p[3] = {2.0}, po[3] = {0.0, 1.0}, out[3] = {};
   int n = 0;
 
-  if(check("grav_undo_planets:pos_app", -1, grav_undo_planets(NULL, po, NOVAS_REDUCED_ACCURACY, &planets, out))) n++;
-  if(check("grav_undo_planets:pos_obs", -1, grav_undo_planets(p, NULL, NOVAS_REDUCED_ACCURACY, &planets, out))) n++;
-  if(check("grav_undo_planets:planets", -1, grav_undo_planets(p, po, NOVAS_REDUCED_ACCURACY, NULL, out))) n++;
-    if(check("grav_undo_planets:pos_src", -1, grav_undo_planets(p, po, NOVAS_REDUCED_ACCURACY, &planets, NULL))) n++;
+  if(check("grav_undo_planets:pos_app", -1, grav_undo_planets(NULL, po, &planets, out))) n++;
+  if(check("grav_undo_planets:pos_obs", -1, grav_undo_planets(p, NULL, &planets, out))) n++;
+  if(check("grav_undo_planets:planets", -1, grav_undo_planets(p, po, NULL, out))) n++;
+    if(check("grav_undo_planets:pos_src", -1, grav_undo_planets(p, po, &planets, NULL))) n++;
 
   planets.mask = 1 << NOVAS_SUN;
-  if(check("grav_undo_planets:converge", -1, grav_undo_planets(p, po, NOVAS_REDUCED_ACCURACY, &planets, out))) {
+  if(check("grav_undo_planets:converge", -1, grav_undo_planets(p, po, &planets, out))) {
     if(check("grav_undo_planets:converge:errno", ECANCELED, errno)) n++;
   }