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

#define MAIN(x) poly##x
#include "test.h"
#include "a/poly.h"

static void print_(a_float const *a, a_size n)
{
    putchar('{');
    a_forenum(a_size, i, n)
    {
        printf(A_FLOAT_PRI("", "g"), a[i]);
        if (i < n - 1) { putchar(','); }
    }
    putchar('}');
}

int main(int argc, char *argv[]) // NOLINT(misc-definitions-in-headers)
{
    (void)argc;
    (void)argv;
    a_float x = 2;
    a_float x4[] = {A_FLOAT_C(1.0), A_FLOAT_C(2.0), A_FLOAT_C(3.0), A_FLOAT_C(4.0)};
    print_(x4, 4);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_eval(x4, 4, x));
    a_poly_swap(x4, 4);
    print_(x4, 4);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_eval(x4, 4, x));
    print_(x4, 4);
    a_poly_swap(x4, 4);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_evar(x4, 4, x));
    print_(x4, 3);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_eval(x4, 3, x));
    a_poly_swap(x4, 3);
    print_(x4, 3);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_eval(x4, 3, x));
    print_(x4, 3);
    printf("=" A_FLOAT_PRI("", "g") "\n", a_poly_evar(x4, 3, x));
    a_poly_swap(x4, 3);
    return 0;
}