dichotomie.ts

import { Debug, Nub, Result } from "base";

export class Dichotomie extends Debug {
    /**  Pas de parcours de l'intervalle pour initialisation dichotomie */
    readonly IDEFINT = 100; 
    /** Nombre d'itérations maximum de la dichotomie */
    readonly IDICMAX = 100;  
   
    
    /**
    * Construction de la classe.
    * @param nub Noeud de calcul contenant la méthode de calcul Equation
    * @param sVarCalc Nom de la variable à calculer
    */
    constructor(private nub: Nub, private sVarCalc: string) {
        super(false);
    }

    /** Valeur inconnue à rechercher */
    get vX() {
        return this.nub.v[this.sVarCalc];
    }

    /** Valeur inconnue à rechercher */
    set vX(vCalc) {
        this.nub.v[this.sVarCalc] = vCalc;
    }
    
    /**
     * Méthode simulant l'opérateur booléen xor
     * @see http://www.howtocreate.co.uk/xor.html
     */
    XOR(a, b) {
        return (a || b) && !(a && b);
    }


    /**
     * Calcul de l'équation analytique.
     * @note Wrapper vers this.nub.Equation pour simplifier le code.
     * On utilise la première variable du tableau des variables pourvant être calculées analytiquement
     * Il faudra s'assurer que cette première variable correspond à la méthode de calcul la plus rapide
     */
    private Calcul() {
        /**
         * @note 
         */
        return this.nub.Equation(this.nub.sVarsEq[0]);
    }


    /**
    * Calcul à l'ouvrage
    * @param sVarCalc Nom de la variable à calculer
    * @param rTarget Valeur cible
    * @param rtol Précision attendue du résultat
    * @param rInit Valeur initiale de l'inconnue à rechercher
    */
    Dichotomie(rTarget: number, rTol: number, rInit: number): Result {
        let res: Result

        let XminInit = 1E-8;
        this.vX = XminInit;
        let v1: number = this.Calcul().vCalc;

        let XmaxInit: number = Math.max(1, rInit) * 100;
        this.vX = XmaxInit;
        let v2 = this.Calcul().vCalc;
        
        let DX = (XmaxInit - XminInit) / this.IDEFINT;