ot_root_find.h

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       ot_root_find.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:54:00 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
#ifndef ot_root_findH
#define ot_root_findH
 
#include <math.h>
#include <stdio.h>
 
 
 
class OT_ROOT_FIND_1D {
public:
 
    typedef double (*Function)(double,void*);
 
    static unsigned ZeroBracOut (Function fn, double &x1, double &x2, void* pvData = 0);
    static unsigned ZeroBracOut (Function fn, double xStart, double xEnd, double &x1, double &x2, void* pvData);
 
    static unsigned ZeroBracIn (Function fn, const double &x1, const double &x2,
                                unsigned n, double xb1[], double xb2[], unsigned nb, void* pvData = 0);
 
    static double RootBisect (Function fn, const double &x1, const double &x2, const double &xacc, void* pvData = 0);
 
    static double RootFlsPos (Function fn, const double &x1, const double &x2,
                              const double &xacc, void* pvData = 0);
 
    /* RootNewton (fn, df, x1, x2, xacc)
       using the newton-raphson method, find the root of fn in [x1, x2] to an
       accuracy +/-xacc. df gives the 1st derivative of fn.                         */
 
    static double RootNewton (Function fn, Function df,
                              const double &x1, const double &x2, const double &xacc, void* pvData = 0);
 
    /* RootSafe (fn, df, x1, x2, xacc)
       using a combination of newton-raphson and bisection, find the root of fn
       in [x1, x2] to an accuracy +/-xacc. df gives the 1st derivative of fn.   */
 
    static double RootSafe (Function fn, Function df,
                            const double &x1, const double &x2, const double &xacc, void* pvData = 0);
};
#endif