// Code starts here-----------
#include <Origin.h>
#include <fft_utils.h>


void deconvolute(string strSignal1, string strSignal2, string strResult)
{
	// Declare datasets and check validity
	Dataset ds1(strSignal1);
	Dataset ds2(strSignal2);
	Dataset ds3(strResult);
	if(!ds1 || !ds2 || !ds3)
	{
		printf("Invalid dataset(s)!\n");
		return;
	}
	
	// Copy datsets to vectors to pass to fft_deconvolute function
	vector vec1, vec2, vec3;
	vec1 = ds1;
	vec2 = ds2;

	int iRet = fft_deconvolute(vec1, vec2, vec3);
	if(0 != iRet)
	{
		printf ("FFT Deconvolution call returned error. Error code: %d\n", iRet);
		return;
	}
	
	// Copy result vector to result dataset
	ds3 = vec3;
}
// Code ends here-----------