/[dtapublic]/projs/trunk/shared_source/c_datd/gmp_rats.h
ViewVC logotype

Contents of /projs/trunk/shared_source/c_datd/gmp_rats.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 56 - (show annotations) (download)
Sat Oct 29 01:53:01 2016 UTC (7 years, 5 months ago) by dashley
File MIME type: text/plain
File size: 10215 byte(s)
License and property (keyword) changes.
1 //$Header$
2 //-------------------------------------------------------------------------------------------------
3 //This file is part of "David T. Ashley's Shared Source Code", a set of shared components
4 //integrated into many of David T. Ashley's projects.
5 //-------------------------------------------------------------------------------------------------
6 //This source code and any program in which it is compiled/used is provided under the MIT License,
7 //reproduced below.
8 //-------------------------------------------------------------------------------------------------
9 //Permission is hereby granted, free of charge, to any person obtaining a copy of
10 //this software and associated documentation files(the "Software"), to deal in the
11 //Software without restriction, including without limitation the rights to use,
12 //copy, modify, merge, publish, distribute, sublicense, and / or sell copies of the
13 //Software, and to permit persons to whom the Software is furnished to do so,
14 //subject to the following conditions :
15 //
16 //The above copyright notice and this permission notice shall be included in all
17 //copies or substantial portions of the Software.
18 //
19 //THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
20 //IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
21 //FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.IN NO EVENT SHALL THE
22 //AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
23 //LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
24 //OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
25 //SOFTWARE.
26 //-------------------------------------------------------------------------------------------------
27 #ifndef GMP_RATS_H_INCLUDED
28 #define GMP_RATS_H_INCLUDED
29
30 #ifdef MODULE_GMP_RATS
31 #define DECMOD_GMP_RATS
32 #else
33 #define DECMOD_GMP_RATS extern
34 #endif
35
36 //The data structure of a rational number. Generally,
37 //it isn't necessary to keep error flags for the rational
38 //number--instead, we can just use the NAN flags of integers,
39 //with the convention that if either of the integers
40 //has any of the NAN flags set, then the rational number
41 //is NAN. This will be good enough. So let's agree
42 //by convention that a rational number is either
43 //valid or NAN--no other possibilities, and nothing more
44 //specific.
45 //
46 //Generally, none of the algorithms in this module make
47 //assumptions about whether the inputs are normalized,
48 //although they may make assurances about whether the
49 //outputs are normalized. Anything except a denominator
50 //of zero is considered acceptable for inputs.
51 //
52 //For the record, "normalized" according to my definition is:
53 // a)Zero has the unique representation 0/1.
54 // b)Numerator and denominator are coprime.
55 // c)Positive rational numbers have positive
56 // components.
57 // d)Negative rational numbers have numerator
58 // negative and denominator positive.
59 typedef struct
60 {
61 GMP_INTS_mpz_struct num;
62 //The numerator.
63 GMP_INTS_mpz_struct den;
64 //The denominator.
65 } GMP_RATS_mpq_struct;
66
67 /******************************************************************/
68 /*** STATUS FUNCTIONS *******************************************/
69 /******************************************************************/
70 //Functions in this category provide information about rational
71 //numbers.
72 //
73 //Returns TRUE if the rational number is NAN, which means
74 //either integer component with flags set or a denominator
75 //of zero.
76 DECMOD_GMP_RATS
77 int GMP_RATS_mpq_is_nan(const GMP_RATS_mpq_struct *rn);
78
79 /******************************************************************/
80 /*** INITIALIZATION, CLEARING, AND SETTING FUNCTIONS ************/
81 /******************************************************************/
82 //Allocates an space for integers, initializes a rational number
83 //to 0/1.
84 DECMOD_GMP_RATS
85 void GMP_RATS_mpq_init(GMP_RATS_mpq_struct *arg);
86 //Destroys a rational number.
87 DECMOD_GMP_RATS
88 void GMP_RATS_mpq_clear(GMP_RATS_mpq_struct *arg);
89 //Sets a rational number to the two integer components
90 //specified. It is not required that these components
91 //be normalized.
92 DECMOD_GMP_RATS
93 void GMP_RATS_mpq_set_si(GMP_RATS_mpq_struct *arg,
94 int num,
95 int den);
96 DECMOD_GMP_RATS
97 void GMP_RATS_mpq_copy( GMP_RATS_mpq_struct *dst,
98 const GMP_RATS_mpq_struct *src);
99
100 //Swaps two rational numbers.
101 DECMOD_GMP_RATS
102 void GMP_RATS_mpq_swap( GMP_RATS_mpq_struct *a,
103 GMP_RATS_mpq_struct *b);
104
105 //Swaps numerator and denominator. Simple swap, no normalization.
106 DECMOD_GMP_RATS
107 void GMP_RATS_mpq_swap_components(GMP_RATS_mpq_struct *arg);
108
109 //Sets a rational number to be a slash-separated integer. The
110 //two components (before and after the slash) may be integers
111 //in scientific notation, integers with commas, etc. The
112 //*failure flag is set TRUE if the number could not be parsed.
113 //The rational number involved must already have been initialized,
114 //so space is allocated. The number is not normalized in any
115 //way--the integer components of the rational number are
116 //exactly the slash-separated components. It is possible to
117 //use this function to simply parse to see if one can form
118 //a rational number, but this is an expensive process because
119 //of the allocation and then destruction involved. It seems
120 //possible in the future that one would introduce a function
121 //which won't take the scientifc notation, and so the word
122 //"complex" is included below. Any subsequent function which
123 //won't accept the scientific notation should be called "simple"
124 //instead.
125 DECMOD_GMP_RATS
126 void GMP_RATS_mpq_set_complex_slash_sepd_rat_num(const char *s,
127 int *failure,
128 GMP_RATS_mpq_struct *rn);
129
130 //Attempts to parse and set a rational number expressed in
131 //pure scientific notation.
132 DECMOD_GMP_RATS
133 void GMP_RATS_mpq_set_sci_not_rat_num(const char *s,
134 int *failure,
135 GMP_RATS_mpq_struct *rn);
136
137 //Attempts to parse and set a rational number which might be in
138 //any accepted format.
139 DECMOD_GMP_RATS
140 void GMP_RATS_mpq_set_all_format_rat_num(const char *s,
141 int *failure,
142 GMP_RATS_mpq_struct *rn);
143
144 /******************************************************************/
145 /*** NORMALIZATION FUNCTIONS ************************************/
146 /******************************************************************/
147 //Normalizes only the signs. Non-negative number expressed as +/+,
148 //negative numbers expressed only as -/+.
149 DECMOD_GMP_RATS
150 void
151 GMP_RATS_mpq_normalize_sign(GMP_RATS_mpq_struct *rn);
152
153
154 DECMOD_GMP_RATS
155 void GMP_RATS_mpq_normalize(GMP_RATS_mpq_struct *rn);
156
157 /******************************************************************/
158 /*** ARITHMETIC FUNCTIONS ***************************************/
159 /******************************************************************/
160 //Adds two rational numbers to produce normalized result.
161 //Args may be same as result. Any NAN inputs result in NAN
162 //output: NAN flagged using 1/0 as return value.
163 DECMOD_GMP_RATS
164 void GMP_RATS_mpq_add( GMP_RATS_mpq_struct *result,
165 const GMP_RATS_mpq_struct *arg1,
166 const GMP_RATS_mpq_struct *arg2);
167 //Subtract, follows the same rules as add.
168 DECMOD_GMP_RATS
169 void GMP_RATS_mpq_sub( GMP_RATS_mpq_struct *result,
170 const GMP_RATS_mpq_struct *arg1,
171 const GMP_RATS_mpq_struct *arg2);
172 //Multiply.
173 DECMOD_GMP_RATS
174 void GMP_RATS_mpq_mul( GMP_RATS_mpq_struct *result,
175 const GMP_RATS_mpq_struct *arg1,
176 const GMP_RATS_mpq_struct *arg2);
177 //Divide.
178 DECMOD_GMP_RATS
179 void GMP_RATS_mpq_div( GMP_RATS_mpq_struct *result,
180 const GMP_RATS_mpq_struct *arg1,
181 const GMP_RATS_mpq_struct *arg2);
182
183 /******************************************************************/
184 /*** COMPARISON FUNCTIONS ***************************************/
185 /******************************************************************/
186 //Compares two rational numbers. The return value is neg if
187 //arg1 < arg2, 0 if arg1 == arg2, and 1 if arg1 > arg2. The
188 //rational numbers passed are not required to be normalized,
189 //either with respect to sign of components or with respect
190 //to coprimality. There are some errors which may occur, either
191 //because a denominator is zero or because intermediate results
192 //overflow as the comparison is made. For that reason, the
193 //"failure" pointer is allowed. If it is NULL, no failure information
194 //is returned and any failure condition results in a return value
195 //of zero and the failure can't be detected by the caller.
196 //If the pointer is not NULL, the failure boolean will be
197 //filled in so the caller can detect it.
198 int GMP_RATS_mpq_cmp(const GMP_RATS_mpq_struct *arg1,
199 const GMP_RATS_mpq_struct *arg2,
200 int *failure);
201
202 /******************************************************************/
203 /*** VERSION CONTROL REPORTING FUNCTIONS ************************/
204 /******************************************************************/
205 DECMOD_GMP_RATS const char *GMP_RATS_cvcinfo(void);
206 DECMOD_GMP_RATS const char *GMP_RATS_hvcinfo(void);
207 #define GMP_RATS_H_VERSION ("$Header$")
208 #endif
209
210 //End of gmp_rats.h.

Properties

Name Value
svn:keywords Header

dashley@gmail.com
ViewVC Help
Powered by ViewVC 1.1.25