Bump version numbers to 6.2.0
[scilab.git] / scilab / modules / statistics / tests / unit_tests / cdft.dia.ref
1 //<-- CLI SHELL MODE -->
2 // =============================================================================
3 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4 // Copyright (C) ????-2008 - INRIA
5 //
6 //  This file is distributed under the same license as the Scilab package.
7 // =============================================================================
8 // =============================================================================
9 // Tests for cdft() function
10 // using a table
11 // =============================================================================
12 // Run with test_run('statistics','cdft',['no_check_error_output']);
13 prec    = 1.e-5;
14 Tab     = [0.9,0.5,0.3,0.20,0.10,0.05,0.02];
15 Df      = [1,2,3,4,5,6,7];
16 Th      = [0.158,0.816,1.250,1.533,2.015,2.447,2.998];
17 [P1,Q1] = cdft("PQ",Th,Df);
18 [P2,Q2] = cdft("PQ",-Th,Df);
19 if norm(Tab-(Q1+P2)) > 0.1 then bugmes();quit;end
20 [P,Q] = cdft("PQ",Th,Df);
21 Th1   = cdft("T",Df,P,Q);
22 Df1   = cdft("Df",P,Q,Th);
23 if norm(Th1-Th) > prec then bugmes();quit;end
24 if norm(Df1-Df) > prec then bugmes();quit;end
25 // IEEE support
26 // See http://bugzilla.scilab.org/show_bug.cgi?id=7296
27 Df       = 1;
28 T        = %inf; // Inf
29 [P,Q]    = cdft("PQ", T, Df);
30 assert_checkequal(P, 1);
31 assert_checkequal(Q, 0);
32 T        = %nan; // NaN
33 [P,Q]    = cdft("PQ", T, Df);
34 assert_checkequal(P, %nan);
35 assert_checkequal(Q, %nan);