source: src/libcfa/gmp @ 2893f6d

aaron-thesisarm-ehcleanup-dtorsdeferred_resndemanglerjacob/cs343-translationjenkins-sandboxnew-astnew-ast-unique-exprnew-envno_listpersistent-indexerresolv-newwith_gc
Last change on this file since 2893f6d was 2893f6d, checked in by Peter A. Buhr <pabuhr@…>, 4 years ago

make gmp.h a CFA standard-header file

  • Property mode set to 100644
File size: 16.4 KB
Line 
1//
2// Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
3//
4// The contents of this file are covered under the licence agreement in the
5// file "LICENCE" distributed with Cforall.
6//
7// gmp --
8//
9// Author           : Peter A. Buhr
10// Created On       : Tue Apr 19 08:43:43 2016
11// Last Modified By : Peter A. Buhr
12// Last Modified On : Sun May 14 23:47:36 2017
13// Update Count     : 9
14//
15
16// https://gmplib.org/gmp-man-6.1.1.pdf
17
18#include <gmp.h>                                                                                // GNU multi-precise integers
19#include <fstream>                                                                              // sout
20
21struct Int { mpz_t mpz; };                                                              // wrap GMP implementation
22
23// constructor
24void ?{}( Int * this ) { mpz_init( this->mpz ); }
25void ?{}( Int * this, Int init ) { mpz_init_set( this->mpz, init.mpz ); }
26void ?{}( Int * this, zero_t ) { mpz_init_set_si( this->mpz, 0 ); }
27void ?{}( Int * this, one_t ) { mpz_init_set_si( this->mpz, 1 ); }
28void ?{}( Int * this, signed long int init ) { mpz_init_set_si( this->mpz, init ); }
29void ?{}( Int * this, unsigned long int init ) { mpz_init_set_ui( this->mpz, init ); }
30void ?{}( Int * this, const char * val ) { if ( mpz_init_set_str( this->mpz, val, 0 ) ) abort(); }
31void ^?{}( Int * this ) { mpz_clear( this->mpz ); }
32
33// assignment
34Int ?=?( Int * lhs, Int rhs ) { mpz_set( lhs->mpz, rhs.mpz ); return *lhs; }
35Int ?=?( Int * lhs, long int rhs ) { mpz_set_si( lhs->mpz, rhs ); return *lhs; }
36Int ?=?( Int * lhs, unsigned long int rhs ) { mpz_set_ui( lhs->mpz, rhs ); return *lhs; }
37//Int ?=?( Int * lhs, const char * rhs ) { if ( mpq_set_str( lhs->mpz, rhs, 0 ) ) abort(); return *lhs; }
38
39char ?=?( char * lhs, Int rhs ) { char val = mpz_get_si( rhs.mpz ); *lhs = val; return val; }
40short int ?=?( short int * lhs, Int rhs ) { short int val = mpz_get_si( rhs.mpz ); *lhs = val; return val; }
41int ?=?( int * lhs, Int rhs ) { int val = mpz_get_si( rhs.mpz ); *lhs = val; return val; }
42long int ?=?( long int * lhs, Int rhs ) { long int val = mpz_get_si( rhs.mpz ); *lhs = val; return val; }
43unsigned char ?=?( unsigned char * lhs, Int rhs ) { unsigned char val = mpz_get_ui( rhs.mpz ); *lhs = val; return val; }
44unsigned short int ?=?( unsigned short int * lhs, Int rhs ) { unsigned short int val = mpz_get_ui( rhs.mpz ); *lhs = val; return val; }
45unsigned int ?=?( unsigned int * lhs, Int rhs ) { unsigned int val = mpz_get_ui( rhs.mpz ); *lhs = val; return val; }
46unsigned long int ?=?( unsigned long int * lhs, Int rhs ) { unsigned long int val = mpz_get_ui( rhs.mpz ); *lhs = val; return val; }
47
48// conversions
49long int narrow( Int val ) { return mpz_get_si( val.mpz ); }
50unsigned long int narrow( Int val ) { return mpz_get_ui( val.mpz ); }
51
52// comparison
53int ?==?( Int oper1, Int oper2 ) { return mpz_cmp( oper1.mpz, oper2.mpz ) == 0; }
54int ?==?( Int oper1, long int oper2 ) { return mpz_cmp_si( oper1.mpz, oper2 ) == 0; }
55int ?==?( long int oper2, Int oper1 ) { return mpz_cmp_si( oper1.mpz, oper2 ) == 0; }
56int ?==?( Int oper1, unsigned long int oper2 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) == 0; }
57int ?==?( unsigned long int oper2, Int oper1 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) == 0; }
58
59int ?!=?( Int oper1, Int oper2 ) { return ! ( oper1 == oper2 ); }
60int ?!=?( Int oper1, long int oper2 ) { return ! ( oper1 == oper2 ); }
61int ?!=?( long int oper1, Int oper2 ) { return ! ( oper1 == oper2 ); }
62int ?!=?( Int oper1, unsigned long int oper2 ) { return ! ( oper1 == oper2 ); }
63int ?!=?( unsigned long int oper1, Int oper2 ) { return ! ( oper1 == oper2 ); }
64
65int ?<?( Int oper1, Int oper2 ) { return mpz_cmp( oper1.mpz, oper2.mpz ) < 0; }
66int ?<?( Int oper1, long int oper2 ) { return mpz_cmp_si( oper1.mpz, oper2 ) < 0; }
67int ?<?( long int oper2, Int oper1 ) { return mpz_cmp_si( oper1.mpz, oper2 ) < 0; }
68int ?<?( Int oper1, unsigned long int oper2 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) < 0; }
69int ?<?( unsigned long int oper2, Int oper1 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) < 0; }
70
71int ?<=?( Int oper1, Int oper2 ) { return mpz_cmp( oper1.mpz, oper2.mpz ) <= 0; }
72int ?<=?( Int oper1, long int oper2 ) { return mpz_cmp_si( oper1.mpz, oper2 ) <= 0; }
73int ?<=?( long int oper2, Int oper1 ) { return mpz_cmp_si( oper1.mpz, oper2 ) <= 0; }
74int ?<=?( Int oper1, unsigned long int oper2 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) <= 0; }
75int ?<=?( unsigned long int oper2, Int oper1 ) { return mpz_cmp_ui( oper1.mpz, oper2 ) <= 0; }
76
77int ?>?( Int oper1, Int oper2 ) { return ! ( oper1 <= oper2 ); }
78int ?>?( Int oper1, long int oper2 ) { return ! ( oper1 <= oper2 ); }
79int ?>?( long int oper1, Int oper2 ) { return ! ( oper1 <= oper2 ); }
80int ?>?( Int oper1, unsigned long int oper2 ) { return ! ( oper1 <= oper2 ); }
81int ?>?( unsigned long int oper1, Int oper2 ) { return ! ( oper1 <= oper2 ); }
82
83int ?>=?( Int oper1, Int oper2 ) { return ! ( oper1 < oper2 ); }
84int ?>=?( Int oper1, long int oper2 ) { return ! ( oper1 < oper2 ); }
85int ?>=?( long int oper1, Int oper2 ) { return ! ( oper1 < oper2 ); }
86int ?>=?( Int oper1, unsigned long int oper2 ) { return ! ( oper1 < oper2 ); }
87int ?>=?( unsigned long int oper1, Int oper2 ) { return ! ( oper1 < oper2 ); }
88
89// arithmetic
90Int +?( Int oper ) { Int pos; mpz_set( pos.mpz, oper.mpz ); return pos; }
91Int -?( Int oper ) { Int neg; mpz_neg( neg.mpz, oper.mpz ); return neg; }
92Int ~?( Int oper ) { Int comp; mpz_com( comp.mpz, oper.mpz ); return comp; }
93
94Int ?&?( Int oper1, Int oper2 ) { Int conjunction; mpz_and( conjunction.mpz, oper1.mpz, oper2.mpz ); return conjunction; }
95Int ?&?( Int oper1, long int oper2 ) { Int conjunction, temp; mpz_set_si( temp.mpz, oper2 ); mpz_and( conjunction.mpz, oper1.mpz, temp.mpz ); return conjunction; }
96Int ?&?( long int oper1, Int oper2 ) { Int conjunction, temp; mpz_set_si( temp.mpz, oper1 ); mpz_and( conjunction.mpz, temp.mpz, oper2.mpz ); return conjunction; }
97Int ?&?( Int oper1, unsigned long int oper2 ) { Int conjunction, temp; mpz_set_ui( temp.mpz, oper2 ); mpz_and( conjunction.mpz, oper1.mpz, temp.mpz ); return conjunction; }
98Int ?&?( unsigned long int oper1, Int oper2 ) { Int conjunction, temp; mpz_set_ui( temp.mpz, oper1 ); mpz_and( conjunction.mpz, temp.mpz, oper2.mpz ); return conjunction; }
99Int ?&=?( Int * lhs, Int rhs ) { return *lhs = *lhs & rhs; }
100
101Int ?|?( Int oper1, Int oper2 ) { Int disjunction; mpz_ior( disjunction.mpz, oper1.mpz, oper2.mpz ); return disjunction; }
102Int ?|?( Int oper1, long int oper2 ) { Int disjunction, temp; mpz_set_si( temp.mpz, oper2 ); mpz_ior( disjunction.mpz, oper1.mpz, temp.mpz ); return disjunction; }
103Int ?|?( long int oper1, Int oper2 ) { Int disjunction, temp; mpz_set_si( temp.mpz, oper1 ); mpz_ior( disjunction.mpz, temp.mpz, oper2.mpz ); return disjunction; }
104Int ?|?( Int oper1, unsigned long int oper2 ) { Int disjunction, temp; mpz_set_ui( temp.mpz, oper2 ); mpz_ior( disjunction.mpz, oper1.mpz, temp.mpz ); return disjunction; }
105Int ?|?( unsigned long int oper1, Int oper2 ) { Int disjunction, temp; mpz_set_ui( temp.mpz, oper1 ); mpz_ior( disjunction.mpz, temp.mpz, oper2.mpz ); return disjunction; }
106Int ?|=?( Int * lhs, Int rhs ) { return *lhs = *lhs | rhs; }
107
108Int ?^?( Int oper1, Int oper2 ) { Int disjunction; mpz_xor( disjunction.mpz, oper1.mpz, oper2.mpz ); return disjunction; }
109Int ?^?( Int oper1, long int oper2 ) { Int disjunction, temp; mpz_set_si( temp.mpz, oper2 ); mpz_ior( disjunction.mpz, oper1.mpz, temp.mpz ); return disjunction; }
110Int ?^?( long int oper1, Int oper2 ) { Int disjunction, temp; mpz_set_si( temp.mpz, oper1 ); mpz_ior( disjunction.mpz, temp.mpz, oper2.mpz ); return disjunction; }
111Int ?^?( Int oper1, unsigned long int oper2 ) { Int disjunction, temp; mpz_set_ui( temp.mpz, oper2 ); mpz_ior( disjunction.mpz, oper1.mpz, temp.mpz ); return disjunction; }
112Int ?^?( unsigned long int oper1, Int oper2 ) { Int disjunction, temp; mpz_set_ui( temp.mpz, oper1 ); mpz_ior( disjunction.mpz, temp.mpz, oper2.mpz ); return disjunction; }
113Int ?^=?( Int * lhs, Int rhs ) { return *lhs = *lhs ^ rhs; }
114
115Int ?+?( Int addend1, Int addend2 ) { Int sum; mpz_add( sum.mpz, addend1.mpz, addend2.mpz ); return sum; }
116Int ?+?( Int addend1, long int addend2 ) { Int sum; if ( addend2 >= 0 ) mpz_add_ui( sum.mpz, addend1.mpz, addend2 ); else mpz_sub_ui( sum.mpz, addend1.mpz, -addend2 ); return sum; }
117Int ?+?( long int addend2, Int addend1 ) { Int sum; if ( addend2 >= 0 ) mpz_add_ui( sum.mpz, addend1.mpz, addend2 ); else mpz_sub_ui( sum.mpz, addend1.mpz, -addend2 ); return sum; }
118Int ?+?( Int addend1, unsigned long int addend2 ) { Int sum; mpz_add_ui( sum.mpz, addend1.mpz, addend2 ); return sum; }
119Int ?+?( unsigned long int addend2, Int addend1 ) { Int sum; mpz_add_ui( sum.mpz, addend1.mpz, addend2 ); return sum; }
120Int ?+=?( Int * lhs, Int rhs ) { return *lhs = *lhs + rhs; }
121Int ?+=?( Int * lhs, long int rhs ) { return *lhs = *lhs + rhs; }
122Int ?+=?( Int * lhs, unsigned long int rhs ) { return *lhs = *lhs + rhs; }
123Int ++?( Int * lhs ) { return *lhs += 1; }
124Int ?++( Int * lhs ) { Int ret = *lhs; *lhs += 1; return ret; }
125
126Int ?-?( Int minuend, Int subtrahend ) { Int diff; mpz_sub( diff.mpz, minuend.mpz, subtrahend.mpz ); return diff; }
127Int ?-?( Int minuend, long int subtrahend ) { Int diff; if ( subtrahend >= 0 ) mpz_sub_ui( diff.mpz, minuend.mpz, subtrahend ); else mpz_add_ui( diff.mpz, minuend.mpz, -subtrahend ); return diff; }
128Int ?-?( long int minuend, Int subtrahend ) { Int diff; if ( subtrahend >= 0 ) mpz_ui_sub( diff.mpz, minuend, subtrahend.mpz ); else { mpz_add_ui( diff.mpz, subtrahend.mpz, -minuend ); mpz_neg( diff.mpz, diff.mpz ); } return diff; }
129Int ?-?( Int minuend, unsigned long int subtrahend ) { Int diff; mpz_sub_ui( diff.mpz, minuend.mpz, subtrahend ); return diff; }
130Int ?-?( unsigned long int minuend, Int subtrahend ) { Int diff; mpz_ui_sub( diff.mpz, minuend, subtrahend.mpz ); return diff; }
131Int ?-=?( Int * lhs, Int rhs ) { return *lhs = *lhs - rhs; }
132Int ?-=?( Int * lhs, long int rhs ) { return *lhs = *lhs - rhs; }
133Int ?-=?( Int * lhs, unsigned long int rhs ) { return *lhs = *lhs - rhs; }
134Int --?( Int * lhs ) { return *lhs -= 1; }
135Int ?--( Int * lhs ) { Int ret = *lhs; *lhs -= 1; return ret; }
136
137Int ?*?( Int multiplicator, Int multiplicand ) { Int product; mpz_mul( product.mpz, multiplicator.mpz, multiplicand.mpz ); return product; }
138Int ?*?( Int multiplicator, long int multiplicand ) { Int product; mpz_mul_si( product.mpz, multiplicator.mpz, multiplicand ); return product; }
139Int ?*?( long int multiplicand, Int multiplicator ) { Int product; mpz_mul_si( product.mpz, multiplicator.mpz, multiplicand ); return product; }
140Int ?*?( Int multiplicator, unsigned long int multiplicand ) { Int product; mpz_mul_ui( product.mpz, multiplicator.mpz, multiplicand ); return product; }
141Int ?*?( unsigned long int multiplicand, Int multiplicator ) { Int product; mpz_mul_ui( product.mpz, multiplicator.mpz, multiplicand ); return product; }
142Int ?*=?( Int * lhs, Int rhs ) { return *lhs = *lhs * rhs; }
143Int ?*=?( Int * lhs, long int rhs ) { return *lhs = *lhs * rhs; }
144Int ?*=?( Int * lhs, unsigned long int rhs ) { return *lhs = *lhs * rhs; }
145
146// some code for operators "/" and "%" taken from g++ gmpxx.h
147Int ?/?( Int dividend, Int divisor ) { Int quotient; mpz_tdiv_q( quotient.mpz, dividend.mpz, divisor.mpz ); return quotient; }
148Int ?/?( Int dividend, unsigned long int divisor ) { Int quotient; mpz_tdiv_q_ui( quotient.mpz, dividend.mpz, divisor ); return quotient; }
149Int ?/?( unsigned long int dividend, Int divisor ) {
150        Int quotient;
151    if ( mpz_sgn( divisor.mpz ) >= 0 ) {
152                if ( mpz_fits_ulong_p( divisor.mpz ) )
153                        mpz_set_ui( quotient.mpz, dividend / mpz_get_ui( divisor.mpz ) );
154                else
155                        mpz_set_ui( quotient.mpz, 0 );
156        } else {
157                mpz_neg( quotient.mpz, divisor.mpz );
158                if ( mpz_fits_ulong_p( quotient.mpz ) ) {
159                        mpz_set_ui( quotient.mpz, dividend / mpz_get_ui( quotient.mpz ) );
160                        mpz_neg( quotient.mpz, quotient.mpz );
161                } else
162                        mpz_set_ui( quotient.mpz, 0 );
163        } // if
164        return quotient;
165} // ?/?
166Int ?/?( Int dividend, long int divisor ) {
167        Int quotient;
168    if ( divisor >= 0 )
169                mpz_tdiv_q_ui( quotient.mpz, dividend.mpz, divisor );
170    else {
171                mpz_tdiv_q_ui( quotient.mpz, dividend.mpz, -divisor );
172                mpz_neg( quotient.mpz, quotient.mpz );
173        } // if
174        return quotient;
175} // ?/?
176Int ?/?( long int dividend, Int divisor ) {
177        Int quotient;
178    if ( mpz_fits_slong_p( divisor.mpz ) )
179                mpz_set_si( quotient.mpz, dividend / mpz_get_si( divisor.mpz ) );
180    else {
181        // if divisor is bigger than a long then the quotient must be zero, unless dividend==LONG_MIN and
182        // dividend==-LONG_MIN in which case the quotient is -1
183        mpz_set_si( quotient.mpz, mpz_cmpabs_ui( divisor.mpz, (dividend >= 0 ? dividend : -dividend)) == 0 ? -1 : 0 );
184        } // if
185        return quotient;
186} // ?/?
187Int ?/=?( Int * lhs, Int rhs ) { return *lhs = *lhs / rhs; }
188Int ?/=?( Int * lhs, long int rhs ) { return *lhs = *lhs / rhs; }
189Int ?/=?( Int * lhs, unsigned long int rhs ) { return *lhs = *lhs / rhs; }
190
191[ Int, Int ] div( Int dividend, Int divisor ) { Int quotient, remainder; mpz_fdiv_qr( quotient.mpz, remainder.mpz, dividend.mpz, divisor.mpz ); return [ quotient, remainder ]; }
192[ Int, Int ] div( Int dividend, unsigned long int divisor ) { Int quotient, remainder; mpz_fdiv_qr_ui( quotient.mpz, remainder.mpz, dividend.mpz, divisor ); return [ quotient, remainder ]; }
193
194Int ?%?( Int dividend, Int divisor ) { Int remainder; mpz_tdiv_r( remainder.mpz, dividend.mpz, divisor.mpz ); return remainder; }
195Int ?%?( Int dividend, unsigned long int divisor ) { Int remainder; mpz_tdiv_r_ui( remainder.mpz, dividend.mpz, divisor ); return remainder; }
196Int ?%?( unsigned long int dividend, Int divisor ) {
197        Int remainder;
198    if ( mpz_sgn( divisor.mpz ) >= 0 ) {
199                if ( mpz_fits_ulong_p( divisor.mpz ) )
200                        mpz_set_ui( remainder.mpz, dividend % mpz_get_ui( divisor.mpz ) );
201                else
202                        mpz_set_ui( remainder.mpz, dividend );
203        } else {
204                mpz_neg( remainder.mpz, divisor.mpz );
205                if ( mpz_fits_ulong_p( remainder.mpz ) )
206                        mpz_set_ui( remainder.mpz, dividend % mpz_get_ui( remainder.mpz ) );
207                else
208                        mpz_set_ui( remainder.mpz, dividend );
209        } // if
210        return remainder;
211} // ?%?
212Int ?%?( Int dividend, long int divisor ) {
213        Int remainder;
214    mpz_tdiv_r_ui( remainder.mpz, dividend.mpz, (divisor >= 0 ? divisor : -divisor));
215        return remainder;
216} // ?%?
217Int ?%?( long int dividend, Int divisor ) {
218        Int remainder;
219    if ( mpz_fits_slong_p( divisor.mpz ) )
220                mpz_set_si( remainder.mpz, dividend % mpz_get_si( divisor.mpz ) );
221        else {
222                // if divisor is bigger than a long then the remainder is dividend unchanged, unless dividend==LONG_MIN and
223                // dividend==-LONG_MIN in which case the remainder is 0
224        mpz_set_si( remainder.mpz, mpz_cmpabs_ui( divisor.mpz, (dividend >= 0 ? dividend : -dividend)) == 0 ? 0 : dividend);
225        } // if
226        return remainder;
227} // ?%?
228Int ?%=?( Int * lhs, Int rhs ) { return *lhs = *lhs % rhs; }
229Int ?%=?( Int * lhs, long int rhs ) { return *lhs = *lhs % rhs; }
230Int ?%=?( Int * lhs, unsigned long int rhs ) { return *lhs = *lhs % rhs; }
231
232Int ?<<?( Int shiften, mp_bitcnt_t shift ) { Int shifted; mpz_mul_2exp( shifted.mpz, shiften.mpz, shift ); return shifted; }
233Int ?<<=?( Int * lhs, mp_bitcnt_t shift ) { return *lhs = *lhs << shift; }
234Int ?>>?( Int shiften, mp_bitcnt_t shift ) { Int shifted; mpz_fdiv_q_2exp( shifted.mpz, shiften.mpz, shift ); return shifted; }
235Int ?>>=?( Int * lhs, mp_bitcnt_t shift ) { return *lhs = *lhs >> shift; }
236
237// number functions
238Int abs( Int oper ) { Int positive; mpz_abs( positive.mpz, oper.mpz ); return positive; }
239Int fact( unsigned long int N ) { Int factorial; mpz_fac_ui( factorial.mpz, N ); return factorial; }
240Int gcd( Int oper1, Int oper2 ) { Int gcdret; mpz_gcd( gcdret.mpz, oper1.mpz, oper2.mpz ); return gcdret; }
241Int pow( Int base, unsigned long int exponent ) { Int power; mpz_pow_ui( power.mpz, base.mpz, exponent ); return power; }
242Int pow( unsigned long int base, unsigned long int exponent ) { Int power; mpz_ui_pow_ui( power.mpz, base, exponent ); return power; }
243void srandom( gmp_randstate_t state ) { gmp_randinit_default( state ); }
244Int random( gmp_randstate_t state, mp_bitcnt_t n ) { Int rand; mpz_urandomb( rand.mpz, state, n ); return rand; }
245Int random( gmp_randstate_t state, Int n ) { Int rand; mpz_urandomm( rand.mpz, state, n.mpz ); return rand; }
246Int random( gmp_randstate_t state, mp_size_t max_size ) { Int rand; mpz_random( rand.mpz, max_size ); return rand; }
247int sgn( Int oper ) { return mpz_sgn( oper.mpz ); }
248Int sqrt( Int oper ) { Int root; mpz_sqrt( root.mpz, oper.mpz ); return root; }
249
250// I/O
251forall( dtype istype | istream( istype ) )
252istype * ?|?( istype * is, Int * mp ) {
253        gmp_scanf( "%Zd", mp );
254        return is;
255} // ?|?
256
257forall( dtype ostype | ostream( ostype ) )
258ostype * ?|?( ostype * os, Int mp ) {
259        if ( sepPrt( os ) ) fmt( os, "%s", sepGetCur( os ) );
260        gmp_printf( "%Zd", mp.mpz );
261        sepOn( os );
262        return os;
263} // ?|?
264
265// Local Variables: //
266// mode: c //
267// tab-width: 4 //
268// End: //
Note: See TracBrowser for help on using the repository browser.