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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
|
// file kernel/x/h/div.h: division of extensible integers
/*-----------------------------------------------------------------------+
| Copyright 2005-2006, Michel Quercia (michel.quercia@prepas.org) |
| |
| This file is part of Numerix. Numerix is free software; you can |
| redistribute it and/or modify it under the terms of the GNU Lesser |
| General Public License as published by the Free Software Foundation; |
| either version 2.1 of the License, or (at your option) any later |
| version. |
| |
| The Numerix Library is distributed in the hope that it will be |
| useful, but WITHOUT ANY WARRANTY; without even the implied warranty |
| of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with the GNU MP Library; see the file COPYING. If not, |
| write to the Free Software Foundation, Inc., 59 Temple Place - |
| Suite 330, Boston, MA 02111-1307, USA. |
+-----------------------------------------------------------------------+
| |
| Division |
| |
+-----------------------------------------------------------------------*/
/* +---------------------+
| Division gnrale |
+---------------------+ */
/*
entre :
a,b = entiers extensibles
_c = NULL ou adresse d'un entier extensible
_d = NULL ou adresse d'un entier extensible
mode = long
contraintes :
b != 0
si _c et _d ne sont pas gaux NULL alors ils sont distincts
sortie :
si mode & 3 = 0 : c <- floor(a/b), d <- a - b*c (division tronque )
si mode & 3 = 1 : c <- floor(a/b+1/2), d <- a - b*c (division centre up )
si mode & 3 = 2 : c <- ceil(a/b), d <- a - b*c (division majore )
si mode & 3 = 3 : c <- ceil(a/b-1/2), d <- a - b*c (division centre down)
si _c != NULL : *_c <- c
si _d != NULL : *_d <- d
valeur de retour : c_api ml_api
si mode & 12 = 0 : NULL unit
si mode & 12 = 4 : c c
si mode & 12 = 8 : d d
si mode & 12 = 12 : NULL (c,d)
erreur :
ZERO_DIVISOR si b = 0
MULTIPLE_RESULT si _c == _d != NULL
*/
#if defined(caml_api) || defined(ocaml_api)
value xx(private_quomod)(xint *_c, xint *_d, xint a, xint b, long mode);
value xx(hquomod) (value mode, xint *_c, xint *_d, xint a, xint b);
value xx(hquo) (value mode, xint *_c, xint a, xint b);
value xx(hmod) (value mode, xint *_d, xint a, xint b);
value xx(f_hquomod)(value mode, xint a, xint b);
value xx(f_hquo) (value mode, xint a, xint b);
value xx(f_hmod) (value mode, xint a, xint b);
value xx(quomod) ( xint *_c, xint *_d, xint a, xint b);
value xx(quo) ( xint *_c, xint a, xint b);
value xx(mod) ( xint *_d, xint a, xint b);
value xx(f_quomod) ( xint a, xint b);
value xx(f_quo) ( xint a, xint b);
value xx(f_mod) ( xint a, xint b);
#elif defined(c_api)
xint xx(private_quomod)(xint *_c, xint *_d, xint a, xint b, long mode);
extern inline void xx(gquomod)(xint *_c, xint *_d, xint a, xint b, long mode) { xx(private_quomod)(_c, _d, a,b, 0|(mode & 3));}
extern inline xint xx(gquo) (xint *_c, xint a, xint b, long mode) {return xx(private_quomod)(_c, NULL, a,b, 4|(mode & 3));}
extern inline xint xx(gmod) ( xint *_d, xint a, xint b, long mode) {return xx(private_quomod)(NULL,_d, a,b, 8|(mode & 3));}
extern inline xint xx(f_gquo) ( xint a, xint b, long mode) {return xx(private_quomod)(NULL,NULL, a,b, 4|(mode & 3));}
extern inline xint xx(f_gmod) ( xint a, xint b, long mode) {return xx(private_quomod)(NULL,NULL, a,b, 8|(mode & 3));}
extern inline void xx(quomod) (xint *_c, xint *_d, xint a, xint b ) { xx(private_quomod)(_c, _d, a,b, 0|0);}
extern inline xint xx(quo) (xint *_c, xint a, xint b ) {return xx(private_quomod)(_c, NULL, a,b, 4|0);}
extern inline xint xx(mod) ( xint *_d, xint a, xint b ) {return xx(private_quomod)(NULL,_d, a,b, 8|0);}
extern inline xint xx(f_quo) ( xint a, xint b ) {return xx(private_quomod)(NULL,NULL, a,b, 4|0);}
extern inline xint xx(f_mod) ( xint a, xint b ) {return xx(private_quomod)(NULL,NULL, a,b, 8|0);}
#endif /* api */
/* +--------------------------+
| Division par un entier |
+--------------------------+ */
/*
entre :
a = entier extensible
b = long ou Caml/Ocaml int
_c = NULL ou adresse d'un entier extensible
mode = long
sortie :
si mode & 3 = 0 : c <- floor(a/b), d <- a - b*c (division tronque)
si mode & 3 = 1 : c <- floor(a/b+1/2), d <- a - b*c (division centre up)
si mode & 3 = 2 : c <- ceil(a/b), d <- a - b*c (division majore)
si mode & 3 = 3 : c <- ceil(a/b-1/2), d <- a - b*c (division centre down)
si _c != NULL : *_c <- c
valeur de retour : c_api ml_api
si mode & 12 = 0 : d unit
si mode & 12 = 4 : c c
si mode & 12 = 8 : d d
si mode & 12 = 12 : c (c,d)
erreur :
ZERO_DIVISOR si b = 0
*/
#if defined(caml_api) || defined(ocaml_api)
value xx(private_quomod_1)(xint *_c, xint a, long b, long mode);
value xx(gquomod_1) (value mode, xint *_c, xint a, long b);
value xx(gquo_1) (value mode, xint *_c, xint a, long b);
value xx(f_gquomod_1)(value mode, xint a, long b);
value xx(f_gquo_1) (value mode, xint a, long b);
value xx(f_gmod_1) (value mode, xint a, long b);
value xx(quomod_1) ( xint *_c, xint a, long b);
value xx(quo_1) ( xint *_c, xint a, long b);
value xx(f_quomod_1) ( xint a, long b);
value xx(f_quo_1) ( xint a, long b);
value xx(f_mod_1) ( xint a, long b);
#elif defined(c_api)
long xx(private_quomod_1)(xint *_c, xint a, long b, long mode);
extern inline long xx(gquomod_1)(xint *_c, xint a, long b, long mode) {return xx(private_quomod_1)(_c, a,b, 8|(mode & 3));}
extern inline xint xx(gquo_1) (xint *_c, xint a, long b, long mode) {return (xint)xx(private_quomod_1)(_c, a,b, 4|(mode & 3));}
extern inline long xx(gmod_1) ( xint a, long b, long mode) {return xx(private_quomod_1)(NULL, a,b, 8|(mode & 3));}
extern inline xint xx(f_gquo_1) ( xint a, long b, long mode) {return (xint)xx(private_quomod_1)(NULL, a,b, 4|(mode & 3));}
extern inline long xx(f_gmod_1) ( xint a, long b, long mode) {return xx(private_quomod_1)(NULL, a,b, 8|(mode & 3));}
extern inline long xx(quomod_1) (xint *_c, xint a, long b ) {return xx(private_quomod_1)(_c, a,b, 8|0);}
extern inline xint xx(quo_1) (xint *_c, xint a, long b ) {return (xint)xx(private_quomod_1)(_c, a,b, 4|0);}
extern inline long xx(mod_1) ( xint a, long b ) {return xx(private_quomod_1)(NULL, a,b, 8|0);}
extern inline xint xx(f_quo_1) ( xint a, long b ) {return (xint)xx(private_quomod_1)(NULL, a,b, 4|0);}
extern inline long xx(f_mod_1) ( xint a, long b ) {return xx(private_quomod_1)(NULL, a,b, 8|0);}
#endif /* api */
|
