1 files changed, 134 insertions, 0 deletions
diff --git a/2.3-1/src/c/signalProcessing/conv/testFloatConv.c b/2.3-1/src/c/signalProcessing/conv/testFloatConv.c
new file mode 100644
index 00000000..932f3752
--- /dev/null
+++ b/2.3-1/src/c/signalProcessing/conv/testFloatConv.c
@@ -0,0 +1,134 @@
+/*
+ *  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"
+#include <stdio.h>
+#include <assert.h>
+
+static void sconvaTest(void){
+	int i;
+	float x1[2]={0.5253563085570931434631f,0.5376229803077876567841f};
+	float h1[1]={0.2256303490139544010162f};
+	float conv1[2]={0.1185363272564196462922f,0.1213040606847684693692f};	
+	float out1[2]; 
+
+	float x2[8]={0.5878720157779753208160f,0.4829179299995303153992f,0.2232865034602582454681f,\
+0.8400885667651891708374f,0.1205995907075703144074f,0.2855364168062806129456f,\
+0.8607514644972980022430f,0.8494101651012897491455f};
+	float h2[5]={0.9931209897622466087341f,0.6488562873564660549164f,0.9923190940171480178833f,\
+0.0500419777818024158478f,0.7485506581142544746399f};
+	float conv2[12]={0.5838280381629499249385f,0.8610403862135279640100f,1.1184514745692031301161f,\
+1.4878174016996970063786f,1.3506563468112671166921f,1.568122169583095804057f,\
+1.3689566747706436444787f,2.3202981934690010135114f,1.5098489496125344633271f,\
+1.099698103860223152850f,0.6868222398318966082087f,0.6358265380955081536030f};	
+	float out2[12]; 
+
+	sconva(x1,2,h1,1,out1);
+	for (i=0;i<2;i++){
+		assert( ( fabs(out1[i]-conv1[i]) / fabs(out1[i]) ) <1e-7);
+	}
+
+	sconva(h1,1,x1,2,out1);
+	for (i=0;i<2;i++){
+		assert( ( fabs(out1[i]-conv1[i]) / fabs(out1[i]) ) <1e-7);
+	}
+
+
+	sconva(x2,8,h2,5,out2);
+	for (i=0;i<12;i++){
+		assert( ( fabs(out2[i]-conv2[i]) / fabs(out2[i]) ) <1e-6);
+	}
+
+	sconva(h2,5,x2,8,out2);
+	for (i=0;i<12;i++){
+		assert( ( fabs(out2[i]-conv2[i]) / fabs(out2[i]) ) <1e-6);
+	}
+}
+
+
+static void cconvaTest(void){
+	int i;
+	float rx1[2]={0.6274093082174658775330f,0.7608432597480714321137f};
+	float ix1[2]={0.0485566202551126480103f,0.6723949727602303028107f};
+	float rh1[1]={0.2017172696068882942200f};
+	float ih1[1]={0.3911573919467628002167f};
+	float rconv1[2]={0.1075660116388346271288f,-0.1095370389478208250633f};
+	float iconv1[2]={0.2552104975446664592020f,0.4332431431659207188289f};
+	floatComplex *x1, *h1, out1[2];
+
+	float rx2[4]={0.9262343775480985641480f,0.5667211269028484821320f,0.5711638936772942543030f,\
+0.8160110483877360820770f};
+	float ix2[4]={0.0568927936255931854248f,0.5595936686731874942780f,0.1249340316280722618103f,\
+0.7279222286306321620941f};
+	float rh2[7]={0.6783956284634768962860f,0.3320095278322696685791f,0.0258709923364222049713f,\
+0.5174468224868178367615f,0.3916873238049447536469f,0.2413538414984941482544f,\
+0.5064434898085892200470f};
+	float ih2[7]={0.4236102001741528511047f,0.2893727528862655162811f,0.0887932181358337402344f,\
+0.6212881752289831638336f,0.3454984454438090324402f,0.7064867629669606685638f,\
+0.5211472492665052413940f};
+	float rconv2[10]={0.6042529849650153384033f,0.4384669630868825995762f,0.3796883200627478305300f,\
+0.8076061278620030314812f,0.3526828755641281398781f,0.3863983221928849465243f,\
+0.3314166433494881713884f,0.113095069738258141578f,-0.0931668091195109737690f,\
+0.0339088159370524999670f};
+	float iconv2[10]={0.4309581525680090119401f,0.9066106897516487128641f,0.7602048387189073164194f,\
+1.7159441161321140612017f,1.5157099506355884432907f,1.5938830446486993430710f,\
+2.1768689123487110492761f,1.5794692240126795290678f,1.1131203493567218831828f,\
+0.7939133870152889072003f};
+	floatComplex *x2, *h2, out2[10];
+	
+
+	x1=FloatComplexMatrix(rx1,ix1,2);
+	h1=FloatComplexMatrix(rh1,ih1,1);
+	x2=FloatComplexMatrix(rx2,ix2,4);
+	h2=FloatComplexMatrix(rh2,ih2,7);
+
+
+	cconva(x1,2,h1,1,out1);
+	for (i=0;i<2;i++){
+		assert( ( fabs(creals(out1[i])-rconv1[i]) / fabs(creals(out1[i])) ) <1e-6);
+		assert( ( fabs(cimags(out1[i])-iconv1[i]) / fabs(cimags(out1[i])) ) <1e-6);
+	}
+
+	cconva(h1,1,x1,2,out1);
+	for (i=0;i<2;i++){
+		assert( ( fabs(creals(out1[i])-rconv1[i]) / fabs(creals(out1[i])) ) <1e-6);
+		assert( ( fabs(cimags(out1[i])-iconv1[i]) / fabs(cimags(out1[i])) ) <1e-6);
+	}
+
+	cconva(x2,4,h2,7,out2);
+	for (i=0;i<10;i++){
+		assert( ( fabs(creals(out2[i])-rconv2[i]) / fabs(creals(out2[i])) ) <3e-6);
+		assert( ( fabs(cimags(out2[i])-iconv2[i]) / fabs(cimags(out2[i])) ) <1e-6);
+	}
+
+	cconva(h2,7,x2,4,out2);
+	for (i=0;i<10;i++){
+		assert( ( fabs(creals(out2[i])-rconv2[i]) / fabs(creals(out2[i])) ) <3e-6);
+		assert( ( fabs(cimags(out2[i])-iconv2[i]) / fabs(cimags(out2[i])) ) <1e-6);
+	}
+
+}
+
+
+
+static int convTest(void){
+	printf("\n>>>> Conv Float Tests\n");
+	sconvaTest();
+	cconvaTest();
+	return 0;
+}
+
+int main (void){
+	assert(convTest()==0);
+	return 0;
+}