SWIFTSIM · AlexKurek · Jul 20, 2024 · Aug 7, 2024
diff --git a/examples/Cooling/CoolingHalo/density_profile.py b/examples/Cooling/CoolingHalo/density_profile.py
@@ -1,7 +1,6 @@
 import numpy as np
 import h5py as h5
 import matplotlib.pyplot as plt
-import sys
 
 n_snaps = 11
 

diff --git a/examples/Cooling/CoolingHalo/makeIC.py b/examples/Cooling/CoolingHalo/makeIC.py
@@ -21,7 +21,6 @@
 import sys
 import numpy as np
 import math
-import random
 
 # Generates N particles in a spherically symmetric distribution with density profile ~r^(-2)
 # usage: python3 makeIC.py 1000: generate 1000 particles

diff --git a/examples/Cooling/CoolingHalo/makeIC_random_box.py b/examples/Cooling/CoolingHalo/makeIC_random_box.py
@@ -20,8 +20,6 @@
 import h5py
 import sys
 import numpy as np
-import math
-import random
 
 # Generates N particles in a spherically symmetric distribution with density profile ~r^(-2)
 # usage: python3 makeIC.py 1000: generate 1000 particles

diff --git a/examples/Cooling/CoolingHalo/test_energy_conservation.py b/examples/Cooling/CoolingHalo/test_energy_conservation.py
@@ -1,7 +1,6 @@
 import numpy as np
 import h5py as h5
 import matplotlib.pyplot as plt
-import sys
 
 n_snaps = 41
 

diff --git a/examples/Cooling/CoolingHaloWithSpin/makeIC.py b/examples/Cooling/CoolingHaloWithSpin/makeIC.py
@@ -21,7 +21,6 @@
 import sys
 import numpy as np
 import math
-import random
 
 # Generates N particles in a spherically symmetric distribution with density profile ~r^(-2)
 # usage: python3 makeIC.py 1000: generate 1000 particles

diff --git a/examples/Cooling/CoolingHaloWithSpin/makeIC_random_box.py b/examples/Cooling/CoolingHaloWithSpin/makeIC_random_box.py
@@ -20,8 +20,6 @@
 import h5py
 import sys
 import numpy as np
-import math
-import random
 
 # Generates N particles in a spherically symmetric distribution with density profile ~r^(-2)
 # usage: python3 makeIC.py 1000: generate 1000 particles

diff --git a/examples/Cooling/CoolingHaloWithSpin/test_energy_conservation.py b/examples/Cooling/CoolingHaloWithSpin/test_energy_conservation.py
@@ -1,7 +1,6 @@
 import numpy as np
 import h5py as h5
 import matplotlib.pyplot as plt
-import sys
 import glob
 
 # Get the total number of snapshots

diff --git a/examples/Cooling/CoolingRates/plot_cooling_rates.py b/examples/Cooling/CoolingRates/plot_cooling_rates.py
@@ -3,7 +3,6 @@
 # by testCooling.
 
 import matplotlib.pyplot as plt
-import numpy as np
 
 # Number of metals tracked by EAGLE cooling
 elements = 11

diff --git a/examples/Cooling/FeedbackEvent_3D/makeIC.py b/examples/Cooling/FeedbackEvent_3D/makeIC.py
@@ -20,7 +20,7 @@
 from swiftsimio import Writer
 from swiftsimio.units import cosmo_units
 
-from unyt import cm, kpc, mh, msun, K, s, kb
+from unyt import cm, mh, msun, K, s, kb
 
 import h5py
 import numpy as np

