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
|
/*
*
* Operations sur les scalaires ( rationnels )
*
*/
#include "scalaires.h"
#include <math.h>
#define PRECISION 1E6
static int pgcd(int a, int b)
{
if (!a)
return b;
if (a < b)
return pgcd(b, a);
return pgcd(b, a % b);
}
rationnel rat_constr_zero(void)
{ /* renvoie 0 */
rationnel temp;
temp.num = 0;
temp.denom = 1;
return temp;
}
rationnel rat_constr(int num, int denom)
{ /* cree une fraction */
rationnel temp;
if (denom < 0) {
denom = -denom;
num = -num;
}
temp.num = num / pgcd(num, denom);
temp.denom = denom / pgcd(num, denom);
return temp;
}
rationnel rat_constr_from_double(double flt)
{ /* cree une fraction a partir d un double */
return rat_constr(floor(flt * PRECISION), PRECISION);
}
void rat_destruct(rationnel rat)
{ /* destructeur */
}
double rat_to_double(rationnel rat)
{ /* obtention du double correspondant a un rationnel */
return ((double) rat.num / (double) rat.denom);
}
rationnel rat_addition(rationnel rat1, rationnel rat2)
{ /* addition */
return rat_constr(rat1.num * rat2.denom + rat2.num * rat1.denom, rat1.denom * rat2.denom);
}
rationnel rat_soustraction(rationnel rat1, rationnel rat2)
{ /* soustraction */
return rat_constr(rat1.num * rat2.denom - rat2.num * rat1.denom, rat1.denom * rat2.denom);
}
rationnel rat_moinsunaire(rationnel rat1)
{ /* moins unaire */
return rat_constr(-rat1.num, rat1.denom);
}
rationnel rat_multiplication(rationnel rat1, rationnel rat2)
{ /* multiplication */
return rat_constr(rat1.num * rat2.num, rat1.denom * rat2.denom);
}
rationnel rat_division(rationnel rat1, rationnel rat2)
{ /* division */
return rat_constr(rat1.num * rat2.denom, rat1.denom * rat2.num);
}
|
