Commit 0225d787 authored by Gaetano Scandariato's avatar Gaetano Scandariato
Browse files

fix: cos i of Pi Men c from the transit impact parameter; add density-prior, dynesty and NUTS tests 

The comparison with allesfitter converted b to cos i as b/(a/R*), missing jaxoplanet's (1 + e sin w)/(1 - e^2)
for the eccentric model: the 3-sigma cos i_c discrepancy disappears (0.01 sigma). New: TOI-216 reweighted with
allesfitter's density prior, dynesty on the known-evidence Gaussian, NUTS runs, round-3 rerun queues.

Co-Authored-By: default avatarClaude Opus 5.5 <noreply@anthropic.com>
parent 1c53d79d

paper_benchmarks/k2140_juliet/runs/nuts_cpu24/k2140.yaml

0 → 100644
+140 −0
Changes for paper_benchmarks/k2140_juliet/runs/nuts_cpu24/k2140.yaml: 140 added lines, 0 removed lines.
Original line number Diff line number Diff line 
data:

  file_path: k2140c.dat

  timeZeroPoint: 2457588.284

  working_dir: ./

instruments:

- K2

- LCOGT

- CORALIE

- FIES

- HARPS

data_type:

- LC

- LC

- RV

- RV

- RV

sampling:

  chain_energy_threshold: 5.0

  checkpoint_interval: 500

  internal_warmup_if_nuts: 500

  max_tree_depth: 10

  num_chains_if_nuts: 4

  num_samples_if_ensemble: 2000

  num_samples_if_nuts: 2000

  num_warmup: 2000

  production_sampler: auto

  random_seed: 1

  reset: true

  target_accept_prob: 0.9

  use_float64: false

  walkers_per_free_parameter: 3

  warmup_sampler: ess

detrending:

  normalization:

  - per_instrument

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.22

    - 0.05

  - - normal

    - 1.215

    - 0.05

  - - normal

    - 1.24

    - 0.05

  normalization_mode:

  - per_instrument

  - multiplicative

  - multiplicative

  - additive

  - additive

  - additive

log_jitter:

- per_instrument

- - uniform

  - -13.82

  - -7.6

- - uniform

  - -13.82

  - -5.3

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

planet_b:

  amp_beam: 0.0

  amp_day_minus_amp_night: 0.0

  amp_ellip: 0.0

  amp_night: 0.0

  ecos: 0.0

  esin: 0.0

  impact:

  - global

  - uniform

  - 0.0

  - 1.2

  kRV:

  - global

  - uniform

  - 0.0

  - 1000.0

  period:

  - global

  - normal

  - 6.5693

  - 0.0001

  phoffset: 0.0

  rp_by_rs:

  - global

  - uniform

  - 0.01

  - 0.5

  tau:

  - global

  - normal

  - 0.0

  - 0.001

  tidallag: 0.0

reports:

  discard_extra_warmup: 0

  run: false

  trace_max_points: 500

star:

  q1:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.5

  - 0.5

  - 0.5

  q2:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.1

  - 0.1

  - 0.1

  rho_star:

  - global

  - uniform

  - 0.01

  - 12.0

paper_benchmarks/k2140_juliet/runs/nuts_epyc60/k2140.yaml

0 → 100644
+140 −0
Changes for paper_benchmarks/k2140_juliet/runs/nuts_epyc60/k2140.yaml: 140 added lines, 0 removed lines.
Original line number Diff line number Diff line 
data:

  file_path: k2140c.dat

  timeZeroPoint: 2457588.284

  working_dir: ./

instruments:

- K2

- LCOGT

- CORALIE

- FIES

- HARPS

data_type:

- LC

- LC

- RV

- RV

- RV

sampling:

  chain_energy_threshold: 5.0

  checkpoint_interval: 500

  internal_warmup_if_nuts: 500

  max_tree_depth: 10

  num_chains_if_nuts: 4

  num_samples_if_ensemble: 2000

  num_samples_if_nuts: 2000

  num_warmup: 2000

  production_sampler: auto

  random_seed: 1

  reset: true

  target_accept_prob: 0.9

  use_float64: false

  walkers_per_free_parameter: 3

  warmup_sampler: ess

detrending:

  normalization:

  - per_instrument

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.22

    - 0.05

  - - normal

    - 1.215

    - 0.05

  - - normal

    - 1.24

    - 0.05

  normalization_mode:

  - per_instrument

  - multiplicative

  - multiplicative

  - additive

  - additive

  - additive

log_jitter:

- per_instrument

- - uniform

  - -13.82

  - -7.6

- - uniform

  - -13.82

  - -5.3

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

