os/ossrv/genericopenlibs/liboil/tsrc/testsuite/dct/src/dct.c
author sl@SLION-WIN7.fritz.box
Fri, 15 Jun 2012 03:10:57 +0200
changeset 0 bde4ae8d615e
permissions -rw-r--r--
First public contribution.
sl@0
     1
/*
sl@0
     2
* Copyright (c) 2010 Nokia Corporation and/or its subsidiary(-ies).
sl@0
     3
* All rights reserved.
sl@0
     4
* This component and the accompanying materials are made available
sl@0
     5
* under the terms of "Eclipse Public License v1.0"
sl@0
     6
* which accompanies this distribution, and is available
sl@0
     7
* at the URL "http://www.eclipse.org/legal/epl-v10.html".
sl@0
     8
*
sl@0
     9
* Initial Contributors:
sl@0
    10
* Nokia Corporation - initial contribution.
sl@0
    11
*
sl@0
    12
* Contributors:
sl@0
    13
*
sl@0
    14
* Description: 
sl@0
    15
*
sl@0
    16
*/
sl@0
    17
sl@0
    18
sl@0
    19
#include <liboil/liboil.h>
sl@0
    20
#include <liboil/liboildebug.h>
sl@0
    21
#include <liboil/liboilfunction.h>
sl@0
    22
#include <stdio.h>
sl@0
    23
#include <stdlib.h>
sl@0
    24
sl@0
    25
#include <liboil/globals.h>
sl@0
    26
sl@0
    27
#define LOG_FILE "c:\\logs\\testsuite_dct_log.txt"
sl@0
    28
#include "std_log_result.h"
sl@0
    29
#include "utils.h"
sl@0
    30
#define LOG_FILENAME_LINE __FILE__, __LINE__
sl@0
    31
sl@0
    32
void create_xml(int result)
sl@0
    33
{
sl@0
    34
    if(result)
sl@0
    35
        assert_failed = 1;
sl@0
    36
    
sl@0
    37
    testResultXml("testsuite_dct");
sl@0
    38
    close_log_file();
sl@0
    39
}
sl@0
    40
sl@0
    41
void test_oil_mdct12_f64()
sl@0
    42
    {
sl@0
    43
    //double * d_6, const double * s_12
sl@0
    44
    double output[6];
sl@0
    45
    double input[12];
sl@0
    46
    double linux_output[] = {-55.78726058370923368557,-9.77999142525906961509,8.30695505006328360764,5.45142094924720943538,-5.13339728852725407648,-4.58890662311381625216};
sl@0
    47
    int i = 0;
sl@0
    48
    
sl@0
    49
    for(i=0; i<12; i++)
sl@0
    50
        {
sl@0
    51
        input[i] = i+1;
sl@0
    52
        
sl@0
    53
        if(i<6)
sl@0
    54
            output[i] = 0;
sl@0
    55
        }
sl@0
    56
    
sl@0
    57
    oil_mdct12_f64(output, input);
sl@0
    58
    
sl@0
    59
    for(i=0; i<6; i++)
sl@0
    60
        {
sl@0
    61
        if(comparefloats(output[i], linux_output[i]))
sl@0
    62
            {
sl@0
    63
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %15.14f, actual value - %15.14f", i,linux_output[i],output[i]);
sl@0
    64
            assert_failed = 1; 
sl@0
    65
            }
sl@0
    66
        }
sl@0
    67
    }
sl@0
    68
sl@0
    69
void test_oil_imdct12_f64()
sl@0
    70
    {
sl@0
    71
    //double * d_12, const double * s_6
sl@0
    72
    double output[12];
sl@0
    73
    double input[6];
sl@0
    74
    double linux_output[] = {-4.33833397075162174872,3.40568926865829446626,-3.31124125233533384005,3.31124125233533384005,-3.40568926865828158768,4.33833397075162441325,-4.80351335916817223648,10.06982028664590167466,-10.35088397845068541869,-10.35088397845079377646,10.06982028664594963629,-4.80351335916817401284};
sl@0
    75
    int i = 0;
sl@0
    76
    
sl@0
    77
    for(i=0; i<12; i++)
sl@0
    78
        {
sl@0
    79
        output[i] = 0;
sl@0
    80
        
sl@0
    81
        if(i<6)
sl@0
    82
            input[i] = i + 1;
sl@0
    83
        }
sl@0
    84
    
sl@0
    85
    oil_imdct12_f64(output, input);
sl@0
    86
    
sl@0
    87
    for(i=0; i<12; i++)
sl@0
    88
        {
sl@0
    89
        if(comparefloats(output[i], linux_output[i]))
sl@0
    90
            {
sl@0
    91
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %15.14f, actual value - %15.14f", i,linux_output[i],output[i]);
sl@0
    92
            assert_failed = 1; 
sl@0
    93
            }
sl@0
    94
        }
sl@0
    95
    }
