macrorugo.spec.ts

import { ParamCalculability, ParamValueMode, Session } from "../../src/index";
import { MacroRugo, MacroRugoFlowType } from "../../src/macrorugo/macrorugo";
import { MacrorugoParams } from "../../src/macrorugo/macrorugo_params";
import { checkResult } from "../test_func";

/*
*** Emergent conditions Cd=1.5***
ks=0.010000
D=0.500000
k=0.700000
Cd0=1.500000
S=0.050000
B=1.000000
h=0.600000
C=0.130000
Q=0.493722  fVal=0.000000
cote_bas = 12.200000
Vdeb = 0.822870
Fr = 0.530417
P = 242.170537
flowcond = emergent
Vmax = 1.799472
V_technique = 1.991299
q_technique = 0.561860
 */

function macroRugoInstanceEmergentCd15(): MacroRugo {
    return new MacroRugo(
        new MacrorugoParams(
            12.5,   // ZF1
            6,      // L
            1,      // B
            0.05,   // If
            0.493722,    // Q
            0.6,    // h
            0.01,   // Ks
            0.13,   // C
            0.5,    // D
            0.7,    // k
            1.5     // Cd0
        )
    );
}

const macroRugoExtraResultEmergentCd15: { [key: string]: number|MacroRugoFlowType } = {
    ENUM_MacroRugoFlowType: MacroRugoFlowType.EMERGENT,
    Fr: 0.530425,
    PV: 242.174003,
    Q_GuideTech: 0.561860,
    V_GuideTech: 1.991299,
    Vdeb: 0.822881,
    Vmax: 1.799502,
    ZF2: 12.2
};

/*
*** Emergent conditions Cd=2***
ks=0.010000
D=0.500000
k=0.700000
Cd0=2.000000
S=0.050000
B=1.000000
h=0.600000
C=0.130000
Q=0.376808  fVal=0.000000
cote_bas = 12.200000
Vdeb = 0.628013
Fr = 0.404814
P = 184.824088
flowcond = emergent
Vmax = 1.536115
V_technique = 1.592932
q_technique = 0.414154
 */

function macroRugoInstanceEmergentCd2(): MacroRugo {
    const nub: MacroRugo = macroRugoInstanceEmergentCd15();
    nub.prms.Cd0.singleValue = 2;
    nub.prms.Q.singleValue = 0.376808;
    return nub;
}

const macroRugoExtraResultEmergentCd2: { [key: string]: number|MacroRugoFlowType } = {
    ENUM_MacroRugoFlowType: MacroRugoFlowType.EMERGENT,
    Fr: 0.404817,
    PV: 184.825734,
    Q_GuideTech: 0.414154,
    V_GuideTech: 1.592932,
    Vdeb: 0.628018,
    Vmax: 1.536129,
    ZF2: 12.2
};

/*
*** Submerged conditions ***
ks=0.010000
D=0.500000
k=0.700000
Cd0=1.500000
S=0.050000
B=1.000000
h=0.800000
C=0.130000
Q=0.908068  fVal=0.000000
cote_bas = 12.200000
Vdeb = 1.135085
Fr = 0.633645
P = 445.407266
flowcond = immerge
q_technique = 0.940450
 */

function macroRugoInstanceSubmerged(): MacroRugo {
    const nub: MacroRugo = macroRugoInstanceEmergentCd15();
    nub.prms.Y.singleValue = 0.8;
    nub.prms.Q.singleValue = 0.908068;
    return nub;
}

const macroRugoExtraResultSubmerged: { [key: string]: number|MacroRugoFlowType } = {
    ENUM_MacroRugoFlowType: MacroRugoFlowType.SUBMERGED,
    Fr: 0.633645,
    PV: 445.407266,
    Q_GuideTech: 0.940450,
    Vdeb: 1.135085,
    ZF2: 12.2
};

