machine   Linux x86_64, 4 CPUs, python 3.12.3, numpy 2.5.4
rietx     1.9.0.dev0 from /home/user/rietx/.claude/worktrees/wp1325-parametric-series/src/rietx ([dev] venv), compiled kernels on
load      1.02 0.48 0.22 (1/5/15 min) at start; 'python -m pytest' processes running: 0
ramp      10 patterns, T = 300..1200 K; a(T) = 4.1566*(1 + 0.0005*(T - 300)/100); T_ref = 750 K, T_half = 450 K; seeds 0..19
units     a1 in 1e-09 Å/K (truth 20783.00); every esd as reported, Bérar-Lelann factor included; |d|/sj = |joint - line(w)| / joint esd

R1: strong signal; phase + instrument per pattern, only a0/a1 shared
  seed  chain s  joint s  ratio |   line(w) ±       esd chi2r |   line(u) ±       esd |     joint ±       esd |d|/sj
     0     1.29     2.84   2.19 | 20780.81 ±     5.45  0.68 | 20780.65 ±     4.49 | 20780.66 ±     5.49  0.027
     1     1.43     2.24   1.56 | 20773.21 ±     5.46  0.64 | 20773.21 ±     4.37 | 20773.19 ±     5.45  0.004
     2     1.50     2.36   1.57 | 20780.90 ±     5.47  0.42 | 20780.84 ±     3.52 | 20780.84 ±     5.46  0.012
     3     1.44     2.38   1.65 | 20782.93 ±     5.43  0.35 | 20782.99 ±     3.23 | 20782.99 ±     5.42  0.011
     4     1.46     2.45   1.68 | 20777.54 ±     5.44  0.81 | 20777.91 ±     4.95 | 20777.91 ±     5.49  0.068
     5     1.32     2.26   1.71 | 20780.50 ±     5.49  0.73 | 20780.33 ±     4.62 | 20780.34 ±     5.46  0.030
     6     1.91     2.68   1.41 | 20785.12 ±     5.42  0.38 | 20785.09 ±     3.34 | 20785.08 ±     5.43  0.008
     7     1.36     2.55   1.87 | 20780.06 ±     5.50  0.29 | 20780.15 ±     2.97 | 20780.15 ±     5.50  0.017
     8     1.37     2.57   1.87 | 20778.21 ±     5.45  0.86 | 20778.10 ±     5.03 | 20778.09 ±     5.45  0.021
     9     1.46     2.49   1.71 | 20782.25 ±     5.51  0.72 | 20782.43 ±     4.60 | 20782.42 ±     5.48  0.032
    10     1.39     2.52   1.81 | 20784.00 ±     5.51  0.64 | 20784.16 ±     4.38 | 20784.15 ±     5.50  0.027
    11     1.50     2.60   1.73 | 20785.63 ±     5.53  0.39 | 20785.79 ±     3.40 | 20785.79 ±     5.50  0.029
    12     1.56     2.58   1.66 | 20783.47 ±     5.55  0.49 | 20783.07 ±     3.89 | 20783.07 ±     5.45  0.074
    13     1.51     2.49   1.66 | 20790.87 ±     5.42  0.86 | 20790.65 ±     5.09 | 20790.66 ±     5.47  0.038
    14     1.42     2.43   1.71 | 20789.51 ±     5.48  0.49 | 20789.31 ±     3.82 | 20789.32 ±     5.50  0.033
    15     1.45     2.63   1.81 | 20780.39 ±     5.46  0.61 | 20780.35 ±     4.28 | 20780.34 ±     5.49  0.009
    16     1.47     2.43   1.65 | 20784.93 ±     5.47  0.62 | 20784.94 ±     4.27 | 20784.94 ±     5.43  0.002
    17     1.38     2.26   1.64 | 20787.82 ±     5.56  0.73 | 20787.86 ±     4.65 | 20787.87 ±     5.50  0.010
    18     1.37     2.72   1.99 | 20786.43 ±     5.50  0.24 | 20786.53 ±     2.70 | 20786.52 ±     5.52  0.017
    19     1.45     2.40   1.66 | 20776.38 ±     5.54  0.27 | 20776.48 ±     2.81 | 20776.48 ±     5.54  0.018
  joint θ: 102 columns (2 shared: vars.a0, vars.a1) x 42000 rows; cell stage freed ['vars.a0', 'vars.a1']
  summary R1, 20 seeds; a1 in 1e-09 Å/K, truth 20783.00
    estimator                 mean      bias   scatter  mean esd  scatter/esd  rms z   min esd   max esd
    chain line, weighted  20782.55     -0.45      4.45      5.48         0.81   0.80      5.42      5.56
    chain line, OLS       20782.54     -0.46      4.41      4.02         1.10   1.09      2.70      5.09
    joint a0 + a1*t       20782.54     -0.46      4.41      5.48         0.80   0.79      5.42      5.54
    a(T_ref) in 1e-06 Å against truth 4.165952 Å: estimator, mean bias, mean esd, rms z
      chain line, weighted    0.223    1.572   0.77
      joint a0 + a1*t         0.207    1.573   0.76
    joint esd / weighted-line esd: mean 0.999, range 0.983-1.008
    |joint - weighted line| / joint esd: mean 0.024, max 0.074
    weighted line chi2_red: mean 0.56, range 0.24-0.86
    per-pattern cell esd (chain), mean over patterns and seeds: 4.976e-06 Å
    Bérar-Lelann factor: chain (mean per pattern) 1.263, joint 1.263
    wall clock chain  :   1.29-  1.91 s
    wall clock joint  :   2.24-  2.84 s
    wall clock joint / chain: 1.41-2.19
    not converged: none; chain cold restarts after pattern 0: 0
    absent-phase / ridge holds fired in a joint fit: none