sl@0
    96
sl@0
    97
void test_oil_mdct36_f64()
sl@0
    98
    {
sl@0
    99
    //double * d_18, const double * s_36
sl@0
   100
    double output[18];
sl@0
   101
    double input[36];
sl@0
   102
    double linux_output[] = {-485.54509125873983066413,-79.64662612510757355722,67.80760448776190685294,39.77358393387979162981,-36.41398520030525531865,-26.85943249683056777144,25.37958877703877291765,20.73179262231989383736,-19.93414220762220523397,-17.30353399124829749667,16.83079383688368935168,15.24349952210559955290,-14.95531193252001500582,-14.00077432972704016834,13.83319274912650698184,13.31880337261293334450,-13.24130282007020475987,-13.07848679902590305346};
sl@0
   103
    int i = 0;
sl@0
   104
    
sl@0
   105
    for(i=0; i<36; i++)
sl@0
   106
        {
sl@0
   107
        input[i] = i+1;
sl@0
   108
        
sl@0
   109
        if(i<18)
sl@0
   110
            output[i] = 0;
sl@0
   111
        }
sl@0
   112
    
sl@0
   113
    oil_mdct36_f64(output, input);
sl@0
   114
    
sl@0
   115
    for(i=0; i<18; i++)
sl@0
   116
        {
sl@0
   117
        if(comparefloats(output[i], linux_output[i]))
sl@0
   118
            {
sl@0
   119
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %15.14f, actual value - %15.14f", i,linux_output[i],output[i]);
sl@0
   120
            assert_failed = 1; 
sl@0
   121
            }
sl@0
   122
        }
sl@0
   123
    }
sl@0
   124
sl@0
   125
void test_oil_imdct36_f64()
sl@0
   126
    {
sl@0
   127
    //double * d_36, const double * s_18
sl@0
   128
    double output[36];
sl@0
   129
    double input[18];
sl@0
   130
    double linux_output[] = {-12.85687531711769970855,11.41757120437667261115,-11.15646549643737905910,10.28156353478235018883,-10.12421323499659031597,9.61624441673711061185,-9.53135483275377914936,9.29662075491665795823,-9.26973797227096696361,9.26973797227097051632,-9.29662075491650519155,9.53135483275377559664,-9.61624441673730778746,10.12421323499673242452,-10.28156353478229512177,11.15646549643748208780,-11.41757120437673655999,12.85687531711769793219,-13.28789667691522780046,15.72998203713914833202,-16.48536211104631377111,21.07260369137263467110,-22.59425383526346919894,33.39775527158744949929,-37.52696988457606863676,85.41533571580391992484,-80.79110748109037842823,-80.79110748109157214003,85.41533571580409045509,-37.52696988457621074531,33.39775527158775503267,-22.59425383526349762064,21.07260369137274835794,-16.48536211104587678733,15.72998203713892451105,-13.28789667691528109117};
sl@0
   131
    int i = 0;
sl@0
   132
    
sl@0
   133
    for(i=0; i<36; i++)
sl@0
   134
        {
sl@0
   135
        output[i] = 0;
sl@0
   136
        
sl@0
   137
        if(i<18)
sl@0
   138
            input[i] = i+1;
sl@0
   139
        }
sl@0
   140
    
sl@0
   141
    oil_imdct36_f64(output, input);
sl@0
   142
    
sl@0
   143
    for(i=0; i<36; i++)
sl@0
   144
        {
sl@0
   145
        if(comparefloats(output[i], linux_output[i]))
sl@0
   146
            {
sl@0
   147
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %15.14f, actual value - %15.14f", i,linux_output[i],output[i]);
sl@0
   148
            assert_failed = 1; 
sl@0
   149
            }
sl@0
   150
        }
sl@0
   151
    }
sl@0
   152
sl@0
   153