function MacroRugoFactory(sInstance: string): MacroRugo {
    let nub: MacroRugo;
    switch (sInstance) {
        case "EmergentCd15": {
            nub = macroRugoInstanceEmergentCd15();
            break;
        }
        case "EmergentCd2": {
            nub = macroRugoInstanceEmergentCd2();
            break;
        }
        case "Submerged": {
            nub = macroRugoInstanceSubmerged();
            break;
        }
        default: {
            throw new Error("Instance name error");
        }
    }
    for (const p of nub.parameterIterator) {
        p.v = p.singleValue;
    }
    return nub;
}

function testMacroRugo(sInstance: string, varTest: string, valRef: number) {
        it(`Calc(${varTest}) should be ${valRef}`, () => {
            const nub = MacroRugoFactory(sInstance);
            nub.prms.Q.v = nub.Calc("Q").vCalc;
            const p = nub.getParameter(varTest);
            p.v = undefined;
            p.valueMode = ParamValueMode.CALCUL;
            const r = nub.Calc(varTest, 0.1);
            checkResult(r, valRef);
        });
}

function testMacroRugoConfig(sInstance: string, Q0: number, fVal0: number, mrExtraRes: { [key: string]: number }) {
    describe(`Condition ${sInstance}` , () => {
        it(`resolveQ(${Q0}) should be ${fVal0}`, () => {
            const nubit = MacroRugoFactory(sInstance);
            nubit.prms.Q.v = Q0;
            // tslint:disable-next-line:no-string-literal
            expect(nubit["resolveQ"]()).toBeCloseTo(fVal0, 5);
        });
        const nub = MacroRugoFactory(sInstance);
        for (const prm of nub.prms) {
            if ([ParamCalculability.DICHO, ParamCalculability.EQUATION].includes(prm.calculability)) {
                testMacroRugo(sInstance, prm.symbol, prm.v);
            }
        }
        for (const sExtraRes in mrExtraRes) {
            if (mrExtraRes.hasOwnProperty(sExtraRes)) {
                it(`${sExtraRes} should be ${mrExtraRes[sExtraRes]}`, () => {
                    expect(MacroRugoFactory(sInstance).Calc("Q").values[sExtraRes])
                    .toBeCloseTo(mrExtraRes[sExtraRes], 5);
                });
            }
        }
    });
}

describe("Class MacroRugo: ", () => {

    testMacroRugoConfig("EmergentCd15", 0.901710, 0.379539, macroRugoExtraResultEmergentCd15);
    testMacroRugoConfig("EmergentCd2", 0.768677, 0.418029, macroRugoExtraResultEmergentCd2);
    // Le test passe en debug mais pas sous Jasmine !??
    // describe(`Condition Submerged` , () => {
    //     it(`resolveAlpha_t(0.07) should be 0.074752`, () => {
    //         // tslint:disable-next-line:no-string-literal
    //         expect(macroRugoInstanceSubmerged()["resolveAlpha_t"](0.07)).toBeCloseTo(0.074752, 5);
    //     });
    // });
    testMacroRugoConfig("Submerged", 1.202280, 0.145051, macroRugoExtraResultSubmerged);

    describe("JalHyd #85", () => {
        it("CalcSerie Q should return the goog result :)", () => {
            const nubMR = macroRugoInstanceEmergentCd2();
            nubMR.prms.C.singleValue = 0.2;
            nubMR.prms.PBH.singleValue = 0.8;
            nubMR.calculatedParam = nubMR.prms.Q;

            nubMR.prms.Y.setValues(0.7, 1.2, 0.1);
            const aQ: number[] = [ 0.363679, 0.433841, 0.595482, 0.768310, 0.991123, 1.441876 ];

            nubMR.CalcSerie();
            for (let i = 0; i < aQ.length; i++) {
                expect(nubMR.result.resultElements[i].vCalc).toBeCloseTo(aQ[i], 4);
            }
        });
    });
});