another test

Author: Jammy2211

autoarray/inversion/pixelization/border_relocator.py

@@ -287,34 +287,31 @@ def relocated_grid_via_ellipse_border_from(grid, origin, a, b, phi, xp=np, eps=1
     # shift to origin
     dy = grid[:, 0] - origin[0]
     dx = grid[:, 1] - origin[1]
-
-    c = xp.cos(phi)
-    s = xp.sin(phi)
-
-    # rotate into ellipse-aligned frame
-    xprime = c * dx + s * dy
-    yprime = -s * dx + c * dy
-
-    # ellipse radius in normalized coords
-    q = (xprime / a) ** 2 + (yprime / b) ** 2
+    #
+    # c = xp.cos(phi)
+    # s = xp.sin(phi)
+    #
+    # # rotate into ellipse-aligned frame
+    # xprime = c * dx + s * dy
+    # yprime = -s * dx + c * dy
+    #
+    # # ellipse radius in normalized coords
+    # q = (xprime / a) ** 2 + (yprime / b) ** 2
 
     outside = q > 1.0
-    scale = 1.0 / xp.sqrt(xp.maximum(q, 1.0 + eps))
-
-    # scale back to boundary
-    xprime2 = xprime * scale
-    yprime2 = yprime * scale
-
-    # rotate back to original frame
-    dx2 = c * xprime2 - s * yprime2
-    dy2 = s * xprime2 + c * yprime2
-
-    moved = xp.stack([origin[0] + dy2, origin[1] + dx2], axis=1)
-
-#    return xp.where(outside[:, None], moved, grid)
-
-    return grid + (moved - grid) * outside.astype(grid.dtype)[:, None]
-
+    # scale = 1.0 / xp.sqrt(xp.maximum(q, 1.0 + eps))
+    #
+    # # scale back to boundary
+    # xprime2 = xprime * scale
+    # yprime2 = yprime * scale
+    #
+    # # rotate back to original frame
+    # dx2 = c * xprime2 - s * yprime2
+    # dy2 = s * xprime2 + c * yprime2
+
+    moved = xp.stack([origin[0] + dy, origin[1] + dx], axis=1)
+
+    return xp.where(outside[:, None], moved, grid)
 
 
 class BorderRelocator: