Zum Inhalt springen

Templates aufrufen!


hunter1670a

Empfohlene Beiträge

Moin an allen,

ich schreibe gerade ein Template mit dem ich numerische Operationen durchführe, leider weiß ich nicht wie man ein Template von einem "Testerprogramm" aufruf. Kann mir da jemand helfen, weil ich sonst meine Algorithmen nicht testen kann.

hier ein Auszug aus meinem Template:

typedef (*func)(const T&, f) function;

Template <class T> bool simpson(T a, T b, function f, unsigned long max_iter){

T s,s1,s2,s4,h,q,I,n;

n=1;

h=(b-a)/2;

s=f(a)+f(B);

s2=0;

q=0;

s4=f(a+h);

s1=h*(s+4*s4);

//if Abfrage --> mit isequal ...

do{

I=s1-q;

if (I <= numeric_limits <T>::epsilon())

max_iter=0;

else

{

h=h/2;

n=2*n;

q=S1;

s4=0;

s2+=s4;

x=a+h;

for(i=1;i=n;i++) {

s4+=f(x);

x+=2*h;

}

s1=(s+2*s2+4*s4)*h;

}

max_iter--;

}while(max_iter != 0);

//Ende der If-Abfrage

s=s1/3;

return s;

Wie kann ich die den simpson von eine Testerprogramm aufrufen ?? Problem ist auch wie versteht das Templet das er die Wert aus der Funktion wie f(x)=x^2 nehmen soll ??

Link zu diesem Kommentar
Auf anderen Seiten teilen

Original geschrieben von hunter1670a

Wie kann ich die den simpson von eine Testerprogramm aufrufen ??

Es gibt keine spezielle Syntax für den Aufruf von Template-Funktionen. Ruf sie einfach auf wie eine normale Funktion.

Problem ist auch wie versteht das Templet das er die Wert aus der Funktion wie f(x)=x^2 nehmen soll ??
Die Funktion übergibst Du doch offenbar als Parameter. Kannst Du ein wenig verständlicher formulieren, wo das Problem liegt?
Link zu diesem Kommentar
Auf anderen Seiten teilen

Das Template soll wie alle möglichen Funktionen wie sin(x) oder Polynome den Flächeninhalt mit den verschiedenen Verfahren ausrechnen.

Das ist mein bisheriges Testprogramm:

#include <iostream.h>

#include <stdio.h>

#include "Integral.h"

void main()

{

double a =-1;

double b =5;

unsigned long max = 15;

printf ("%.15f\n",Numerik::trapez(a,b,max));

printf ("%.15f\n",Numerik::simpson(a,b,max));

max = 4;

printf ("%.15f\n",Numerik::romberg(a,b,max));

}

Diese Funktion habe ich jetzt in mein Template hingeschrieben um meine Algorithmen wenigstens zu testen.

double f (double x)

{

double g = ((x*x)+(2*x)-1);

return g;

}

Diese Funktion darf aber gar nicht in dem Template stehen sondern im Testprogramm, und da ist mein Problem wie findet mein Template die Funktion aus dem Test Programm

Link zu diesem Kommentar
Auf anderen Seiten teilen

Ich habe vor langer Zeit auch mal so ein Programm geschrieben. Ist nicht besonders hübsch, aber wen es stört, kann es ja ändern. Heute morgen hatte ich jedenfalls keine Lust dazu. Ich habe meine Funktion, die berechnet werden soll als Klasse implementiert und dann an das Template übergeben.


# include <iostream>

# include <math>


class TIntervall

{

    private:

        double FIntervall[2];

        int FIntervallSwitch;

    public:

        TIntervall ( void );

        TIntervall ( double, double );

        void setIntervall ( double, double );

        double getLowerBorder ( void );

        double getHigherBorder ( void );

};


TIntervall :: TIntervall ( void )

{

    FIntervallSwitch = 0;


    FIntervall[0] = 0.0;

    FIntervall[1] = 0.0;

}


TIntervall :: TIntervall ( double a, double b )

{

    setIntervall ( a, b );

}


void TIntervall :: setIntervall ( double a, double b )

{

    if ( a < b )

        FIntervallSwitch = 0;

    else

        FIntervallSwitch = 1;


    FIntervall[FIntervallSwitch%2] = a;

    FIntervall[FIntervallSwitch%2+1] = b;

}


double TIntervall :: getLowerBorder ( void )

{

    return FIntervall[FIntervallSwitch%2];

}


double TIntervall :: getHigherBorder ( void )

{

    return FIntervall[FIntervallSwitch%2+1];

}


//---------------------------------------------------------------------------


template <class Function> class TTrapez

{

    public:

        double integrate ( TIntervall, int );

};


template <class Function>

double TTrapez <Function> :: integrate ( TIntervall inter, int n )

{

    Function f;

    double ret = f ( inter.getLowerBorder () ) + f ( inter.getHigherBorder () );

    double step = ( inter.getHigherBorder () - inter.getLowerBorder () ) / n;


    for ( int i = 1; i < n; i++ )

        ret += 2 * f ( inter.getLowerBorder () + i * step );


    return ret * ( step / 2 );

}


//---------------------------------------------------------------------------


template <class Function> class TSimpson

{

    public:

        double integrate ( TIntervall, int );

};


template <class Function>

double TSimpson <Function> :: integrate ( TIntervall inter, int n )

{

    Function f;

    double ret = f ( inter.getLowerBorder () ) + f ( inter.getHigherBorder () );

    double step = ( inter.getHigherBorder () - inter.getLowerBorder () ) /

                    ( 2 * n );


    for ( int i = 1; i < ( 2 * n ); i++ )

        ret += ( 2 * ( i % 2 ) + 2 ) * f ( inter.getLowerBorder () + i * step );


    return ret * ( step / 3 );

}


//---------------------------------------------------------------------------


class TSinus

{

    public:

        double operator() ( double );

};


double TSinus :: operator () ( double x )

{

    return std::sin ( x );

}


//---------------------------------------------------------------------------


int main ( void )

{

    TTrapez <TSinus> TrapezIntegral;

    TSimpson <TSinus> SimpsonIntegral;

    TIntervall intervall ( 0, 1 );


    std::cout.precision ( 20 );

    std::cout << "Ergebnisse fuer f(x) = sin(x) mit Trapezregel: " << std::endl;


    for ( int i = 1; i < 11; i++ )

    {

	std::cout << "Fuer n = ";

        std::cout.width ( 2 );

        std::cout << i << ": " << TrapezIntegral.integrate ( intervall, i ) << std::endl;

    }


	std::cout << std::endl << "Ergebnisse fuer f(x) = sin(x) mit Simpson-Regel: " << std::endl;


    for ( int i = 1; i < 11; i++ )

    {

	std::cout << "Fuer n = ";

        std::cout.width ( 2 );

        std::cout << i << ": " << SimpsonIntegral.integrate ( intervall, i ) << std::endl;

    }


    return 0;

}

HTH

Jan

Link zu diesem Kommentar
Auf anderen Seiten teilen

Dein Kommentar

Du kannst jetzt schreiben und Dich später registrieren. Wenn Du ein Konto hast, melde Dich jetzt an, um unter Deinem Benutzernamen zu schreiben.

Gast
Auf dieses Thema antworten...

×   Du hast formatierten Text eingefügt.   Formatierung wiederherstellen

  Nur 75 Emojis sind erlaubt.

×   Dein Link wurde automatisch eingebettet.   Einbetten rückgängig machen und als Link darstellen

×   Dein vorheriger Inhalt wurde wiederhergestellt.   Editor leeren

×   Du kannst Bilder nicht direkt einfügen. Lade Bilder hoch oder lade sie von einer URL.

Fachinformatiker.de, 2024 by SE Internet Services

fidelogo_small.png

Schicke uns eine Nachricht!

Fachinformatiker.de ist die größte IT-Community
rund um Ausbildung, Job, Weiterbildung für IT-Fachkräfte.

Fachinformatiker.de App

Download on the App Store
Get it on Google Play

Kontakt

Hier werben?
Oder sende eine E-Mail an

Social media u. feeds

Jobboard für Fachinformatiker und IT-Fachkräfte

×
×
  • Neu erstellen...