29:null 2a:["$","title",null,{"children":"Metalog Distribution - Range Slider | I001962 | Observable"}] 2f:T734,function newtonRaphson(f, fp, x0, options) { var x1, y, yp, tol, maxIter, iter, yph, ymh, yp2h, ym2h, h, hr, verbose, eps; // Iterpret variadic forms: // fp is first derivative of - currently not using or tested with metalog if (typeof fp !== 'function') { console.log("fp is not a function"); options = x0; x0 = fp; fp = null; } options = options || {}; tol = options.tolerance === undefined ? 1e-7 : options.tolerance; eps = options.epsilon === undefined ? 2.220446049250313e-16 : options.epsilon; maxIter = options.maxIterations === undefined ? 20 : options.maxIterations; h = options.h === undefined ? 1e-4 : options.h; verbose = options.verbose === undefined ? false : options.verbose; hr = 1 / h; iter = 0; while (iter++ < maxIter) { // Compute the value of the function: y = f(x0); if (fp) { yp = fp(x0); } else { // Needs numerical derivatives: yph = f(x0 + h); ymh = f(x0 - h); yp2h = f(x0 + 2 * h); ym2h = f(x0 - 2 * h); yp = ((ym2h - yp2h) + 8 * (yph - ymh)) * hr / 12; } // Check for badly conditioned update (extremely small first deriv relative to function): if (Math.abs(yp) <= eps * Math.abs(y)) { if (verbose) { console.log('Newton-Raphson: failed to converged due to nearly zero first derivative'); } return false; } // Update the guess: x1 = x0 - y / yp; // Check for convergence: if (Math.abs(x1 - x0) <= tol * Math.abs(x1)) { if (verbose) { console.log('Newton-Raphson: converged to x = ' + x1 + ' after ' + iter + ' iterations'); } return x1; } // Transfer update to the new guess: x0 = x1; } if (verbose) { console.log('Newton-Raphson: Maximum iterations reached (' + maxIter + ')'); } return false; } 30:T5ec,function revMetalogUpper(y) { function convert_to_float(a) { var floatValue = parseFloat(a); return floatValue; } let vector = []; let np_a = a; for (let index = 1; index < a.length + 1; ++index) { //cant start with 0 coeeffs for each aCoeff vector.push(basis(y, index)); } let wrappedVector = [vector]; let wrappedNp_x = [np_a]; let wrappedNp_a = wrappedNp_x[0].map((e) => [e]); let mky = multiply(wrappedVector, wrappedNp_a); // Unbounded if (typeof bl == String && typeof bu == String) { return mky; } if ( (typeof bl === "string" || bl instanceof String) && (typeof bu === "string" || bu instanceof String) ) { console.log("mky ", mky, mky - targetValueUpper); return mky - targetValueUpper; } // Bounded lower else if (typeof bl !== "string" && typeof bu == "string") { convert_to_float(bl); return bl + Math.exp(mky); } // Bounded upper else if (typeof bl == "string" && typeof bu != "string") { convert_to_float(bu); return bu - Math.exp(-mky); } // Bounded else if (typeof bl != "string" && typeof bu != "string") { return bl + (bu * Math.exp(mky)) / (1 + Math.exp(mky)); } } // Test cases for newtonRaphson //console.log('Using the derivative:', newtonRaphson(f, fp, 42, {verbose: true})); // console.log('Using numerical first derivative:', newtonRaphson(f, 94, {verbose: true})); // Metalog workzone - Enter aCoefs and targetValue then set bl and bu and finally run metalog with an initial guess31:T608,//function metalog(y, a, bl = "", bu = "") function revMetalog(y) { function convert_to_float(a) { var floatValue = parseFloat(a); return floatValue; } let vector = []; let np_a = a; for (let index = 1; index < a.length + 1; ++index) { //cant start with 0 coeeffs for each aCoeff vector.push(basis(y, index)); } let wrappedVector = [vector]; let wrappedNp_x = [np_a]; let wrappedNp_a = wrappedNp_x[0].map((e) => [e]); let mky = multiply(wrappedVector, wrappedNp_a); // Unbounded if (typeof bl == String && typeof bu == String) { return mky; } if ( (typeof bl === "string" || bl instanceof String) && (typeof bu === "string" || bu instanceof String) ) { console.log("mky ", mky, mky - targetValue); return mky - targetValue; } // Bounded lower else if (typeof bl !== "string" && typeof bu == "string") { convert_to_float(bl); return bl + Math.exp(mky); } // Bounded upper else if (typeof bl == "string" && typeof bu != "string") { convert_to_float(bu); return bu - Math.exp(-mky); } // Bounded else if (typeof bl != "string" && typeof bu != "string") { return bl + (bu * Math.exp(mky)) / (1 + Math.exp(mky)); } } // Test cases for newtonRaphson //console.log('Using the derivative:', newtonRaphson(f, fp, 42, {verbose: true})); // console.log('Using numerical first derivative:', newtonRaphson(f, 94, {verbose: true})); // Metalog workzone - Enter aCoefs and targetValue then set bl and bu and finally run metalog with an initial guess32:T46d,function metalog(y, a, bl = "", bu = "") { function convert_to_float(a) { var floatValue = parseFloat(a); return floatValue; } let vector = []; let np_a = a; for (let index = 1; index < a.length + 1; ++index) { //cant start with 0 coeeffs for each aCoeff vector.push(basis(y, index)); } let wrappedVector = [vector]; let wrappedNp_x = [np_a]; let wrappedNp_a = wrappedNp_x[0].map((e) => [e]); let mky = multiply(wrappedVector, wrappedNp_a); // Unbounded if (typeof bl == String && typeof bu == String) { return mky; } if ( (typeof bl === "string" || bl instanceof String) && (typeof bu === "string" || bu instanceof String) ) { return mky; } // Bounded lower else if (typeof bl !== "string" && typeof bu == "string") { convert_to_float(bl); return bl + Math.exp(mky); } // Bounded upper else if (typeof bl == "string" && typeof bu != "string") { convert_to_float(bu); return bu - Math.exp(-mky); } // Bounded else if (typeof bl != "string" && typeof bu != "string") { return bl + (bu * Math.exp(mky)) / (1 + Math.exp(mky)); } }33:Tea0,function rangeInput(options = {}) { const { min = 0, max = 100, step = 1, value = [min, max], theme = theme_Flat, } = options; const controls = {}; const scope = randomScope(); // Will be used to sanitize values while avoiding floating point issues. const input = DOM.element('input', {type: 'range', min, max, step}); const dom = html`

${controls.track = html`

${controls.zone = html`

${controls.range = html`

${controls.min = html`

`} ${controls.max = html`

`} `} `} `} ${html`