tpl_fonction.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       tpl_fonction.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:53:59 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef _TPLFONCTION_
#define _TPLFONCTION_
 
#include <math.h>
 
template <class A,class B,class C>
class TPL_FONCTION1
{
public:
TPL_FONCTION1(B& bb,A (*fonc)(B&,C)):f(fonc),b(bb)
{
};
 
virtual ~TPL_FONCTION1()
{
};
 
 
 
A integrer_gauss_2(C t1,C t2,int nb_pas=32)
{
double integrale=0.;
for (int i=0;i<nb_pas;i++)
        {
        C ti=t1+(t2-t1)*i/nb_pas;
        C tii=t1+(t2-t1)*(i+1)/nb_pas;
        double t=0.7886751345*tii+0.2113248654*ti;
        integrale=integrale+0.5*(tii-ti)*f(b,t);
        t=0.7886751345*ti+0.2113248654*tii;
        integrale=integrale+0.5*(tii-ti)*f(b,t);
        }
return integrale;
};
 
 
void integrer_jusqua_gauss_2 (C x1, C dx, C valeur_integrale_cible, C * x2)
{
    C xi, xii;
    A integrale,integrale_pas_precedent;
 
    integrale=0;
    xii = x1;
 
    do
        {
        xi = xii;
        xii = xi+dx;
        integrale_pas_precedent=integrale;
        double x=0.7886751345*xii+0.2113248654*xi;
        integrale=integrale+0.5*dx*f(b,x);
        x=0.7886751345*xi+0.2113248654*xii;
        integrale=integrale+0.5*dx*f(b,x);
        }
    while ( fabs (integrale) < fabs (valeur_integrale_cible) );
    *x2=xi+(xii-xi)*(valeur_integrale_cible-integrale_pas_precedent)/(integrale-integrale_pas_precedent);
 
    if ( fabs((*x2-x1)/dx) < 32)
        integrer_jusqua_gauss_2 (x1, (*x2-x1)/40, valeur_integrale_cible, x2);
 
}
 
 
 
private:
A (*f)(B&,C);
B& b;
};
 
 
 
 
#endif