doc/mpc.texi: Add documentation for new functions mpc_set_str and mpc_strtoc.

doc/version.texi: Update to february 2009.
NEWS src/mpc.h src/Makefile.am: Add new functions mpc_set_str and mpc_strtoc.
configure.ac: Check if locale.h exists.
src/strtoc.c: new function.
src/set_str.c: new function.
tests/mpc-tests.h: make public some helper functions from read_data.c.
tests/read_data.c: make public some helper functions.
tests/Makefile.am: Add tstrtoc.c and tstrtoc.dat.
tests/tstrtoc.c: test file for mpc_strtoc.
tests/strtoc.dat: data file for mpc_strtoc.


git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@414 211d60ee-9f03-0410-a15a-8952a2c7a4e4

1 files changed, 254 insertions, 0 deletions

diff --git a/tests/tstrtoc.c b/tests/tstrtoc.c
new file mode 100644
index 0000000..326665c
--- /dev/null
+++ b/tests/tstrtoc.c
@@ -0,0 +1,254 @@
+/* test file for mpc_strtoc.
+
+Copyright (C) 2009 Philippe Th\'eveny
+
+This file is part of the MPC Library.
+
+The MPC Library is free software; you can redistribute it and/or modify
+it under the terms of the GNU Lesser General Public License as published by
+the Free Software Foundation; either version 2.1 of the License, or (at your
+option) any later version.
+
+The MPC Library is distributed in the hope that it will be useful, but
+WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
+or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
+License for more details.
+
+You should have received a copy of the GNU Lesser General Public License
+along with the MPC Library; see the file COPYING.LIB.  If not, write to
+the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
+MA 02111-1307, USA. */
+
+#include <ctype.h>
+#include <string.h>
+#include <stdio.h>
+#include <stdlib.h>
+
+#include "mpc.h"
+#include "mpc-impl.h"
+#include "mpc-tests.h"
+
+#ifndef __SRCDIR
+#define __SRCDIR .
+#endif
+
+extern unsigned long line_number;
+extern int nextchar;
+extern char *pathname;
+
+static const char *rnd_mode[] =
+  { "MPC_RNDNN", "MPC_RNDZN", "MPC_RNDUN", "MPC_RNDDN",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined",
+    "MPC_RNDNZ", "MPC_RNDZZ", "MPC_RNDUZ", "MPC_RNDDZ",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined",
+    "MPC_RNDNU", "MPC_RNDZU", "MPC_RNDUU", "MPC_RNDDU",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined",
+    "MPC_RNDND", "MPC_RNDZD", "MPC_RNDUD", "MPC_RNDDD",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined", "undefined", "undefined", "undefined",
+    "undefined", "undefined",
+  };
+
+static char file_name[] = "strtoc.dat";
+
+static void
+read_int (FILE *fp, int *nread, const char *name)
+{
+  int n = 0;
+
+  if (nextchar == EOF)
+    {
+      printf ("Error: Unexpected EOF when reading mpfr precision "
+              "in file '%s' line %ld\n",
+              pathname, line_number);
+      exit (1);
+    }
+  ungetc (nextchar, fp);
+  n = fscanf (fp, "%i", nread);
+  if (ferror (fp) || n == 0 || n == EOF)
+    {
+      printf ("Error: Cannot read %s in file '%s' line %ld\n",
+              name, pathname, line_number);
+      exit (1);
+    }
+  nextchar = getc (fp);
+  skip_whitespace_comments (fp);
+}
+
+static void
+read_string (FILE *fp, char *str, const char *name)
+{
+  size_t len;
+  size_t str_len;
+
+  len = 0;
+  str_len = str == NULL ? 0 : sizeof (str);
+
+  if (nextchar == '"')
+    nextchar = getc (fp);
+  else
+    goto error;
+
+  while (nextchar != EOF && nextchar != '"')
+    {
+      if (nextchar == '\n')
+        line_number++;
+      if (len + 1 > str_len)
+        {
+          str = realloc (str, 2 * str_len);
+          if (str == NULL)
+            {
+              printf ("Cannot allocate memory\n");
+              exit (1);
+            }
+          str_len *= 2;
+        }
+      str[len++] = nextchar;
+      nextchar = getc (fp);
+    }
+
+  if (nextchar != '"')
+    goto error;
+
+  if (len + 1 > str_len)
+    {
+      str = realloc (str, str_len + 1);
+      if (str == NULL)
+        {
+          printf ("Cannot allocate memory\n");
+          exit (1);
+        }
+      str_len *= 2;
+    }
+  str[len] = '\0';
+
+  nextchar = getc (fp);
+  skip_whitespace_comments (fp);
+  return;
+
+ error:
+  printf ("Error: Unable to read %s in file '%s' line '%ld'\n",
+          name, pathname, line_number);
+  exit (1);
+}
+
+int
+main (void)
+{
+  FILE *fp;
+
+  char *str = NULL;
+  char *rstr = NULL;
+  char *end = NULL;
+
+  int base;
+  int inex_re;
+  int inex_im;
+  mpc_t expected, got;
+  mpc_rnd_t rnd = MPC_RNDNN;
+  int inex = 0;
+  known_signs_t ks = {1, 1};
+
+
+  /* 1. open data file */
+  pathname = (char *)malloc ((strlen (QUOTE (__SRCDIR)) + strlen (file_name)
+                              + 2) * sizeof (char));
+  if (pathname == NULL)
+    {
+      printf ("Cannot allocate memory\n");
+      exit (1);
+    }
+  sprintf (pathname, QUOTE (__SRCDIR)"/%s", file_name);
+  fp = fopen (pathname, "r");
+  if (fp == NULL)
+    {
+      fprintf (stderr, "Unable to open %s\n", pathname);
+      exit (1);
+    }
+
+  /* 2. init needed variables */
+  str = (char *) malloc (255 * sizeof (char));
+  if (str == NULL)
+    {
+      printf ("Cannot allocate memory\n");
+      exit (1);
+    }
+  rstr = (char *) malloc (255 * sizeof (char));
+  if (str == NULL)
+    {
+      printf ("Cannot allocate memory\n");
+      exit (1);
+    }
+  mpc_init2 (expected, 53);
+  mpc_init2 (got, 53);
+
+  /* 3. read data file */
+  line_number = 1;
+  nextchar = getc (fp);
+  while (nextchar != EOF)
+    {
+      skip_whitespace_comments (fp);
+
+      /* 3.1 read a line of data: expected result, base, rounding mode */
+      read_ternary (fp, &inex_re);
+      read_ternary (fp, &inex_im);
+      read_mpc (fp, expected, NULL);
+
+      read_string (fp, str, "number string");
+      read_string (fp, rstr, "string remainder");
+      read_int (fp, &base, "base");
+      read_mpc_rounding_mode (fp, &rnd);
+
+      /* 3.2 convert string at the same precision as the expected result */
+      mpfr_set_prec (MPC_RE (got), MPC_PREC_RE (expected));
+      mpfr_set_prec (MPC_IM (got), MPC_PREC_IM (expected));
+      inex = mpc_strtoc (got, str, &end, base, rnd);
+
+      /* 3.3 compare this result with the expected one */
+      if (inex != MPC_INEX (inex_re, inex_im)
+          || !same_mpc_value (got, expected, ks)
+          || strcmp (end, rstr) != 0)
+        {
+          printf ("mpc_strtoc(str) failed (line %lu)\nwith base=%d and "
+                  "rounding mode %s\n", line_number-1, base, rnd_mode[rnd]);
+          if (inex != MPC_INEX (inex_re, inex_im))
+            printf ("ternary value: got %s, expected (%s, %s)\n",
+                    MPC_INEX_STR (inex),
+                    (inex_re == +1 ? "+1" : (inex_re == -1 ? "-1" : "0")),
+                    (inex_im == +1 ? "+1" : (inex_im == -1 ? "-1" : "0")));
+          printf ("str = \"%s\"\n", str);
+          if (strcmp (end, rstr) != 0)
+            printf ("string remainder expected \"%s\"\n"
+                    "                 got      \"%s\"\n",
+                    rstr, end);
+          else
+            {
+              printf ("     ");
+              OUT (got);
+              OUT (expected);
+            }
+          exit (1);
+        }
+
+      end = NULL;
+    }
+
+  fclose (fp);
+
+  /* 4. Clear used variables */
+  mpc_clear (expected);
+  mpc_clear (got);
+  if (str != NULL)
+    free (str);
+  if (rstr != NULL)
+    free (rstr);
+  free (pathname);
+
+  return 0;
+}