Commit 5b319203 authored by Gaetano Scandariato's avatar Gaetano Scandariato
Browse files

chore: keep PyORBIT's run.out of the K2-140 benchmark (the ln Z source) and ignore its pickles 

Co-Authored-By: default avatarClaude Sonnet 5.5 <noreply@anthropic.com>
parent cfc0ee3b

paper_benchmarks/.gitignore

+5 −0
Changes for paper_benchmarks/.gitignore: 5 added lines, 0 removed lines.
Original line number Diff line number Diff line
@@ -7,3 +7,8 @@ wasp18_phase_curve/runs/aies_de*/PyAETNAouts/*_samples.npz 
# Posterior chains and checkpoints of the TTV-scaling NUTS runs (up to 33 MB each); summaries are kept.

ttv_scaling/runs/nuts*/PyAETNAouts/*_samples.npz

ttv_scaling/runs/*/PyAETNAouts/*.pckl

# PyORBIT's outputs of the K2-140 runs (pickles, dynesty checkpoint; ~25 MB each) and the per-run copies of the inputs;

# run.out, time.out and the exported posterior (results/pyorbit_posterior.npz) are kept.

k2140_pyorbit/runs/*/k2140_pyorbit/

k2140_pyorbit/runs/*/k2140_*.dat

k2140_pyorbit/runs/*/run_pyorbit.py

paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run1/run.out

0 → 100644
+104 −0
Changes for paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run1/run.out: 104 added lines, 0 removed lines.
Original line number Diff line number Diff line 


Internal reformatting for Nested Sampling compatibility of priors



Opening:  LC_K2

Opening:  LC_LCOGT

Opening:  RV_CORALIE

Opening:  RV_FIES

Opening:  RV_HARPS



Adding common model of planet:  b

  LC model: transit_K2_b is using ld_K2 LD parameters

  LC model: transit_LCOGT_b is using ld_LCOGT LD parameters

*** planet b global parameters:

    orbital model: circular

    orbital parametrization:  Standard

    time of inferior conjunction replacing mean longitude as a free parameter:  True

    longitude of ascending node as a free parameter:  False



Using limb darkening coefficient parametrization:  Kipping

Using limb darkening coefficient parametrization:  Kipping



====================================================================================================

     Ids, spaces (s), boundaries (b) and priors (p) of the sampler parameters     

====================================================================================================



----- dataset:  LC_K2

jitter_0      id:  10  s:Logarithmic b:[    -19.9374,     -10.9658]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_K2                     ----- model: normalization_K2              

n_factor      id:  11  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  LC_LCOGT

jitter_0      id:  12  s:Logarithmic b:[    -19.9374,      -7.6439]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_LCOGT                  ----- model: normalization_LCOGT           

n_factor      id:  13  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  RV_CORALIE

jitter_0      id:  14  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  15  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1220.   50.]



----- dataset:  RV_FIES

jitter_0      id:  16  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  17  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1215.   50.]



----- dataset:  RV_HARPS

jitter_0      id:  18  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  19  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1240.   50.]



----- common model:  b

P             id:   0  s:Linear      b:[      6.5000,       6.7000]   p:Gaussian   [6.5693e+00 1.0000e-04]

K             id:   1  s:Linear      b:[      0.0000,    1000.0000]   p:Uniform   []

Tc            id:   2  s:Linear      b:[2457588.0000, 2457588.6000]   p:Gaussian   [2.45758828e+06 1.00000000e-03]

R_Rs          id:   3  s:Linear      b:[      0.0100,       0.5000]   p:Uniform   []

b             id:   4  s:Linear      b:[      0.0000,       1.2000]   p:Uniform   []

e             fixed (no id, space, bound)                             p:Fixed   [0. 0.]

omega         fixed (no id, space, bound)                             p:Fixed   [90.  0.]



----- common model:  star_parameters

density       id:   7  s:Linear      b:[      0.0071,       8.5107]   p:Uniform   []



----- common model:  ld_K2

ld_q1         id:   5  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   6  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  ld_LCOGT

ld_q1         id:   8  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   9  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  normalization_K2



----- common model:  normalization_LCOGT



Number of live points: 500



number of multiprocessing threads: 24



Reference Time Tref:  2457588.284



*************************************************************



Using threading pool for dynesty: True

Setting up the Dynamic Nested Sampling, number of live points =    500

                                        posterior/evidence split = 0.000

                                        inizial stopping criterion = 0.1000



Running Dynamic Nested Sampling





Getting the results



Summary

=======

niter: 27720

ncall: 1080538

eff(%):  2.612

logz: 14666.585 +/-  0.513

paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run10/run.out

0 → 100644
+104 −0
Changes for paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run10/run.out: 104 added lines, 0 removed lines.
Original line number Diff line number Diff line 


Internal reformatting for Nested Sampling compatibility of priors



Opening:  LC_K2

Opening:  LC_LCOGT

Opening:  RV_CORALIE

Opening:  RV_FIES

Opening:  RV_HARPS



