Check for NULL jd_tdb[2] in solarsystem_hp() variants.

Author: attipaci

src/solsys1.c

@@ -172,6 +172,10 @@ short solarsystem(double jd_tdb, short body, short origin, double *position, dou
 }
 
 short solarsystem_hp(const double jd_tdb[2], short body, short origin, double *position, double *velocity) {
+  if(!jd_tdb) {
+    errno = EINVAL;
+    return -1;
+  }
   return planet_eph_manager_hp(jd_tdb, body, origin, position, velocity);
 }
 #endif

src/solsys2.c

@@ -176,6 +176,11 @@ short planet_jplint_hp(const double jd_tdb[2], enum novas_planet body, enum nova
   double posvel[6] = { };
   int i;
 
+  if(!jd_tdb) {
+    errno = EINVAL;
+    return -1;
+  }
+
   // Perform sanity checks on the input body and origin.
   if((body < NOVAS_MERCURY) || (body > NOVAS_MOON)) {
     errno = EINVAL;

src/solsys3.c

@@ -292,6 +292,11 @@ short earth_sun_calc(double jd_tdb, enum novas_planet body, enum novas_origin or
 short earth_sun_calc_hp(const double jd_tdb[2], enum novas_planet body, enum novas_origin origin, double *position,
         double *velocity) {
 
+  if(!jd_tdb) {
+    errno = EINVAL;
+    return -1;
+  }
+
   int error = earth_sun_calc(jd_tdb[0] + jd_tdb[1], body, origin, position, velocity);
   if(error)
     return error;

test/src/test-errors.c

@@ -775,6 +775,22 @@ static int test_earth_sun_calc() {
   return n;
 }
 
+static int test_earth_sun_calc_hp() {
+  double p[3], v[3], tdb2[2] = { NOVAS_JD_J2000 };
+  int n = 0;
+
+  enable_earth_sun_hp(1);
+
+  if(check("earth_sun_calc_hp:tdb", -1, earth_sun_calc_hp(NULL, NOVAS_SUN, NOVAS_BARYCENTER, p, v))) n++;
+  if(check("earth_sun_calc_hp:pos", -1, earth_sun_calc_hp(tdb2, NOVAS_SUN, NOVAS_BARYCENTER, NULL, v))) n++;
+  if(check("earth_sun_calc_hp:vel", -1, earth_sun_calc_hp(tdb2, NOVAS_SUN, NOVAS_BARYCENTER, p, NULL))) n++;
+  if(check("earth_sun_calc_hp:number", 2, earth_sun_calc_hp(tdb2, NOVAS_JUPITER, NOVAS_BARYCENTER, p, v))) n++;
+
+  enable_earth_sun_hp(0);
+
+  return n;
+}
+
 static int test_sun_eph() {
   extern int sun_eph(double jd, double *ra, double *dec, double *dis);
 
@@ -873,6 +889,7 @@ int main() {
   if(test_grav_def()) n++;
 
   if(test_earth_sun_calc()) n++;
+  if(test_earth_sun_calc_hp()) n++;
   if(test_sun_eph()) n++;
 
   if(n) fprintf(stderr, " -- FAILED %d tests\n", n);