/[gsl]/gsl/cdf/test.c
ViewVC logotype

Diff of /gsl/cdf/test.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.2.12 by bjg, Mon Jun 23 18:19:31 2003 UTC revision 1.1.2.13 by bjg, Tue Jun 24 14:13:15 2003 UTC
# Line 243  main (void) Line 243  main (void)
243    
244    /* Tests for F distribution */    /* Tests for F distribution */
245    
246      /* p(x, nu1, nu2) := betaI(1 / (1 + nu2 / nu1 * x), nu1 / 2, nu2 / 2) */
247    
248    TEST(gsl_cdf_Fdist_P, (0.0, 1.2, 1.3), 0.0, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (0.0, 1.2, 1.3), 0.0, TEST_TOL6);
249    TEST(gsl_cdf_Fdist_P, (1e-100, 1.2, 1.3), 6.98194275525039002e-61, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (1e-100, 1.2, 1.3), 6.98194275525039002e-61, TEST_TOL6);
250    TEST(gsl_cdf_Fdist_P, (0.001, 1.2, 1.3), 1.10608485860238564e-2, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (0.001, 1.2, 1.3), 1.10608485860238564e-2, TEST_TOL6);
# Line 289  main (void) Line 291  main (void)
291    TEST(gsl_cdf_Fdist_P, (1000.0, 1.2, 500), 1.0, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (1000.0, 1.2, 500), 1.0, TEST_TOL6);
292    TEST(gsl_cdf_Fdist_P, (10000.0, 1.2, 500), 1.0, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (10000.0, 1.2, 500), 1.0, TEST_TOL6);
293    
 /* [(p(x, nu1, nu2) := betaI(1 / (1 + nu2 / nu1 * x), nu1 / 2, nu2 / 2))] */  
