hauteur.ts

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

/**
 * Calcul de la hauteur critique
 */
export class cHautCritique extends acNewton {
        /**
         * Section sur laquuelle porte le calcul
         */
        private Sn: acSection;

        /**
         * Constructeur de la classe
         * @param Sn Section sur laquelle on fait le calcul
         */
        constructor(Sn: acSection, maxIter: number, dbg: boolean = false) {
                super(Sn.prms, maxIter, dbg);
                this.Sn = Sn.clone();
        }

        /**
         * Calcul de la fonction dont on cherche le zéro
         * @param rX Variable dont dépend la fonction
         */
        CalcFn(rX: number) {
                // Calcul de la fonction
                if (this.Sn.Calc("S", rX) != 0)
                        return (Math.pow(this.Sn.prms.Q.v, 2) * this.Sn.Calc("B", rX) / Math.pow(this.Sn.Calc("S", rX), 3) / cParamsCanal.G - 1);

                return Infinity;
        }

        /**
         * 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: number) {
                let S = this.Sn.Calc("S");
                if (S != 0) {
                        let B = this.Sn.Calc("B");
                        // L'initialisation à partir de rX a été faite lors de l'appel à CalcFn
                        let Der = (this.Sn.Calc("dB") * S - 3 * B * B);
                        return Math.pow(this.Sn.prms.Q.v, 2) / cParamsCanal.G * Der / Math.pow(S, 4);
                }

                return Infinity;
        }
}

/**
 * Calcul de la hauteur normale
 */
export class cHautNormale extends acNewton {
        /**
       * Section sur laquuelle porte le calcul
       */
        private Sn: acSection;

        /**
         * Débit connu
         */
        private Q: number;

        /**
         * Coefficient de Strickler
         */
        private Ks: number;