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

/** @file pobyso.h

 * Integration of Sollya to C programs

 * @author S.T.

 * @date 2011-10-11

 * Note: pobyso stands for POwered BY SOllya.

 *                         --      -- --

 */

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

/*

 * Add below all the headers needed to get this header work.

 */

/* <stdio.h> is needed *before* <mpfr.h> for all MPFR input/output functions

 * prototypes be defined. */

#include <string.h>

#include <stdio.h>

#include <sollya.h>

#include <mpfr.h>

#ifndef POBYSO_h

/* Typedefs to make code more readable. */

typedef sollya_obj_t pobyso_error_t;

typedef sollya_obj_t pobyso_func_exp_t;

typedef sollya_obj_t pobyso_precision_t;

typedef sollya_obj_t pobyso_range_t;

typedef sollya_obj_t sollya_verbosity_t;

#define POBYSO_ABSOLUTE 1

#define POBYSO_RELATIVE 2

/* Mimic the default behavior of interactive Sollya. */

#define POBYSO_DEFAULT_POINTS 501

#define POBYSO_INF_NORM_NUM_POINTS (POBYSO_DEFAULT_POINTS)

#define POBYSO_GUESS_DEGREE_BOUND 1024

/* Very thin wrappers around a lot of Sollya functions.

 */

static inline pobyso_error_t pobyso_error(void)

{return(sollya_lib_error());}

static inline int pobyso_is_error(sollya_obj_t errorCandidate)

{return(sollya_lib_obj_is_error(errorCandidate));}

static inline int pobyso_is_function(sollya_obj_t functionCandidate)

{return(sollya_lib_obj_is_function(functionCandidate));}

/**

 * Print object(s) to stdout.

 * A very thin wrapper around the lib_sollya_v_autoprint() function.

 * Should be called with NULL as the final argument.

 */

void

pobyso_autoprint(sollya_obj_t soObj, ...);

/**

 * Print object(s) to stdout: the va_list companion function.

 * A very thin wrapper around the lib_sollya_v_autoprint() function.

 * The last argument in the va_list should be NULL.

 */

void

pobyso_autoprint_v(sollya_obj_t soObj, va_list va);

/**

 * Parse a string to create a Sollya object.

 * A very thin wrapper around the sollya_lib_parse_string() function.

 * If the final ";" is forgotten in the expression, it is added by the

 * function.

 * @return a Sollya expression if successes, a Sollya error if fails.

 */

pobyso_error_t

pobyso_parse(const char* expression);

/**

 * A wrapper around the Sollya Remez function.

 */

/**

 * Parse a string to create a Sollya object.

 * A very thin wrapper around the sollya_lib_parse_string() function.

 * If the final ";" is forgotten in the expression, it is added by the

 * function.

 * @return a Sollya functional expression if successes, a Sollya error if fails.

 */

pobyso_error_t

pobyso_parse_string(const char* expression);

/**

 * A wrapper around the Sollya Remez function with the canonical monomials

 * base.

 */

pobyso_func_exp_t

pobyso_remez_canonical_monomials_base(pobyso_func_exp_t function,

                                      long int degree,

                                      pobyso_range_t interval,

                                      pobyso_func_exp_t weight,

                                      double quality,

                                      pobyso_range_t bounds);

/**

 * A wrapper around the Sollya Remez function with the a sparse monomials

 * base.

 */

pobyso_func_exp_t

pobyso_remez_sparse_monomials_base(pobyso_func_exp_t);

/**

 * A wrapper around the Sollya Remez function with the canonical monomials

 * base.

 */

pobyso_func_exp_t

pobyso_remez_arbitrary_base(pobyso_func_exp_t);

/**

 * Set the verbosity mode off (level 0).

 * The current level of vebositiy is returned.

 */

sollya_verbosity_t

pobyso_set_verbosity_off(void);

/**

 * Set the verbosity level to newVerbosityLevel.

 * @param newVerbosityLevel must be a Sollay object corresponding to an

 *        integer constant.

 */

sollya_verbosity_t

pobyso_set_verbosity_to(sollya_obj_t newVerbosityLevel);

#if 0

/**

 * Create the canonical (non sparse) base of monomials for a given degree.

 */

chain*

pobyso_create_canonical_monomials_base(const unsigned int degree);

