Merge remote-tracking branch 'origin/main'

Tests/alpha_conversion_test.ipynb

@@ -1,30 +1,21 @@
 {
  "cells": [
   {
+   "metadata": {},
    "cell_type": "code",
-   "id": "initial_id",
-   "metadata": {
-    "collapsed": true,
-    "ExecuteTime": {
-     "end_time": "2026-01-20T15:30:31.978159Z",
-     "start_time": "2026-01-20T15:30:31.835157Z"
-    }
-   },
+   "outputs": [],
+   "execution_count": null,
    "source": [
     "%load_ext autoreload\n",
     "%autoreload 2"
    ],
-   "outputs": [],
-   "execution_count": 1
+   "id": "a78faf7f4883772f"
   },
   {
-   "metadata": {
-    "ExecuteTime": {
-     "end_time": "2026-01-20T15:30:33.910807Z",
-     "start_time": "2026-01-20T15:30:32.803089Z"
-    }
-   },
+   "metadata": {},
    "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
    "source": [
     "%reload_ext autoreload\n",
     "%autoreload 2\n",
@@ -33,9 +24,56 @@
     "from GHA_triaxial.utils import alpha_para2ell, alpha_ell2para\n",
     "import numpy as np"
    ],
-   "id": "9ad815aea55574e3",
+   "id": "46aa84a937fea491"
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
    "outputs": [],
-   "execution_count": 2
+   "execution_count": null,
+   "source": [
+    "ell = EllipsoidTriaxial.init_name(\"KarneyTest2024\")\n",
+    "diffs = []\n",
+    "for beta_deg in range(-180, 181, 45):\n",
+    "    for lamb_deg in range(-90, 91, 45):\n",
+    "        for alpha_deg in range(0, 360, 45):\n",
+    "            beta = wu.deg2rad(beta_deg)\n",
+    "            lamb = wu.deg2rad(lamb_deg)\n",
+    "            u, v = ell.ell2para(beta, lamb)\n",
+    "            alpha = wu.deg2rad(alpha_deg)\n",
+    "\n",
+    "            alpha_para_1, *_ = alpha_ell2para(ell, beta, lamb, alpha)\n",
+    "            alpha_ell_1, *_ = alpha_para2ell(ell, u, v, alpha_para_1)\n",
+    "            diff_1 = wu.deg2rad(abs(alpha_ell_1 - alpha))/3600\n",
+    "\n",
+    "            alpha_ell_2, *_ = alpha_para2ell(ell, u, v, alpha)\n",
+    "            alpha_para_2, *_ = alpha_ell2para(ell, beta, lamb, alpha_ell_2)\n",
+    "            diff_2 = wu.deg2rad(abs(alpha_para_2 - alpha))/3600\n",
+    "\n",
+    "            diffs.append((beta_deg, lamb_deg, alpha_deg, diff_1, diff_2))\n",
+    "diffs = np.array(diffs)"
+   ],
+   "id": "82fc6cbbe7d5abcb"
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
+   "source": [
+    "i_max_ell = np.argmax(diffs[:, 3])\n",
+    "max_ell = diffs[i_max_ell, 3]\n",
+    "point_max_ell = diffs[i_max_ell, :3]\n",
+    "\n",
+    "i_max_para = np.argmax(diffs[:, 4])\n",
+    "max_para = diffs[i_max_para, 4]\n",
+    "point_max_para = diffs[i_max_para, :4]\n",
+    "\n",
+    "print(f'Für elliptisches Alpha = {point_max_ell[2]}° und beta = {point_max_ell[0]}°, lamb = {point_max_ell[1]}°: diff = {max_ell}\"')\n",
+    "print(f'Für parametrisches Alpha = {point_max_para[2]}° und beta = {point_max_para[0]}°, lamb = {point_max_para[1]}°: diff = {max_ell}\"')\n",
+    "pass"
+   ],
+   "id": "97b5b8c9ca5377ab"
   },
   {
    "metadata": {
@@ -70,40 +108,6 @@
    "id": "98b9b220118deb3f",
    "outputs": [],
    "execution_count": 6
-  },
-  {
-   "metadata": {
-    "ExecuteTime": {
-     "end_time": "2026-01-20T15:33:50.497990Z",
-     "start_time": "2026-01-20T15:33:50.261115Z"
-    }
-   },
-   "cell_type": "code",
-   "source": [
-    "i_max_ell = np.argmax(diffs[:, 3])\n",
-    "max_ell = diffs[i_max_ell, 3]\n",
-    "point_max_ell = diffs[i_max_ell, :3]\n",
-    "\n",
-    "i_max_para = np.argmax(diffs[:, 4])\n",
-    "max_para = diffs[i_max_para, 4]\n",
-    "point_max_para = diffs[i_max_para, :4]\n",
-    "\n",
-    "print(f'Für elliptisches Alpha = {point_max_ell[2]}° und beta = {point_max_ell[0]}°, lamb = {point_max_ell[1]}°: diff = {max_ell}\"')\n",
-    "print(f'Für parametrisches Alpha = {point_max_para[2]}° und beta = {point_max_para[0]}°, lamb = {point_max_para[1]}°: diff = {max_ell}\"')\n",
-    "pass"
-   ],
-   "id": "3c74b65b0e85e3c2",
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Für elliptisches Alpha = 315.0° und beta = -90.0°, lamb = -90.0°: diff = 3.426945967752335e-05\"\n",
-      "Für parametrisches Alpha = 315.0° und beta = -90.0°, lamb = -90.0°: diff = 3.426945967752335e-05\"\n"
-     ]
-    }
-   ],
-   "execution_count": 7
   }
  ],
  "metadata": {