/*
Copyright (c) 2013 Daniel Stahlke

Permission is hereby granted, free of charge, to any person obtaining a copy
of this software and associated documentation files (the "Software"), to deal
in the Software without restriction, including without limitation the rights
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
copies of the Software, and to permit persons to whom the Software is
furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in
all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
THE SOFTWARE.
*/

#include <fstream>
#include <vector>
#include <math.h>

#include "gnuplot-iostream.h"

// Yes, I'm including a *.cc file.  It contains main().
#include "examples-framework.cc"

// http://stackoverflow.com/a/1658429
#ifdef _WIN32
	#include <windows.h>
	inline void mysleep(unsigned millis) {
		::Sleep(millis);
	}
#else
	#include <unistd.h>
	inline void mysleep(unsigned millis) {
		::usleep(millis * 1000);
	}
#endif

void demo_basic() {
	// -persist option makes the window not disappear when your program exits
	Gnuplot gp("gnuplot -persist");
	// For debugging or manual editing of commands:
	//Gnuplot gp(fopen("plot.gnu"));
	// or
	//Gnuplot gp("tee plot.gnu | gnuplot -persist");

	// NOTE: we can use map here because the X values are intended to be
	// sorted.  If this was not the case, vector<pair<double,double>> could be
	// used instead.

	std::map<double, double> xy_pts_A;
	for(double x=-2; x<2; x+=0.01) {
		double y = x*x*x;
		xy_pts_A[x] = y;
	}

	std::map<double, double> xy_pts_B;
	for(double alpha=0; alpha<1; alpha+=1.0/24.0) {
		double theta = alpha*2.0*3.14159;
		xy_pts_B[cos(theta)] = sin(theta);
	}

	gp << "set xrange [-2:2]\nset yrange [-2:2]\n";
	gp << "plot '-' with lines title 'cubic', '-' with points title 'circle'\n";
	gp.send(xy_pts_A).send(xy_pts_B);
}

void demo_array() {
	// -persist option makes the window not disappear when your program exits
	Gnuplot gp("gnuplot -persist");

	double arr[] = { 1, 3, 2 };

	gp << "plot '-' with lines\n";
	gp.send(arr);
}

namespace xyz {
template <typename T>
struct MyTriple {
	MyTriple(
		T _x,
		T _y,
		T _z
	) : x(_x), y(_y), z(_z) { }

	T x, y, z;
};
}
using xyz::MyTriple;

namespace gnuplotio {
	template<typename T>
	struct BinfmtSender<MyTriple<T> > {
		static void send(std::ostream &stream) {
			BinfmtSender<T>::send(stream);
			BinfmtSender<T>::send(stream);
			BinfmtSender<T>::send(stream);
		}
	};

	template <typename T>
	struct BinarySender<MyTriple<T> > {
		static void send(std::ostream &stream, const MyTriple<T> &v) {
			BinarySender<T>::send(stream, v.x);
			BinarySender<T>::send(stream, v.y);
			BinarySender<T>::send(stream, v.z);
		}
	};

	template<typename T>
	struct TextSender<MyTriple<T> > {
		static void send(std::ostream &stream, const MyTriple<T> &v) {
			TextSender<T>::send(stream, v.x);
			stream << " ";
			TextSender<T>::send(stream, v.y);
			stream << " ";
			TextSender<T>::send(stream, v.z);
		}
	};
}

void demo_tuple() {
	// -persist option makes the window not disappear when your program exits
	Gnuplot gp("gnuplot -persist");

	std::vector<MyTriple<double> > pts;
	for(double alpha=0; alpha<1; alpha+=1.0/120.0) {
		double theta = alpha*2.0*3.14159;
		double x = (2+cos(3*theta))*cos(2*theta);
		double y = (2+cos(3*theta))*sin(2*theta);
		double z = sin(3*theta);
		pts.push_back(MyTriple<double> (x, y, z));
	}

	gp << "splot '-' binary" << gp.binfmt(pts, "record") << "with lines notitle\n";
	gp.sendBinary(pts);
}

void demo_tmpfile() {
	// -persist option makes the window not disappear when your program exits
	Gnuplot gp("gnuplot -persist");

	// NOTE: we can use map here because the X values are intended to be
	// sorted.  If this was not the case, vector<pair<double,double>> could be
	// used instead.

	std::map<double, double> xy_pts_A;
	for(double x=-2; x<2; x+=0.01) {
		double y = x*x*x;
		xy_pts_A[x] = y;
	}

	std::map<double, double> xy_pts_B;
	for(double alpha=0; alpha<1; alpha+=1.0/24.0) {
		double theta = alpha*2.0*3.14159;
		xy_pts_B[cos(theta)] = sin(theta);
	}

	gp << "set xrange [-2:2]\nset yrange [-2:2]\n";
	// Data will be sent via a temporary file.  These are erased when you call
	// gp.clearTmpfiles() or when gp goes out of scope.  If you pass a filename
	// (i.e. "gp.file(pts, 'mydata.dat')"), then the named file will be created
	// and won't be deleted.
	gp << "plot" << gp.file(xy_pts_A) << "with lines title 'cubic',"
		<< gp.file(xy_pts_B) << "with points title 'circle'\n";
}

