--- a/interface/tests/python/check_mixed_mesh.py
+++ b/interface/tests/python/check_mixed_mesh.py
@@ -96,12 +96,12 @@
         #print 'index of HEXA: ', i
         listHexa.append(i)
     else:
-        print 'Geometric transformation: ', gt[0]
+        print( 'Geometric transformation: ', gt[0])
 
-print 'num Tetra: ', len(listTetra)
-print 'num Hexa: ', len(listHexa)
-print 'num Hexa: ', len(listHexa)
-print 'num Hexa: ', len(listHexa)
+print( 'num Tetra: ', len(listTetra))
+print( 'num Hexa: ', len(listHexa))
+print( 'num Hexa: ', len(listHexa))
+print( 'num Hexa: ', len(listHexa))
 
 
 mfu = gf.MeshFem(m, 3)
@@ -172,7 +172,7 @@
 
 ETA3 = ETA3tmp [ ETA3tmp.size - mfer.nbdof() : ETA3tmp.size ]
 
-print 'sum(ETA3): ', sum(ETA3)
+print( 'sum(ETA3): ', sum(ETA3))
 
 # 1d) jump at NEUMANN_BOUNDARY_NO_LOAD
 g = "[0; 0; 0]"
@@ -185,7 +185,7 @@
 
 ETA4 = ETA4tmp [ ETA4tmp.size - mfer.nbdof() : ETA4tmp.size ]
 
-print 'sum(ETA4): ', sum(ETA4)
+print( 'sum(ETA4): ', sum(ETA4))
 
 
 ETA_square = ETA1 + ETA2 + ETA3 + ETA4 # element wise