Fix coverage from not including starlink in test suite anymore

Author: rmjarvis

galsim/table.py

@@ -951,7 +951,7 @@ def __call__(self, x, y, grid=False):
         Returns:
             a scalar value if x and y are scalar, or a numpy array if x and y are arrays.
         """
-        if np.__version__ >= "2.0":
+        if np.__version__ >= "2.0":  # pragma: no branch
             # I'm not sure if there is a simpler way to do this in 2.0.
             # We want a copy when edge_mode == wrap. Otherwise, only copy if dtype changes or
             # x,y aren't already arrays. That used to be simple...
@@ -1071,7 +1071,7 @@ def gradient(self, x, y, grid=False):
             A tuple of (dfdx, dfdy) where dfdx, dfdy are single values (if x,y were single
             values) or numpy arrays.
         """
-        if np.__version__ >= "2.0":
+        if np.__version__ >= "2.0":  # pragma: no branch
             copy = True if self.edge_mode=='wrap' else None
         else:
             copy = self.edge_mode=='wrap'

tests/test_fitsheader.py

@@ -75,6 +75,11 @@ def check_tpv(header):
     header = galsim.FitsHeader(file_name=file_name, dir=dir, hdu=0)
     check_tpv(header)
     check_pickle(header)
+    # Can also pass the hdu itself.
+    with pyfits.open(os.path.join(dir,file_name)) as hdu_list:
+        header = galsim.FitsHeader(hdu_list=hdu_list[0])
+    check_tpv(header)
+    check_pickle(header)
     # If you pass in a pyfits Header object, that should also work
     with pyfits.open(os.path.join(dir,file_name)) as hdu_list:
         header = galsim.FitsHeader(header=hdu_list[0].header)