---

output_to_GNUPLOT.h

Go to the documentation of this file.

#ifndef __OUTPUT_TO_GNUPLOT__
#define __OUTPUT_TO_GNUPLOT__

#include <fstream>
#include <iostream>

#include "function_2D_types.h"

#include "msg_stream.h"


namespace Function_2D{

  template <typename Function>
  void output_to_GNUPLOT(const char*     file_name,
                         const Function& f,
                         const real_type x_begin = -16,
                         const real_type x_end = 16,
                         const real_type x_step = 1,
                         const real_type y_begin = -16,
                         const real_type y_end = 16,
                         const real_type y_step = 1);



  
/*
  ###########################################################
  ###########################################################
  ###########################################################
  ##############                               ##############
  ##############  I M P L E M E N T A T I O N  ##############
  ##############                               ##############
  ###########################################################
  ###########################################################
  ###########################################################
*/

  
  
  template <typename Function>
  void output_to_GNUPLOT(const char*     file_name,
                         const Function& f,
                         const real_type x_begin,
                         const real_type x_end,
                         const real_type x_step,
                         const real_type y_begin,
                         const real_type y_end,
                         const real_type y_step){

    using namespace std;
    
    ofstream out(file_name);

    for(real_type x=x_begin;x<=x_end;x+=x_step){
      for(real_type y=y_begin;y<=y_end;y+=y_step){
        out<<x<<' '<<y<<' '<<f(x,y)<<'\n';
      }
      out<<'\n';
    }
    
    out.flush();
    out.close();
  }
  
}

#endif

---