294    TEST(gsl_cdf_Fdist_P, (0.0, 200, 500), 0.0, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (0.0, 200, 500), 0.0, TEST_TOL6);
295    TEST(gsl_cdf_Fdist_P, (1e-100, 200, 500), 0.0, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (1e-100, 200, 500), 0.0, TEST_TOL6);
296    TEST(gsl_cdf_Fdist_P, (0.001, 200, 500), 4.09325080403669893e-251, TEST_TOL6);    TEST(gsl_cdf_Fdist_P, (0.001, 200, 500), 4.09325080403669893e-251, TEST_TOL6);
# Line 307  main (void) Line 308  main (void)
308    
309    /* Tests for gamma distribution */    /* Tests for gamma distribution */
310    
311      /* p(x, a, b) := gammaP(b, x / a) */
312    
313    TEST(gsl_cdf_gamma_P, (0.0, 1, 1), 0.0, TEST_TOL6);    TEST(gsl_cdf_gamma_P, (0.0, 1, 1), 0.0, TEST_TOL6);
314    TEST(gsl_cdf_gamma_P, (1e-100, 1, 1), 1e-100, TEST_TOL6);    TEST(gsl_cdf_gamma_P, (1e-100, 1, 1), 1e-100, TEST_TOL6);
315    TEST(gsl_cdf_gamma_P, (0.001, 1, 1), 9.99500166625008332e-4 , TEST_TOL6);    TEST(gsl_cdf_gamma_P, (0.001, 1, 1), 9.99500166625008332e-4 , TEST_TOL6);
# Line 322  main (void) Line 325  main (void)
325    TEST(gsl_cdf_gamma_P, (1000.0, 1, 1), 1e0, TEST_TOL6);    TEST(gsl_cdf_gamma_P, (1000.0, 1, 1), 1e0, TEST_TOL6);
326    TEST(gsl_cdf_gamma_P, (10000.0, 1, 1), 1e0, TEST_TOL6);    TEST(gsl_cdf_gamma_P, (10000.0, 1, 1), 1e0, TEST_TOL6);
327    
328      TEST(gsl_cdf_gamma_P, (0.0, 10, 1), 0.0, TEST_TOL6);
329      TEST(gsl_cdf_gamma_P, (1e-100, 10, 1), 1e-101, TEST_TOL6);
330      TEST(gsl_cdf_gamma_P, (0.001, 10, 1), 9.99950001666625001e-5 , TEST_TOL6);
331      TEST(gsl_cdf_gamma_P, (0.01, 10, 1), 9.99500166625008332e-4, TEST_TOL6);
332      TEST(gsl_cdf_gamma_P, (0.1, 10, 1), 9.95016625083194643e-3, TEST_TOL6);
333      TEST(gsl_cdf_gamma_P, (0.325, 10, 1), 3.19775501686939529e-2, TEST_TOL6);
334      TEST(gsl_cdf_gamma_P, (1.0, 10, 1), 9.51625819640404268e-2, TEST_TOL6);
335      TEST(gsl_cdf_gamma_P, (1.5, 10, 1), 1.39292023574942193e-1, TEST_TOL6);
336      TEST(gsl_cdf_gamma_P, (2.0, 10, 1), 1.81269246922018141e-1, TEST_TOL6);
337      TEST(gsl_cdf_gamma_P, (10.0, 10, 1), 6.32120558828557678e-1, TEST_TOL6);
338      TEST(gsl_cdf_gamma_P, (20.0, 10, 1), 8.64664716763387308e-1, TEST_TOL6);
339      TEST(gsl_cdf_gamma_P, (100.0, 10, 1), 9.99954600070237515e-1, TEST_TOL6);
340      TEST(gsl_cdf_gamma_P, (1000.0, 10, 1), 1e0, TEST_TOL6);
341      TEST(gsl_cdf_gamma_P, (10000.0, 10, 1), 1e0, TEST_TOL6);
342    
343      TEST(gsl_cdf_gamma_P, (0.0, 10, 17), 0e0, TEST_TOL6);
344      TEST(gsl_cdf_gamma_P, (1e-100, 10, 17), 0e0, TEST_TOL6);
345      TEST(gsl_cdf_gamma_P, (0.001, 10, 17), 2.81119174040422844e-83, TEST_TOL6);
346      TEST(gsl_cdf_gamma_P, (0.01, 10, 17), 2.80880324651985887e-66, TEST_TOL6);
347      TEST(gsl_cdf_gamma_P, (0.1, 10, 17), 2.78502998087492130e-49, TEST_TOL6);
348      TEST(gsl_cdf_gamma_P, (0.325, 10, 17), 1.37283653245125844e-40, TEST_TOL6);
349      TEST(gsl_cdf_gamma_P, (1.0, 10, 17), 2.55811932292544243e-32, TEST_TOL6);
350      TEST(gsl_cdf_gamma_P, (1.5, 10, 17), 2.40420441175422372e-29, TEST_TOL6);
351      TEST(gsl_cdf_gamma_P, (2.0, 10, 17), 3.05092926217898577e-27, TEST_TOL6);
352      TEST(gsl_cdf_gamma_P, (10.0, 10, 17), 1.094920130378183e-15, TEST_TOL6);
353      TEST(gsl_cdf_gamma_P, (20.0, 10, 17), 5.60605096173161688e-11, TEST_TOL6);
354      TEST(gsl_cdf_gamma_P, (100.0, 10, 17), 2.70416097848011280e-2, TEST_TOL6);
355      TEST(gsl_cdf_gamma_P, (1000.0, 10, 17), 1.000000000000000000e0, TEST_TOL6);
356      TEST(gsl_cdf_gamma_P, (10000.0, 10, 17), 1.000000000000000000e0, TEST_TOL6);
357    
358    exit (gsl_test_summary());    exit (gsl_test_summary());
359  }  }

Legend:
Removed from v.1.1.2.12  
changed lines
  Added in v.1.1.2.13

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26