void demo_png() {
	Gnuplot gp;

	gp << "set terminal png\n";

	std::vector<double> y_pts;
	for(int i=0; i<1000; i++) {
		double y = (i/500.0-1) * (i/500.0-1);
		y_pts.push_back(y);
	}

	std::cout << "Creating my_graph_1.png" << std::endl;
	gp << "set output 'my_graph_1.png'\n";
	gp << "plot '-' with lines, sin(x/200) with lines\n";
	gp.send(y_pts);

	// NOTE: we can use map here because the X values are intended to be
	// sorted.  If this was not the case, vector<pair<double,double>> could be
	// used instead.

	std::map<double, double> xy_pts_A;
	for(double x=-2; x<2; x+=0.01) {
		double y = x*x*x;
		xy_pts_A[x] = y;
	}

	std::map<double, double> xy_pts_B;
	for(double alpha=0; alpha<1; alpha+=1.0/24.0) {
		double theta = alpha*2.0*3.14159;
		xy_pts_B[cos(theta)] = sin(theta);
	}

	std::cout << "Creating my_graph_2.png" << std::endl;
	gp << "set output 'my_graph_2.png'\n";
	gp << "set xrange [-2:2]\nset yrange [-2:2]\n";
	gp << "plot '-' with lines title 'cubic', '-' with points title 'circle'\n";
	gp.send(xy_pts_A).send(xy_pts_B);
}

void demo_vectors() {
	// -persist option makes the window not disappear when your program exits
	Gnuplot gp("gnuplot -persist");

	std::vector<std::vector<double> > vecs(4);
	for(double alpha=0; alpha<1; alpha+=1.0/24.0) {
		double theta = alpha*2.0*3.14159;
		vecs[0].push_back( cos(theta));
		vecs[1].push_back( sin(theta));
		vecs[2].push_back(-cos(theta)*0.1);
		vecs[3].push_back(-sin(theta)*0.1);
	}

	gp << "set xrange [-2:2]\nset yrange [-2:2]\n";
	gp << "plot '-' with vectors title 'circle'\n";
	gp.send(vecs);
}

std::vector<std::vector<double> > get_trefoil() {
	std::vector<std::vector<double> > vecs(3);
	for(double alpha=0; alpha<1; alpha+=1.0/120.0) {
		double theta = alpha*2.0*3.14159;
		vecs[0].push_back((2+cos(3*theta))*cos(2*theta));
		vecs[1].push_back((2+cos(3*theta))*sin(2*theta));
		vecs[2].push_back(sin(3*theta));
	}
	return vecs;
}

void demo_inline_text() {
	std::cout << "Creating inline_text.gnu" << std::endl;
	// This file handle will be closed automatically when gp goes out of scope.
	Gnuplot gp(fopen("inline_text.gnu", "w"));

	std::vector<std::vector<double> > vecs = get_trefoil();

	gp << "splot '-' with lines notitle\n";
	gp.send(vecs);
}

void demo_inline_binary() {
	std::cout << "Creating inline_binary.gnu" << std::endl;
	// This file handle will be closed automatically when gp goes out of scope.
	Gnuplot gp(fopen("inline_binary.gnu", "wb"));

	std::vector<std::vector<double> > vecs = get_trefoil();

	gp << "splot '-' binary" << gp.binfmt(vecs, "record") << "with lines notitle\n";
	gp.sendBinary(vecs);
}

void demo_external_text() {
	std::cout << "Creating external_text.gnu" << std::endl;
	// This file handle will be closed automatically when gp goes out of scope.
	Gnuplot gp(fopen("external_text.gnu", "w"));

	std::vector<std::vector<double> > vecs = get_trefoil();

	std::cout << "Creating external_text.dat" << std::endl;
	gp << "splot" << gp.file(vecs, "external_text.dat") << "with lines notitle\n";
}

void demo_external_binary() {
	std::cout << "Creating external_binary.gnu" << std::endl;
	// This file handle will be closed automatically when gp goes out of scope.
	Gnuplot gp(fopen("external_binary.gnu", "w"));

	std::vector<std::vector<double> > vecs = get_trefoil();

	std::cout << "Creating external_binary.dat" << std::endl;
	gp << "splot" << gp.binaryFile(vecs, "external_binary.dat", "record") << "with lines notitle\n";
}

void demo_animation() {
	Gnuplot gp;

	std::cout << "Press Ctrl-C to quit (closing gnuplot window doesn't quit)." << std::endl;

	gp << "set yrange [-1:1]\n";

	const int N = 1000;
	std::vector<double> pts(N);

	double theta = 0;
	while(1) {
		for(int i=0; i<N; i++) {
			double alpha = (double(i)/N-0.5) * 10;
			pts[i] = sin(alpha*8.0 + theta) * exp(-alpha*alpha/2.0);
		}

		gp << "plot '-' binary" << gp.binfmt(pts) << "with lines notitle\n";
		gp.sendBinary(pts);
		gp.flush();

		theta += 0.2;
		mysleep(100);
	}
}

void register_demos() {
	register_demo("basic",                  demo_basic);
	register_demo("array",                  demo_array);
	register_demo("tuple",                  demo_tuple);
	register_demo("tmpfile",                demo_tmpfile);
	register_demo("png",                    demo_png);
	register_demo("vectors",                demo_vectors);
	register_demo("script_inline_text",     demo_inline_text);
	register_demo("script_inline_binary",   demo_inline_binary);
	register_demo("script_external_text",   demo_external_text);
	register_demo("script_external_binary", demo_external_binary);
	register_demo("animation",              demo_animation);
}