/********************************************************
 *  ██████╗  ██████╗████████╗██╗
 * ██╔════╝ ██╔════╝╚══██╔══╝██║
 * ██║  ███╗██║        ██║   ██║
 * ██║   ██║██║        ██║   ██║
 * ╚██████╔╝╚██████╗   ██║   ███████╗
 *  ╚═════╝  ╚═════╝   ╚═╝   ╚══════╝
 * Geophysical Computational Tools & Library (GCTL)
 *
 * Copyright (c) 2022  Yi Zhang (yizhang-geo@zju.edu.cn)
 *
 * GCTL is distributed under a dual licensing scheme. You can redistribute 
 * it and/or modify it under the terms of the GNU Lesser General Public 
 * License as published by the Free Software Foundation, either version 2 
 * of the License, or (at your option) any later version. You should have 
 * received a copy of the GNU Lesser General Public License along with this 
 * program. If not, see <http://www.gnu.org/licenses/>.
 * 
 * If the terms and conditions of the LGPL v.2. would prevent you from using 
 * the GCTL, please consider the option to obtain a commercial license for a 
 * fee. These licenses are offered by the GCTL's original author. As a rule, 
 * licenses are provided "as-is", unlimited in time for a one time fee. Please 
 * send corresponding requests to: yizhang-geo@zju.edu.cn. Please do not forget 
 * to include some description of your company and the realm of its activities. 
 * Also add information on how to contact you by electronic and paper mail.
 ******************************************************/

#include "FractalTopo.h"

