librobotcontrol-sys 0.4.0

Rust port of librobotcontrol
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109

/**
 * @example    rc_test_polynomial.c
 *
 * @brief      Tests the functions in rc_polynomial.h
 *
 * @author     James Strawson
 * @date       1/29/2018
 */

#include <stdio.h>
#include <rc/math.h>

#define LEN 4 // length of polynomial to test

int main()
{
	int i;
	rc_vector_t a = RC_VECTOR_INITIALIZER;
	rc_vector_t b = RC_VECTOR_INITIALIZER;
	rc_vector_t c = RC_VECTOR_INITIALIZER;
	rc_vector_t d = RC_VECTOR_INITIALIZER;
	rc_vector_t e = RC_VECTOR_INITIALIZER;

	printf("\nRandom polynomials of orders 0-9\n");
	for(i=1;i<=10;i++){
		rc_vector_random(&a,i);
		rc_poly_print(a);
	}

	printf("\npolynomial of ones a:\n");
	rc_vector_ones(&a, 2);
	rc_poly_print(a);

	printf("polynomial of ones b:\n");
	rc_vector_ones(&b, 3);
	rc_poly_print(b);

	printf("Polynomial convolution a*b:\n");
	rc_poly_conv(a,b,&c);
	rc_poly_print(c);

	// test powers
	printf("\na to the power of 0-5:\n");
	for(i=0;i<=5;i++){
		rc_poly_power(a,i,&c);
		rc_poly_print(c);
	}

	// test subtract
	printf("\nrc_poly_subtract c=b-a\n");
	rc_poly_subtract(b,a,&c);
	rc_poly_print(c);

	// test subtract
	printf("\nrc_poly_subtract_inplace c=c-a\n");
	rc_poly_subtract_inplace(&c,a);
	rc_poly_print(c);

	// test add
	printf("\nrc_poly_add a+b=c\n");
	rc_poly_add(a,b,&c);
	rc_poly_print(c);

	// test add
	printf("\nrc_poly_add_inplace c=c+a\n");
	rc_poly_add_inplace(&c,a);
	rc_poly_print(c);

	// differentiate
	printf("0st-3rd derivative of b\n");
	for(i=0;i<=3;i++){
		rc_poly_differentiate(b,i,&c);
		rc_poly_print(c);
	}

	// divide
	rc_vector_fibonnaci(&a, 2);
	rc_vector_fibonnaci(&b, 5);
	printf("\nNew polynomial a:\n");
	rc_poly_print(a);
	printf("\nNew polynomial b:\n");
	rc_poly_print(b);
	printf("\ndivide a into b\n");
	rc_poly_divide(b,a,&c,&d);
	rc_poly_print(c);
	printf("remainder:\n");
	rc_poly_print(d);

	// confirm
	printf("\nConfirm by multiplying the result and adding remainder\n");
	rc_poly_conv(c,a,&e);
	rc_poly_add_inplace(&e,d);
	rc_poly_print(e);

	// butterworth
	printf("\nfirst 4 butterworth polynomials\n");
	for(i=1;i<=4;i++){
		rc_poly_butter(i,1,&a);
		rc_poly_print(a);
	}

	rc_vector_free(&a);
	rc_vector_free(&b);
	rc_vector_free(&c);
	rc_vector_free(&d);
	rc_vector_free(&e);
	printf("\nDONE\n");
	return 0;
}