// This file is released into Public Domain.
use crate::common::*;

use gnuplot::*;
use std::vec::Vec;

mod common;

fn example(c: Common)
{
	let w = 61;
	let h = 61;
	let mut z1 = Vec::with_capacity((w * h) as usize);
	for i in 0..h
	{
		for j in 0..w
		{
			let y = 8.0 * (i as f64) / h as f64 - 4.0;
			let x = 8.0 * (j as f64) / w as f64 - 4.0;
			z1.push(x.cos() * y.cos() / ((x * x + y * y).sqrt() + 1.0));
		}
	}

	let mut fg = Figure::new();

	fg.axes3d()
		.set_title("Surface fg3.1", &[])
		.surface(z1.iter(), w, h, Some((-4.0, -4.0, 4.0, 4.0)), &[])
		.set_x_label("X", &[])
		.set_y_label("Y", &[])
		.set_z_label("Z", &[])
		.set_z_range(Fix(-1.0), Fix(1.0))
		.set_z_ticks(Some((Fix(1.0), 1)), &[Mirror(false)], &[])
		.set_view(45.0, 45.0);

	c.show(&mut fg, "example3_1");

	let mut fg = Figure::new();

	fg.axes3d()
		.set_title("Map fg3.2", &[])
		.surface(z1.iter(), w, h, None, &[])
		.set_x_label("X", &[])
		.set_y_label("Y", &[])
		.set_view_map();

	c.show(&mut fg, "example3_2");

	let mut fg = Figure::new();

	fg.axes3d()
		.set_pos_grid(2, 2, 0)
		.set_title("Base fg3.3", &[])
		.show_contours(true, false, Cubic(10), Fix(""), Auto)
		.surface(z1.iter(), w, h, Some((-4.0, -4.0, 4.0, 4.0)), &[])
		.set_view(45.0, 45.0);

	fg.axes3d()
		.set_pos_grid(2, 2, 1)
		.set_title("Surface", &[])
		.show_contours(false, true, Linear, Fix(""), Auto)
		.surface(z1.iter(), w, h, Some((-4.0, -4.0, 4.0, 4.0)), &[])
		.set_view(45.0, 45.0);

	fg.axes3d()
		.set_pos_grid(2, 2, 2)
		.set_title("Both + Fix Levels", &[])
		.show_contours(true, true, Linear, Fix("%f"), Fix(1))
		.surface(z1.iter(), w, h, Some((-4.0, -4.0, 4.0, 4.0)), &[])
		.set_view(45.0, 45.0);

	fg.axes3d()
		.set_pos_grid(2, 2, 3)
		.set_title("Custom Levels", &[])
		.show_contours_custom(true, false, Linear, Fix(""), Some(0f32).iter())
		.surface(z1.iter(), w, h, Some((-4.0, -4.0, 4.0, 4.0)), &[])
		.set_view(45.0, 45.0);

	c.show(&mut fg, "example3_3");
}

fn main()
{
	Common::new().map(|c| example(c));
}