Added output for Ex 04, 06, 09. Renamed directory: figures --> output…

… (so including text output)
rpoleski · Nov 30, 2023 · fb91742 · fb91742
1 parent 2e6280d
commit fb91742
Show file tree

Hide file tree

Showing 26 changed files with 69 additions and 0 deletions.
diff --git a/...mple_figures/example_01/Ex01_Figure_1.png → ...ample_output/example_01/Ex01_Figure_1.png b/...mple_figures/example_01/Ex01_Figure_1.png → ...ample_output/example_01/Ex01_Figure_1.png
diff --git a/...mple_figures/example_01/Ex01_Figure_2.png → ...ample_output/example_01/Ex01_Figure_2.png b/...mple_figures/example_01/Ex01_Figure_2.png → ...ample_output/example_01/Ex01_Figure_2.png
diff --git a/...mple_figures/example_01/Ex01_Figure_3.png → ...ample_output/example_01/Ex01_Figure_3.png b/...mple_figures/example_01/Ex01_Figure_3.png → ...ample_output/example_01/Ex01_Figure_3.png
diff --git a/...mple_figures/example_01/Ex01_Figure_4.png → ...ample_output/example_01/Ex01_Figure_4.png b/...mple_figures/example_01/Ex01_Figure_4.png → ...ample_output/example_01/Ex01_Figure_4.png
diff --git a/...mple_figures/example_01/Ex01_Figure_5.png → ...ample_output/example_01/Ex01_Figure_5.png b/...mple_figures/example_01/Ex01_Figure_5.png → ...ample_output/example_01/Ex01_Figure_5.png
diff --git a/...mple_figures/example_02/Ex02_Figure_1.png → ...ample_output/example_02/Ex02_Figure_1.png b/...mple_figures/example_02/Ex02_Figure_1.png → ...ample_output/example_02/Ex02_Figure_1.png
diff --git a/...mple_figures/example_02/Ex02_Figure_2.png → ...ample_output/example_02/Ex02_Figure_2.png b/...mple_figures/example_02/Ex02_Figure_2.png → ...ample_output/example_02/Ex02_Figure_2.png
diff --git a/...mple_figures/example_03/Ex03_Figure_1.png → ...ample_output/example_03/Ex03_Figure_1.png b/...mple_figures/example_03/Ex03_Figure_1.png → ...ample_output/example_03/Ex03_Figure_1.png
diff --git a/...mple_figures/example_03/Ex03_Figure_2.png → ...ample_output/example_03/Ex03_Figure_2.png b/...mple_figures/example_03/Ex03_Figure_2.png → ...ample_output/example_03/Ex03_Figure_2.png
diff --git a/data/expected_example_output/example_04/rE.tab b/data/expected_example_output/example_04/rE.tab
@@ -0,0 +1,10 @@
+# Physical Einstein ring radius in AU (Lens plane).
+# columns = distance in kpc
+# rows = mass in solMass
+        1.00  2.00  3.00  4.00  5.00  6.00  7.00
+  10.0  8.44 11.05 12.36 12.76 12.36 11.05  8.44
+   1.0  2.67  3.50  3.91  4.04  3.91  3.50  2.67
+   0.3  1.46  1.91  2.14  2.21  2.14  1.91  1.46
+   0.1  0.84  1.11  1.24  1.28  1.24  1.11  0.84
+  0.01  0.27  0.35  0.39  0.40  0.39  0.35  0.27
+ 0.001  0.08  0.11  0.12  0.13  0.12  0.11  0.08
diff --git a/data/expected_example_output/example_04/rEtilde.tab b/data/expected_example_output/example_04/rEtilde.tab
@@ -0,0 +1,10 @@
+# Projection of Einstein ring radius in AU (Observer plane).
+# columns = distance in kpc
+# rows = mass in solMass
+        1.00  2.00  3.00  4.00  5.00  6.00  7.00
+  10.0  9.65 14.74 19.77 25.52 32.95 44.21 67.53
+   1.0  3.05  4.66  6.25  8.07 10.42 13.98 21.36
+   0.3  1.67  2.55  3.42  4.42  5.71  7.66 11.70
+   0.1  0.96  1.47  1.98  2.55  3.30  4.42  6.75
+  0.01  0.31  0.47  0.63  0.81  1.04  1.40  2.14
+ 0.001  0.10  0.15  0.20  0.26  0.33  0.44  0.68
diff --git a/data/expected_example_output/example_04/tE.tab b/data/expected_example_output/example_04/tE.tab
@@ -0,0 +1,10 @@
+# Einstein timescale in days for mu_rel =  4.0 mas/yr.
+# columns = distance in kpc
+# rows = mass in solMass
+        1.00  2.00  3.00  4.00  5.00  6.00  7.00
+  10.0 770.8 504.6 376.1 291.3 225.7 168.2 110.1
+   1.0 243.8 159.6 118.9  92.1  71.4  53.2  34.8
+   0.3 133.5  87.4  65.1  50.5  39.1  29.1  19.1
+   0.1  77.1  50.5  37.6  29.1  22.6  16.8  11.0
+  0.01  24.4  16.0  11.9   9.2   7.1   5.3   3.5
+ 0.001   7.7   5.0   3.8   2.9   2.3   1.7   1.1
diff --git a/data/expected_example_output/example_04/thetaE.tab b/data/expected_example_output/example_04/thetaE.tab
@@ -0,0 +1,10 @@
+# Angular Einstein ring radius in milliarcseconds.
+# columns = distance in kpc
+# rows = mass in solMass
+        1.00  2.00  3.00  4.00  5.00  6.00  7.00
+  10.0  8.44  5.53  4.12  3.19  2.47  1.84  1.21
+   1.0  2.67  1.75  1.30  1.01  0.78  0.58  0.38
+   0.3  1.46  0.96  0.71  0.55  0.43  0.32  0.21
+   0.1  0.84  0.55  0.41  0.32  0.25  0.18  0.12
+  0.01  0.27  0.17  0.13  0.10  0.08  0.06  0.04
+ 0.001  0.08  0.06  0.04  0.03  0.02  0.02  0.01
diff --git a/...mple_figures/example_05/Ex05_Figure_1.png → ...ample_output/example_05/Ex05_Figure_1.png b/...mple_figures/example_05/Ex05_Figure_1.png → ...ample_output/example_05/Ex05_Figure_1.png
diff --git a/...figures/example_05/Ex05_Figure_1_zoom.png → ..._output/example_05/Ex05_Figure_1_zoom.png b/...figures/example_05/Ex05_Figure_1_zoom.png → ..._output/example_05/Ex05_Figure_1_zoom.png
diff --git a/...mple_figures/example_05/Ex05_Figure_2.png → ...ample_output/example_05/Ex05_Figure_2.png b/...mple_figures/example_05/Ex05_Figure_2.png → ...ample_output/example_05/Ex05_Figure_2.png
diff --git a/...mple_figures/example_05/Ex05_Figure_3.png → ...ample_output/example_05/Ex05_Figure_3.png b/...mple_figures/example_05/Ex05_Figure_3.png → ...ample_output/example_05/Ex05_Figure_3.png
diff --git a/...figures/example_05/Ex05_Figure_3_zoom.png → ..._output/example_05/Ex05_Figure_3_zoom.png b/...figures/example_05/Ex05_Figure_3_zoom.png → ..._output/example_05/Ex05_Figure_3_zoom.png
diff --git a/...mple_figures/example_05/Ex05_Figure_4.png → ...ample_output/example_05/Ex05_Figure_4.png b/...mple_figures/example_05/Ex05_Figure_4.png → ...ample_output/example_05/Ex05_Figure_4.png
diff --git a/...mple_figures/example_05/Ex05_Figure_5.png → ...ample_output/example_05/Ex05_Figure_5.png b/...mple_figures/example_05/Ex05_Figure_5.png → ...ample_output/example_05/Ex05_Figure_5.png
diff --git a/...mple_figures/example_05/Ex05_Figure_6.png → ...ample_output/example_05/Ex05_Figure_6.png b/...mple_figures/example_05/Ex05_Figure_6.png → ...ample_output/example_05/Ex05_Figure_6.png
diff --git a/data/expected_example_output/example_06/Ex06_Figure_1.png b/data/expected_example_output/example_06/Ex06_Figure_1.png
diff --git a/data/expected_example_output/example_06/Ex06_output.txt b/data/expected_example_output/example_06/Ex06_output.txt
@@ -0,0 +1,10 @@
+Fitted parameters:
+2453630.39088 0.22670 0.23167
+0.49048 0.02013 0.01869
+93.45656 2.19216 2.15444
+0.27419 0.02800 0.03255
+0.10357 0.00732 0.00715
+
+Smallest chi2 model:
+2453630.376707949 0.48814776622836326 93.89323919007633 0.2722480223300758 0.10093639484208665
+946.7504141947664
diff --git a/data/expected_example_output/example_09/Ex09_Figure_1.png b/data/expected_example_output/example_09/Ex09_Figure_1.png
diff --git a/data/expected_example_output/example_09/Ex09_Figure_2.png b/data/expected_example_output/example_09/Ex09_Figure_2.png
diff --git a/data/expected_example_output/example_09/Ex09_output.txt b/data/expected_example_output/example_09/Ex09_output.txt
@@ -0,0 +1,19 @@
+Initial Trial
+    t_0 (HJD)       u_0    t_E (d)
+   5380.00000  0.100000    18.0000
+Chi2 = 32241.84124710082
+
+Best Fit: t_0 =   5379.57091, u_0 = 0.5230, t_E =   17.940
+Chi2 =       427.20
+
+scipy.optimize.minimize result:
+     fun: 427.2038219878052
+     jac: array([-0.03365149, -3.52971227, -0.1673083 ])
+ message: 'Optimization terminated successfully.'
+    nfev: 11
+    nhev: 0
+     nit: 9
+    njev: 28
+  status: 0
+ success: True
+       x: array([5.37957091e+03, 5.22985608e-01, 1.79398771e+01])