diff --git a/CHANGELOG.md b/CHANGELOG.md
index ca7fc97981db34eb884545f780f752ef3f6c6f3c..92b3554c23ae2e13bd4d2fde1453208d64323eb6 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -16,6 +16,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 ### Fixed
 
 + Avoided using deprecated functionality in latest scipy, numpy, and sympy
++ Fix syntax for python 3.9
 
 ### Removed
 
diff --git a/environment-latest.yml b/environment-latest.yml
index 4ed09bbedf78c25e5ea2c6211abe8210e517c369..33dff4c2f684e1f74808318a0806f920c958d686 100644
--- a/environment-latest.yml
+++ b/environment-latest.yml
@@ -2,7 +2,7 @@ name: qsymm-latest
 channels:
     - conda-forge
 dependencies:
-    - python=3.8
+    - python=3.9
     - numpy
     - scipy
     - sympy
diff --git a/pytest.ini b/pytest.ini
index 8e76f5a2b2fce71f20c31b9ff2944b232aa9fdd4..64b068597c69f9c8ca36e3d199e20884c0187070 100644
--- a/pytest.ini
+++ b/pytest.ini
@@ -6,3 +6,5 @@ filterwarnings =
     ignore:the matrix subclass is not the recommended way to represent matrices:PendingDeprecationWarning
     # Remove once we rely on sympy v1.3
     ignore:np\.asscalar\(a\) is deprecated since NumPy v1\.16, use a\.item\(\) instead:DeprecationWarning
+    # Remove once setuptools figures it out
+    ignore:lib2to3:PendingDeprecationWarning
diff --git a/qsymm/groups.py b/qsymm/groups.py
index 1ebeae1291e6415fc11dc78f34cfad3c8efc7e22..b34de1315528a7963ee7defa87449b8d338150c9 100644
--- a/qsymm/groups.py
+++ b/qsymm/groups.py
@@ -5,7 +5,6 @@ from functools import lru_cache
 from fractions import Fraction
 from numbers import Number
 from collections import OrderedDict
-from copy import deepcopy
 
 import numpy as np
 import tinyarray as ta
diff --git a/qsymm/tests/test_symmetry_finder.py b/qsymm/tests/test_symmetry_finder.py
index 3a57d3a0f8ef6e7f6653ce16f8b9b08d73cb857c..36aa2a01ab0d1dd84ce85ed277f8037da0bb58d2 100644
--- a/qsymm/tests/test_symmetry_finder.py
+++ b/qsymm/tests/test_symmetry_finder.py
@@ -29,7 +29,7 @@ def sumrep(*args):
 
 def test_cont_finder():
     # Test symmetry adapted basis
-    gens = sumrep((*2*[0.5*sigma[[3, 1, 2]]]))
+    gens = sumrep(*2*[0.5*sigma[[3, 1, 2]]])
     U2 = kwant_rmt.circular(len(gens[0]))
     gens2 = np.einsum('ij,ajk,kl->ail',(U2.T).conjugate(),gens,U2)
     U = symmetry_adapted_sun(gens, check=True)