INNER CODE UNIT · JavaScript

tmp

lh3/psmc · utils/split-time.js:158

			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];
}

function main(args) {
	let opt = { round:20, d0:1e-5, scale:100 };

View source record →