Adding common model of planet:  b

  LC model: transit_K2_b is using ld_K2 LD parameters

  LC model: transit_LCOGT_b is using ld_LCOGT LD parameters

*** planet b global parameters:

    orbital model: circular

    orbital parametrization:  Standard

    time of inferior conjunction replacing mean longitude as a free parameter:  True

    longitude of ascending node as a free parameter:  False



Using limb darkening coefficient parametrization:  Kipping

Using limb darkening coefficient parametrization:  Kipping



====================================================================================================

     Ids, spaces (s), boundaries (b) and priors (p) of the sampler parameters     

====================================================================================================



----- dataset:  LC_K2

jitter_0      id:  10  s:Logarithmic b:[    -19.9374,     -10.9658]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_K2                     ----- model: normalization_K2              

n_factor      id:  11  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  LC_LCOGT

jitter_0      id:  12  s:Logarithmic b:[    -19.9374,      -7.6439]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_LCOGT                  ----- model: normalization_LCOGT           

n_factor      id:  13  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  RV_CORALIE

jitter_0      id:  14  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  15  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1220.   50.]



----- dataset:  RV_FIES

jitter_0      id:  16  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  17  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1215.   50.]



----- dataset:  RV_HARPS

jitter_0      id:  18  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  19  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1240.   50.]



----- common model:  b

P             id:   0  s:Linear      b:[      6.5000,       6.7000]   p:Gaussian   [6.5693e+00 1.0000e-04]

K             id:   1  s:Linear      b:[      0.0000,    1000.0000]   p:Uniform   []

Tc            id:   2  s:Linear      b:[2457588.0000, 2457588.6000]   p:Gaussian   [2.45758828e+06 1.00000000e-03]

R_Rs          id:   3  s:Linear      b:[      0.0100,       0.5000]   p:Uniform   []

b             id:   4  s:Linear      b:[      0.0000,       1.2000]   p:Uniform   []

e             fixed (no id, space, bound)                             p:Fixed   [0. 0.]

omega         fixed (no id, space, bound)                             p:Fixed   [90.  0.]



----- common model:  star_parameters

density       id:   7  s:Linear      b:[      0.0071,       8.5107]   p:Uniform   []



----- common model:  ld_K2

ld_q1         id:   5  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   6  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  ld_LCOGT

ld_q1         id:   8  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   9  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  normalization_K2



----- common model:  normalization_LCOGT



Number of live points: 500



number of multiprocessing threads: 24



Reference Time Tref:  2457588.284



*************************************************************



Using threading pool for dynesty: True

Setting up the Dynamic Nested Sampling, number of live points =    500

                                        posterior/evidence split = 0.000

                                        inizial stopping criterion = 0.1000



Running Dynamic Nested Sampling





Getting the results



Summary

=======

niter: 27959

ncall: 1089935

eff(%):  2.611

logz: 14665.979 +/-  0.539

paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run2/run.out

0 → 100644
+104 −0
Changes for paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run2/run.out: 104 added lines, 0 removed lines.
Original line number Diff line number Diff line 


Internal reformatting for Nested Sampling compatibility of priors



Opening:  LC_K2

Opening:  LC_LCOGT

Opening:  RV_CORALIE

Opening:  RV_FIES

Opening:  RV_HARPS



Adding common model of planet:  b

  LC model: transit_K2_b is using ld_K2 LD parameters

  LC model: transit_LCOGT_b is using ld_LCOGT LD parameters

*** planet b global parameters:

    orbital model: circular

    orbital parametrization:  Standard

    time of inferior conjunction replacing mean longitude as a free parameter:  True

    longitude of ascending node as a free parameter:  False



Using limb darkening coefficient parametrization:  Kipping

Using limb darkening coefficient parametrization:  Kipping



====================================================================================================

     Ids, spaces (s), boundaries (b) and priors (p) of the sampler parameters     

====================================================================================================



----- dataset:  LC_K2

jitter_0      id:  10  s:Logarithmic b:[    -19.9374,     -10.9658]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_K2                     ----- model: normalization_K2              

n_factor      id:  11  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  LC_LCOGT

jitter_0      id:  12  s:Logarithmic b:[    -19.9374,      -7.6439]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_LCOGT                  ----- model: normalization_LCOGT           

n_factor      id:  13  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  RV_CORALIE

jitter_0      id:  14  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  15  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1220.   50.]



----- dataset:  RV_FIES

jitter_0      id:  16  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  17  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1215.   50.]



----- dataset:  RV_HARPS

jitter_0      id:  18  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  19  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1240.   50.]



----- common model:  b

P             id:   0  s:Linear      b:[      6.5000,       6.7000]   p:Gaussian   [6.5693e+00 1.0000e-04]

K             id:   1  s:Linear      b:[      0.0000,    1000.0000]   p:Uniform   []

Tc            id:   2  s:Linear      b:[2457588.0000, 2457588.6000]   p:Gaussian   [2.45758828e+06 1.00000000e-03]

R_Rs          id:   3  s:Linear      b:[      0.0100,       0.5000]   p:Uniform   []

