liba 0.1.15

An algorithm library based on C/C++
Documentation 
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

#define MAIN(x) mf##x
#include "test.h"
#include "a/mf.h"
#include <string.h>

int main(int argc, char *argv[]) // NOLINT(misc-definitions-in-headers)
{
    main_init(argc, argv, 1);
    unsigned int e = A_MF_NUL;
    if (argc > 1)
    {
        if (strstr(argv[1], "gauss2.")) { e = A_MF_GAUSS2; }
        else if (strstr(argv[1], "gauss.")) { e = A_MF_GAUSS; }
        else if (strstr(argv[1], "gbell.")) { e = A_MF_GBELL; }
        else if (strstr(argv[1], "dsig.")) { e = A_MF_DSIG; }
        else if (strstr(argv[1], "psig.")) { e = A_MF_PSIG; }
        else if (strstr(argv[1], "sig.")) { e = A_MF_SIG; }
        else if (strstr(argv[1], "trap.")) { e = A_MF_TRAP; }
        else if (strstr(argv[1], "tri.")) { e = A_MF_TRI; }
        else if (strstr(argv[1], "lins.")) { e = A_MF_LINS; }
        else if (strstr(argv[1], "linz.")) { e = A_MF_LINZ; }
        else if (strstr(argv[1], "s.")) { e = A_MF_S; }
        else if (strstr(argv[1], "z.")) { e = A_MF_Z; }
        else if (strstr(argv[1], "pi.")) { e = A_MF_PI; }
    }
    a_float params[6] = {0, 0, 0, 0, 0, 0};
    for (int i = 2; i < argc && i < 8; ++i)
    {
        params[i - 2] = strtonum(argv[i], A_NULL);
    }
    if (e != A_MF_NUL)
    {
        a_float delta = (params[1] - params[0]) / A_FLOAT_C(100.0);
        for (int i = 0; i < 100; ++i)
        {
            a_float x = params[0] + delta * a_float_c(i);
            a_float y = a_mf(e, x, params + 2);
            debug(A_FLOAT_PRI(, "g,") A_FLOAT_PRI(, "g\n"), x, y);
        }
    }
    else
    {
        TEST_BUG(a_mf_gauss(0, 1, 0) >= 1);
        TEST_BUG(a_mf_gauss2(0, 1, -1, +1, 1) >= 1);
        TEST_BUG(a_mf_gbell(0, 2, 4, 0) >= 1);
        TEST_BUG(a_mf_sig(0, 2, 0) >= 0.5);
        TEST_BUG(a_mf_dsig(0, 5, -2, +5, +2) > 0.99);
        TEST_BUG(a_mf_psig(0, 5, -2, -5, +2) > 0.99);
        TEST_BUG(a_mf_trap(0, -2, -1, +1, 2) >= 1);
        TEST_BUG(a_mf_tri(0, -1, 0, +1) >= 1);
        TEST_BUG(a_mf_lins(0, -1, 0) >= 1);
        TEST_BUG(a_mf_linz(0, 0, +1) >= 1);
        TEST_BUG(a_mf_s(0, -1, 0) >= 1);
        TEST_BUG(a_mf_z(0, 0, +1) >= 1);
        TEST_BUG(a_mf_pi(0, -2, -1, +1, +2) >= 1);
    }
    return 0;
}