Example2D.h

#ifndef __EXAMPLE2D__
#define __EXAMPLE2D__

#include <MooNMD_Io.h>
#include <Constants.h>
#include <vector>


class Example2D 
{
  protected:
    Example2D();
  
  public:
    Example2D(std::vector <DoubleFunct2D*> exact,
              std::vector <BoundCondFunct2D*> bc,
              std::vector <BoundValueFunct2D*> bd, CoeffFct2D *coeffs);

    ~Example2D();

    /* functions representing the exact solution */
    std::vector <DoubleFunct2D*> exact_solution;
    /* functions representing the boundary conditions */
    std::vector <BoundCondFunct2D*> boundary_conditions;
    /* functions representing the boundary data */
    std::vector <BoundValueFunct2D*> boundary_data;
    /* functions representing the coefficients of the pde */
    CoeffFct2D *problem_coefficients;
    
    //void *example_info();

    const std::vector <DoubleFunct2D*> & get_exact() const 
    { return exact_solution; }

    DoubleFunct2D* get_exact(unsigned int i) const
    { return exact_solution.at(i); }

    BoundCondFunct2D** get_bc()
    { return &boundary_conditions[0]; }

    BoundCondFunct2D* get_bc(unsigned int i) const
    { return boundary_conditions.at(i); }

    BoundValueFunct2D** get_bd()
    { return &boundary_data[0]; }

    BoundValueFunct2D* get_bd(unsigned int i) const
    { return boundary_data.at(i); }

    CoeffFct2D* get_coeffs() const
    { return problem_coefficients; }
};

#endif // __EXAMPLE2D__