diff options
Diffstat (limited to 'src/statisticsFunctions/variance/testFloatVariance.c')
| -rw-r--r-- | src/statisticsFunctions/variance/testFloatVariance.c | 329 |
1 files changed, 0 insertions, 329 deletions
diff --git a/src/statisticsFunctions/variance/testFloatVariance.c b/src/statisticsFunctions/variance/testFloatVariance.c deleted file mode 100644 index 97b56d4d..00000000 --- a/src/statisticsFunctions/variance/testFloatVariance.c +++ /dev/null @@ -1,329 +0,0 @@ -/* - * 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-06 ) ; - } - -} - - - - -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; -} |