int main(int argc, char* argv[]) try
{
	const char* pro_name = "fractopo";
	const char* pro_info = "1.5 - Fractal topography generater. \
		This program is a toolkit of the GCTL package. The GCTL comes with ABSOLUTE NO WARRANTY. \
		Please see instructions or contact the author for more information.";

	const char* p_str = "Roughness and maximal changing range of the output topography data. \
		Small roughness generates topography with more high frequent changes. Big roughness results \
		in long wavelength terrain undulations. The defaults are 1.2/600.";
	const char* r_str = "Range of the output topography data. The program can output topography data \
		under a rectangular grid or at random positions. The defaults are 0/1000/0/1000/10/10.";
	const char* c_str = "Altitudes at the four corner positions. 'dl' represents 'down left' and \
		'ur' is 'up right'. The defaults are 0/0/0/0.";
	const char* g_str = "Filter parameters. 'cx/cy' are center coordinates of the filter. \
		'sigma' is radius of the shorter axis of the 0.5 filter amplitude. \
		'ang' is angle between the shorter axis of the filter and the x-axis. \
		'lenx/leny' represents length factor the long and short axises of the filter. \
		usually we take lenx = 1 to ensure the radius of 0.5 filter amplitude equals \
		the distance given by sigma in the direction of the shorter axis of the filter. \
		a bigger 'len' value indicates faster changing rate in the corresponding \
		direction. For instance 1/0.5 yields an oval which flattening equals 2.0. \
		'magnify' is filter amplitude at the center. 'rho' is changing rate of the arctan function. \
		This argument is optional if the filer type is set to 'exp'.";
	const char* t_str = "Choose weighting function, input 'pow' for power function (default) and 'atan' for arctg function.";


	double Para[2] = {1.2, 600};
	double Range[6] = {0, 1000, 0, 1000, 10, 10};
	double Topo[4] = {0, 0, 0, 0};
	double Gauss[8] = {1e+30, 1e+30, 1e+30, 1e+30, 1e+30, 1e+30, 1e+30, 1e+30};
	char flType[1024] = "atan";
	bool bilinear = false;

	// option format
	// char* info, char* format, bool manda
	static struct gctl::option_info long_opts_help[] = 
	{
		{p_str, "<roughness>/<range>", false},
		{r_str, "<xmin>/<xmax>/<ymin>/<ymax>/[<dx>/<dy>|<point-num>]", false},
		{c_str, "<dl>/<ul>/<ur>/<dr>", false},
		{g_str, "<cx>/<cy>/<sigma>/<ang>/<lenx>/<leny>/<magnify>[/rho]", false},
		{t_str, "<type>", false},
		{"Instead of a random topography, create a bilinear topography which four corners are define by the -c option.", 0, false},
		{"Show help information.", 0, false},
		{0, 0, 0}
	};

	static struct option long_opts[] = 
	{
		{"parameter", required_argument, NULL, 'p'},
		{"range", required_argument, NULL, 'r'},
		{"corner-value", required_argument, NULL, 'c'},
		{"filter-para", required_argument, NULL, 'f'},
		{"filter-type", required_argument, NULL, 't'},
		{"bilinear", no_argument, NULL, 'b'},
		{"help", no_argument, NULL, 'h'},
		{0, 0, 0, 0}
	};

	int curr;
	while(1)
	{
		int optIndex = 0;

		curr = getopt_long(argc, argv, "hbp:r:c:f:t:", long_opts, &optIndex);

		if (curr == -1) break;

		switch (curr)
		{
			case 'h': //显示帮助信息
				gctl::display_logo(std::clog);
				gctl::getopt_long_help(long_opts, long_opts_help, pro_name, pro_info, std::clog, " > <outfile>");
				return 0;
			case 'b':
				bilinear = true;
				break;
			case 'p':
				if (2 != sscanf(optarg, "%lf/%lf", &Para[0],&Para[1])) //格式化读入参数
				{
					GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
					gctl::getopt_long_option_info('p', long_opts, long_opts_help);
					return 0;
				}
				break;
			case 'r':
				if (6 != sscanf(optarg, "%lf/%lf/%lf/%lf/%lf/%lf", &Range[0], &Range[1], &Range[2], 
					&Range[3], &Range[4], &Range[5]))
				{
					if (5 != sscanf(optarg, "%lf/%lf/%lf/%lf/%lf", &Range[0], &Range[1], &Range[2], 
						&Range[3], &Range[4]))
					{
						GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
						gctl::getopt_long_option_info('r', long_opts, long_opts_help);
						return 0;
					}
					else
					{
						//确定只有5个范围参数 将最后一个range参数设置为1e+30
						Range[5] = 1e+30;

						if (Range[0] > Range[1] || Range[2] > Range[3] || Range[4] <= 0)
						{
							GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
							gctl::getopt_long_option_info('r', long_opts, long_opts_help);
							return 0;
						}
					}
				}
				else
				{
					if (Range[0]>Range[1] || (Range[0]+2*Range[4]) > Range[1] || 
						Range[2]>Range[3] || (Range[2]+2*Range[5]) > Range[3])
					{
						GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
						gctl::getopt_long_option_info('r', long_opts, long_opts_help);
						return 0;
					}
				}
				break;
			case 'c':
				if (4 != sscanf(optarg,"%lf/%lf/%lf/%lf",&Topo[0],&Topo[1],&Topo[2],&Topo[3]))
				{
					GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
					gctl::getopt_long_option_info('c', long_opts, long_opts_help);
					return 0;
				}
				break;
			case 'f':
				if (8 != sscanf(optarg,"%lf/%lf/%lf/%lf/%lf/%lf/%lf/%lf", &Gauss[0], &Gauss[1], 
					&Gauss[2], &Gauss[3], &Gauss[4], &Gauss[5], &Gauss[6], &Gauss[7]))
				{
					if (7 != sscanf(optarg,"%lf/%lf/%lf/%lf/%lf/%lf/%lf", &Gauss[0], &Gauss[1], 
						&Gauss[2], &Gauss[3], &Gauss[4], &Gauss[5], &Gauss[6]))
					{
						GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
						gctl::getopt_long_option_info('f', long_opts, long_opts_help);
						return 0;
					}
				}
				break;
			case 't':
				if (1 != sscanf(optarg,"%s",flType)) //格式化读入参数
				{
					GCTL_ShowWhatError("Bad syntax:", GCTL_ERROR_ERROR, optarg, 0, 0);
					gctl::getopt_long_option_info('t', long_opts, long_opts_help);
					return 0;
				}
				break;
			case '?':
				gctl::getopt_long_help(long_opts, long_opts_help, pro_name, " - Unknown options. Please see the help information below.");
				return 0;
			default:
				abort();
		}
	}

	FracTopo FT;
	FT.run(Range,Topo,Para,Gauss,flType,bilinear);
	return 0;
}
catch(std::exception &e)
{
	GCTL_ShowWhatError(e.what(), GCTL_ERROR_ERROR, 0, 0, 0);
}