diff options
Diffstat (limited to '2.3-1/src/c/signalProcessing/conv/dconva.c')
-rw-r--r-- | 2.3-1/src/c/signalProcessing/conv/dconva.c | 39 |
1 files changed, 39 insertions, 0 deletions
diff --git a/2.3-1/src/c/signalProcessing/conv/dconva.c b/2.3-1/src/c/signalProcessing/conv/dconva.c new file mode 100644 index 00000000..975a5b19 --- /dev/null +++ b/2.3-1/src/c/signalProcessing/conv/dconva.c @@ -0,0 +1,39 @@ + +/* + * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab + * Copyright (C) 2008 - INRIA - Arnaud TORSET + * + * This file must be used under the terms of the CeCILL. + * This source file is licensed as described in the file COPYING, which + * you should have received as part of this distribution. The terms + * are also available at + * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + * + */ + +#include "conv.h" + +void dconva(double *in1, int size1, double *in2,int size2, double *out){ + int i; + doubleComplex *in1Cpx, *in2Cpx, *result; + + in1Cpx=(doubleComplex*)malloc(2*(unsigned int)size1*sizeof(double)); + in2Cpx=(doubleComplex*)malloc(2*(unsigned int)size2*sizeof(double)); + for (i=0;i<size1;i++){ + in1Cpx[i]=DoubleComplex(in1[i],0); + } + for (i=0;i<size2;i++){ + in2Cpx[i]=DoubleComplex(in2[i],0); + } + result=(doubleComplex*)malloc(2*(unsigned int)(size1+size2-1)*sizeof(double)); + + zconva(in1Cpx,size1,in2Cpx,size2,result); + zreala(result,size1+size2-1,out); + + free(result); + free(in2Cpx); + free(in1Cpx); + +} + + |