home *** CD-ROM | disk | FTP | other *** search
/ GEMini Atari / GEMini_Atari_CD-ROM_Walnut_Creek_December_1993.iso / files / gnu / gplibt02 / trationa.cc < prev    next >
C/C++ Source or Header  |  1993-07-24  |  4KB  |  186 lines
  1. /*
  2.  a test file for Rational class
  3.  */
  4.  
  5. #include <xrationa.h>
  6.  
  7. #include <assert.h>
  8.  
  9. #define tassert(ex) {if ((ex)) cerr << #ex << "\n"; \
  10.                        else _assert(#ex, __FILE__,__LINE__); }
  11.  
  12. #ifndef PI
  13. #  define PI (3.14159265358979323846)
  14. #endif
  15.  
  16. void identtest(Rational& a, Rational& b, Rational& c)
  17. {
  18.   Rational one(1, 1);
  19.  
  20.   assert(-(-a) == a);
  21.   assert((a + b) ==  (b + a));
  22.   assert((a + (-b)) ==  (a - b));
  23.   assert((a * b) ==  (b * a));
  24.   assert((a * (-b)) == -(a * b));
  25.   assert((a / (-b)) == -(a / b));
  26.   assert((a / b) == (a * (one / b)));
  27.   assert((a / b) == (one / (b / a)));
  28.   assert((a - b) ==  -(b - a));
  29.   assert((a + (b + c)) == ((a + b) + c));
  30.   assert((a * (b * c)) == ((a * b) * c));
  31.   assert((a * (b + c)) == ((a * b) + (a * c)));
  32.   assert(((a - b) + b) == a);
  33.   assert(((a + b) - b) == a);
  34.   assert(((a * b) / b) == a);
  35.   assert(((a / b) * b) == a);
  36.  
  37.   Rational x = a;
  38.   x *= b;
  39.   assert(x == (a * b));
  40.   x += c;
  41.   assert(x == ((a * b) + c));
  42.   x -= a;
  43.   assert(x == (((a * b) + c) - a));
  44.   x /= b;
  45.   assert(x == ((((a * b) + c) - a) / b));
  46.  
  47.   assert(x.OK());
  48. }
  49.  
  50.  
  51.  
  52. void simpletest()
  53. {
  54.   Rational one = 1;
  55.   assert(one.OK());
  56.   Rational third(1, 3);
  57.   assert(third.OK());
  58.   Rational half(1, 2);
  59.   assert(half.OK());
  60.  
  61.   cout << "one = " << one << "\n";
  62.   cout << "third = " << third << "\n";
  63.   cout << "half = " << half << "\n";
  64.  
  65.   cout << "third + half = " << third + half << "\n";
  66.   cout << "third - half = " << third - half << "\n";
  67.   cout << "third * half = " << third * half << "\n";
  68.   cout << "third / half = " << third / half << "\n";
  69.  
  70.   Rational onePointTwo = 1.2;
  71.   cout << "onePointTwo = " << onePointTwo << "\n";
  72.   cout << "double(onePointTwo) = " << double(onePointTwo) << "\n";
  73.  
  74.   Rational a = one;
  75.   cout << "a = " << a << "\n";
  76.   assert(a.OK());
  77.   a += half;
  78.   cout << "a += half = " << a << "\n";
  79.   assert(a == Rational(3, 2));
  80.   a -= half;
  81.   cout << "a -= half = " << a << "\n";
  82.   assert(a == Rational(1));
  83.   a *= half;
  84.   cout << "a *= half = " << a << "\n";
  85.   assert(a == half);
  86.   a /= half;
  87.   cout << "a /= half = " << a << "\n";
  88.   assert(a == Rational(1));
  89.   assert(a.OK());
  90.  
  91.   identtest(one, one, one);
  92.   identtest(one, third, half);
  93.   identtest(third, half, one);
  94.   identtest(onePointTwo, half, a);
  95. }
  96.  
  97. void pitest()
  98. {
  99.   Rational half(1, 2);
  100.   Rational approxpi(355, 113);
  101.   assert(approxpi.OK());
  102.   cout << "approxpi = " << approxpi << "\n";
  103.   cout << "double(approxpi) = " << double(approxpi) << "\n";
  104.  
  105.   Rational rpi = Rational(PI);
  106.   cout << "rpi = Rational(PI) = " << rpi << "\n";
  107.   assert(rpi.OK());
  108.   cout << "double(rpi) = " << double(rpi) << "\n";
  109.  
  110.   cout << "approxpi + rpi = " << approxpi + rpi << "\n";
  111.   cout << "approxpi - rpi = " << approxpi - rpi << "\n";
  112.   cout << "approxpi * rpi = " << approxpi * rpi << "\n";
  113.   cout << "approxpi / rpi = " << approxpi / rpi << "\n";
  114.  
  115.   Rational negapproxpi = -approxpi;
  116.  
  117.   cout << "-approxpi = " << negapproxpi << "\n";
  118.   assert(sign(negapproxpi) < 0);
  119.   cout << "abs(negapproxpi) = " << abs(negapproxpi) << "\n";
  120.   assert(abs(negapproxpi) == approxpi);
  121.  
  122.   assert(approxpi != rpi);
  123.   assert(approxpi >= rpi);
  124.   assert(approxpi > rpi);
  125.   assert(!(approxpi == rpi));
  126.   assert(!(approxpi <= rpi));
  127.   assert(!(approxpi < rpi));
  128.   assert((approxpi >? rpi) == approxpi);
  129.   assert((approxpi <? rpi) == rpi);
  130.  
  131.   assert(floor(approxpi) == 3);
  132.   assert(ceil(approxpi) == 4);
  133.   assert(trunc(approxpi) == 3);
  134.   assert(round(approxpi) == 3);
  135.  
  136.   assert(floor(negapproxpi + half) == -3);
  137.   assert(ceil(negapproxpi + half) == -2);
  138.   assert(trunc(negapproxpi + half) == -2);
  139.   assert(round(negapproxpi + half) == -3);
  140.  
  141.   identtest(approxpi, rpi, negapproxpi);
  142.   identtest(rpi, approxpi, rpi);
  143.   identtest(negapproxpi, half, rpi);
  144. }
  145.  
  146.  
  147.  
  148. void IOtest()
  149. {
  150.   Rational a;
  151.   cout << "\nenter a Rational in form a/b or a: ";
  152.   cin >> a;
  153.   cout << "number = " << a << "\n";
  154.   assert(a.OK());
  155. }
  156.  
  157. // as a fct just to test Rational fcts
  158. Rational estimate_e(long n)
  159. {
  160.   Rational x = Rational(n + 1, n);
  161.   Rational e = pow(x, n);
  162.   return e;
  163. }
  164.  
  165. void etest(long n)
  166. {
  167.   cout << "approximating e as pow(1+1/n),n) for n =" << n << "\n";
  168.   Rational approxe = estimate_e(n);
  169.   assert(approxe.OK());
  170.   cout << "double(approxe) = " << double(approxe) << "\n";
  171.   cout << "log(approxe) = " << log(approxe) << "\n";
  172.   assert(log(approxe) <= 1.0);
  173.   cout << "approxe = " << approxe << "\n";
  174. }
  175.  
  176. main()
  177. {
  178.   simpletest();
  179.   pitest();
  180.   IOtest();
  181.   etest(10);
  182.   etest(100);
  183.   etest(1000);
  184.   cout << "\nEnd of test\n";
  185. }
  186.