ex08 - expected output

rpoleski · Jun 3, 2024 · a136ed3 · a136ed3
1 parent 6470f56
commit a136ed3
Show file tree

Hide file tree

Showing 15 changed files with 580 additions and 0 deletions.
diff --git a/data/expected_example_output/example_08/Ex08_Figure_01.png b/data/expected_example_output/example_08/Ex08_Figure_01.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_02.png b/data/expected_example_output/example_08/Ex08_Figure_02.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_03.png b/data/expected_example_output/example_08/Ex08_Figure_03.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_04.png b/data/expected_example_output/example_08/Ex08_Figure_04.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_05.png b/data/expected_example_output/example_08/Ex08_Figure_05.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_06.png b/data/expected_example_output/example_08/Ex08_Figure_06.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_07.png b/data/expected_example_output/example_08/Ex08_Figure_07.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_08.png b/data/expected_example_output/example_08/Ex08_Figure_08.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_09.png b/data/expected_example_output/example_08/Ex08_Figure_09.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_10.png b/data/expected_example_output/example_08/Ex08_Figure_10.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_11.png b/data/expected_example_output/example_08/Ex08_Figure_11.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_12.png b/data/expected_example_output/example_08/Ex08_Figure_12.png
diff --git a/data/expected_example_output/example_08/Ex08_Figure_13.png b/data/expected_example_output/example_08/Ex08_Figure_13.png
diff --git a/data/expected_example_output/example_08/Ex08_output.txt b/data/expected_example_output/example_08/Ex08_output.txt
@@ -0,0 +1,248 @@
+Initial Guess
+    t_0 (HJD)       u_0    t_E (d) 
+2460962.36458  0.411823    22.8092 
+default magnification method: point_source
+(array([7.80252182]), 0.7542833943450846)
+[2.58000914]
+Best-fit Point Lens
+    t_0 (HJD)       u_0    t_E (d) 
+2460963.60793  0.445415    24.5259 
+default magnification method: point_source
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00100   1.56450   0.00010000   -28.77442 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+        chi2      s       q   alpha     rho
+         inf 0.6104 0.00001  154.39 deg 1.31513
+   109585.90 0.6104 0.00002  172.79 deg 1.31524
+   126863.76 0.6104 0.00003  155.62 deg 0.00001
+   126755.39 0.6104 0.00006  155.68 deg 0.00000
+   126564.83 0.6104 0.00010  155.77 deg 0.00000
+   126233.20 0.6104 0.00018  155.94 deg 0.00000
+   125667.51 0.6104 0.00032  156.21 deg 0.00001
+   124734.87 0.6104 0.00056  156.66 deg 0.00000
+   123266.77 0.6104 0.00100  157.33 deg 0.00001
+   121049.72 0.6104 0.00178  158.36 deg 0.00001
+   117879.53 0.6104 0.00316  159.99 deg 0.00001
+   113717.79 0.6104 0.00562  162.38 deg 0.00000
+   126757.84 0.6246 0.00001  153.66 deg 0.00001
+   126605.54 0.6246 0.00002  153.85 deg 0.00000
+   126374.03 0.6246 0.00003  154.09 deg 0.00000
+   126035.52 0.6246 0.00006  154.41 deg 0.00000
+   125561.54 0.6246 0.00010  154.84 deg 0.00001
+   123288.75 0.6246 0.00018  154.14 deg 0.00001
+   120876.61 0.6246 0.00032  154.65 deg 0.00001
+   117323.91 0.6246 0.00056  155.38 deg 0.00000
+   112342.06 0.6246 0.00100  156.42 deg 0.00000
+   108694.43 0.6246 0.00178  145.47 deg 0.00000
+   119015.41 0.6246 0.00316  162.45 deg 0.00001
+    87927.11 0.6246 0.00562  162.52 deg 0.00000
+   119486.27 0.6392 0.00001  150.74 deg 0.00074
+   122072.20 0.6392 0.00002  151.90 deg 0.00045
+   118976.86 0.6392 0.00003  152.12 deg 0.00051
+   114304.97 0.6392 0.00006  152.42 deg 0.00065
+   107689.38 0.6392 0.00010  152.82 deg 0.00095
+    98991.65 0.6392 0.00018  153.35 deg 0.00113
+    87935.87 0.6392 0.00032  154.04 deg 0.00121
+    74210.57 0.6392 0.00056  154.97 deg 0.00116
+    57922.12 0.6392 0.00100  156.20 deg 0.00131
+    40879.29 0.6392 0.00178  157.86 deg 0.00138
+    52646.47 0.6392 0.00316  143.13 deg 0.00730
+    37713.51 0.6392 0.00562  162.72 deg 0.00396
+   109587.22 0.6541 0.00001  167.36 deg 1.31518
+   126872.88 0.6541 0.00002  154.68 deg 0.00000
+   126770.48 0.6541 0.00003  154.71 deg 0.00000
+   126589.58 0.6541 0.00006  154.77 deg 0.00000
+   126269.39 0.6541 0.00010  154.88 deg 0.00000
+   125710.03 0.6541 0.00018  155.08 deg 0.00000
+   124758.15 0.6541 0.00032  155.42 deg 0.00000
+   123217.56 0.6541 0.00056  155.99 deg 0.00000
+   120919.43 0.6541 0.00100  156.90 deg 0.00000
+   117903.50 0.6541 0.00178  158.28 deg 0.00000
+   108443.79 0.6541 0.00316  143.43 deg 0.00000
+   111267.68 0.6541 0.00562  165.28 deg 0.06254
+   126797.90 1.4941 0.00001  -28.43 deg 0.00000
+   126637.63 1.4941 0.00002  -28.43 deg 0.00000
+   126353.10 1.4941 0.00003  -28.43 deg 0.00000
+   125848.68 1.4941 0.00006  -28.42 deg 0.00000
+   124956.55 1.4941 0.00010  -28.41 deg 0.00000
+   123385.65 1.4941 0.00018  -28.39 deg 0.00000
+   120642.08 1.4941 0.00032  -28.36 deg 0.00000
+   115924.35 1.4941 0.00056  -28.32 deg 0.00000
+   108051.26 1.4941 0.00100  -28.27 deg 0.00000
+    95658.48 1.4941 0.00178  -28.22 deg 0.00000
+    78488.50 1.4941 0.00316  -28.22 deg 0.00000
+    63051.31 1.4941 0.00562  -28.27 deg 0.00000
+   126328.40 1.5289 0.00001  -28.56 deg 0.00000
+   125804.89 1.5289 0.00002  -28.55 deg 0.00000
+   124877.69 1.5289 0.00003  -28.55 deg 0.00000
+   123236.39 1.5289 0.00006  -28.54 deg 0.00000
+   120329.32 1.5289 0.00010  -28.53 deg 0.00000
+   115182.20 1.5289 0.00018  -28.51 deg 0.00000
+   106176.99 1.5289 0.00032  -28.49 deg 0.00000
+    91217.24 1.5289 0.00056  -28.47 deg 0.00000
+    70956.63 1.5289 0.00100  -28.43 deg 0.00000
+    68154.46 1.5289 0.00178  -28.28 deg 0.00766
+   109581.42 1.5289 0.00316 -226.53 deg 1.31384
+   109587.84 1.5289 0.00562   -6.58 deg 1.31642
+   106275.48 1.5645 0.00001  -28.64 deg 0.00116
+    98230.89 1.5645 0.00002  -28.61 deg 0.00146
+    87557.39 1.5645 0.00003  -28.56 deg 0.00191
+    74292.16 1.5645 0.00006  -28.49 deg 0.00249
+    58455.76 1.5645 0.00010  -28.40 deg 0.00338
+   109585.54 1.5645 0.00018 -279.69 deg 1.31480
+   109571.98 1.5645 0.00032 -1191.35 deg 1.31998
+   109585.67 1.5645 0.00056 -311.33 deg 1.31621
+    75018.23 1.5645 0.00100  -27.61 deg 0.02555
+   141081.50 1.5645 0.00178  -28.44 deg 0.00967
+   220471.70 1.5645 0.00316  -28.60 deg 0.01309
+   317860.09 1.5645 0.00562  -28.84 deg 0.01750
+   125450.61 1.6009 0.00001  -28.13 deg 0.00001
+   124227.52 1.6009 0.00002  -28.13 deg 0.00000
+   122023.03 1.6009 0.00003  -28.13 deg 0.00000
+   118037.01 1.6009 0.00006  -28.14 deg 0.00000
+   110906.28 1.6009 0.00010  -28.15 deg 0.00000
+    98759.83 1.6009 0.00018  -28.17 deg 0.00000
+    81170.52 1.6009 0.00032  -28.22 deg 0.00000
+    66662.37 1.6009 0.00056  -28.36 deg 0.00575
+    73045.70 1.6009 0.00100  -28.56 deg 0.00933
+   109587.04 1.6009 0.00178  -36.87 deg 1.31468
+   101566.33 1.6009 0.00316  -31.00 deg 0.05921
+   221215.63 1.6009 0.00562  -27.81 deg 0.00939
+Best Models
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00396   0.63918   0.00562341   162.72441 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00138   0.63918   0.00177828   157.86424 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00730   0.63918   0.00316228   143.13384 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00131   0.63918   0.00100000   156.20416 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00338   1.56450   0.00010000   -28.40424 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00396   0.63918   0.00562341   162.72441 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+       message: Maximum number of function evaluations has been exceeded.
+       success: False
+        status: 1
+           fun: 32187.282231863494
+             x: [ 2.461e+06  3.763e-01  2.696e+01  2.726e-03  1.618e+02
+                  6.672e-01  2.990e-03]
+           nit: 895
+          nfev: 1400
+ final_simplex: (array([[ 2.461e+06,  3.763e-01, ...,  6.672e-01,
+                         2.990e-03],
+                       [ 2.461e+06,  3.763e-01, ...,  6.672e-01,
+                         2.993e-03],
+                       ...,
+                       [ 2.461e+06,  3.763e-01, ...,  6.672e-01,
+                         2.991e-03],
+                       [ 2.461e+06,  3.762e-01, ...,  6.672e-01,
+                         2.986e-03]]), array([ 3.219e+04,  3.219e+04,  3.219e+04,  3.219e+04,
+                        3.219e+04,  3.219e+04,  3.219e+04,  3.219e+04]))
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00138   0.63918   0.00177828   157.86424 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+       message: Optimization terminated successfully.
+       success: True
+        status: 0
+           fun: 35104.622391109806
+             x: [ 2.461e+06  4.448e-01  2.560e+01  1.137e-03  1.572e+02
+                  6.483e-01  1.924e-03]
+           nit: 553
+          nfev: 970
+ final_simplex: (array([[ 2.461e+06,  4.448e-01, ...,  6.483e-01,
+                         1.924e-03],
+                       [ 2.461e+06,  4.448e-01, ...,  6.483e-01,
+                         1.924e-03],
+                       ...,
+                       [ 2.461e+06,  4.448e-01, ...,  6.483e-01,
+                         1.924e-03],
+                       [ 2.461e+06,  4.448e-01, ...,  6.483e-01,
+                         1.924e-03]]), array([ 3.510e+04,  3.510e+04,  3.510e+04,  3.510e+04,
+                        3.510e+04,  3.510e+04,  3.510e+04,  3.510e+04]))
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00730   0.63918   0.00316228   143.13384 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+       message: Optimization terminated successfully.
+       success: True
+        status: 0
+           fun: 29065.412476542464
+             x: [ 2.461e+06  3.476e-01  2.958e+01  3.689e-03  1.438e+02
+                  6.879e-01  3.328e-03]
+           nit: 420
+          nfev: 781
+ final_simplex: (array([[ 2.461e+06,  3.476e-01, ...,  6.879e-01,
+                         3.328e-03],
+                       [ 2.461e+06,  3.476e-01, ...,  6.879e-01,
+                         3.328e-03],
+                       ...,
+                       [ 2.461e+06,  3.476e-01, ...,  6.879e-01,
+                         3.328e-03],
+                       [ 2.461e+06,  3.476e-01, ...,  6.879e-01,
+                         3.328e-03]]), array([ 2.907e+04,  2.907e+04,  2.907e+04,  2.907e+04,
+                        2.907e+04,  2.907e+04,  2.907e+04,  2.907e+04]))
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00131   0.63918   0.00100000   156.20416 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+       message: Optimization terminated successfully.
+       success: True
+        status: 0
+           fun: 43231.18862833552
+             x: [ 2.461e+06  4.609e-01  2.509e+01  8.803e-04  1.556e+02
+                  6.410e-01  1.320e-03]
+           nit: 627
+          nfev: 1111
+ final_simplex: (array([[ 2.461e+06,  4.609e-01, ...,  6.410e-01,
+                         1.320e-03],
+                       [ 2.461e+06,  4.609e-01, ...,  6.410e-01,
+                         1.320e-03],
+                       ...,
+                       [ 2.461e+06,  4.609e-01, ...,  6.410e-01,
+                         1.320e-03],
+                       [ 2.461e+06,  4.609e-01, ...,  6.410e-01,
+                         1.320e-03]]), array([ 4.323e+04,  4.323e+04,  4.323e+04,  4.323e+04,
+                        4.323e+04,  4.323e+04,  4.323e+04,  4.323e+04]))
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.60793  0.445415    24.5259 0.00338   1.56450   0.00010000   -28.40424 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+       message: Maximum number of function evaluations has been exceeded.
+       success: False
+        status: 1
+           fun: 40312.11049729267
+             x: [ 2.461e+06  3.918e-01  2.999e+01  3.464e-03 -3.037e+01
+                  1.452e+00  1.610e-04]
+           nit: 864
+          nfev: 1400
+ final_simplex: (array([[ 2.461e+06,  3.918e-01, ...,  1.452e+00,
+                         1.610e-04],
+                       [ 2.461e+06,  3.918e-01, ...,  1.452e+00,
+                         1.610e-04],
+                       ...,
+                       [ 2.461e+06,  3.918e-01, ...,  1.452e+00,
+                         1.610e-04],
+                       [ 2.461e+06,  3.918e-01, ...,  1.452e+00,
+                         1.610e-04]]), array([ 4.031e+04,  4.031e+04,  4.031e+04,  4.031e+04,
+                        4.031e+04,  4.031e+04,  4.031e+04,  4.031e+04]))
+    t_0 (HJD)       u_0    t_E (d)     rho         s            q alpha (deg) 
+2460963.64558  0.347621    29.5829 0.00369   0.68785   0.00332758   143.76107 
+default magnification method: point_source
+other magnification methods: [0.0, 'point_source_point_lens', 2460982.0, 'VBBL', 2460985.0, 'point_source_point_lens', 2470000.0]
+chi2: 29065.41247654246
+Total Runtime: 0:32:44.949161