| [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: // | 
|---|