diff --git a/.travis.yml b/.travis.yml
index 07719bba9..6eee7f742 100644
--- a/.travis.yml
+++ b/.travis.yml
@@ -39,11 +39,9 @@ matrix:
       env: OMP=n MPI=n PHDF5=n
     - python: "3.8"
       env: OMP=y MPI=n PHDF5=n
-    - python: "3.8"
-      env: OMP=n MPI=y PHDF5=n
     - python: "3.8"
       env: OMP=n MPI=y PHDF5=n VECTFIT=y
-    - python: "3.7"
+    - python: "3.8"
       env: OMP=n MPI=y PHDF5=y
     - python: "3.8"
       env: OMP=y MPI=y PHDF5=y DAGMC=y
diff --git a/openmc/data/multipole.py b/openmc/data/multipole.py
index e40b0d0d0..2dfc94ca2 100644
--- a/openmc/data/multipole.py
+++ b/openmc/data/multipole.py
@@ -20,8 +20,6 @@ from .data import K_BOLTZMANN
 from .neutron import IncidentNeutron
 from .resonance import ResonanceRange
 
-import vectfit as vf
-
 
 # Constants that determine which value to access
 _MP_EA = 0       # Pole
@@ -353,7 +351,7 @@ def _vectfit_xs(energy, ce_xs, mts, rtol=1e-3, atol=1e-5, orders=None,
     if log:
         print("Found {} real poles and {} conjugate complex pairs.".format(
                len(real_idx), len(conj_idx)))
-    mp_poles = best_poles[real_idx+conj_idx]
+    mp_poles = best_poles[real_idx + conj_idx]
     mp_residues = np.concatenate((best_residues[:, real_idx],
                                   best_residues[:, conj_idx]*2), axis=1)/1j
     if log: