dace-rs 0.2.1

Pure Rust implementation of DACE, the Differential Algebra Computational Toolbox
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

/* Print a few DAs through daceWrite for exact-format comparison. */
#include <stdio.h>
#include "dace/dacebase.h"

int main(void)
{
    daceInitialize(3, 2);
    DACEDA x, y, f, z;
    daceAllocateDA(&x, 0); daceAllocateDA(&y, 0); daceAllocateDA(&f, 0); daceAllocateDA(&z, 0);
    daceCreateVariable(&x, 1, 1.0);
    daceCreateVariable(&y, 2, 1.0);
    daceCreateConstant(&f, 1.0);
    daceAddDouble(&f, 0.0, &f);
    DACEDA t1, t2;
    daceAllocateDA(&t1, 0); daceAllocateDA(&t2, 0);
    daceMultiplyDouble(&x, 2.0, &t1);
    daceAdd(&f, &t1, &f);
    daceMultiply(&y, &y, &t2);
    daceMultiplyDouble(&t2, -0.5, &t2);
    daceAdd(&f, &t2, &f);
    unsigned int n = 0;
    char buf[140 * 64];
    daceWrite(&f, buf, &n);
    for (unsigned int i = 0; i < n; i++)
        printf("[%s]\n", &buf[i * 140]);
    daceWrite(&z, buf, &n);
    for (unsigned int i = 0; i < n; i++)
        printf("[%s]\n", &buf[i * 140]);
    return 0;
}