summaryrefslogtreecommitdiff
path: root/src/util/rational.h
diff options
context:
space:
mode:
authorTim King <taking@cs.nyu.edu>2010-03-26 22:37:12 +0000
committerTim King <taking@cs.nyu.edu>2010-03-26 22:37:12 +0000
commit3b8d92b17c17aadfa920f36a1ab631e36c70e00e (patch)
tree3d19639afd33f5329d5c3d29c6b020acdb229ad6 /src/util/rational.h
parent6a383befdf0fd88ff3c76dc001777475b34cf694 (diff)
Added GMP backed Rational and Integer classes, and white box tests for them. You may have to reconfigure after this update.
Diffstat (limited to 'src/util/rational.h')
-rw-r--r--src/util/rational.h172
1 files changed, 172 insertions, 0 deletions
diff --git a/src/util/rational.h b/src/util/rational.h
new file mode 100644
index 000000000..e3f760031
--- /dev/null
+++ b/src/util/rational.h
@@ -0,0 +1,172 @@
+/********************* */
+/** rational.h
+ ** Original author: taking
+ ** Major contributors: none
+ ** Minor contributors (to current version): none
+ ** This file is part of the CVC4 prototype.
+ ** Copyright (c) 2009, 2010 The Analysis of Computer Systems Group (ACSys)
+ ** Courant Institute of Mathematical Sciences
+ ** New York University
+ ** See the file COPYING in the top-level source directory for licensing
+ ** information.
+ **
+ ** A multiprecision rational constant.
+ **/
+
+#include <gmpxx.h>
+#include <string>
+#include "integer.h"
+
+#ifndef __CVC4__RATIONAL_H
+#define __CVC4__RATIONAL_H
+
+namespace CVC4 {
+
+class Rational {
+private:
+ /**
+ * Stores the value of the rational is stored in a C++ GMP rational class.
+ * Using this instead of mpq_t allows for easier destruction.
+ */
+ mpq_class d_value;
+
+ /**
+ * Constructs a Rational from a mpq_class object.
+ * Does a deep copy.
+ * Assumes that the value is in canonical form, and thus does not
+ * have to call canonicalize() on the value.
+ */
+ Rational(const mpq_class& val) : d_value(val) { }
+
+public:
+
+ /** Constructs a rational with the value 0/1. */
+ Rational() : d_value(0){
+ d_value.canonicalize();
+ }
+
+ /**
+ * Constructs a Rational from a C string.
+ * Throws std::invalid_argument if the stribng is not a valid rational.
+ * For more information about what is a vaid rational string,
+ * see GMP's documentation for mpq_set_str().
+ */
+ Rational(const char * s, int base = 10): d_value(s,base) {
+ d_value.canonicalize();
+ }
+ Rational(const std::string& s, unsigned base = 10) : d_value(s, base) {
+ d_value.canonicalize();
+ }
+
+ Rational(const Rational& q) : d_value(q.d_value) {
+ d_value.canonicalize();
+ }
+
+ Rational( signed int n, signed int d) : d_value(n,d) {
+ d_value.canonicalize();
+ }
+ Rational(unsigned int n, unsigned int d) : d_value(n,d) {
+ d_value.canonicalize();
+ }
+ Rational( signed long int n, signed long int d) : d_value(n,d) {
+ d_value.canonicalize();
+ }
+ Rational(unsigned long int n, unsigned long int d) : d_value(n,d) {
+ d_value.canonicalize();
+ }
+
+ Rational(const Integer& n, const Integer& d):
+ d_value(n.get_mpz(), d.get_mpz())
+ {
+ d_value.canonicalize();
+ }
+
+ ~Rational() {}
+
+
+ Integer getNumerator() const {
+ return Integer(d_value.get_num());
+ }
+
+ Integer getDenominator() const{
+ return Integer(d_value.get_den());
+ }
+
+
+ Rational& operator=(const Rational& x){
+ if(this == &x) return *this;
+ d_value = x.d_value;
+ return *this;
+ }
+
+ Rational operator-() const{
+ return Rational(-(d_value));
+ }
+
+ bool operator==(const Rational& y) const {
+ return d_value == y.d_value;
+ }
+
+ bool operator!=(const Rational& y) const {
+ return d_value != y.d_value;
+ }
+
+ bool operator< (const Rational& y) const {
+ return d_value < y.d_value;
+ }
+
+ bool operator<=(const Rational& y) const {
+ return d_value <= y.d_value;
+ }
+
+ bool operator> (const Rational& y) const {
+ return d_value > y.d_value;
+ }
+
+ bool operator>=(const Rational& y) const {
+ return d_value >= y.d_value;
+ }
+
+
+ Rational operator+(const Rational& y) const{
+ return Rational( d_value + y.d_value );
+ }
+
+ Rational operator-(const Rational& y) const {
+ return Rational( d_value - y.d_value );
+ }
+
+ Rational operator*(const Rational& y) const {
+ return Rational( d_value * y.d_value );
+ }
+ Rational operator/(const Rational& y) const {
+ return Rational( d_value / y.d_value );
+ }
+
+ std::string toString(int base = 10) const {
+ return d_value.get_str(base);
+ }
+
+
+ friend std::ostream& operator<<(std::ostream& os, const Rational& n);
+
+
+ /**
+ * Computes the hash of the node from the first word of the
+ * numerator, the denominator.
+ */
+ size_t hash() const {
+ size_t numeratorHash = gmpz_hash(d_value.get_num_mpz_t());
+ size_t denominatorHash = gmpz_hash(d_value.get_den_mpz_t());
+
+ return numeratorHash xor denominatorHash;
+ }
+
+};/* class Rational */
+
+std::ostream& operator<<(std::ostream& os, const Rational& n);
+
+}/* CVC4 namespace */
+
+#endif /* __CVC4__RATIONAL_H */
+
generated by cgit on debian on lair
contact matthew@masot.net with questions or feedback