#include <vector>
#include "top_insert.h"
#include "geo.h"
#include "standard_contour.h"

double * get_inter_line(int up_num, double * up_data, 
		int down_num, double *down_data,
		double init_step_ratio,
		double ratio,
		int * return_num
		)
{
	std::vector<point<double> > up_line;
	std::vector<point<double> > down_line;
	for (int i=0; i<up_num; ++i){
		int index = i*3;
		point<double> p(*(up_data+index), 
				*(up_data+index+1),
				*(up_data+index+2)
			       );
		up_line.push_back(p);
	}
	for (int i=0; i<up_num; ++i){
		int index = i*3;
		point<double> p(*(down_data+index), 
				*(down_data+index+1),
				*(down_data+index+2)
			       );
		down_line.push_back(p);
	}
	boundary_point bound(up_line, down_line);
	bound.cacu();
	std::vector<point<double> > val = bound.gen_inter_line(ratio);
	int size = (int)val.size();
	int located = size*3;
	double *r_val = new double[located];
	for (int i = 0; i<size; i++){
		int index = i*3;
		*(r_val+index+0) = val[i].coord[0];
		*(r_val+index+1) = val[i].coord[1];
		*(r_val+index+2) = val[i].coord[2];
	}
	*return_num = size;
	return r_val;
}

void del_inter_line(double * r_val)
{
	delete [] r_val;
}