Laboratoire de l'Informatique et du Parallélisme » Pobyso

/** @file test-pobyso-dirty-find-zeros.c

 * Name & purpose

 *

 * Three arguments are expected : an functional expression, a lowerBound

 * and an upper bound. The zeroes of the functional expression are

 * searched for over the [lower bound, upper bound] interval.

 * If only the functional expression is literally "NULL" the

 * pobyso_dirty_find_zeros is tested with a NULL argument as for the function.

 *

 * The others arguments are not looked for.

 *

 * @author S. T.

 * @date 2014-12-05

 *

 * @details Verbosity manipulation are performed  to make sur  that verbosity

 *  suppression is correctly managed in pobyso_is_constant_expression()

 *  function.

 */

/******************************************************************************/

/* Headers, applying the "particular to general" convention.*/

#include "pobyso.h"

/* includes of local headers */

/* includes of project headers */

/* includes of system headers */

#include <stdlib.h>

/* Other declarations */

/* Internal prototypes */

/* Types, constants and macros definitions */

/* Global variables */

/* Functions */

int

main(int argc, char** argv)

{

  pobyso_func_exp_t funcExpSo  = NULL;

  sollya_obj_t lowerBoundSo = NULL;

  sollya_obj_t upperBoundSo = NULL;

  mpfr_t lowerBoundMp;

  mpfr_t upperBoundMp;

  int zerosCount = 0;

  mpfr_t* zerosArrayMp = NULL;

  int i;

  sollya_lib_init();

  /* Two command line arguments expected. */

  if (argc < 4)

  {

    fprintf(stderr,"\n\n");

    fprintf(stderr, "Usage: %s funcExp lowerBound upperBound\n", argv[0]);

    fprintf(stderr,

            "            funcExp: a quoted functional expression for Sollya;\n");

    fprintf(stderr,

            "            lowerBound: a (possibly quoted) Sollya constant expression for the lower bound;\n");

    fprintf(stderr,

            "            upperBound: a (possibly quoted) Sollya constant expression for the upper bound.");

    fprintf(stderr,"\n\n");

    sollya_lib_close();

    return 1;

  }

  if (! strcmp(argv[1], "NULL")) /* First parameter is the "NULL" string) */

  {

    zerosArrayMp = pobyso_dirty_find_zeros_bounds(NULL,

                                                  lowerBoundMp,

                                                  upperBoundMp,

                                                  &zerosCount);

    sollya_lib_close();

    /* No clean up needed. */

    return 1;

  }

  pobyso_set_canonical_on();

  lowerBoundSo = pobyso_parse_string(argv[2]);

  if (lowerBoundSo == NULL)

  {

    fprintf(stderr,

        "\n\n%s: can't parse \"%s\" to a valid expression. Aborting!\n\n",

        argv[0], argv[2]);

    sollya_lib_close();

    return 1;

  }

  upperBoundSo = pobyso_parse_string(argv[3]);

  if (upperBoundSo == NULL)

  {

    fprintf(stderr,

        "\n\n%s: can't parse \"%s\" to a valid expression. Aborting!\n\n",

        argv[0], argv[3]);

    sollya_lib_clear_obj(lowerBoundSo);

    sollya_lib_close();

    return 1;

  }

  if (! pobyso_is_constant_expression(lowerBoundSo) ||

      ! pobyso_is_constant_expression(upperBoundSo))

  {

    fprintf(stderr,

        "\n\n%s: can't evaluate \"%s\" nor  \"%s\" to a constant expression. Aborting!\n\n",

        argv[0], argv[2], argv[3]);

    sollya_lib_clear_obj(lowerBoundSo);

    sollya_lib_clear_obj(upperBoundSo);

    sollya_lib_close();

    return 1;

  }

  mpfr_init2(lowerBoundMp, 512);

  if (! sollya_lib_get_constant(lowerBoundMp, lowerBoundSo))

  {

    fprintf(stderr,

        "\n\n%s: can't recover the value of \"%s\". Aborting!\n\n",

        argv[0], argv[2]);

    mpfr_clear(lowerBoundMp);

    sollya_lib_clear_obj(lowerBoundSo);

    sollya_lib_clear_obj(upperBoundSo);

    sollya_lib_close();

    return 1;

  }

  mpfr_init2(upperBoundMp, 512);

  if (! sollya_lib_get_constant(upperBoundMp, upperBoundSo))

  {

    fprintf(stderr,

        "\n\n%s: can't recover the value of \"%s\". Aborting!\n\n",

        argv[0], argv[3]);

    mpfr_clear(lowerBoundMp);

    mpfr_clear(upperBoundMp);

    sollya_lib_clear_obj(lowerBoundSo);

    sollya_lib_clear_obj(upperBoundSo);

    sollya_lib_close();

    return 1;

  }

  if (argc > 4)

  {

    fprintf(stdout, "Bounds: [");

    mpfr_out_str(stdout, 10, 0, lowerBoundMp, MPFR_RNDN);

    fprintf(stdout, ", ");

    mpfr_out_str(stdout, 10, 0, upperBoundMp, MPFR_RNDN);

    fprintf(stdout, "]\n");

  }

  /* Not needed any more. */

  sollya_lib_clear_obj(lowerBoundSo);

  sollya_lib_clear_obj(upperBoundSo);

  funcExpSo = pobyso_parse_string(argv[1]);

  if (funcExpSo == NULL)

  {

    fprintf(stderr,

        "\n\n%s: could not parse \"%s\" as a valid Sollya expression. Aborting!\n\n",

        argv[0], argv[1]);

    mpfr_clear(lowerBoundMp);

    mpfr_clear(upperBoundMp);

    sollya_lib_close();

    return 1;

  }

  if (argc > 4)

  {

    fprintf(stdout, "Functional expression: ");

    pobyso_autoprint(funcExpSo);

  }

  zerosArrayMp = pobyso_dirty_find_zeros_bounds(funcExpSo,

                                              lowerBoundMp,

                                              upperBoundMp,

                                              &zerosCount);

  sollya_lib_clear_obj(funcExpSo);

  mpfr_clear(lowerBoundMp);

  mpfr_clear(upperBoundMp);

  if (zerosCount < 0)

  {

    return 1;

  }

  if (zerosCount == 0)

  {

    return 0;

  }

  if (zerosCount > 0)

  {

    fprintf(stdout, "Zeroes: ");

    for (i = 0 ; i < zerosCount ; i++)

    {

      mpfr_out_str(stdout, 10, 0, zerosArrayMp[i], MPFR_RNDN);

      mpfr_clear(zerosArrayMp[i]);

      if (i < zerosCount-1) fprintf(stdout, ", ");

    }

    fprintf(stdout, ".\n");

  }

  else

  {

    fprintf(stdout, "No zeroes!\n");

  }

  if (zerosCount > 0) free(zerosArrayMp);

  sollya_lib_close();

  return 0;

} /* End main */