diff --git a/.github/workflows/pythonapp.yml b/.github/workflows/pythonapp.yml
index 8bf19874..9bf7cd87 100644
--- a/.github/workflows/pythonapp.yml
+++ b/.github/workflows/pythonapp.yml
@@ -14,10 +14,10 @@ jobs:
 
     steps:
     - uses: actions/checkout@v1
-    - name: Set up Python 3.7
+    - name: Set up Python 3.10.16
       uses: actions/setup-python@v1
       with:
-        python-version: 3.7
+        python-version: 3.10.16
     - name: Install dependencies
       run: |
         python -m pip install --upgrade pip
diff --git a/cellconstructor/ForceTensor.py b/cellconstructor/ForceTensor.py
index 97f9aea9..97f2ce9b 100644
--- a/cellconstructor/ForceTensor.py
+++ b/cellconstructor/ForceTensor.py
@@ -408,7 +408,8 @@ def Center(self, nneigh=None, Far=2,tol=1.0e-5):
             nq2=self.supercell_size[2]
             
             # by default no maximum distances
-            dist_range=np.ones((self.nat))*np.infty
+            #dist_range=np.ones((self.nat))*np.infty
+            dist_range=np.ones((self.nat))*np.inf
             #
             if nneigh!=None:
             # =======================================================================
@@ -1520,7 +1521,8 @@ def Center(self, nneigh=None, Far=2,tol=1.0e-5):
             nq2=self.supercell_size[2]
             
             # by default no maximum distances
-            dist_range=np.ones((self.nat))*np.infty
+            #dist_range=np.ones((self.nat))*np.infty
+            dist_range=np.ones((self.nat))*np.inf
             #
             if nneigh!=None:
             # =======================================================================
diff --git a/cellconstructor/Spectral.py b/cellconstructor/Spectral.py
index 2304c4a2..488daeae 100644
--- a/cellconstructor/Spectral.py
+++ b/cellconstructor/Spectral.py
@@ -1945,7 +1945,8 @@ def findne(val,e0,de):
                     else:
                         xtriv=findne(wq[ifreq],e0,de)
                         osval=np.real(z[xtriv-1,ism,ifreq])
-                        diff=np.infty
+                        #diff=np.infty
+                        diff=np.inf
                         for x in range(ne):
                             value=np.real(z[x,ism,ifreq])-energies_cm[x]
                             if( abs(value ) < 2*de) :
@@ -2661,7 +2662,8 @@ def findne(val,e0,de):
                     else:
                         xtriv=findne(wq[ifreq],e0,de)
                         osval=np.real(z[xtriv-1,ism,ifreq])
-                        diff=np.infty
+                        #diff=np.infty
+                        diff=np.inf
                         for x in range(ne):
                             value=np.real(z[x,ism,ifreq])-energies_cm[x]
                             if( abs(value ) < 2*de) :