File: incomplete_gamma_inv.cpp

package info (click to toggle)
gcem 1.18.0-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 752 kB
  • sloc: cpp: 5,162; makefile: 279; python: 24
file content (77 lines) | stat: -rw-r--r-- 3,905 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
/*################################################################################
  ##
  ##   Copyright (C) 2016-2024 Keith O'Hara
  ##
  ##   This file is part of the GCE-Math C++ library.
  ##
  ##   Licensed under the Apache License, Version 2.0 (the "License");
  ##   you may not use this file except in compliance with the License.
  ##   You may obtain a copy of the License at
  ##
  ##       http://www.apache.org/licenses/LICENSE-2.0
  ##
  ##   Unless required by applicable law or agreed to in writing, software
  ##   distributed under the License is distributed on an "AS IS" BASIS,
  ##   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  ##   See the License for the specific language governing permissions and
  ##   limitations under the License.
  ##
  ################################################################################*/

#ifndef _WIN32
    #define TEST_ERR_TOL 1e-12
#endif

#define TEST_PRINT_PRECISION_1 3
#define TEST_PRINT_PRECISION_2 18

#include "gcem_tests.hpp"

int main()
{
    print_begin("incomplete_gamma_inv");

    //

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        1.0L,    2.0L, 0.26424111765711527644L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        1.0L,    1.5L, 0.42759329552912134220L);

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        3.0L,    2.0L, 0.80085172652854419439L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        5.0L,    2.0L, 0.95957231800548714595L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        9.0L,    2.0L, 0.99876590195913317327L);

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv, TEST_POSINF,    0.0L,                    1.0L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        0.0L,    2.0L,                    0.0L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,        0.0L,    0.0L,                    0.0L);

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       11.0L,   11.5L, 0.47974821959920432857L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       18.0L,   15.5L, 0.75410627202774938027L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       18.0L,   19.0L, 0.43775501395596266851L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       19.0L,   20.0L, 0.43939261060849132967L);

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       39.0L,   38.0L, 0.58504236243630125536L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       55.0L,   56.0L, 0.46422093592347296598L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,       99.0L,   98.0L, 0.55342727426212001696L);

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      101.0L,  102.0L, 0.47356929040257916830L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      297.0L,  298.0L, 0.48457398488934400049L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      301.0L,  302.0L, 0.48467673854389853316L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      497.0L,  498.0L, 0.48807306199561317772L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      501.0L,  502.0L, 0.48812074555706047585L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      799.0L,  798.0L, 0.51881664512582725823L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      800.0L,  801.0L, 0.49059834200005758564L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,      999.0L,  997.0L, 0.52943655952003709775L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv,     1001.0L, 1005.0L, 0.45389544705967349580L);

    //

    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv, TEST_NAN, TEST_NAN, TEST_NAN);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv, TEST_NAN, TEST_NAN,     0.5L);
    GCEM_TEST_EXPECTED_VAL(gcem::incomplete_gamma_inv, TEST_NAN,     2.0L, TEST_NAN);

    //

    print_final("incomplete_gamma_inv");

    return 0;
}