void test_oil_dct36_f32()
sl@0
   154
    {
sl@0
   155
    //float * d_36, int dstr, const float * s_36, int sstr
sl@0
   156
    float output[36];
sl@0
   157
    float input[36];
sl@0
   158
    float linux_output[] = {2.000000,3.000000,4.000000,5.000000,6.000000,7.000000,8.000000,9.000000,10.000000,11.000000,12.000000,13.000000,14.000000,15.000000,16.000000,17.000000,18.000000,19.000000,20.000000,21.000000,22.000000,23.000000,24.000000,25.000000,26.000000,27.000000,28.000000,29.000000,30.000000,31.000000,32.000000,33.000000,34.000000,35.000000,36.000000,37.000000};
sl@0
   159
    int i = 0;
sl@0
   160
    
sl@0
   161
    for(i=0; i<36; i++)
sl@0
   162
        {
sl@0
   163
        output[i] = i+2;
sl@0
   164
        input[i] = i+1;
sl@0
   165
        }
sl@0
   166
    
sl@0
   167
    oil_dct36_f32(output, 1, input, 2);
sl@0
   168
    
sl@0
   169
    for(i=0; i<36; i++)
sl@0
   170
        {
sl@0
   171
        if(output[i] != linux_output[i])
sl@0
   172
            {
sl@0
   173
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %f, actual value - %f", i,linux_output[i],output[i]);
sl@0
   174
            assert_failed = 1; 
sl@0
   175
            }
sl@0
   176
        }
sl@0
   177
    }
sl@0
   178
sl@0
   179