R2: weak signal (phase scale x0.035); as R1 otherwise
  seed  chain s  joint s  ratio |   line(w) ±       esd chi2r |   line(u) ±       esd |     joint ±       esd |d|/sj
     0     1.50     2.56   1.70 | 20778.25 ±    27.02  1.02 | 20776.45 ±    27.66 | 20776.56 ±    27.31  0.062
     1     1.36     2.88   2.11 | 20801.13 ±    27.15  0.33 | 20801.19 ±    15.46 | 20801.13 ±    27.01  0.000
     2     1.31     2.49   1.90 | 20769.06 ±    27.05  0.49 | 20770.71 ±    18.78 | 20770.50 ±    26.98  0.054
     3     1.43     2.41   1.68 | 20798.57 ±    27.21  2.06 | 20799.23 ±    38.82 | 20799.27 ±    27.07  0.026
     4     1.42     2.46   1.73 | 20756.53 ±    27.31  0.31 | 20756.61 ±    15.12 | 20756.44 ±    27.29  0.003
     5     1.37     2.45   1.78 | 20755.55 ±    26.91  0.90 | 20755.01 ±    25.63 | 20755.04 ±    27.01  0.019
     6     1.45     2.68   1.84 | 20802.25 ±    26.93  0.76 | 20801.15 ±    23.34 | 20801.25 ±    26.85  0.037
     7     1.42     2.37   1.66 | 20783.72 ±    27.27  1.24 | 20783.48 ±    30.17 | 20783.00 ±    27.15  0.027
     8     1.56     2.66   1.70 | 20795.33 ±    27.20  0.72 | 20795.41 ±    23.19 | 20795.35 ±    27.22  0.001
     9     1.37     2.54   1.86 | 20787.06 ±    27.01  0.83 | 20785.88 ±    24.47 | 20785.88 ±    27.05  0.044
    10     1.31     2.58   1.97 | 20805.00 ±    27.16  0.68 | 20805.40 ±    22.53 | 20805.28 ±    27.18  0.010
    11     1.35     2.53   1.87 | 20790.01 ±    27.45  0.19 | 20790.34 ±    11.87 | 20790.36 ±    27.36  0.013
    12     1.45     2.70   1.86 | 20755.69 ±    27.39  0.37 | 20754.24 ±    16.49 | 20754.53 ±    27.08  0.043
    13     1.41     2.68   1.90 | 20780.34 ±    26.98  1.39 | 20779.67 ±    31.93 | 20779.56 ±    27.02  0.029
    14     1.50     2.93   1.95 | 20763.31 ±    27.18  1.22 | 20763.03 ±    29.46 | 20762.93 ±    27.08  0.014
    15     1.31     2.65   2.01 | 20787.13 ±    26.79  0.52 | 20787.09 ±    19.53 | 20787.10 ±    26.94  0.001
    16     1.40     2.74   1.96 | 20788.54 ±    26.83  0.37 | 20788.01 ±    16.24 | 20788.19 ±    26.87  0.013
    17     1.44     2.71   1.88 | 20788.37 ±    27.46  1.21 | 20788.68 ±    29.94 | 20788.62 ±    27.28  0.009
    18     1.47     2.87   1.96 | 20783.08 ±    27.14  0.15 | 20783.15 ±    10.40 | 20783.04 ±    27.23  0.001
    19     1.44     2.73   1.90 | 20742.17 ±    27.09  0.78 | 20742.13 ±    24.17 | 20742.47 ±    27.29  0.011
  joint θ: 102 columns (2 shared: vars.a0, vars.a1) x 42000 rows; cell stage freed ['vars.a0', 'vars.a1']
  summary R2, 20 seeds; a1 in 1e-09 Å/K, truth 20783.00
    estimator                 mean      bias   scatter  mean esd  scatter/esd  rms z   min esd   max esd
    chain line, weighted  20780.55     -2.45     17.89     27.13         0.66   0.65     26.79     27.46
    chain line, OLS       20780.34     -2.66     18.00     22.76         0.79   0.87     10.40     38.82
    joint a0 + a1*t       20780.33     -2.67     17.95     27.11         0.66   0.65     26.85     27.36
    a(T_ref) in 1e-06 Å against truth 4.165952 Å: estimator, mean bias, mean esd, rms z
      chain line, weighted   -0.211    7.786   0.56
      joint a0 + a1*t        -0.277    7.790   0.56
    joint esd / weighted-line esd: mean 1.000, range 0.988-1.011
    |joint - weighted line| / joint esd: mean 0.021, max 0.062
    weighted line chi2_red: mean 0.78, range 0.15-2.06
    per-pattern cell esd (chain), mean over patterns and seeds: 2.463e-05 Å
    Bérar-Lelann factor: chain (mean per pattern) 1.137, joint 1.137
    wall clock chain  :   1.31-  1.56 s
    wall clock joint  :   2.37-  2.93 s
    wall clock joint / chain: 1.66-2.11
    not converged: none; chain cold restarts after pattern 0: 0
    absent-phase / ridge holds fired in a joint fit: none
    per-pattern cell esd R2/R1: 4.95

