diff --git a/examples/fmm-error.py b/examples/fmm-error.py
index ea50d70f66b4c842c0404ab326f3e02c2cf3eccd..4f766ff0be5b408982fe5f345bdf7c210d7579f1 100644
--- a/examples/fmm-error.py
+++ b/examples/fmm-error.py
@@ -61,7 +61,7 @@ def main():
     density_discr = places.get_discretization("unaccel_qbx")
 
     nodes = thaw(actx, density_discr.nodes())
-    angle = actx.np.atan2(nodes[1], nodes[0])
+    angle = actx.np.arctan2(nodes[1], nodes[0])
 
     from pytential import bind, sym
     if k:
diff --git a/examples/layerpot-3d.py b/examples/layerpot-3d.py
index 78112858c16890ef13720d5305839d62e7f6cd8b..75fb3412649f2ec0a32b8b84a6dc6b3d8fec6a41 100644
--- a/examples/layerpot-3d.py
+++ b/examples/layerpot-3d.py
@@ -73,7 +73,7 @@ def main(mesh_name="ellipsoid"):
     density_discr = places.get_discretization("qbx")
 
     nodes = thaw(actx, density_discr.nodes())
-    angle = actx.np.atan2(nodes[1], nodes[0])
+    angle = actx.np.arctan2(nodes[1], nodes[0])
 
     if k:
         kernel = HelmholtzKernel(3)
diff --git a/examples/layerpot.py b/examples/layerpot.py
index 40a060fb4a2e52aeab7bf20cfce3ccf73cf20bc6..605faf177cdeacda07d0d708a1bf672c4fcd5fab 100644
--- a/examples/layerpot.py
+++ b/examples/layerpot.py
@@ -64,7 +64,7 @@ def main(curve_fn=starfish, visualize=True):
 
     from meshmode.dof_array import thaw
     nodes = thaw(actx, density_discr.nodes())
-    angle = actx.np.atan2(nodes[1], nodes[0])
+    angle = actx.np.arctan2(nodes[1], nodes[0])
 
     if k:
         kernel = HelmholtzKernel(2)
diff --git a/examples/scaling-study.py b/examples/scaling-study.py
index a14a5aed0b41085379c79f8cd90684ba4ec65daa..8f3b4639eeeac5a12ad05036940e847710a73803 100644
--- a/examples/scaling-study.py
+++ b/examples/scaling-study.py
@@ -124,7 +124,7 @@ def timing_run(nx, ny, visualize=False):
 
     from meshmode.dof_array import thaw
     nodes = thaw(actx, density_discr.nodes())
-    angle = actx.np.atan2(nodes[1], nodes[0])
+    angle = actx.np.arctan2(nodes[1], nodes[0])
 
     sigma = actx.np.cos(mode_nr*angle)
 
diff --git a/test/test_layer_pot_eigenvalues.py b/test/test_layer_pot_eigenvalues.py
index b579791307cdbc55298618942bfe6977f276c1b5..68baf63396b329914b4e7cba87c4abfabdfb1f6d 100644
--- a/test/test_layer_pot_eigenvalues.py
+++ b/test/test_layer_pot_eigenvalues.py
@@ -138,7 +138,7 @@ def test_ellipse_eigenvalues(ctx_factory, ellipse_aspect, mode_nr, qbx_order,
             pt.gca().set_aspect("equal")
             pt.show()
 
-        angle = actx.np.atan2(nodes[1]*ellipse_aspect, nodes[0])
+        angle = actx.np.arctan2(nodes[1]*ellipse_aspect, nodes[0])
 
         ellipse_fraction = ((1-ellipse_aspect)/(1+ellipse_aspect))**mode_nr
 
@@ -310,7 +310,7 @@ def test_sphere_eigenvalues(ctx_factory, mode_m, mode_n, qbx_order,
         nodes = thaw(actx, density_discr.nodes())
         r = actx.np.sqrt(nodes[0]*nodes[0] + nodes[1]*nodes[1] + nodes[2]*nodes[2])
         phi = actx.np.acos(nodes[2]/r)
-        theta = actx.np.atan2(nodes[0], nodes[1])
+        theta = actx.np.arctan2(nodes[0], nodes[1])
 
         ymn = unflatten(actx, density_discr,
                 actx.from_numpy(
diff --git a/test/test_symbolic.py b/test/test_symbolic.py
index a510596cdc9f9d56c2f974a7a64666fe36e5cea8..c5f8fbc1f16e0192995f44b090531b638a96901b 100644
--- a/test/test_symbolic.py
+++ b/test/test_symbolic.py
@@ -67,7 +67,7 @@ def get_ellipse_with_ref_mean_curvature(actx, nelements, aspect=1):
 
     a = 1
     b = 1/aspect
-    t = actx.np.atan2(nodes[1] * aspect, nodes[0])
+    t = actx.np.arctan2(nodes[1] * aspect, nodes[0])
 
     return discr, a*b / ((a*actx.np.sin(t))**2 + (b*actx.np.cos(t))**2)**(3/2)
 
@@ -90,11 +90,11 @@ def get_torus_with_ref_mean_curvature(actx, h):
     a = r_major
     b = r_minor
 
-    u = actx.np.atan2(nodes[1], nodes[0])
+    u = actx.np.arctan2(nodes[1], nodes[0])
     from pytools.obj_array import flat_obj_array
     rvec = flat_obj_array(actx.np.cos(u), actx.np.sin(u), 0*u)
     rvec = sum(nodes * rvec) - a
-    cosv = actx.np.cos(actx.np.atan2(nodes[2], rvec))
+    cosv = actx.np.cos(actx.np.arctan2(nodes[2], rvec))
 
     return discr, (a + 2.0 * b * cosv) / (2 * b * (a + b * cosv))