void test_oil_fdct8x8s_s16()
sl@0
   180
    {
sl@0
   181
    //int16_t * d_8x8, int ds, const int16_t * s_8x8, int ss
sl@0
   182
    int16_t output[64];
sl@0
   183
    int16_t input[64];
sl@0
   184
    int16_t linux_output[] = {13824,-4608,-1024,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
sl@0
   185
    int i = 0;
sl@0
   186
    
sl@0
   187
    for(i=0; i<64; i++)
sl@0
   188
        {
sl@0
   189
        output[i] = 0;
sl@0
   190
        input[i] = i+1;
sl@0
   191
        }
sl@0
   192
    
sl@0
   193
    oil_fdct8x8s_s16(output, 1, input, 2);
sl@0
   194
    
sl@0
   195
    for(i=0; i<64; i++)
sl@0
   196
        {
sl@0
   197
        if(output[i] != linux_output[i])
sl@0
   198
            {
sl@0
   199
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %d, actual value - %d", i,linux_output[i],output[i]);
sl@0
   200
            assert_failed = 1; 
sl@0
   201
            }
sl@0
   202
        }
sl@0
   203
    }
sl@0
   204
sl@0
   205
void test_oil_fdct8x8theora()
sl@0
   206
    {
sl@0
   207
    //const int16_t * s_8x8, int16_t * d_8x8
sl@0
   208
    int16_t output[64];
sl@0
   209
    int16_t input[64];
sl@0
   210
    int16_t linux_output[] = {1036,-62,0,-5,0,-5,0,-5,-582,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-61,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,-20,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-6,0,0,0,0,0,0,0};
sl@0
   211
    int i = 0;
sl@0
   212
    
sl@0
   213
    for(i=0; i<64; i++)
sl@0
   214
        {
sl@0
   215
        output[i] = 0;
sl@0
   216
        input[i] = i+1;
sl@0
   217
        }
sl@0
   218
    
sl@0
   219
    oil_fdct8x8theora(input, output);
sl@0
   220
    
sl@0
   221
    for(i=0; i<64; i++)
sl@0
   222
        {
sl@0
   223
        if(output[i] != linux_output[i])
sl@0
   224
            {
sl@0
   225
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %d, actual value - %d", i,linux_output[i],output[i]);
sl@0
   226
            assert_failed = 1; 
sl@0
   227
            }
sl@0
   228
        }
sl@0
   229
    }
sl@0
   230
sl@0
   231
void test_oil_idct8x8lim10_s16()
sl@0
   232
    {
sl@0
   233
    //int16_t * d_8x8, int dstr, const int16_t * s_8x8, int sstr
sl@0
   234
    int16_t output[64];
sl@0
   235
    int16_t input[64];
sl@0
   236
    int16_t linux_output[] = {1292,-771,768,-511,-1,255,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
sl@0
   237
    int i = 0;
sl@0
   238
    
sl@0
   239
    for(i=0; i<64; i++)
sl@0
   240
        {
sl@0
   241
        output[i] = 0;
sl@0
   242
        input[i] = i+1;
sl@0
   243
        }
sl@0
   244
    
sl@0
   245
    oil_idct8x8lim10_s16(output, 1, input, 2);
sl@0
   246
    
sl@0
   247
    for(i=0; i<64; i++)
sl@0
   248
        {
sl@0
   249
        if(output[i] != linux_output[i])
sl@0
   250
            {
sl@0
   251
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %d, actual value - %d", i,linux_output[i],output[i]);
sl@0
   252
            assert_failed = 1; 
sl@0
   253
            }
sl@0
   254
        }
sl@0
   255
    }
sl@0
   256
sl@0
   257
void test_oil_idct8theora_s16()
sl@0
   258
    {
sl@0
   259
    //int16_t * d_8, int dstr, const int16_t * s_8, int sstr
sl@0
   260
    int16_t output[8];
sl@0
   261
    int16_t input[8];
sl@0
   262
    int16_t linux_output[] = {16,-18,6,-5,5,-2,2,0};
sl@0
   263
    int i = 0;
sl@0
   264
    
sl@0
   265
    for(i=0; i<8; i++)
sl@0
   266
        {
sl@0
   267
        output[i] = 0;
sl@0
   268
        input[i] = i+1;
sl@0
   269
        }
sl@0
   270
    
sl@0
   271
    oil_idct8theora_s16(output, 1, input, 2);
sl@0
   272
    
sl@0
   273
    for(i=0; i<8; i++)
sl@0
   274
        {
sl@0
   275
        if(output[i] != linux_output[i])
sl@0
   276
            {
sl@0
   277
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %d, actual value - %d", i,linux_output[i],output[i]);
sl@0
   278
            assert_failed = 1; 
sl@0
   279
            }
sl@0
   280
        }
sl@0
   281
    }
sl@0
   282
sl@0
   283
void test_oil_idct8x8theora_s16()
sl@0
   284
    {
sl@0
   285
    //int16_t * d_8x8, int dstr, const int16_t * s_8x8, int sstr
sl@0
   286
    int16_t output[64];
sl@0
   287
    int16_t input[64];
sl@0
   288
    int16_t linux_output[] = {157,-123,50,-36,26,-14,10,1,-110,51,-30,12,-16,2,-10,-2,52,-31,10,-10,10,0,6,0,-34,19,-11,6,-7,1,-3,-2,26,-15,7,-4,5,-1,3,0,-14,5,-2,2,0,-2,0,-2,12,-9,2,0,6,0,4,0,-1,-1,2,-2,4,2,-2,1};
sl@0
   289
    int i = 0;
sl@0
   290
    
sl@0
   291
    for(i=0; i<64; i++)
sl@0
   292
        {
sl@0
   293
        output[i] = 0;
sl@0
   294
        input[i] = i+1;
sl@0
   295
        }
sl@0
   296
    
sl@0
   297
    oil_idct8x8theora_s16(output, 1, input, 2);
sl@0
   298
    
sl@0
   299
    for(i=0; i<64; i++)
sl@0
   300
        {
sl@0
   301
        if(output[i] != linux_output[i])
sl@0
   302
            {
sl@0
   303
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %d, actual value - %d", i,linux_output[i],output[i]);
sl@0
   304
            assert_failed = 1; 
sl@0
   305
            }
sl@0
   306
        }
sl@0
   307
    }
sl@0
   308
sl@0
   309
void test_oil_imdct32_f32()
sl@0
   310
    {
sl@0
   311
    //float * d_32, const float * s_32
sl@0
   312
    float output[32];
sl@0
   313
    float input[32];
sl@0
   314
    float expected_output[] = {528.00000000000000000000,-207.42237854003906250000,-0.00000000000002593932,-22.97223091125488281250,-0.00000000000002312733,-8.21512222290039062500,0.00000000000000710542,-4.14796209335327148438,-0.00000000000003045654,-2.47257304191589355469,0.00000000000000630225,-1.62263202667236328125,0.00000000000008801467,-1.13173198699951171875,-0.00000000000003877801,-0.82146918773651123047,0.00000000000035616302,-0.61161065101623535156,0.00000000000014964419,-0.46168053150177001953,0.00000000000006150462,-0.34939777851104736328,-0.00000000000018143386,-0.26159864664077758789,0.00000000000003763656,-0.19001007080078125000,-0.00000000000012429814,-0.12911282479763031006,-0.00000000000002208823,-0.07497953623533248901,-0.00000000000003485602,-0.02459304779767990112};
sl@0
   315
    int i = 0;
sl@0
   316
    
sl@0
   317
    for(i=0; i<32; i++)
sl@0
   318
        {
sl@0
   319
        input[i] = i+1;
sl@0
   320
        output[i] = 0;
sl@0
   321
        }
sl@0
   322
    
sl@0
   323
    oil_imdct32_f32(output, input);
sl@0
   324
    
sl@0
   325
    for(i=0; i<32; i++)
sl@0
   326
        {
sl@0
   327
        if(comparefloats(output[i], expected_output[i]))
sl@0
   328
            {
sl@0
   329
            std_log(LOG_FILENAME_LINE, "output[%d]: expected value - %15.14f, actual value - %15.14f", i,expected_output[i],output[i]);
sl@0
   330
            assert_failed = 1; 
sl@0
   331
            }
sl@0
   332
        }
sl@0
   333
    }
sl@0
   334
sl@0
   335
int main (int argc, char *argv[])
sl@0
   336
{
sl@0
   337
    oil_init ();
sl@0
   338
 
sl@0
   339
	std_log(LOG_FILENAME_LINE,"START oil_mdct12_f64 TEST");
sl@0
   340
    test_oil_mdct12_f64(); //fails because of float precision diff
sl@0
   341
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   342
sl@0
   343
	std_log(LOG_FILENAME_LINE,"START oil_imdct12_f64 TEST");
sl@0
   344
    test_oil_imdct12_f64(); //fails because of float precision diff
sl@0
   345
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   346
sl@0
   347
	std_log(LOG_FILENAME_LINE,"START oil_mdct36_f64 TEST");
sl@0
   348
    test_oil_mdct36_f64(); //fails because of float precision diff
sl@0
   349
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   350
sl@0
   351
	std_log(LOG_FILENAME_LINE,"START oil_imdct36_f64 TEST");
sl@0
   352
    test_oil_imdct36_f64(); //fails because of float precision diff
sl@0
   353
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   354
    
sl@0
   355
	std_log(LOG_FILENAME_LINE,"START oil_dct36_f32 TEST");
sl@0
   356
    test_oil_dct36_f32();
sl@0
   357
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   358
sl@0
   359
	std_log(LOG_FILENAME_LINE,"START oil_fdct8x8s_s16 TEST");
sl@0
   360
    test_oil_fdct8x8s_s16();
sl@0
   361
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   362
sl@0
   363
	std_log(LOG_FILENAME_LINE,"START oil_fdct8x8theora TEST");
sl@0
   364
    test_oil_fdct8x8theora();
sl@0
   365
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   366
sl@0
   367
	std_log(LOG_FILENAME_LINE,"START oil_idct8x8lim10_s16 TEST");
sl@0
   368
    test_oil_idct8x8lim10_s16();
sl@0
   369
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   370
sl@0
   371
	std_log(LOG_FILENAME_LINE,"START oil_idct8theora_s16 TEST");
sl@0
   372
    test_oil_idct8theora_s16();
sl@0
   373
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   374
sl@0
   375
	std_log(LOG_FILENAME_LINE,"START oil_idct8x8theora_s16 TEST");
sl@0
   376
    test_oil_idct8x8theora_s16();
sl@0
   377
	std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   378
sl@0
   379
	std_log(LOG_FILENAME_LINE,"START oil_imdct32_f32 TEST");
sl@0
   380
    test_oil_imdct32_f32(); //getting different output than linux
sl@0
   381
    std_log(LOG_FILENAME_LINE,"END TEST\n");
sl@0
   382
sl@0
   383
    if(assert_failed)
sl@0
   384
      std_log(LOG_FILENAME_LINE,"Test Failed");
sl@0
   385
    else
sl@0
   386
      std_log(LOG_FILENAME_LINE,"Test Successful");
sl@0
   387
    
sl@0
   388
    create_xml(0);
sl@0
   389
    return 0;
sl@0
   390
}
sl@0
   391