INNER CODE UNIT · JavaScript
tao
lh3/psmc · utils/split-time.js:156
let tao = r.a[i+1].t - r.a[i].t;
let x = tao / r.a[i].n;
let tmp = 1.0 / x - 1.0 / (Math.exp(x) - 1);
if (tmp < 0.0) {
warn(`Warning: approximation failure at interval ${i}: ${tmp} < 0 (x=${x})`);
tmp = 0.0;
} else if (tmp >= 1.0) {
warn(`Warning: approximation failure at interval ${i}: ${tmp} >= 1 (x=${x})`);
tmp = 1.0;
}
t = tao * tmp;
}
f += r.a[i].sigma1;
y += (1.0 - Math.exp(-r.theta0 * (t + (r.a[i].t - r.a[k0].t)) / scale)) * r.a[i].sigma1;
}
return [z, y/f];
}