R3: shared instrument (zero_shift, profile.*) in the joint fit; the chain per pattern
  seed  chain s  joint s  ratio |   line(w) ±       esd chi2r |   line(u) ±       esd |     joint ±       esd |d|/sj |   control ±       esd
     0     1.63     2.63   1.62 | 20780.81 ±     5.45  0.68 | 20780.65 ±     4.49 | 20783.65 ±     2.13  1.334 | 20783.65 ±     2.13
     1     1.55     2.88   1.86 | 20773.21 ±     5.46  0.64 | 20773.21 ±     4.37 | 20780.71 ±     2.11  3.548 | 20780.72 ±     2.12
     2     1.53     2.62   1.71 | 20780.90 ±     5.47  0.42 | 20780.84 ±     3.52 | 20785.22 ±     2.12  2.036 | 20785.22 ±     2.12
     3     1.46     3.00   2.05 | 20782.93 ±     5.43  0.35 | 20782.99 ±     3.23 | 20783.92 ±     2.10  0.472 | 20783.92 ±     2.10
     4     1.47     2.77   1.88 | 20777.54 ±     5.44  0.81 | 20777.91 ±     4.95 | 20784.41 ±     2.12  3.235 | 20784.41 ±     2.12
     5     1.50     2.78   1.85 | 20780.50 ±     5.49  0.73 | 20780.33 ±     4.62 | 20782.25 ±     2.12  0.821 | 20782.25 ±     2.12
     6     1.37     3.23   2.36 | 20785.12 ±     5.42  0.38 | 20785.09 ±     3.34 | 20783.96 ±     2.11  0.550 | 20783.96 ±     2.11
     7     1.59     2.76   1.73 | 20780.06 ±     5.50  0.29 | 20780.15 ±     2.97 | 20782.61 ±     2.13  1.196 | 20782.61 ±     2.13
     8     1.44     2.54   1.76 | 20778.21 ±     5.45  0.86 | 20778.10 ±     5.03 | 20782.92 ±     2.11  2.229 | 20782.92 ±     2.12
     9     1.57     2.48   1.58 | 20782.25 ±     5.51  0.72 | 20782.43 ±     4.60 | 20778.98 ±     2.13  1.534 | 20778.98 ±     2.13
    10     1.63     2.43   1.49 | 20784.00 ±     5.51  0.64 | 20784.16 ±     4.38 | 20783.81 ±     2.14  0.087 | 20783.81 ±     2.13
    11     1.37     2.69   1.96 | 20785.63 ±     5.53  0.39 | 20785.79 ±     3.40 | 20781.67 ±     2.14  1.854 | 20781.67 ±     2.13
    12     1.47     2.26   1.54 | 20783.47 ±     5.55  0.49 | 20783.07 ±     3.89 | 20784.49 ±     2.12  0.482 | 20784.49 ±     2.11
    13     1.44     2.37   1.64 | 20790.87 ±     5.42  0.86 | 20790.65 ±     5.09 | 20782.48 ±     2.12  3.949 | 20782.48 ±     2.12
    14     1.46     2.96   2.03 | 20789.51 ±     5.48  0.49 | 20789.31 ±     3.82 | 20785.85 ±     2.13  1.716 | 20785.85 ±     2.13
    15     1.51     2.53   1.67 | 20780.39 ±     5.46  0.61 | 20780.35 ±     4.28 | 20783.18 ±     2.14  1.307 | 20783.18 ±     2.13
    16     1.41     2.59   1.84 | 20784.93 ±     5.47  0.62 | 20784.94 ±     4.27 | 20783.46 ±     2.10  0.700 | 20783.46 ±     2.10
    17     1.34     2.57   1.92 | 20787.82 ±     5.56  0.73 | 20787.86 ±     4.65 | 20781.91 ±     2.14  2.768 | 20781.90 ±     2.13
    18     1.67     2.41   1.44 | 20786.43 ±     5.50  0.24 | 20786.53 ±     2.70 | 20784.02 ±     2.14  1.125 | 20784.02 ±     2.14
    19     1.45     2.57   1.78 | 20776.38 ±     5.54  0.27 | 20776.48 ±     2.81 | 20782.56 ±     2.16  2.872 | 20782.56 ±     2.15
  joint θ: 48 columns (8 shared: instrument.zero_shift, instrument.profile.u, instrument.profile.v, instrument.profile.w, instrument.profile.x, instrument.profile.y, vars.a0, vars.a1) x 42000 rows; cell stage freed ['vars.a0', 'vars.a1']
  summary R3, 20 seeds; a1 in 1e-09 Å/K, truth 20783.00
    estimator                 mean      bias   scatter  mean esd  scatter/esd  rms z   min esd   max esd
    chain line, weighted  20782.55     -0.45      4.45      5.48         0.81   0.80      5.42      5.56
    chain line, OLS       20782.54     -0.46      4.41      4.02         1.10   1.09      2.70      5.09
    joint a0 + a1*t       20783.10      0.10      1.57      2.13         0.74   0.72      2.10      2.16
    control (no tie)      20783.10      0.10      1.57      2.12         0.74   0.72      2.10      2.15
    a(T_ref) in 1e-06 Å against truth 4.165952 Å: estimator, mean bias, mean esd, rms z
      chain line, weighted    0.223    1.572   0.77
      joint a0 + a1*t         0.207    1.574   0.76
      control (no tie)        0.207    1.574   0.76
    joint esd / weighted-line esd: mean 0.388, range 0.381-0.392
    |joint - weighted line| / joint esd: mean 1.691, max 3.949
    weighted line chi2_red: mean 0.56, range 0.24-0.86
    per-pattern cell esd (chain), mean over patterns and seeds: 4.976e-06 Å
    Bérar-Lelann factor: chain (mean per pattern) 1.263, joint 1.264
    |joint - control| / joint esd: mean 0.000, max 0.001; control esd / joint esd 0.997-1.001
    wall clock chain  :   1.34-  1.67 s
    wall clock joint  :   2.26-  3.23 s
    wall clock control:   2.28-  3.15 s
    wall clock joint / chain: 1.44-2.36
    not converged: none; chain cold restarts after pattern 0: 0
    absent-phase / ridge holds fired in a joint fit: none
    per-pattern cell esd R3/R1: 1.00

