refactored newton method

doc/api.rst

@@ -42,9 +42,7 @@ $\partial$SGP4 API
    util.from_datetime_to_fractional_day
    util.from_datetime_to_mjd
    util.from_datetime_to_jd
-   util.from_cartesian_to_tle_elements
    util.from_cartesian_to_keplerian
-   util.from_cartesian_to_keplerian_torch
    sgp4
    sgp4_batched
    sgp4init.sgp4init

doc/conf.py

@@ -85,6 +85,13 @@
 
 html_logo = "_static/logo_dsgp4.png"
 
+
+linkcheck_ignore = [
+     r'https://www\.esa\.int/gsp/ACT/team/giacomo_acciarini/',
+     r'https://www.esa.int/gsp/ACT/team/dario_izzo/',
+ ]
+
+
 html_theme_options = {
     "repository_url": "https://github.com/esa/dSGP4/",
     "repository_branch": "master",

doc/notebooks/gradient_based_optimization.ipynb

@@ -48,7 +48,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 1,
    "id": "95627789",
    "metadata": {},
    "outputs": [],
@@ -59,7 +59,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 15,
+   "execution_count": 2,
    "id": "91d2d6c5",
    "metadata": {},
    "outputs": [],
@@ -72,16 +72,15 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 16,
+   "execution_count": 8,
    "id": "c376c651",
    "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "F(y): 1.0751821836135118e-13\n",
-      "Solution found, at iter: 6\n",
+      "Solution F(y) = 2.4867357777951112e-14, converged in 5 iterations\n",
       "TLE(\n",
       "0 COSMOS 2251 DEB\n",
       "1 34454U 93036SX  22068.91971155  .00000319  00000-0  11812-3 0  9996\n",
@@ -95,23 +94,22 @@
    ],
    "source": [
     "#when I propagate to zero, I expect the returned TLE to be identical to my_tle:\n",
-    "found_tle, y=dsgp4.newton_method(tle_0=my_tle,time_mjd=my_tle.date_mjd,new_tol=1e-12,max_iter=100)\n",
+    "found_tle, y=dsgp4.newton_method(tle0=my_tle,time_mjd=my_tle.date_mjd,new_tol=1e-12,max_iter=10,verbose=True)\n",
     "\n",
     "print(my_tle,found_tle)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 7,
    "id": "e0f85268",
    "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "F(y): 4.547475677268469e-13\n",
-      "Solution found, at iter: 7\n",
+      "Solution F(y) = 4.547517012268725e-13, converged in 5 iterations\n",
       "TLE(\n",
       "0 COSMOS 2251 DEB\n",
       "1 34454U 93036SX  22068.91971155  .00000319  00000-0  11812-3 0  9996\n",
@@ -125,18 +123,19 @@
    ],
    "source": [
     "#I now propagate until 1000 minutes after\n",
-    "found_tle, y=dsgp4.newton_method(tle_0=my_tle,\n",
+    "found_tle, y=dsgp4.newton_method(tle0=my_tle,\n",
     "                                 time_mjd=my_tle.date_mjd+800.,\n",
     "                                 new_tol=1e-12,\n",
-    "                                 max_iter=20)\n",
+    "                                 max_iter=10,\n",
+    "                                 verbose=True)\n",
     "#Newton still converges, and the TLE is of course now different:\n",
     "print(my_tle, found_tle)"
    ]
   }
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "Python 3 (ipykernel)",
+   "display_name": "dsgp4",
    "language": "python",
    "name": "python3"
   },

dsgp4/__init__.py

@@ -1,4 +1,4 @@
-__version__ = '1.1.4'
+__version__ = '1.1.5'
 
 import torch
 torch.set_default_dtype(torch.float64)