RungeKuttaAdaptive.h

gehe zur Dokumentation dieser Datei

/*
Autor: $Author: kunkel $ State: $State: Exp $
Datum: $Date: 2005/05/30 12:35:25 $
Version: $Revision: 1.1 $
*/

#ifndef c_H
#define RUNGEKUTTAADAPTIVE_H

#include <iterator>
#include <list>

#include "Integrator.h"
#include "Object.h"

// numerical recipes Makros
#define SAFETY 0.9
#define PGROW -0.2
#define PSHRNK -0.25
#define ERRCON 1.89e-4
#define NR_END 1
#define FREE_ARG char*

class RungeKuttaAdaptive : public Integrator {
private:

    double* m;

    double* radius;

    double t;

    int n;

    double eps;

    double nextsteplength;

    double actualsteplength;

    int objnumber;

    list<Collisionspair> collisions;

    double distance(const int i, const int j, const double x[]);
    double acceleration(const int i, const int komponente,
                        const double x[] );
    void derivs(double x, double y[], double dydx[]);
    void rkck(double y[], double dydx[], int n, double x, double h,
              double yout[], double yerr[]);
    bool rkqs(double y[], double dydx[], int n, double *x,
              double htry, double eps,double yscal[], double *hdid,
              double *hnext);
    double *vec(long nl, long nh);
    void free_vec(double *v, long nl, long nh);

protected:

    virtual void saveAttributes();
    virtual void loadAttributes();
    virtual void deleteAdditionalInformation(void * additionalInformation);

public:

    virtual double integrate(double time);

    virtual string getName() {
        return "RungeKuttaAdaptive";
    }

    virtual double getTime() {
        return t;
    }

    virtual void save(ostream & FILE);
    virtual void load(istream & FILE);
    virtual void createAttributeWindow();

    virtual list<Collisionspair>* getCollisionslist();

    RungeKuttaAdaptive();
    virtual ~RungeKuttaAdaptive() {}
    ;
};


#endif

---