famous_test.c

   1 // Test famous continued fraction expansions are correct.
   2 // as well as expansions based on Taylor series.
   3
   4 #include <stdio.h>
   5 #include <stdlib.h>
   6 #include <string.h>
   7 #include <gmp.h>
   8 #include "cf.h"
   9 #include "test.h"
  10
  11 int main() {
  12   CF_NEW_EXPECT_DEC(cf_new_sqrt2, "1.41421356237309504880");
  13   CF_NEW_EXPECT_DEC(cf_new_e, "2.71828182845904523536");
  14   CF_NEW_EXPECT_DEC(cf_new_pi, "3.1415926535897932384");
  15   CF_NEW_EXPECT_DEC(cf_new_tan1, "1.5574077246549022305");
  16   CF_NEW_EXPECT_DEC(cf_new_sin1, "0.8414709848078965066");
  17   CF_NEW_EXPECT_DEC(cf_new_cos1, "0.5403023058681397174");
  18
  19   mpz_t z;
  20   mpz_init(z);
  21   cf_t x;
  22
  23   mpz_set_si(z, 3);
  24   x = cf_new_tanh(z);
  25   CF_EXPECT_DEC(x, "0.99505475368673045133188018525548847509781385470028");
  26   cf_free(x);
  27
  28   mpz_set_ui(z, 2);
  29   x = cf_new_epow(z);
  30   CF_EXPECT_DEC(x, "7.38905609893065022723042746057500781318031557055184");
  31   cf_free(x);
  32
  33   mpz_clear(z);
  34   return 0;
  35 }