b             id:   4  s:Linear      b:[      0.0000,       1.2000]   p:Uniform   []

e             fixed (no id, space, bound)                             p:Fixed   [0. 0.]

omega         fixed (no id, space, bound)                             p:Fixed   [90.  0.]



----- common model:  star_parameters

density       id:   7  s:Linear      b:[      0.0071,       8.5107]   p:Uniform   []



----- common model:  ld_K2

ld_q1         id:   5  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   6  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  ld_LCOGT

ld_q1         id:   8  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   9  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  normalization_K2



----- common model:  normalization_LCOGT



Number of live points: 500



number of multiprocessing threads: 24



Reference Time Tref:  2457588.284



*************************************************************



Using threading pool for dynesty: True

Setting up the Dynamic Nested Sampling, number of live points =    500

                                        posterior/evidence split = 0.000

                                        inizial stopping criterion = 0.1000



Running Dynamic Nested Sampling





Getting the results



Summary

=======

niter: 27941

ncall: 1086489

eff(%):  2.618

logz: 14666.346 +/-  0.527

paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run3/run.out

0 → 100644
+104 −0
Changes for paper_benchmarks/k2140_pyorbit/runs/laptop_cpu24_run3/run.out: 104 added lines, 0 removed lines.
Original line number Diff line number Diff line 


Internal reformatting for Nested Sampling compatibility of priors



Opening:  LC_K2

Opening:  LC_LCOGT

Opening:  RV_CORALIE

Opening:  RV_FIES

Opening:  RV_HARPS



Adding common model of planet:  b

  LC model: transit_K2_b is using ld_K2 LD parameters

  LC model: transit_LCOGT_b is using ld_LCOGT LD parameters

*** planet b global parameters:

    orbital model: circular

    orbital parametrization:  Standard

    time of inferior conjunction replacing mean longitude as a free parameter:  True

    longitude of ascending node as a free parameter:  False



Using limb darkening coefficient parametrization:  Kipping

Using limb darkening coefficient parametrization:  Kipping



====================================================================================================

     Ids, spaces (s), boundaries (b) and priors (p) of the sampler parameters     

====================================================================================================



----- dataset:  LC_K2

jitter_0      id:  10  s:Logarithmic b:[    -19.9374,     -10.9658]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_K2                     ----- model: normalization_K2              

n_factor      id:  11  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  LC_LCOGT

jitter_0      id:  12  s:Logarithmic b:[    -19.9374,      -7.6439]   p:Uniform   []

offset_0      fixed (no id, space, bound)                             p:Fixed   [0. 0.]



----- dataset: LC_LCOGT                  ----- model: normalization_LCOGT           

n_factor      id:  13  s:Linear      b:[      0.5000,       1.5000]   p:Gaussian   [1.  0.1]



----- dataset:  RV_CORALIE

jitter_0      id:  14  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  15  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1220.   50.]



----- dataset:  RV_FIES

jitter_0      id:  16  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  17  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1215.   50.]



----- dataset:  RV_HARPS

jitter_0      id:  18  s:Logarithmic b:[     -3.3219,       6.6439]   p:Uniform   []

offset_0      id:  19  s:Linear      b:[   1000.0000,    1500.0000]   p:Gaussian   [1240.   50.]



----- common model:  b

P             id:   0  s:Linear      b:[      6.5000,       6.7000]   p:Gaussian   [6.5693e+00 1.0000e-04]

K             id:   1  s:Linear      b:[      0.0000,    1000.0000]   p:Uniform   []

Tc            id:   2  s:Linear      b:[2457588.0000, 2457588.6000]   p:Gaussian   [2.45758828e+06 1.00000000e-03]

R_Rs          id:   3  s:Linear      b:[      0.0100,       0.5000]   p:Uniform   []

b             id:   4  s:Linear      b:[      0.0000,       1.2000]   p:Uniform   []

e             fixed (no id, space, bound)                             p:Fixed   [0. 0.]

omega         fixed (no id, space, bound)                             p:Fixed   [90.  0.]



----- common model:  star_parameters

density       id:   7  s:Linear      b:[      0.0071,       8.5107]   p:Uniform   []



----- common model:  ld_K2

ld_q1         id:   5  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   6  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  ld_LCOGT

ld_q1         id:   8  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_q2         id:   9  s:Linear      b:[      0.0000,       1.0000]   p:Uniform   []

ld_c1         derived (no id, space, bound)                           p:Uniform   [0.0, 1.0]

ld_c2         derived (no id, space, bound)                           p:Uniform   [-1.0, 1.0]



----- common model:  normalization_K2



----- common model:  normalization_LCOGT



Number of live points: 500



number of multiprocessing threads: 24



Reference Time Tref:  2457588.284



*************************************************************



Using threading pool for dynesty: True

Setting up the Dynamic Nested Sampling, number of live points =    500

                                        posterior/evidence split = 0.000

                                        inizial stopping criterion = 0.1000



Running Dynamic Nested Sampling





Getting the results



Summary

=======

niter: 27499

ncall: 1069693

eff(%):  2.617

logz: 14666.404 +/-  0.549