diff --git a/setup.py b/setup.py
index 3288524c0e794a29b7316e6e453c4d7fe5e1a284..20d63fff244a5afb2a52ce5fa7896d163925c3ba 100644
--- a/setup.py
+++ b/setup.py
@@ -45,7 +45,7 @@ setup(name="pymbolic",
         'Topic :: Software Development :: Libraries',
         'Topic :: Utilities',
         ],
-      author=u"Andreas Kloeckner",
+      author="Andreas Kloeckner",
       author_email="inform@tiker.net",
       license = "MIT",
       url="http://mathema.tician.de/software/pymbolic",
diff --git a/test/test_maxima.py b/test/test_maxima.py
index 67a6e23175fed6933a67b1424e6ab74bb1614dcd..04a670c84efed3ded82002f82b8db7c39bfbcc37 100644
--- a/test/test_maxima.py
+++ b/test/test_maxima.py
@@ -75,13 +75,13 @@ def test_strict_round_trip(knl):
         result = knl.eval_expr(expr)
         round_trips_correctly = result == expr
         if not round_trips_correctly:
-            print "ORIGINAL:"
-            print
-            print expr
-            print
-            print "POST-MAXIMA:"
-            print
-            print result
+            print("ORIGINAL:")
+            print("")
+            print(expr)
+            print("")
+            print("POST-MAXIMA:")
+            print("")
+            print(result)
         assert round_trips_correctly
 
 def test_lax_round_trip(knl):
diff --git a/test/test_pymbolic.py b/test/test_pymbolic.py
index b9b9a80f1ef4667f24a61c87a4e8c407b84a5dbd..71e384822894a96a7989870857447ebbfd51cd9f 100644
--- a/test/test_pymbolic.py
+++ b/test/test_pymbolic.py
@@ -90,9 +90,9 @@ def test_fft():
     from pymbolic.algorithm import fft, sym_fft
 
     vars = numpy.array([var(chr(97+i)) for i in range(16)], dtype=object)
-    print vars
+    print(vars)
 
-    print fft(vars)
+    print(fft(vars))
     traced_fft = sym_fft(vars)
 
     from pymbolic.mapper.stringifier import PREC_NONE
@@ -102,10 +102,10 @@ def test_fft():
     code = [ccm(tfi, PREC_NONE) for tfi in traced_fft]
 
     for cse_name, cse_str in enumerate(ccm.cse_name_list):
-        print "%s = %s" % (cse_name, cse_str)
+        print("%s = %s" % (cse_name, cse_str))
 
     for i, line in enumerate(code):
-        print "result[%d] = %s" % (i, line)
+        print("result[%d] = %s" % (i, line))
 
 
 
@@ -157,26 +157,26 @@ def test_parser():
     parse("(2*a[1]*b[1]+2*a[0]*b[0])*(hankel_1(-1,sqrt(a[1]**2+a[0]**2)*k) "
             "-hankel_1(1,sqrt(a[1]**2+a[0]**2)*k))*k /(4*sqrt(a[1]**2+a[0]**2)) "
             "+hankel_1(0,sqrt(a[1]**2+a[0]**2)*k)")
-    print repr(parse("d4knl0"))
-    print repr(parse("0."))
-    print repr(parse("0.e1"))
-    print repr(parse("0.e1"))
-    print repr(parse("a >= 1"))
-    print repr(parse("a <= 1"))
-
-    print repr(parse(":"))
-    print repr(parse("1:"))
-    print repr(parse(":2"))
-    print repr(parse("1:2"))
-    print repr(parse("::"))
-    print repr(parse("1::"))
-    print repr(parse(":1:"))
-    print repr(parse("::1"))
-    print repr(parse("3::1"))
-    print repr(parse(":5:1"))
-    print repr(parse("3:5:1"))
-
-    print parse("3::1")
+    print(repr(parse("d4knl0")))
+    print(repr(parse("0.")))
+    print(repr(parse("0.e1")))
+    print(repr(parse("0.e1")))
+    print(repr(parse("a >= 1")))
+    print(repr(parse("a <= 1")))
+
+    print(repr(parse(":")))
+    print(repr(parse("1:")))
+    print(repr(parse(":2")))
+    print(repr(parse("1:2")))
+    print(repr(parse("::")))
+    print(repr(parse("1::")))
+    print(repr(parse(":1:")))
+    print(repr(parse("::1")))
+    print(repr(parse("3::1")))
+    print(repr(parse(":5:1")))
+    print(repr(parse("3:5:1")))
+
+    print(parse("3::1"))
 
     assert parse("e1") == prim.Variable("e1")
     assert parse("d1") == prim.Variable("d1")
@@ -312,11 +312,11 @@ def playground():
     vec5 = MV(np.random.randn(dims))
 
     a = vec3^vec4
-    print (a.rev()*a).project(0)
-    print a.scalar_product(a.rev())
+    print((a.rev()*a).project(0))
+    print(a.scalar_product(a.rev()))
 
-    #print a.norm_squared()
-    #print ((a.rev()*a).project(0) ).close_to( a.norm_squared() )
+    #print(a.norm_squared())
+    #print(((a.rev()*a).project(0) ).close_to( a.norm_squared() ))