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