diff --git a/examples/Cooling/FeedbackEvent_3D/plotEnergy.py b/examples/Cooling/FeedbackEvent_3D/plotEnergy.py
@@ -11,7 +11,6 @@
 
 from makeIC import (
     gamma,
-    initial_density,
     initial_temperature,
     inject_temperature,
     mu,

diff --git a/examples/Cooling/FeedbackEvent_3D/plotEnergyAll.py b/examples/Cooling/FeedbackEvent_3D/plotEnergyAll.py
@@ -7,12 +7,11 @@
 
 from swiftsimio import load
 
-from unyt import Gyr, erg, mh, kb, Myr
+from unyt import erg, mh, kb, Myr
 from scipy.interpolate import interp1d
 
 from makeIC import (
     gamma,
-    initial_density,
     initial_temperature,
     inject_temperature,
     mu,

diff --git a/examples/GEAR/AgoraCosmo/convert_ic.py b/examples/GEAR/AgoraCosmo/convert_ic.py
@@ -3,11 +3,9 @@
 from h5py import File
 import shutil
 import numpy as np
-import sys
 from swiftsimio.visualisation.smoothing_length_generation import (
     generate_smoothing_lengths,
 )
-import time
 from unyt import kpc
 
 filename = "agora_arepo.hdf5"

diff --git a/examples/GravityTests/DiscPatch/HydroStatic/makeIC.py b/examples/GravityTests/DiscPatch/HydroStatic/makeIC.py
@@ -21,10 +21,8 @@
 ##############################################################################
 
 import h5py
-import sys
 import numpy as np
 import math
-import random
 
 # Generates a disc-patch in hydrostatic equilibrium
 #

diff --git a/examples/GravityTests/DiscPatch/HydroStatic_1D/makeIC.py b/examples/GravityTests/DiscPatch/HydroStatic_1D/makeIC.py
@@ -21,10 +21,8 @@
 ##############################################################################
 
 import h5py
-import sys
 import numpy as np
 import math
-import random
 
 # Generates a disc-patch in hydrostatic equilibrium
 #

diff --git a/examples/GravityTests/ExternalPointMass/energy_plot.py b/examples/GravityTests/ExternalPointMass/energy_plot.py
@@ -2,7 +2,6 @@
 
 matplotlib.use("Agg")
 from pylab import *
-import h5py
 
 # Plot parameters
 params = {

diff --git a/examples/GravityTests/ExternalPointMass/makeIC.py b/examples/GravityTests/ExternalPointMass/makeIC.py
@@ -22,7 +22,6 @@
 import sys
 import numpy
 import math
-import random
 
 # Generates a random distriution of particles, for motion in an external potential centred at (0,0,0)
 

diff --git a/examples/GravityTests/Hernquist_circularorbit/makeIC.py b/examples/GravityTests/Hernquist_circularorbit/makeIC.py
@@ -17,8 +17,6 @@
 #
 ################################################################################
 import numpy as np
-import matplotlib.pyplot as plt
-from astropy import units
 import h5py as h5
 
 C = 8.0

diff --git a/examples/GravityTests/Hernquist_circularorbit/plotprog.py b/examples/GravityTests/Hernquist_circularorbit/plotprog.py
@@ -20,7 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-from scipy.integrate import odeint
 
 t = np.linspace(0, 40, 100000)
 y0 = [0, 10]

diff --git a/examples/GravityTests/Hernquist_radialinfall/makeIC.py b/examples/GravityTests/Hernquist_radialinfall/makeIC.py
@@ -18,10 +18,8 @@
 ##############################################################################
 
 import h5py
-import sys
 import numpy
 import math
-import random
 import numpy as np
 
 # Generates N particles in a spherical distribution centred on [0,0,0], to be moved in an isothermal potential

diff --git a/examples/GravityTests/Hernquist_radialinfall/plotprog.py b/examples/GravityTests/Hernquist_radialinfall/plotprog.py
@@ -20,7 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-from scipy.integrate import odeint
 
 
 lengthrun = 2001

diff --git a/examples/GravityTests/HydrostaticHalo/makeIC.py b/examples/GravityTests/HydrostaticHalo/makeIC.py
@@ -21,7 +21,6 @@
 import sys
 import numpy as np
 import math
-import random
 
 # Generates N particles in a spherically symmetric distribution with density profile ~r^(-2)
 # usage: python3 makeIC.py 1000: generate 1000 particles

diff --git a/examples/GravityTests/IsothermalPotential/energy_plot.py b/examples/GravityTests/IsothermalPotential/energy_plot.py
@@ -2,7 +2,6 @@
 
 matplotlib.use("Agg")
 from pylab import *
-import h5py
 
 # Plot parameters
 params = {

diff --git a/examples/GravityTests/IsothermalPotential/makeIC.py b/examples/GravityTests/IsothermalPotential/makeIC.py
@@ -22,7 +22,6 @@
 import sys
 import numpy
 import math
-import random
 
 # Generates N particles in a spherical distribution centred on [0,0,0], to be moved in an isothermal potential
 # usage: python3 makeIC.py 1000 0 : generate 1000 particles on circular orbits

diff --git a/examples/GravityTests/NFW_Halo/makeIC.py b/examples/GravityTests/NFW_Halo/makeIC.py
@@ -19,8 +19,6 @@
 ################################################################################
 
 import numpy as np
-import matplotlib.pyplot as plt
-from astropy import units
 import h5py as h5
 
 C = 8.0

diff --git a/examples/GravityTests/Plummer_Selfgravitating/plotdensity.py b/examples/GravityTests/Plummer_Selfgravitating/plotdensity.py
@@ -1,7 +1,6 @@
 import numpy as np
 from matplotlib import pyplot as plt
 import argparse
-from scipy.optimize import curve_fit
 import h5py
 
 # Parse user input

diff --git a/examples/HydroTests/BlobTest_3D/makeIC.py b/examples/HydroTests/BlobTest_3D/makeIC.py
@@ -3,7 +3,6 @@
 """
 
 import numpy as np
-import h5py
 
 from unyt import cm, g, s, erg
 from unyt.unit_systems import cgs_unit_system

diff --git a/examples/HydroTests/BlobTest_Braspenning2022/plot_cloud_evolution.py b/examples/HydroTests/BlobTest_Braspenning2022/plot_cloud_evolution.py
@@ -3,9 +3,7 @@
 import numpy as np
 import matplotlib
 import matplotlib.pyplot as plt
-import os
-from unyt import s, pc, Myr, kb, mp, cm, K, erg, g, G, km
-import h5py
+from unyt import s, kb, mp
 import glob
 
 

diff --git a/examples/HydroTests/EvrardCollapse_3D/plotEnergy.py b/examples/HydroTests/EvrardCollapse_3D/plotEnergy.py
@@ -22,8 +22,6 @@
 matplotlib.use("Agg")
 import matplotlib.pyplot as plt
 import numpy as np
-from scipy import stats
-import sys
 
 plt.style.use("../../../tools/stylesheets/mnras.mplstyle")
 

diff --git a/examples/HydroTests/Gradients/makeICs.py b/examples/HydroTests/Gradients/makeICs.py
@@ -18,7 +18,6 @@
 ################################################################################
 
 import h5py
-import random
 import numpy as np
 import sys
 

diff --git a/examples/HydroTests/Gradients/plot.py b/examples/HydroTests/Gradients/plot.py
@@ -17,7 +17,6 @@
 #
 ################################################################################
 
-import scipy as sp
 import pylab as pl
 import numpy as np
 import h5py

diff --git a/examples/HydroTests/KelvinHelmholtz_2D/makeIC.py b/examples/HydroTests/KelvinHelmholtz_2D/makeIC.py
@@ -19,7 +19,6 @@
 
 import h5py
 from numpy import *
-import sys
 
 # Generates a swift IC file for the Kelvin-Helmholtz vortex in a periodic box
 

diff --git a/examples/HydroTests/KelvinHelmholtz_2D/makeMovieSwiftsimIO.py b/examples/HydroTests/KelvinHelmholtz_2D/makeMovieSwiftsimIO.py
@@ -10,9 +10,7 @@
 """
 
 import os
-import h5py as h5
 import numpy as np
-import scipy.interpolate as si
 
 from swiftsimio import load
 from swiftsimio.visualisation import project_gas_pixel_grid

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_NFW_MN/generate_MN.py b/examples/IsolatedGalaxy/IsolatedGalaxy_NFW_MN/generate_MN.py
@@ -19,7 +19,6 @@
 
 import numpy as np
 from scipy import special
-import matplotlib.pyplot as plt
 import h5py
 
 G = 4.299e-6

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_dmparticles/angularmomentum.py b/examples/IsolatedGalaxy/IsolatedGalaxy_dmparticles/angularmomentum.py
@@ -20,7 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-import scipy.optimize as sco
 
 
 Nmax = 2001

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_dmparticles/profilefit.py b/examples/IsolatedGalaxy/IsolatedGalaxy_dmparticles/profilefit.py
@@ -20,8 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-from matplotlib.colors import BoundaryNorm
-from matplotlib.ticker import MaxNLocator
 import scipy.optimize as sco
 import os
 

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_feedback/SFH.py b/examples/IsolatedGalaxy/IsolatedGalaxy_feedback/SFH.py
@@ -9,7 +9,6 @@
 Written by James Willis (james.s.willis@durham.ac.uk)
 """
 
-import glob
 import h5py as h5
 import numpy as np
 import matplotlib.pyplot as plt

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_feedback/plot_box_evolution.py b/examples/IsolatedGalaxy/IsolatedGalaxy_feedback/plot_box_evolution.py
@@ -22,7 +22,6 @@
 
 matplotlib.use("Agg")
 from pylab import *
-from scipy import stats
 import h5py
 import numpy as np
 import glob

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_potential/angularmomentum.py b/examples/IsolatedGalaxy/IsolatedGalaxy_potential/angularmomentum.py
@@ -20,7 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-import scipy.optimize as sco
 
 
 Nmax = 2001

diff --git a/examples/IsolatedGalaxy/IsolatedGalaxy_potential/profilefit.py b/examples/IsolatedGalaxy/IsolatedGalaxy_potential/profilefit.py
@@ -20,8 +20,6 @@
 import numpy as np
 import h5py
 import matplotlib.pyplot as plt
-from matplotlib.colors import BoundaryNorm
-from matplotlib.ticker import MaxNLocator
 import scipy.optimize as sco
 import os
 

diff --git a/tools/analyse_dump_cells.py b/tools/analyse_dump_cells.py
@@ -31,7 +31,6 @@
 import pylab as pl
 import numpy as np
 import sys
-import pandas
 
 xcol = 0
 ycol = 1

diff --git a/tools/convert_snapshot_to_ICs.py b/tools/convert_snapshot_to_ICs.py
@@ -10,7 +10,6 @@
 """
 
 from swiftsimio import Writer, load
-import numpy as np
 import unyt
 import os
 

diff --git a/tools/parallel_replicate_ICs.py b/tools/parallel_replicate_ICs.py
@@ -39,18 +39,13 @@
 """
 
 import h5py as h
-import numpy as np
 import matplotlib
 
 matplotlib.use("Agg")
 from pylab import *
 import os.path
-from tqdm import tqdm
-from tqdm import trange
 import time
 from numba import jit, prange
-from swiftsimio import Writer
-from swiftsimio.units import cosmo_units
 
 replicate = 1
 box_size = 1

diff --git a/tools/plot_gravity_checks.py b/tools/plot_gravity_checks.py
@@ -2,7 +2,6 @@
 
 import sys
 import glob
-import re
 import numpy as np
 import matplotlib.pyplot as plt
 

diff --git a/tools/plot_scaling_results.py b/tools/plot_scaling_results.py
@@ -16,7 +16,6 @@
 import re
 import numpy as np
 import matplotlib.pyplot as plt
-import scipy.stats
 import ntpath
 
 params = {