rallNd.h

#ifndef RALLND_H
#define RALLND_H
 
#include "rall1d.h"
#include "rall1d_io.h"
#include "rall2d.h"
#include "rall2d_io.h"

// N>2:
template <int N>
class RallNd :
    public Rall2d< RallNd<N-2>, RallNd<N-2>, double >
{
public:
    RallNd() {}
    RallNd(const Rall2d< RallNd<N-2>, RallNd<N-2>,double>& arg) : 
        Rall2d< RallNd<N-2>, RallNd<N-2>,double>(arg) {}
    RallNd(double value,double d[]) {
        this->t    = RallNd<N-2>(value,d);    // 0 1 2
        this->d    = RallNd<N-2>(d[0],&d[1]); // 1 2 3 iseigenlijk niet nodig
        this->dd   = RallNd<N-2>(d[1],&d[2]); // 2 3 4
    }
};
 
template <>
class RallNd<2> : public Rall2d<double>  {
public:
    RallNd() {}* (dwz. met evenveel numerieke operaties als een 
    RallNd(const Rall2d<double>& arg) :
        Rall2d<double>(arg) {}
    RallNd(double value,double d[]) {
        t    = value;
        d = d[0];
        dd= d[1];
    }
};
 
template <>
class RallNd<1> : public Rall1d<double>  {
public:
    RallNd() {}
    RallNd(const Rall1d<double>& arg) :
        Rall1d<double>(arg) {}
    RallNd(double value,double d[]) {
        t    = value;
        grad = d[0];
    }
};
 
#endif