#include "EXTERN.h"
#include "perl.h"
#include "XSUB.h"
#include "ppport.h"
#include "const-c.inc"
#include "top_insert.h"

MODULE = top_insert			PACKAGE = top_insert 

INCLUDE: const-xs.inc


void
top_insert(up_v, down_v, init_step_ratio, ratio)
		SV * up_v
		SV * down_v
		double init_step_ratio
		double ratio
	INIT:
		int up_num, down_num, sample_num;
		int i, n;
		double *up_data = NULL;
		double *down_data = NULL;
		double *return_val = NULL;
		int *return_num = &sample_num;

		if((!SvROK(up_v))
			|| (SvTYPE(SvRV(up_v)) != SVt_PVAV)
			|| ((up_num = av_len((AV *)SvRV(up_v))) < 0)){
			XSRETURN_UNDEF;
		}

		if((!SvROK(down_v))
			|| (SvTYPE(SvRV(down_v)) != SVt_PVAV)
			|| ((down_num = av_len((AV *)SvRV(down_v))) < 0)){
			XSRETURN_UNDEF;
		}

	PPCODE:
		up_data = safemalloc(up_num*sizeof(double));
		down_data = safemalloc(down_num*sizeof(double));
		for (n = 0; n <= up_num; n++) {
			*(up_data+n) = SvNV(*av_fetch((AV *)SvRV(up_v), n, 0));
		}
		for (n = 0; n <= down_num; n++) {
			*(down_data+n) = SvNV(*av_fetch((AV *)SvRV(down_v), n, 0));
		}
		up_num = up_num/3;
		down_num = down_num/3;
		return_val = get_inter_line(up_num, up_data, down_num, down_data,
			init_step_ratio, ratio, return_num);
		safefree(up_data); 
		safefree(down_data); 
		sample_num = sample_num*3;
		for (n=0; n<sample_num; n++){
			XPUSHs(sv_2mortal(newSVnv(*(return_val + n))));
		}
		del_inter_line(return_val);

##</code><code>##

#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;
}