//One universal basic required here to get things going once loaded
window.onload = function () {
    //restoreDefaultValues(); //Un-comment this if you want to start with defaults
    Main();
};
//Any global variables go here


//Main is hard wired as THE place to start calculating when input changes
//It does no calculations itself, it merely sets them up, sends off variables, gets results and, if necessary, plots them.
function Main() {
    saveSettings();

    //Send all the inputs as a structured object
    //If you need to convert to, say, SI units, do it here!
    const inputs = {
        A: sliders.SlideA.value / 100,
        Trel: sliders.SlideTrel.value,
        Width: (6 - sliders.SlideWidth.value) * 3,
        xStart: sliders.SlidexStart.value,
    }

    //Send inputs off to CalcIt where the names are instantly available
    //Get all the resonses as an object, result
    const result = CalcIt(inputs)

    //Set all the text box outputs
    document.getElementById('Results').value = result.Results
    //Do all relevant plots by calling plotIt - if there's no plot, nothing happens
    //plotIt is part of the app infrastructure in app.new.js
    if (result.plots) {
        for (let i = 0; i < result.plots.length; i++) {
            plotIt(result.plots[i], result.canvas[i]);
        }
    }

    //You might have some other stuff to do here, but for most apps that's it for CalcIt!
}

//Here's the app calculation
//The inputs are just the names provided - their order in the curly brackets is unimportant!
//By convention the input values are provided with the correct units within Main
function CalcIt({ A, Trel, Width, xStart }) {

    let plotData = [], pPts = [], tPts = [], xPts = [], dPts1 = [], dPts2 = [], i, j
    let x1 = 0, x2 = 1
    const iStep = 0.01
    for (i = 0; i <= 1.01; i += iStep) {
        y=getY(i,A,Width)
        pPts.push({ x: i, y: y })
        //Inelegant code to get more accurate x1, x2 values
        if (i>0 && y >= Trel && x1 == 0) {
            for (j = i - iStep; j <= i; j += 0.00001) {
                y=getY(j,A,Width)
                if (y >= Trel) { x1 = j;break }
            }
         }
        if (i>xStart && y <= Trel && x2 ==1){
            for (j = i - iStep; j <= i; j += 0.00001) {
                y=getY(j,A,Width)
                if (y <= Trel) { x2 = j;break }
            }
        }
    }
    if (x1<=0.015) x1=0
    if (x2>=0.990) x2=1
    const xg1=x1 //x1>0?x1:-0.5
    const xg2=x2//<1?x2:1.5

    tPts.push({ x: 0, y: Trel }, { x: 1, y: Trel })
    xPts.push({ x: xStart, y: 0 }, { x: xStart, y: 100 })
    dPts1.push({ x: xg1, y: 0 }, { x: xg1, y: Trel })
    dPts2.push({ x: xg2, y: 0 }, { x: xg2, y: Trel })
    x1=Math.round(x1 * 100+Number.EPSILON) / 100
    x2=Math.round(x2 * 100+Number.EPSILON) / 100
    console.log(x1,x2,xStart)
    if (xStart < x1 || xStart > x2 || Trel>75) {
        Results = "Single phase"
        dPts1=[];dPts2=[]
    } else {
        const L1 = 100 * (x2 - xStart) / (x2 - x1)
        const L2 = 100 * (xStart - x1) / (x2 - x1)
        const N1 = (x2 - xStart) / (xStart - x1)
        const N2 = (xStart - x1) / (x2 - xStart)
        const N12 = N1 + N2
        const N1N12 = N1 / N12
        const N2N12 = 100 * N2 / N12
        Results = "Separates into:\nSolvent Rich Phase x=" + x1.toFixed(2) + " & Solute Rich Phase x= " + x2.toFixed(2) + "\nat ratio\nSolvent Rich : Solute Rich = " + L1.toFixed(0) + " : " + L2.toFixed(0) + "\n"
        Results += "\nSolvent Rich/Solute Rich = (" + x2.toFixed(2) + "-" + xStart.toFixed(2) + ")/(" + xStart.toFixed(2) + "-" + x1.toFixed(2) + ") =" + N1.toFixed(2)
        Results += "\nSolute Rich/Solvent Rich = (" + xStart.toFixed(2) + "-" + x1.toFixed(2) + ")/(" + x2.toFixed(2) + "-" + xStart.toFixed(2) + ") =" + N2.toFixed(2)
        Results+="\nSolvent Rich / (Solute Rich + Solvent Rich) = (" + x2.toFixed(2) + "-" + xStart.toFixed(2) + ")/(" + x2.toFixed(2) + "-" + x1.toFixed(2) + ") =" + (L1/100).toFixed(2)
    }

    plotData.push(pPts, tPts, xPts, dPts1, dPts2)
    const prmap = {
        plotData: plotData, //An array of 1 or more datasets
        lineLabels: ["Phase Boundary", "T", "Init.", "Solv Rich", "Solute Rich"], //An array of labels for each dataset
        //hideLegend: true,
        dottedLine: [false, true, true, true, true],
        // colors: theColors,
        // showLines: showLines,
        // showPoints: showPoints,
        xLabel: "x_solute& ", //Label for the x axis, with an & to separate the units
        yLabel: "T (relative)& ", //Label for the y axis, with an & to separate the units
        y2Label: null, //Label for the y2 axis, null if not needed
        yAxisL1R2: [], //Array to say which axis each dataset goes on. Blank=Left=1
        logX: false, //Is the x-axis in log form?
        xTicks: undefined, //We can define a tick function if we're being fancy
        logY: false, //Is the y-axis in log form?
        yTicks: undefined, //We can define a tick function if we're being fancy
        legendPosition: 'top', //Where we want the legend - top, bottom, left, right
        xMinMax: [0, 1], //Set min and max, e.g. [-10,100], leave one or both blank for auto
        yMinMax: [0, 100], //Set min and max, e.g. [-10,100], leave one or both blank for auto
        y2MinMax: [,], //Set min and max, e.g. [-10,100], leave one or both blank for auto
        xSigFigs: 'F2', //These are the sig figs for the Tooltip readout. A wide choice!
        ySigFigs: 'F0', //F for Fixed, P for Precision, E for exponential
    };
    return {
        Results: Results,
        plots: [prmap],
        canvas: ['canvas'],
    };
}
function getY(x,A,Width){
    let y=0
    if (x <= A) {
        y = 75 * (1 - Math.pow(Width * Math.pow(A, 1.5) * (A - x), 2))
   } else {
       y = 75 * (1 - Math.pow(Width * Math.pow(1 - A, 1.5) * (A - x), 2))
   }
return y
}