/*
* 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.323312455
#define LINES 7
#define COLUMNS 10
#define MATRIX {0.102326436434,0.939044147730,0.638607255649,0.059883427806,0.311199211050,\
0.377510951832,0.909709410276,0.928138985299,0.855155082885,0.308869514614,\
0.729239261709,0.127006734256,0.242402311414,0.807912770193,0.621117659844,\
0.081452373881,0.088844371960,0.587072817143,0.068080495112,0.995457543526,\
0.835322400089,0.656235548668,0.767905956600,0.706818781793,0.187699332368,\
0.568810169585,0.401039466728,0.397196868435,0.483212330844,0.161224133335,\
0.932864250615,0.719208867755,0.262409000657,0.931377100293,0.596724558622,\
0.310327019077,0.350161732174,0.104841318447,0.830732351169,0.804621329065,\
0.050784428138,0.551920620725,0.125952405389,0.083156193607,0.420623748098,\
0.168616071343,0.707545555662,0.104633864947,0.465604968835,0.489758284763,\
0.167727907188,0.760689914692,0.454900974408,0.063926273491,0.87152507063,0.978668818716,\
0.135898929089,0.484713065904,0.785042577423,0.049412418623,\
0.947593610268,0.884073690511,0.212796610314,0.654832502827,0.007690255996,\
0.031808692031,0.706323316786,0.008731670678,0.270014822017,0.153652121313}
#define MATRIXROW {0.168824942350,0.072908122527,0.057987977092,0.149998975703,0.085086319404,\
0.108508638451,0.096937985957,0.107281660264,0.091363871817,0.128577818421}
#define MATRIXCOL { 0.123307827428,0.123213594199,0.054037725959,0.092551230301,0.055505564488,\
0.127629921958,0.138516743375 }
#define IMATRIX {0.05826747464016080,0.99946373142302036,0.85600351681932807,0.62329693790525198,\
0.49792320514097810,0.44210509583353996,0.55833499459549785,0.09633230511099100,\
0.80100952507928014,0.34466254524886608,\
0.48831309471279383,0.95087061496451497,0.01332767866551876,0.76261807000264525,\
0.74346329551190138,0.65157829830422997,0.57006288319826126,0.70580983115360141,\
0.51323400903493166,0.64977857517078519,\
0.80352442665025592,0.69910932797938585,0.69258948462083936,0.91169391758739948,\
0.92454590043053031,0.27637310232967138,0.31692579621449113,0.86305770650506020,\
0.20109100220724940,0.01225362811237574,\
0.38199013099074364,0.96006405679509044,0.16440964583307505,0.59608811559155583,\
0.72101737372577190,0.58661046391353011,0.99326277803629637,0.00761850038543344,\
0.78608208894729614,0.89965870184823871,\
0.43115562805905938,0.81857266277074814,0.83410377753898501,0.55516970623284578,\
0.51345925079658628,0.04770902730524540,0.80747798969969153,0.80489510949701071,\
0.78650354826822877,0.30791273340582848,\
0.75731822755187750,0.53423820668831468,0.42618893459439278,0.85448804078623652,\
0.92601215932518244,0.94169309409335256,0.85547966323792934,0.59637623047456145,\
0.69515300076454878,0.18359116325154901,\
0.01200280850753188,0.2255702270194888,0.34853330114856362,0.55313225090503693,\
0.63695094687864184,0.96132039744406939,0.50314606027677655,0.11768362112343311,\
0.42481321236118674,0.52963322307914495}
#define RMATRIX {0.53868199465796351,0.65326874051243067,0.37497402401641011,0.60422550700604916,\
0.33622304117307067,0.03531436901539564,0.33213760564103723,0.39493087679147720,\
0.04170337272807956,0.87363853026181459,\
0.65950810909271240,0.31221040291711688,0.15333442110568285,0.66478573577478528,\
0.44319023378193378,0.40884594758972526,0.94536898937076330,0.67733758920803666,\
0.18673646822571754,0.92331133363768458,\
0.55921846115961671,0.75012728199362755,0.23257926432415843,0.57102064136415720,\
0.60969385923817754,0.67336730472743511,0.07624059682711959,0.10437540244311094,\
0.95672677317634225,0.19771346449851990,\
0.05278092902153730,0.75808868417516351,0.17633479088544846,0.01146994484588504,\
0.36732212174683809,0.19948644982650876,0.52714426256716251,0.21445603063330054,\
0.94751045759767294,0.04272260749712586,\
0.03257346292957664,0.53217577841132879,0.87635089689865708,0.15471007302403450,\
0.25448470888659358,0.91815057490020990,0.72509902389720082,0.87545845471322536,\
0.05635281419381499,0.01432112138718367,\
0.57406943850219250,0.80984140699729323,0.01664119493216276,0.93888836959376931,\
0.11264799535274506,0.66036546928808093,0.76042845565825701,0.25908330874517560,\
0.88014078326523304,0.70332178613170981,\
0.94614937948063016,0.78529163636267185,0.20241560926660895,0.10981105919927359,\
0.75415370846167207,0.15320260450243950,0.59521253732964396,0.71318271104246378,\
0.51107599260285497,0.52125945501029491}
#define IMATRIXROW {-0.0419293790250184,-0.05416405376783551,+0.11775191182077104,+0.08015791690354153,+0.00011032748454003,\
-0.10124210045965254,0.06720614014223115,+0.02470873836336094,-0.07082037296251543,-0.00621987869853545}
#define RMATRIXROW {0.01439621349781126,-0.04697393442028135,-0.03018604392150223,0.09765990436118642,0.01606884471097247,\
-0.00604250542333471,0.02870576510687836,-0.04985742926427023,0.12356001472881707,0.05909180726174274}
#define IMATRIXCOL {-0.01904810719805961,+0.04121453132770223,+0.00378553967089328,\
+0.10307670507342410,+0.03497427750676810,+0.01761849804019567,-0.10182321471060848}
#define RMATRIXCOL {-0.02739893277693160,0.01792852626358409,-0.02205084134302381,\
-0.01012610864032573,0.07317682977907164,0.04762062457313784,0.00849011627117621}
static void dvariancesTest (void ) {
double toTest = dvariances (SINGLE) ;
assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
}
static void drowvariancesTest (void ) {
double toTest = drowvariances (SINGLE) ;
assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
}
static void dcolumnvariancesTest (void) {
double toTest = dcolumnvariances (SINGLE) ;
assert ( fabs ( toTest - SINGLE ) / fabs( toTest) < 1e-06 ) ;
}
static void zvariancesTest (void ) {
doubleComplex toTest = zvariances ( DoubleComplex ( SINGLE , SINGLE ) );
assert ( fabs ( zreals( toTest) - SINGLE ) / fabs( zreals( toTest)) < 1e-06 ) ;
assert ( fabs ( zimags( toTest) - SINGLE ) / fabs( zimags( toTest)) < 1e-06 ) ;
}
static void zrowvariancesTest (void ) {
doubleComplex toTest = zrowvariances ( DoubleComplex ( SINGLE , SINGLE ) );
assert ( fabs ( zreals( toTest) - SINGLE ) / fabs( zreals( toTest)) < 1e-06 ) ;
assert ( fabs ( zimags( toTest) - SINGLE ) / fabs( zimags( toTest)) < 1e-06 ) ;
}
static void zcolumnvariancesTest (void ) {
doubleComplex toTest = zrowvariances ( DoubleComplex ( SINGLE , SINGLE ) );
assert ( fabs ( zreals( toTest) - SINGLE ) / fabs( zreals( toTest)) < 1e-06 ) ;
assert ( fabs ( zimags( toTest) - SINGLE ) / fabs( zimags( toTest)) < 1e-06 ) ;
}
static void dvarianceaTest ( void){
double mtoTest[] = MATRIX ;
double result = 0.100753297434;
double out = dvariancea ( mtoTest , LINES*COLUMNS ); ;
printf ( "out : %e\t result : %e\t\n" , out , result ) ;
assert ( fabs ( out - result ) / fabs( out ) < 1e-06 ) ;
}
static void drowvarianceaTest ( void )
{
int i = 0 ;
double mtoTest[] = MATRIX ;
double result[COLUMNS] = MATRIXROW ;
double out[COLUMNS] ;
drowvariancea ( 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 dcolumnvarianceaTest ( void )
{
int i = 0 ;
double mtoTest[] = MATRIX ;
double result[LINES] = MATRIXCOL ;
double out[ LINES] ;
dcolumnvariancea ( 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 zvarianceaTest ( void){
double rmtoTest [] = RMATRIX ;
double imtoTest [] = IMATRIX ;
doubleComplex Result = DoubleComplex ( 0.01225832057910708 , 0.01071188208687752 );
doubleComplex* mtoTest = DoubleComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS ) ;
doubleComplex out = zvariancea ( mtoTest , LINES*COLUMNS ); ;
assert ( fabs( zreals(out) - zreals (Result) ) / fabs (zreals (out)) < 1e-06 );
assert( fabs( zimags(out) - zimags (Result )) / fabs (zimags (out)) < 1e-06 );
}
static void zrowvarianceaTest (void ) {
int i = 0 ;
double rmtoTest [] = RMATRIX ;
double imtoTest [] = IMATRIX ;
double rResult [] = RMATRIXROW;
double iResult [] = IMATRIXROW ;
doubleComplex out[COLUMNS];
doubleComplex* mtoTest = DoubleComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS ) ;
doubleComplex* Result = DoubleComplexMatrix ( rResult , iResult , COLUMNS );
zrowvariancea ( 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 , zreals(out[i]) , zimags(out[i]), zreals(Result[i]) ,zimags( Result[i]));
assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-04 );
assert( fabs( zimags(out[i]) - zimags (Result[i])) / fabs (zimags (out[i])) < 1e-04);
}
}
static void zcolumnvarianceaTest ( void) {
int i = 0 ;
double rmtoTest [] = RMATRIX ;
double imtoTest [] = IMATRIX ;
double rResult [] = RMATRIXCOL;
double iResult [] = IMATRIXCOL ;
doubleComplex out[LINES];
doubleComplex* mtoTest = DoubleComplexMatrix ( rmtoTest , imtoTest , LINES*COLUMNS ) ;
doubleComplex* Result = DoubleComplexMatrix ( rResult , iResult , LINES );
zcolumnvariancea ( mtoTest , LINES , COLUMNS , out ) ;
for ( i = 0 ; i < LINES ; i++)
{
printf ( "%d out : %e %e \t result %e %e \n" , i , zreals(out[i]) , zimags(out[i]), zreals(Result[i]) ,zimags( Result[i]));
assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-06 );
assert( fabs( zimags(out[i]) - zimags (Result[i])) / fabs (zimags (out[i])) < 1e-06);
}
}
static int testDoubleVariance (void) {
printf("\n\n\n\n*********************\n");
printf("***** Double Tests ****\n");
printf("*********************\n");
printf("\n\t>>>>singleton real Tests\n");
dcolumnvariancesTest () ;
dvariancesTest () ;
drowvariancesTest () ;
printf("\n\t>>>>singleton complex Tests\n");
zrowvariancesTest () ;
zcolumnvariancesTest () ;
zvariancesTest () ;
printf("\n\t>>>>array real Tests\n");
dvarianceaTest () ;
dcolumnvarianceaTest () ;
drowvarianceaTest ( );
printf("\n\t>>>>array complex Tests\n");
zvarianceaTest () ;
zrowvarianceaTest ( );
zcolumnvarianceaTest () ;
return 0;
}
int main(void) {
assert(testDoubleVariance () == 0);
return 0;
}