[bb82c03] | 1 | //
|
---|
[53ba273] | 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.
|
---|
[bb82c03] | 6 | //
|
---|
[630a82a] | 7 | // rational -- Rational numbers are numbers written as a ratio, i.e., as a fraction, where the numerator (top number)
|
---|
| 8 | // and the denominator (bottom number) are whole numbers. When creating and computing with rational numbers, results
|
---|
| 9 | // are constantly reduced to keep the numerator and denominator as small as possible.
|
---|
[bb82c03] | 10 | //
|
---|
[53ba273] | 11 | // Author : Peter A. Buhr
|
---|
| 12 | // Created On : Wed Apr 6 17:56:25 2016
|
---|
| 13 | // Last Modified By : Peter A. Buhr
|
---|
[6c6455f] | 14 | // Last Modified On : Mon May 15 21:30:12 2017
|
---|
| 15 | // Update Count : 90
|
---|
[bb82c03] | 16 | //
|
---|
[f621a148] | 17 |
|
---|
[17e5e2b] | 18 | #ifndef RATIONAL_H
|
---|
| 19 | #define RATIONAL_H
|
---|
[53ba273] | 20 |
|
---|
[3d9b5da] | 21 | #include "iostream"
|
---|
[53ba273] | 22 |
|
---|
[561f730] | 23 | trait scalar( otype T ) {
|
---|
| 24 | };
|
---|
| 25 |
|
---|
| 26 | trait arithmetic( otype T | scalar( T ) ) {
|
---|
| 27 | int !?( T );
|
---|
| 28 | int ?==?( T, T );
|
---|
| 29 | int ?!=?( T, T );
|
---|
| 30 | int ?<?( T, T );
|
---|
| 31 | int ?<=?( T, T );
|
---|
| 32 | int ?>?( T, T );
|
---|
| 33 | int ?>=?( T, T );
|
---|
| 34 | void ?{}( T *, zero_t );
|
---|
| 35 | void ?{}( T *, one_t );
|
---|
| 36 | T +?( T );
|
---|
| 37 | T -?( T );
|
---|
| 38 | T ?+?( T, T );
|
---|
| 39 | T ?-?( T, T );
|
---|
| 40 | T ?*?( T, T );
|
---|
| 41 | T ?/?( T, T );
|
---|
| 42 | T ?%?( T, T );
|
---|
| 43 | T ?/=?( T *, T );
|
---|
| 44 | T abs( T );
|
---|
| 45 | };
|
---|
| 46 |
|
---|
[630a82a] | 47 | // implementation
|
---|
[561f730] | 48 |
|
---|
| 49 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
[53ba273] | 50 | struct Rational {
|
---|
[561f730] | 51 | RationalImpl numerator, denominator; // invariant: denominator > 0
|
---|
[53ba273] | 52 | }; // Rational
|
---|
| 53 |
|
---|
[630a82a] | 54 | // constructors
|
---|
[561f730] | 55 |
|
---|
| 56 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 57 | void ?{}( Rational(RationalImpl) * r );
|
---|
| 58 |
|
---|
| 59 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 60 | void ?{}( Rational(RationalImpl) * r, RationalImpl n );
|
---|
| 61 |
|
---|
| 62 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 63 | void ?{}( Rational(RationalImpl) * r, RationalImpl n, RationalImpl d );
|
---|
| 64 |
|
---|
| 65 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 66 | void ?{}( Rational(RationalImpl) * r, zero_t );
|
---|
| 67 |
|
---|
| 68 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 69 | void ?{}( Rational(RationalImpl) * r, one_t );
|
---|
[630a82a] | 70 |
|
---|
[f621a148] | 71 | // getter for numerator/denominator
|
---|
[561f730] | 72 |
|
---|
| 73 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 74 | RationalImpl numerator( Rational(RationalImpl) r );
|
---|
| 75 |
|
---|
| 76 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 77 | RationalImpl denominator( Rational(RationalImpl) r );
|
---|
| 78 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 79 | [ RationalImpl, RationalImpl ] ?=?( * [ RationalImpl, RationalImpl ] dest, Rational(RationalImpl) src );
|
---|
| 80 |
|
---|
[f621a148] | 81 | // setter for numerator/denominator
|
---|
[561f730] | 82 |
|
---|
| 83 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 84 | RationalImpl numerator( Rational(RationalImpl) r, RationalImpl n );
|
---|
| 85 |
|
---|
| 86 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 87 | RationalImpl denominator( Rational(RationalImpl) r, RationalImpl d );
|
---|
[630a82a] | 88 |
|
---|
| 89 | // comparison
|
---|
[561f730] | 90 |
|
---|
| 91 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 92 | int ?==?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 93 |
|
---|
| 94 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 95 | int ?!=?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 96 |
|
---|
| 97 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 98 | int ?<?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 99 |
|
---|
| 100 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 101 | int ?<=?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 102 |
|
---|
| 103 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 104 | int ?>?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 105 |
|
---|
| 106 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 107 | int ?>=?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
[630a82a] | 108 |
|
---|
| 109 | // arithmetic
|
---|
[561f730] | 110 |
|
---|
| 111 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 112 | Rational(RationalImpl) +?( Rational(RationalImpl) r );
|
---|
| 113 |
|
---|
| 114 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 115 | Rational(RationalImpl) -?( Rational(RationalImpl) r );
|
---|
| 116 |
|
---|
| 117 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 118 | Rational(RationalImpl) ?+?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 119 |
|
---|
| 120 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 121 | Rational(RationalImpl) ?-?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 122 |
|
---|
| 123 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 124 | Rational(RationalImpl) ?*?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
| 125 |
|
---|
| 126 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 127 | Rational(RationalImpl) ?/?( Rational(RationalImpl) l, Rational(RationalImpl) r );
|
---|
[630a82a] | 128 |
|
---|
| 129 | // conversion
|
---|
[6c6455f] | 130 | forall ( otype RationalImpl | arithmetic( RationalImpl ) | { double convert( RationalImpl ); } )
|
---|
| 131 | double widen( Rational(RationalImpl) r );
|
---|
| 132 | forall ( otype RationalImpl | arithmetic( RationalImpl ) | { double convert( RationalImpl ); RationalImpl convert( double );} )
|
---|
| 133 | Rational(RationalImpl) narrow( double f, RationalImpl md );
|
---|
[630a82a] | 134 |
|
---|
| 135 | // I/O
|
---|
[561f730] | 136 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 137 | forall( dtype istype | istream( istype ) | { istype * ?|?( istype *, RationalImpl * ); } )
|
---|
| 138 | istype * ?|?( istype *, Rational(RationalImpl) * );
|
---|
| 139 |
|
---|
| 140 | forall ( otype RationalImpl | arithmetic( RationalImpl ) )
|
---|
| 141 | forall( dtype ostype | ostream( ostype ) | { ostype * ?|?( ostype *, RationalImpl ); } )
|
---|
| 142 | ostype * ?|?( ostype *, Rational(RationalImpl ) );
|
---|
[53ba273] | 143 |
|
---|
[17e5e2b] | 144 | #endif // RATIONAL_H
|
---|
| 145 |
|
---|
[53ba273] | 146 | // Local Variables: //
|
---|
| 147 | // mode: c //
|
---|
| 148 | // tab-width: 4 //
|
---|
| 149 | // End: //
|
---|