hauteur.ts

import { IParamsEquation, Debug } from "../base"
import { acSection, cParamsCanal } from "./section_type";
import { acNewton } from "./newton";
import { cLog } from "./log";

/**
 * Calcul de la hauteur critique
 */
export class cHautCritique extends acNewton {
        private Sn: acSection;
        private oP: cParamsCanal;
        /**
         * Constructeur de la classe
         * @param oSn Section sur laquelle on fait le calcul
         * @param oP Paramètres supplémentaires (Débit, précision...)
         */
        constructor(Sn: acSection, oP: cParamsCanal) {
                super(oP);
                this.Sn = Sn;
                this.oP = oP;

        }
        /**
         * Calcul de la fonction dont on cherche le zéro
         * @param rX Variable dont dépend la fonction
         */
        CalcFn(rX) {
                // Calcul de la fonction
                if (this.Sn.Calc('S', rX) != 0) {
                        var Fn = (Math.pow(this.oP.v.Q, 2) * this.Sn.Calc('B', rX) / Math.pow(this.Sn.Calc('S', rX), 3) / cParamsCanal.G - 1);
                }
                else {
                        Fn = Infinity;
                }
                return Fn;
        }
        /**
         * Calcul analytique de la dérivée de la fonction dont on cherche le zéro
         * @param rX Variable dont dépend la fonction
         */
        CalcDer(rX) {
                if (this.Sn.Calc('S') != 0) {
                        // L'initialisation à partir de rX a été faite lors de l'appel à CalcFn
                        var Der = (this.Sn.Calc('dB') * this.Sn.Calc('S') - 3 * this.Sn.Calc('B') * this.Sn.Calc('B'));
                        Der = Math.pow(this.oP.v.Q, 2) / cParamsCanal.G * Der / Math.pow(this.Sn.Calc('S'), 4);
                }
                else {
                        Der = Infinity;
                }
                //spip_log('cHautCritique:CalcDer('.rX.')='.rDer,'hydraulic.'._LOG_DEBUG);
                return Der;
        }
}
/**
 * Calcul de la hauteur normale
 */
export class cHautNormale extends acNewton {
        private Sn: acSection;
        private Q: number;
        private Ks: number;
        private If: number;
        /**
         * Constructeur de la classe
         * @param oSn Section sur laquelle on fait le calcul
         * @param oP Paramètres supplémentaires (Débit, précision...)
         */
        constructor(Sn: acSection, oP: cParamsCanal) {
                super(oP);
                this.Sn = Sn;
                this.Q = oP.v.Q;