weno.h

Go to the documentation of this file.

#ifndef WENO_H
#define WENO_H

/*
    This header provides and API for any order WENO reconstruction. The
    coefficients are copmuted using Ian's weno-coefficients.ipynb, which
    autogenerates the code for wenoUpwinds.h and wenoUpwinds.cc.
*/


#include "simData.h"
#include <stdexcept>


class WenoBase
{
  public:

    Data * data; 

    int order;    
    int shift;    


    WenoBase(Data * data, int order);


    virtual ~WenoBase() { }     


    virtual double upwindX(double * arr, int var, int i, int j, int k) = 0;


    virtual double upwindY(double * arr, int var, int i, int j, int k) = 0;


    virtual double upwindZ(double * arr, int var, int i, int j, int k) = 0;


    virtual double downwindX(double * arr, int var, int i, int j, int k) = 0;


    virtual double downwindY(double * arr, int var, int i, int j, int k) = 0;


    virtual double downwindZ(double * arr, int var, int i, int j, int k) = 0;


    virtual void reconstructUpwind(double * arr, double * recon, int nvars, int dir);


    virtual void reconstructDownwind(double * arr, double * recon, int nvars, int dir);


    virtual void checkSufficientGhostZones();

};


class Weno3 : public WenoBase
{
  public:
    Weno3(Data * data) : WenoBase(data, 3) { }

    virtual ~Weno3() { }    


    virtual double upwindX(double * arr, int var, int i, int j, int k);


    virtual double upwindY(double * arr, int var, int i, int j, int k);


    virtual double upwindZ(double * arr, int var, int i, int j, int k);


    virtual double downwindX(double * arr, int var, int i, int j, int k);


    virtual double downwindY(double * arr, int var, int i, int j, int k);


    virtual double downwindZ(double * arr, int var, int i, int j, int k);

};


class Weno5 : public WenoBase
{
  public:
    Weno5(Data * data) : WenoBase(data, 5) { }

    virtual ~Weno5() { }    


    virtual double upwindX(double * arr, int var, int i, int j, int k);


    virtual double upwindY(double * arr, int var, int i, int j, int k);


    virtual double upwindZ(double * arr, int var, int i, int j, int k);


    virtual double downwindX(double * arr, int var, int i, int j, int k);


    virtual double downwindY(double * arr, int var, int i, int j, int k);


    virtual double downwindZ(double * arr, int var, int i, int j, int k);

};

class Weno7 : public WenoBase
{
  public:
    Weno7(Data * data) : WenoBase(data, 7) { }

    virtual ~Weno7() { }    


    virtual double upwindX(double * arr, int var, int i, int j, int k);


    virtual double upwindY(double * arr, int var, int i, int j, int k);


    virtual double upwindZ(double * arr, int var, int i, int j, int k);


    virtual double downwindX(double * arr, int var, int i, int j, int k);


    virtual double downwindY(double * arr, int var, int i, int j, int k);


    virtual double downwindZ(double * arr, int var, int i, int j, int k);

};

class Weno9 : public WenoBase
{
  public:
    Weno9(Data * data) : WenoBase(data, 9) { }

    virtual ~Weno9() { }    


    virtual double upwindX(double * arr, int var, int i, int j, int k);


    virtual double upwindY(double * arr, int var, int i, int j, int k);


    virtual double upwindZ(double * arr, int var, int i, int j, int k);


    virtual double downwindX(double * arr, int var, int i, int j, int k);


    virtual double downwindY(double * arr, int var, int i, int j, int k);


    virtual double downwindZ(double * arr, int var, int i, int j, int k);

};

class Weno11 : public WenoBase
{
  public:
    Weno11(Data * data) : WenoBase(data, 11) { }

    virtual ~Weno11() { }    


    virtual double upwindX(double * arr, int var, int i, int j, int k);


    virtual double upwindY(double * arr, int var, int i, int j, int k);


    virtual double upwindZ(double * arr, int var, int i, int j, int k);


    virtual double downwindX(double * arr, int var, int i, int j, int k);


    virtual double downwindY(double * arr, int var, int i, int j, int k);


    virtual double downwindZ(double * arr, int var, int i, int j, int k);

};

#endif