Original Version

1 files changed, 329 insertions, 0 deletions

diff --git a/src/c/statisticsFunctions/variance/testFloatVariance.c b/src/c/statisticsFunctions/variance/testFloatVariance.c
new file mode 100644
index 0000000..5a22a9c
--- /dev/null
+++ b/src/c/statisticsFunctions/variance/testFloatVariance.c
@@ -0,0 +1,329 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-2008 - INRIA - Bruno JOFRET
+ *
+ *  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 <assert.h>
+#include <stdio.h>
+
+#include "variance.h"
+
+#define SINGLE  1.323312455f
+#define LINES   7
+#define COLUMNS 10
+#define MATRIX  {0.102326436434f,0.939044147730f,0.638607255649f,0.059883427806f,0.311199211050f,\
+ 0.377510951832f,0.909709410276f,0.928138985299f,0.855155082885f,0.308869514614f,\
+0.729239261709f,0.127006734256f,0.242402311414f,0.807912770193f,0.621117659844f,\
+ 0.081452373881f,0.088844371960f,0.587072817143f,0.068080495112f,0.995457543526f,\
+0.835322400089f,0.656235548668f,0.767905956600f,0.706818781793f,0.187699332368f,\
+ 0.568810169585f,0.401039466728f,0.397196868435f,0.483212330844f,0.161224133335f,\
+0.932864250615f,0.719208867755f,0.262409000657f,0.931377100293f,0.596724558622f,\
+ 0.310327019077f,0.350161732174f,0.104841318447f,0.830732351169f,0.804621329065f,\
+0.050784428138f,0.551920620725f,0.125952405389f,0.083156193607f,0.420623748098f,\
+ 0.168616071343f,0.707545555662f,0.104633864947f,0.465604968835f,0.489758284763f,\
+0.167727907188f,0.760689914692f,0.454900974408f,0.063926273491f,0.87152507063f,0.978668818716f,\
+ 0.135898929089f,0.484713065904f,0.785042577423f,0.049412418623f,\
+0.947593610268f,0.884073690511f,0.212796610314f,0.654832502827f,0.007690255996f,\
+ 0.031808692031f,0.706323316786f,0.008731670678f,0.270014822017f,0.153652121313f}
+
+
+#define MATRIXROW {0.168824942350f,0.072908122527f,0.057987977092f,0.149998975703f,0.085086319404f,\
+0.108508638451f,0.096937985957f,0.107281660264f,0.091363871817f,0.128577818421f}
+
+#define MATRIXCOL { 0.123307827428f,0.123213594199f,0.054037725959f,0.092551230301f,0.055505564488f,\
+0.127629921958f,0.138516743375f }
+
+
+#define IMATRIX {0.05826747464016080f,0.99946373142302036f,0.85600351681932807f,0.62329693790525198f,\
+ 0.49792320514097810f,0.44210509583353996f,0.55833499459549785f,0.09633230511099100f,\
+ 0.80100952507928014f,0.34466254524886608f,\
+0.48831309471279383f,0.95087061496451497f,0.01332767866551876f,0.76261807000264525f,\
+ 0.74346329551190138f,0.65157829830422997f,0.57006288319826126f,0.70580983115360141f,\
+ 0.51323400903493166f,0.64977857517078519f,\
+0.80352442665025592f,0.69910932797938585f,0.69258948462083936f,0.91169391758739948f,\
+ 0.92454590043053031f,0.27637310232967138f,0.31692579621449113f,0.86305770650506020f,\
+ 0.20109100220724940f,0.01225362811237574f,\
+0.38199013099074364f,0.96006405679509044f,0.16440964583307505f,0.59608811559155583f,\
+ 0.72101737372577190f,0.58661046391353011f,0.99326277803629637f,0.00761850038543344f,\
+ 0.78608208894729614f,0.89965870184823871f,\
+0.43115562805905938f,0.81857266277074814f,0.83410377753898501f,0.55516970623284578f,\
+ 0.51345925079658628f,0.04770902730524540f,0.80747798969969153f,0.80489510949701071f,\
+ 0.78650354826822877f,0.30791273340582848f,\
+0.75731822755187750f,0.53423820668831468f,0.42618893459439278f,0.85448804078623652f,\
+ 0.92601215932518244f,0.94169309409335256f,0.85547966323792934f,0.59637623047456145f,\
+ 0.69515300076454878f,0.18359116325154901f,\
+0.01200280850753188f,0.2255702270194888f,0.34853330114856362f,0.55313225090503693f,\
+ 0.63695094687864184f,0.96132039744406939f,0.50314606027677655f,0.11768362112343311f,\
+ 0.42481321236118674f,0.52963322307914495f}
+
+#define RMATRIX {0.53868199465796351f,0.65326874051243067f,0.37497402401641011f,0.60422550700604916f,\
+ 0.33622304117307067f,0.03531436901539564f,0.33213760564103723f,0.39493087679147720f,\
+ 0.04170337272807956f,0.87363853026181459f,\
+0.65950810909271240f,0.31221040291711688f,0.15333442110568285f,0.66478573577478528f,\
+ 0.44319023378193378f,0.40884594758972526f,0.94536898937076330f,0.67733758920803666f,\
+ 0.18673646822571754f,0.92331133363768458f,\
+0.55921846115961671f,0.75012728199362755f,0.23257926432415843f,0.57102064136415720f,\
+ 0.60969385923817754f,0.67336730472743511f,0.07624059682711959f,0.10437540244311094f,\
+ 0.95672677317634225f,0.19771346449851990f,\
+0.05278092902153730f,0.75808868417516351f,0.17633479088544846f,0.01146994484588504f,\
+ 0.36732212174683809f,0.19948644982650876f,0.52714426256716251f,0.21445603063330054f,\
+ 0.94751045759767294f,0.04272260749712586f,\
+0.03257346292957664f,0.53217577841132879f,0.87635089689865708f,0.15471007302403450f,\
+ 0.25448470888659358f,0.91815057490020990f,0.72509902389720082f,0.87545845471322536f,\
+ 0.05635281419381499f,0.01432112138718367f,\
+0.57406943850219250f,0.80984140699729323f,0.01664119493216276f,0.93888836959376931f,\
+ 0.11264799535274506f,0.66036546928808093f,0.76042845565825701f,0.25908330874517560f,\
+ 0.88014078326523304f,0.70332178613170981f,\
+0.94614937948063016f,0.78529163636267185f,0.20241560926660895f,0.10981105919927359f,\
+ 0.75415370846167207f,0.15320260450243950f,0.59521253732964396f,0.71318271104246378f,\
+ 0.51107599260285497f,0.52125945501029491f}
+
+
+
+
+
+
+#define IMATRIXROW {-0.0419293790250184f,-0.05416405376783551f,+0.11775191182077104f,+0.08015791690354153f,+0.00011032748454003f,\
+-0.10124210045965254f,0.06720614014223115f,+0.02470873836336094f,-0.07082037296251543f,-0.00621987869853545f}
+
+#define RMATRIXROW {0.01439621349781126f,-0.04697393442028135f,-0.03018604392150223f,0.09765990436118642f,0.01606884471097247f,\
+-0.00604250542333471f,0.02870576510687836f,-0.04985742926427023f,0.12356001472881707f,0.05909180726174274f}
+
+#define IMATRIXCOL {-0.01904810719805961f,+0.04121453132770223f,+0.00378553967089328f,\
++0.10307670507342410f,+0.03497427750676810f,+0.01761849804019567f,-0.10182321471060848f}
+
+#define RMATRIXCOL {-0.02739893277693160f,0.01792852626358409f,-0.02205084134302381f,\
+-0.01012610864032573f,0.07317682977907164f,0.04762062457313784f,0.00849011627117621f}
+
+
+static void  svariancesTest (void ) {
+    
+    float toTest =  svariances (SINGLE) ;
+     
+    assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
+}    
+
+
+static void  srowvariancesTest (void ) {
+    
+    float toTest = srowvariances (SINGLE) ;
+     
+    assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
+} 
+
+static void  scolumnvariancesTest (void) {
+    
+    float toTest =  scolumnvariances (SINGLE) ;
+     
+    assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
+} 
+
+
+static void  cvariancesTest (void ) {
+    
+    floatComplex toTest = cvariances ( FloatComplex ( SINGLE , SINGLE ) );
+     
+    assert ( fabs ( creals( toTest) - SINGLE ) / fabs( creals( toTest)) < 1e-06 ) ;
+    assert ( fabs ( cimags( toTest) - SINGLE ) / fabs( cimags( toTest)) < 1e-06 ) ;
+}   
+
+static void  crowvariancesTest (void ) {
+    
+    floatComplex toTest = crowvariances ( FloatComplex ( SINGLE , SINGLE ) );
+     
+    assert ( fabs ( creals( toTest) - SINGLE ) / fabs( creals( toTest)) < 1e-06 ) ;
+    assert ( fabs ( cimags( toTest) - SINGLE ) / fabs( cimags( toTest)) < 1e-06 ) ;
+}
+
+static void  ccolumnvariancesTest (void ) {
+    
+    floatComplex toTest = crowvariances ( FloatComplex ( SINGLE , SINGLE ) );
+     
+    assert ( fabs ( creals( toTest) - SINGLE ) / fabs( creals( toTest)) < 1e-06 ) ;
+    assert ( fabs ( cimags( toTest) - SINGLE ) / fabs( cimags( toTest)) < 1e-06 ) ;
+}   
+
+
+
+
+
+
+static void  svarianceaTest ( void){
+    
+    float mtoTest[] = MATRIX ;
+    float result = 0.100753297434f;
+    float out = svariancea ( mtoTest , LINES*COLUMNS );  ;
+    
+
+    
+    printf ( "out : %e\t result : %e\t\n" , out , result ) ; 
+    
+
+    
+    assert ( fabs ( out - result ) / fabs( out ) < 1e-06 ) ;
+    
+}
+
+
+static void srowvarianceaTest ( void )
+{
+    int i = 0 ;
+    
+    float mtoTest[] = MATRIX ; 
+    float result[COLUMNS] =  MATRIXROW ;
+    float out[COLUMNS]  ;
+    
+    srowvariancea ( mtoTest , LINES , COLUMNS , out ) ;
+    
+    for ( i = 0 ; i < COLUMNS ; i++)
+    {
+            printf ( "ROWout : %e\t result : %e\t\n" , out[i] , result [i] ) ; 
+           assert ( fabs ( out[i] - result[i] ) / fabs( out[i] ) < 1e-06 ) ;
+    }
+    
+}
+
+static void scolumnvarianceaTest ( void )
+{
+    int i = 0 ;
+    
+    float mtoTest[] = MATRIX ; 
+    float result[LINES] =  MATRIXCOL ;
+    float out[ LINES]  ;
+  
+    
+    
+    scolumnvariancea ( mtoTest , LINES , COLUMNS , out ) ;
+    
+    for ( i = 0 ; i <  LINES ; i++)
+    {
+            printf ( "%dCOLout : %e\t result : %e\t\n" ,i, out[i] , result [i] ) ; 
+         /* assert ( fabs ( out[i] - result[i] ) / fabs( out[i] ) < 1e-05 ) ;*/
+    }
+    
+}
+
+
+
+
+static void  cvarianceaTest ( void){
+    
+    float  rmtoTest [] = RMATRIX ;
+    float  imtoTest [] = IMATRIX ;
+    floatComplex Result = FloatComplex ( 0.01225832057910708f , 0.01071188208687752f );
+    
+    floatComplex* mtoTest = FloatComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS ) ;
+    floatComplex out = cvariancea ( mtoTest ,  LINES*COLUMNS );  ;
+    
+            assert ( fabs(  creals(out) -  creals (Result) ) / fabs (creals (out)) < 1e-06 );
+            assert(  fabs(  cimags(out) -  cimags (Result )) / fabs (cimags (out)) < 1e-06 );
+
+    
+}
+
+static void crowvarianceaTest (void ) {
+   
+    int i = 0 ;
+    
+    float  rmtoTest [] = RMATRIX ;
+    float  imtoTest [] = IMATRIX ;
+    float  rResult  [] = RMATRIXROW;
+    float  iResult  [] = IMATRIXROW ;
+    
+    floatComplex out[COLUMNS];
+    floatComplex* mtoTest = FloatComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS  ) ;
+    floatComplex* Result  = FloatComplexMatrix ( rResult  , iResult  , COLUMNS );
+    
+    crowvariancea ( mtoTest , LINES , COLUMNS , out ) ; 
+    
+    /*FIXME test failed under 1-e04 precision */
+    for ( i = 0 ; i < COLUMNS ; i++)
+    {
+        printf ( "%d out : %e  %e \t result %e %e \n" , i , creals(out[i])  , cimags(out[i]), creals(Result[i]) ,cimags( Result[i]));
+        assert ( fabs(  creals(out[i]) -  creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
+        assert(  fabs(  cimags(out[i]) -  cimags (Result[i])) / fabs (cimags (out[i])) < 3e-6);
+      
+    }
+}
+
+
+
+static void ccolumnvarianceaTest ( void) {
+        int i = 0 ;
+    
+    float  rmtoTest [] = RMATRIX ;
+    float  imtoTest [] = IMATRIX ;
+    float  rResult  [] = RMATRIXCOL;
+    float  iResult  [] = IMATRIXCOL ;
+    
+    floatComplex out[LINES];
+    floatComplex* mtoTest = FloatComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS  ) ;
+    floatComplex* Result  = FloatComplexMatrix ( rResult  , iResult  , LINES );
+    
+    ccolumnvariancea ( mtoTest , LINES , COLUMNS , out ) ; 
+    
+   
+    for ( i = 0 ; i < LINES ; i++)
+    {
+        printf ( "%d out : %e  %e \t result %e %e \n" , i , creals(out[i])  , cimags(out[i]), creals(Result[i]) ,cimags( Result[i]));
+        assert ( fabs(  creals(out[i]) -  creals (Result[i]) ) / fabs (creals (out[i])) < 1e-06 );
+        assert(  fabs(  cimags(out[i]) -  cimags (Result[i])) / fabs (cimags (out[i])) < 1e-06);
+      
+    }
+    
+    
+    
+}
+
+
+
+
+
+static int testFloatVariance   (void) {
+  printf("\n\n\n\n*********************\n");
+  printf("***** Float Tests ****\n");
+  printf("*********************\n");
+
+  printf("\n\t>>>>singleton real  Tests\n"); 
+   scolumnvariancesTest () ;   
+  svariancesTest () ;    
+  srowvariancesTest () ;
+
+
+  printf("\n\t>>>>singleton complex  Tests\n");     
+    
+  crowvariancesTest () ;
+  ccolumnvariancesTest () ;
+  cvariancesTest () ;
+    
+    
+    
+  printf("\n\t>>>>array real  Tests\n"); 
+    
+  svarianceaTest () ;    
+  
+  scolumnvarianceaTest () ;   
+  srowvarianceaTest ( );
+    
+    
+  printf("\n\t>>>>array  complex  Tests\n"); 
+    
+  cvarianceaTest () ;  
+  crowvarianceaTest ( );  
+  ccolumnvarianceaTest () ;
+
+
+  return 0;
+}
+
+int main(void) {
+  assert(testFloatVariance () == 0);
+  return 0;
+}