/**

 * Create a chain from an array of int.

 */

sollya_int_list

pobyso_create_int_list_from_int_Array(int* intArray,

                                    const unsigned int arrayLength);

/**

 * Create a chain from an array of int.

 */

sollya_int_list_t

pobyso_create_int_list_from_unsigned_int_array(unsigned int* intArray,

                                            const unsigned int arrayLength);

/**

 * Differentiation of a function.

 * A slim wrapper around the Sollya differentiate function.

 * @param functionNode - the Sollya node to differentiate;

 * @return a node representing the function differentiated or NULL, if

 *         something goes wrong.

 */

sollya_obj_t

pobyso_diff(sollya_obj_t function);

/**

 * A match to the Sollya dirtyinfnorm.

 * A slim wrapper around the Sollya function.

 * @param infnorm - out parameter to return the result, must be "inited"

 *                  and "cleared" by the caller;

 * @param functionNode - the Sollya node to compute the infinite norm of;

 * @param lowerBound - the lower bound of the interval;

 * @param upperBound - the upper bound of the interval;

 * @param precision  - the internal precision Sollya must use.

 * @return 0 if everything is OK, != 0 if something goes wrong.

 */

int

pobyso_dirty_infnorm(mpfr_t infNorm,

                      node *functionNode,

                      mpfr_t lowerBound,

                      mpfr_t upperBound,

                      mp_prec_t precision);

/**

 * Faithful evaluation of an expression.

 *

 * @param faitufulEvaluation - holds the result, must be "inited" by the

 *                             caller;

 */

int

pobyso_evaluate_faithful(mpfr_t faithfulEvaluation,

                          node *nodeToEvaluate,

                          mpfr_t argument,

                          mpfr_prec_t precision);

/**

 * Find the zeros of a function on a given interval.

 */

chain*

pobyso_find_zeros(node *function,

                  mpfr_t *lowerBound,

                  mpfr_t *upperBound);

/**

 * Free a chain of node.

 * All elements of the chain have to be nodes.

 */

void

pobyso_free_chain_of_nodes(chain *theChain);

/**

 * Free a range.

 * It involves clearing the mpfr_t elements and deallocating the

 * pointers.

 */

void

pobyso_free_range(rangetype range);

/**

 * Computes a good polynomial approximation with fixed-point or floating-point

 * coefficients.

 */

node*

pobyso_fp_minimax_canonical_monomials_base(node *function,

                                          int degree,

                                          chain *formats,

                                          chain *points,

                                          mpfr_t lowerBound,

                                          mpfr_t upperBound,

                                          int fpFixedArg,

                                          int absRel,

                                          node *constPart,

                                          node *minimax);

/**

 * Parses a string to build the node representing the function.

 * In fact, does nothing for the moment: the string must be a correct function

 * definition. No error correction takes place here.

 */

node*

pobyso_parse_function(char *functionString,

                      char *freeVariableNameString);

/** Compute a polynomial approximation in the canonical monomials basis for

 *  a function, for a given precision. The returned polynomial has the minimal

 *  degree to achieve the required precision.

 */

node*

pobyso_remez_approx_canonical_monomials_base_for_error(node *functionNode,

                                                      unsigned int mode,

                                                      mpfr_t lowerBound,

                                                      mpfr_t upperBound,

                                                      mpfr_t eps);

/**

 * Computes a the remez approximation of a function.

 * @param function   - the node holding the function to approximate;

 * @param weight     - the node holding the weight function, can be NULL. In

 *                     this case a default weight of "1" will be provided and

 *                     the approximation is related is with respect to the

 *                     absolute error.

 * @param degree     - the degree of the approximation polynomial;

 * @param lowerBound - the lower bound of the approximation interval;

 * @param upperBound - the upper bound of the approximation interval;

 * @param quality    - quality = (eps - eps*) / eps*; the search stop when the required

 *                     quality is achieved.

 *

 * @return a node holding the approximation polynomial in Horner form.

 */

node*

pobyso_remez_canonical_monomials_base(node *function,

                                       node *weight,

                                       unsigned int degree,

                                       mpfr_t lowerBound,

                                       mpfr_t upperBound,

                                       mpfr_t quality);

#endif

#define POBYSO_h

#endif