planet_b:

  amp_beam: 0.0

  amp_day_minus_amp_night: 0.0

  amp_ellip: 0.0

  amp_night: 0.0

  ecos: 0.0

  esin: 0.0

  impact:

  - global

  - uniform

  - 0.0

  - 1.2

  kRV:

  - global

  - uniform

  - 0.0

  - 1000.0

  period:

  - global

  - normal

  - 6.5693

  - 0.0001

  phoffset: 0.0

  rp_by_rs:

  - global

  - uniform

  - 0.01

  - 0.5

  tau:

  - global

  - normal

  - 0.0

  - 0.001

  tidallag: 0.0

reports:

  discard_extra_warmup: 0

  run: false

  trace_max_points: 500

star:

  q1:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.5

  - 0.5

  - 0.5

  q2:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.1

  - 0.1

  - 0.1

  rho_star:

  - global

  - uniform

  - 0.01

  - 12.0

paper_benchmarks/k2140_juliet/runs/nuts_h100/k2140.yaml

0 → 100644
+140 −0
Changes for paper_benchmarks/k2140_juliet/runs/nuts_h100/k2140.yaml: 140 added lines, 0 removed lines.
Original line number Diff line number Diff line 
data:

  file_path: k2140c.dat

  timeZeroPoint: 2457588.284

  working_dir: ./

instruments:

- K2

- LCOGT

- CORALIE

- FIES

- HARPS

data_type:

- LC

- LC

- RV

- RV

- RV

sampling:

  chain_energy_threshold: 5.0

  checkpoint_interval: 500

  internal_warmup_if_nuts: 500

  max_tree_depth: 10

  num_chains_if_nuts: 4

  num_samples_if_ensemble: 2000

  num_samples_if_nuts: 2000

  num_warmup: 2000

  production_sampler: auto

  random_seed: 1

  reset: true

  target_accept_prob: 0.9

  use_float64: false

  walkers_per_free_parameter: 3

  warmup_sampler: ess

detrending:

  normalization:

  - per_instrument

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.22

    - 0.05

  - - normal

    - 1.215

    - 0.05

  - - normal

    - 1.24

    - 0.05

  normalization_mode:

  - per_instrument

  - multiplicative

  - multiplicative

  - additive

  - additive

  - additive

log_jitter:

- per_instrument

- - uniform

  - -13.82

  - -7.6

- - uniform

  - -13.82

  - -5.3

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

planet_b:

  amp_beam: 0.0

  amp_day_minus_amp_night: 0.0

  amp_ellip: 0.0

  amp_night: 0.0

  ecos: 0.0

  esin: 0.0

  impact:

  - global

  - uniform

  - 0.0

  - 1.2

  kRV:

  - global

  - uniform

  - 0.0

  - 1000.0

  period:

  - global

  - normal

  - 6.5693

  - 0.0001

  phoffset: 0.0

  rp_by_rs:

  - global

  - uniform

  - 0.01

  - 0.5

  tau:

  - global

  - normal

  - 0.0

  - 0.001

  tidallag: 0.0

reports:

  discard_extra_warmup: 0

  run: false

  trace_max_points: 500

star:

  q1:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.5

  - 0.5

  - 0.5

  q2:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.1

  - 0.1

  - 0.1

  rho_star:

  - global

  - uniform

  - 0.01

  - 12.0

paper_benchmarks/k2140_juliet/runs/nuts_rtx/k2140.yaml

0 → 100644
+140 −0
Changes for paper_benchmarks/k2140_juliet/runs/nuts_rtx/k2140.yaml: 140 added lines, 0 removed lines.
Original line number Diff line number Diff line 
data:

  file_path: k2140c.dat

  timeZeroPoint: 2457588.284

  working_dir: ./

instruments:

- K2

- LCOGT

- CORALIE

- FIES

- HARPS

data_type:

- LC

- LC

- RV

- RV

- RV

sampling:

  chain_energy_threshold: 5.0

  checkpoint_interval: 500

  internal_warmup_if_nuts: 500

  max_tree_depth: 10

  num_chains_if_nuts: 4

  num_samples_if_ensemble: 2000

  num_samples_if_nuts: 2000

  num_warmup: 2000

  production_sampler: auto

  random_seed: 1

  reset: true

  target_accept_prob: 0.9

  use_float64: false

  walkers_per_free_parameter: 3

  warmup_sampler: ess

detrending:

  normalization:

  - per_instrument

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.0

    - 0.1

  - - normal

    - 1.22

    - 0.05

  - - normal

    - 1.215

    - 0.05

  - - normal

    - 1.24

    - 0.05

  normalization_mode:

  - per_instrument

  - multiplicative

  - multiplicative

  - additive

  - additive

  - additive

log_jitter:

- per_instrument

- - uniform

  - -13.82

  - -7.6

- - uniform

  - -13.82

  - -5.3

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

- - uniform

  - -2.3

  - 4.61

