Original Version

1 files changed, 139 insertions, 0 deletions

diff --git a/2.3-1/src/c/auxiliaryFunctions/conj/testConj.c b/2.3-1/src/c/auxiliaryFunctions/conj/testConj.c
new file mode 100644
index 00000000..bdb443fc
--- /dev/null
+++ b/2.3-1/src/c/auxiliaryFunctions/conj/testConj.c
@@ -0,0 +1,139 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ *  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 "conj.h"
+#include <assert.h>
+#include <stdio.h>
+#include <math.h>
+#define SIZE 100
+#define FMATRIX  {0.0373321f,0.4293466f,0.3157233f,0.3682477f,0.1458774f,0.6768379f,0.5261979f,0.4003626f,\
+0.0029108f,0.3068181f,0.7902694f,0.9577950f,0.6689271f,0.2929616f,0.8223899f,0.0179846f,\
+0.8710701f,0.3181024f,0.5724473f,0.5738658f,0.3492018f,0.7067298f,0.3599201f,0.4052311f,\
+0.6140811f,0.6524047f,0.0295171f,0.5668896f,0.7114017f,0.4588491f,0.4720517f,0.0623731f,\
+0.0854401f,0.0134564f,0.3543002f,0.6719395f,0.1360619f,0.2119744f,0.4015942f,0.4036219f,\
+0.5628382f,0.5531093f,0.0768984f,0.5360758f,0.0900212f,0.6218026f,0.0001215f,0.7219727f,\
+0.6661293f,0.6799288f,0.9514125f,0.1206901f,0.7318910f,0.0831662f,0.8021909f,0.0290492f,\
+0.6502991f,0.7427882f,0.2176611f,0.9788486f,0.8925237f,0.3511275f,0.4745619f,0.7699462f,\
+0.0734342f,0.5894619f,0.2378993f,0.4076361f,0.4304818f,0.0370717f,0.6404148f,0.4130743f,\
+0.8396902f,0.0002195f,0.5085374f,0.4754997f,0.9625065f,0.9798011f,0.0785698f,0.9150032f,\
+0.6172610f,0.9964195f,0.0468599f,0.9245926f,0.9689865f,0.4916290f,0.2977053f,0.0603055f,\
+0.6313472f,0.0241030f,0.1527438f,0.9481178f,0.2744266f,0.4794727f,0.4855768f,0.6764004f,\
+0.8321250f,0.0125876f,0.5453780f,0.8426717f}
+#define DMATRIX {0.70964269945397973061,0.6610225578770041465759,0.7252350724302232265472,\
+0.510143497958779335022,0.2824866441078484058380,0.2214915128424763679504,\
+0.5767575385980308055878,0.7048253491520881652832,0.5607945644296705722809,\
+0.3628880484029650688171,0.3088326244615018367767,0.9813414867967367172241,\
+0.7428482254035770893097,0.1616391474381089210510,0.3652521963231265544891,\
+0.3835070468485355377197,0.1169181005097925662994,0.8911246126517653465271,\
+0.0634272000752389431000,0.6263941843062639236450,0.6807702076621353626251,\
+0.2588596837595105171204,0.4394804346375167369843,0.8598818778991699218750,\
+0.5482823201455175876617,0.8135995296761393547058,0.4898237003944814205170,\
+0.0240026097744703292847,0.7369195525534451007843,0.2411156157031655311584,\
+0.1529693226329982280731,0.2648359201848506927490,0.4298193217255175113678,\
+0.7673939457163214683533,0.8753260499797761440277,0.3796988371759653091431,\
+0.3062356826849281787872,0.3880051793530583381653,0.1047293278388679027557,\
+0.7603831812739372253418,0.0340930395759642124176,0.1423966242000460624695,\
+0.5554559468291699886322,0.8031897451728582382202,0.5183992316015064716339,\
+0.1518561029806733131409,0.5884730662219226360321,0.2545093484222888946533,\
+0.9992679939605295658112,0.6398976957425475120544,0.5506716123782098293304,\
+0.4607007671147584915161,0.5933007937856018543243,0.6538193570449948310852,\
+0.4168340521864593029022,0.9910155385732650756836,0.3720780410803854465485,\
+0.7061824081465601921082,0.5780865414999425411224,0.602319641038775444031,\
+0.5715096746571362018585,0.0549629041925072669983,0.1205854485742747784,\
+0.0143620483577251434326,0.0257951230742037296295,0.8411248764023184776306,\
+0.4381882525049149990082,0.4643401596695184707642,0.4197426405735313892365,\
+0.8023654492571949958801,0.5228588166646659374237,0.5095121040940284729004,\
+0.4531980180181562900543,0.7817818326875567436218,0.776557037141174077988,\
+0.0466059204190969467163,0.5113326688297092914581,0.7502101892605423927307,\
+0.6883628661744296550751,0.0167756117880344390869,0.9345410899259150028229,\
+0.5609863763675093650818,0.2205847105942666530609,0.7648540753871202468872,\
+0.4816976976580917835236,0.1558785634115338325500,0.2476022052578628063202,\
+0.5216529071331024169922,0.0805230387486517429352,0.4409417239949107170105,\
+0.8911933614872395992279,0.9780590813606977462769,0.4588093762286007404327,\
+0.5510440031066536903381,0.0960960905067622661591,0.4196785055100917816162,\
+0.7502556503750383853912,0.7218149593099951744080,0.2895541018806397914886,\
+0.3907764498144388198853}
+
+
+static void cconjsTest ( void )
+{
+    float realin = 0.126244857f;
+    float imagin = -0.47864959f;
+    floatComplex in = FloatComplex ( realin , imagin  );
+
+    assert ( fabs ( creals ( in  ) - realin ) / fabs( creals (in)) < 1e-06 );
+    assert ( fabs ( cimags ( in  ) - imagin ) / fabs( cimags (in)) < 1e-06 );
+}
+
+
+static void cconjaTest ( void )
+{
+    int i = 0 ;
+    float tin [] = FMATRIX ;
+    floatComplex* in = FloatComplexMatrix ( tin , tin , SIZE );
+    floatComplex out[100] ;
+
+    cconja ( in , SIZE , out ) ;
+
+    for ( i = 0 ; i < SIZE ; i++ )
+    {
+    assert ( fabs ( creals ( in[i] ) - tin[i] ) / fabs( creals (in[i])) < 1e-06 );
+    assert ( fabs ( cimags ( in[i] ) - tin[i] ) / fabs( cimags (in[i])) < 1e-06 );
+    }
+
+}
+
+static void zconjsTest ( void )
+{
+    double realin = 0.126244857788445;
+    double imagin = -0.47864959451267;
+    doubleComplex in = DoubleComplex ( realin , imagin  );
+
+    assert ( fabs ( zreals ( in  ) - realin ) / fabs( zreals (in)) < 1e-16 );
+    assert ( fabs ( zimags ( in  ) - imagin ) / fabs( zimags (in)) < 1e-16 );
+}
+
+static void zconjaTest ( void )
+{
+    int i = 0 ;
+    double tin [] = DMATRIX ;
+    doubleComplex* in = DoubleComplexMatrix ( tin , tin , SIZE );
+    doubleComplex out[100] ;
+
+    zconja ( in , SIZE , out ) ;
+
+    for ( i = 0 ; i < SIZE ; i++ )
+    {
+    assert ( fabs ( zreals ( in[i] ) - tin[i] ) / fabs( zreals (in[i])) < 1e-16 );
+    assert ( fabs ( zimags ( in[i] ) - tin[i] ) / fabs( zimags (in[i])) < 1e-16 );
+    }
+
+}
+
+static int testConj   (void) {
+
+  printf("\n>>>> Float a Tests\n");
+  cconjsTest();
+  cconjaTest();
+printf("\t>>>> Double a Tests\n");
+  zconjsTest();
+  zconjaTest();
+
+
+  return 0;
+}
+
+
+int main(void) {
+  assert(testConj() == 0);
+  return 0;
+}