#include"mex.h"

void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[]){
	int i,j,ln,sn;
	double *lv,*sv,*v1,*v2;
	if(nlhs>2)
		mexErrMsgTxt("Two return values");
	if(nrhs!=2)
		mexErrMsgTxt("Require 2 inputs");
	for(i=0;i<2;++i){
		if(!mxIsDouble(prhs[i]))
			mexErrMsgTxt("Only for double value.");
		if(mxIsComplex(prhs[i]))
			mexErrMsgTxt("Not for complex value.");
		if(mxGetM(prhs[i])!=1)
			mexErrMsgTxt("Only for row vector.");
	}
	ln=mxGetN(prhs[0]);
	sn=mxGetN(prhs[1]);
	if(ln<sn*2-1)
		mexErrMsgTxt("Empty result.");
	plhs[0]=mxCreateDoubleMatrix(1,(ln-2*sn+3)/2,mxREAL);
	plhs[1]=mxCreateDoubleMatrix(1,(ln-2*sn+2)/2,mxREAL);
	lv=mxGetPr(prhs[0]);
	sv=mxGetPr(prhs[1]);
	v1=mxGetPr(plhs[0]);
	v2=mxGetPr(plhs[1]);
	for(i=0;i<(ln-2*sn+3);i+=2){
		for(j=0;j<sn;++j)
			*(v1+i/2)+=(*(lv+i+2*j))*(*(sv+j));
	}
	for(i=0;i<(ln-2*sn+2);i+=2){
		for(j=0;j<sn;++j)
			*(v2+i/2)+=(*(lv+i+2*j+1))*(*(sv+j));
	}
};