planet_b:

  amp_beam: 0.0

  amp_day_minus_amp_night: 0.0

  amp_ellip: 0.0

  amp_night: 0.0

  ecos: 0.0

  esin: 0.0

  impact:

  - global

  - uniform

  - 0.0

  - 1.2

  kRV:

  - global

  - uniform

  - 0.0

  - 1000.0

  period:

  - global

  - normal

  - 6.5693

  - 0.0001

  phoffset: 0.0

  rp_by_rs:

  - global

  - uniform

  - 0.01

  - 0.5

  tau:

  - global

  - normal

  - 0.0

  - 0.001

  tidallag: 0.0

reports:

  discard_extra_warmup: 0

  run: false

  trace_max_points: 500

star:

  q1:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.5

  - 0.5

  - 0.5

  q2:

  - per_instrument

  - - uniform

    - 0

    - 1

  - - uniform

    - 0

    - 1

  - 0.1

  - 0.1

  - 0.1

  rho_star:

  - global

  - uniform

  - 0.01

  - 12.0

paper_benchmarks/machines/rerun_round3_laptop.sh

0 → 100755
+44 −0
Changes for paper_benchmarks/machines/rerun_round3_laptop.sh: 44 added lines, 0 removed lines.
Original line number Diff line number Diff line 
#!/bin/bash

# 2026-10-02: laptop reruns after the second code round (4 slice steps per parameter, RM anomaly averaged over the

# exposures, NUTS chains from distinct walkers; commit 1c53d79 or later) and the new NUTS benchmarks: every

# nested-sampling run and the NUTS runs. Ensemble runs are unchanged by this round. One job at a time; each run

# refuses to start with capped clocks or outside the performance profile, and records both in power_profile.txt.

# Takes about 3 h; start it on mains, and do not suspend the laptop meanwhile.

set -u

B=$(realpath "$(dirname "$0")/..")

LOG=/dati/software/pyLib/PyAETNA_bench_work/logs/rerun_round3_laptop.log

mkdir -p "$(dirname "$LOG")"

export PATH=/dati/pyenvs/pyaetnaEnv/bin:$PATH

say() { echo "$(date '+%F %T') $*" >> "$LOG"; }

check() {

  local f p ac

  f=$(cat /sys/devices/system/cpu/cpu0/cpufreq/scaling_max_freq); p=$(powerprofilesctl get)

  ac=$(cat /sys/class/power_supply/A*/online 2>/dev/null | head -1)

  if [ "$f" -lt 4000000 ] || [ "$p" != performance ] || [ "$ac" != 1 ]; then

    say "REFUSED: max_freq=$f profile=$p ac=$ac"; exit 1

  fi

  echo "profile=$p max_freq_khz=$f ac=$ac $(date '+%F %T')" > "$B/$1/power_profile.txt"

}

run() {  # folder platforms cpus

  check "$1"

  say "start $1 ($2, cpus $3)"

  (cd "$B" && taskset -c "$3" nice -n 19 ./run.sh "$1" "$2" 24G > /dev/null)

  check "$1"  # cap after the run invalidates it: logged by check

  say "end $1 $(cat "$B/$1/done.txt") $(grep 'Elapsed (wall' "$B/$1/time.out")"

}

say "code: $(git -C "$B/.." rev-parse --short HEAD) + diff sha256 $(git -C "$B/.." diff HEAD -- src | sha256sum | cut -c1-16)"

for r in k2140_juliet/runs/k2c_e8_rtx pimen_allesfitter/runs/circ_sqrt_rtx pimen_allesfitter/runs/ecc_sqrt_rtx \

         pimen_allesfitter/runs/circ_sqrt pimen_allesfitter/runs/ecc_sqrt pimen_allesfitter/runs/circ pimen_allesfitter/runs/ecc \

         wasp18_phase_curve/runs/ns_gpu_rtx wasp18_phase_curve/runs/ns_gpu k2100_pyaneti/runs/ns_rtx \

         toi5398_rm_pyorbit/runs/ccf_rtx toi5398_rm_pyorbit/runs/template_rtx toi5398_rm_pyorbit/runs/ccf \

         toi5398_rm_pyorbit/runs/template toi5398_rm_pyorbit/runs/ns toi5398_rm_pyorbit/runs/linear toi216_ttv/runs/ns_rtx \

         k2140_juliet/runs/nuts_rtx wasp18_phase_curve/runs/nuts_auto_rtx; do

  run "$r" cuda,cpu 0-23

done

for r in k2140_juliet/runs/k2c_e8_cpu24 pimen_allesfitter/runs/circ_sqrt_cpu24 pimen_allesfitter/runs/ecc_sqrt_cpu24 \

         wasp18_phase_curve/runs/ns_cpu24 wasp18_phase_curve/runs/ns_cpu k2100_pyaneti/runs/ns_cpu24 \

         toi5398_rm_pyorbit/runs/ccf_cpu24 k2140_juliet/runs/nuts_cpu24 wasp18_phase_curve/runs/nuts_auto_cpu24; do

  run "$r" cpu 0-23

done

run toi216_ttv/runs/ns_cpu12 cpu 0-11

say "done"