From ba3c0614f869038791ad021887f4794d39c5b72a Mon Sep 17 00:00:00 2001 From: Andreas Kloeckner Date: Mon, 13 Jul 2020 15:38:36 -0500 Subject: [PATCH] Replace actx.np.atan2 with arctan2 --- examples/fmm-error.py | 2 +- examples/layerpot-3d.py | 2 +- examples/layerpot.py | 2 +- examples/scaling-study.py | 2 +- test/test_layer_pot_eigenvalues.py | 4 ++-- test/test_symbolic.py | 6 +++--- 6 files changed, 9 insertions(+), 9 deletions(-) diff --git a/examples/fmm-error.py b/examples/fmm-error.py index ea50d70f..4f766ff0 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 78112858..75fb3412 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 40a060fb..605faf17 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 a14a5aed..8f3b4639 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 b5797913..68baf633 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 a510596c..c5f8fbc1 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)) -- GitLab