#ifndef M_GRF_PDE_1D_H
#define M_GRF_PDE_1D_H
#include <fstream>
#include <string>
#include <pde_1d.h>

class pde_1d;

class grf_pde_1d {
public:
  grf_pde_1d(pde_1d *pde,std::string file_name);

  virtual ~grf_pde_1d() {};

  virtual void write_1_grf(const char *name,int index,double *p, 
                                       int vsize,double t);
  virtual std::string make_file_name(std::string name,int index,
                                     std::string suffix);

protected:
  int nx_;
  double xmin_;
  double xmax_;
  double dx_;

  std::string file_name_;
};

/****************************** PERIODIC GRF *************************/
class grf_pde_1d_periodic : grf_pde_1d {
public:
  grf_pde_1d_periodic(pde_1d *pde,std::string file_name) : 
               grf_pde_1d(pde,file_name) {};

  virtual void write_1_grf(const char *name,int index,double *p, 
                                    int vsize,double t);

};

#endif