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