monomial-poly.h

Go to the documentation of this file.

/*
 * This file is part of MPSolve 3.2.2
 *
 * Copyright (C) 2001-2020, Dipartimento di Matematica "L. Tonelli", Pisa.
 * License: http://www.gnu.org/licenses/gpl.html GPL version 3 or higher
 *
 * Authors:
 *   Leonardo Robol <leonardo.robol@unipi.it>
 */
 
#ifndef MPS_MONOMIAL_POLY_H_
#define MPS_MONOMIAL_POLY_H_

 
 
  #include <mps/polynomial.h>
  #include <mps/mps.h>
  #include <gmp.h>
  #include <pthread.h>
 
#define MPS_MONOMIAL_POLY(t) (MPS_POLYNOMIAL_CAST (mps_monomial_poly, t))
#define MPS_IS_MONOMIAL_POLY(t) (mps_polynomial_check_type (t, "mps_monomial_poly"))
 
MPS_BEGIN_DECLS
 
#ifdef _MPS_PRIVATE
 
struct mps_monomial_poly_double_buffer {
  char active;
  mpc_t *mfpc1;
  mpc_t *mfpc2;
};
 

struct mps_monomial_poly {
  struct mps_polynomial methods;
 
  struct mps_monomial_poly_double_buffer db;

  mps_boolean *spar;

  double *fpr;

  cplx_t *fpc;

  cplx_t *fppc;

  rdpe_t *dpr;

  cdpe_t *dpc;

  mpf_t *mfpr;

  mpc_t *mfpc;

  pthread_mutex_t * mfpc_mutex;

  mpc_t *mfppc;

  double *fap;

  rdpe_t *dap;

  mpq_t *initial_mqp_r;

  mpq_t *initial_mqp_i;

  pthread_mutex_t regenerating;

  long int prec;
};
#endif /* #ifdef _MPS_PRIVATE */
 
/* These routines are thought for polynomial handling, i.e. allocating and
 * setting coefficients of the polynomials, and setting the precision of the
 * floating point coefficients that are in there */
 
mps_monomial_poly * mps_monomial_poly_new (mps_context * s, long int degree);
 
void mps_monomial_poly_free (mps_context * s, mps_polynomial * mp);
 
long int mps_monomial_poly_get_precision (mps_context * s, mps_monomial_poly * mp);
 
long int mps_monomial_poly_raise_precision (mps_context * s, mps_polynomial * mp, long int prec);
 
void mps_monomial_poly_set_coefficient_q (mps_context * s, mps_monomial_poly * mp, long int i,
                                          mpq_t real_part, mpq_t imag_part);
void mps_monomial_poly_set_coefficient_d (mps_context * s, mps_monomial_poly * mp, long int i,
                                          double real_part, double imag_part);
void mps_monomial_poly_set_coefficient_f (mps_context * s, mps_monomial_poly * p, long int i,
                                          mpc_t coeff);
void mps_monomial_poly_set_coefficient_int (mps_context * s, mps_monomial_poly * mp, long int i,
                                            long long real_part, long long imag_part);
 
void mps_monomial_poly_set_coefficient_s (mps_context * s, mps_monomial_poly * p,
                                          int i, const char * real_coeff,
                                          const char * imag_coeff);
 
void mps_monomial_poly_get_coefficient_d (mps_context * s, mps_monomial_poly * p,
                                          int i, cplx_t output);
 
void mps_monomial_poly_get_coefficient_q (mps_context * s, mps_monomial_poly * p,
                                          int i, mpq_t real_output, mpq_t imag_output);
 
mps_monomial_poly * mps_monomial_poly_derive (mps_context * s, mps_monomial_poly * p, int k, long int wp);
 
mps_boolean mps_monomial_poly_feval (mps_context * ctx, mps_polynomial *p, cplx_t x, cplx_t value, double * error);
 
mps_boolean mps_monomial_poly_deval (mps_context * ctx, mps_polynomial *p, cdpe_t x, cdpe_t value, rdpe_t error);
 
mps_boolean mps_monomial_poly_meval (mps_context * ctx, mps_polynomial *p, mpc_t x, mpc_t value, rdpe_t error);
 
void mps_monomial_poly_fstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
void mps_monomial_poly_dstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
void mps_monomial_poly_mstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
void mps_monomial_poly_fnewton (mps_context * ctx, mps_polynomial * p,
                                mps_approximation * root, cplx_t corr);
 
void mps_monomial_poly_dnewton (mps_context * ctx, mps_polynomial * p,
                                mps_approximation * root, cdpe_t corr);
 
void mps_monomial_poly_mnewton (mps_context * ctx, mps_polynomial * p,
                                mps_approximation * root, mpc_t corr, long int wp);
 
void mps_monomial_poly_get_leading_coefficient (mps_context * ctx, mps_polynomial * p,
                                                mpc_t leading_coefficient);
 
void mps_monomial_poly_deflate (mps_context * ctx, mps_polynomial * p);
 
mps_monomial_poly * mps_monomial_poly_read_from_stream (mps_context * s, mps_input_buffer * buffer,
                                                        mps_structure structure, mps_density density,
                                                        long int precision);
 
MPS_END_DECLS
 
#endif