monomial-matrix-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_MATRIX_POLY_H_
#define MPS_MONOMIAL_MATRIX_POLY_H_

 
#include <mps/polynomial.h>
 
#define MPS_MONOMIAL_MATRIX_POLY(t) (MPS_POLYNOMIAL_CAST (mps_monomial_matrix_poly, t))
#define MPS_IS_MONOMIAL_MATRIX_POLY(t) \
  (mps_polynomial_check_type (t, "mps_monomial_matrix_poly"))
 
#ifdef  __cplusplus
extern "C"
{
#endif
 
#define MPS_MONOMIAL_MATRIX_POLY_HESSENBERG 0x0001
 
#ifdef _MPS_PRIVATE

struct mps_monomial_matrix_poly {
  mps_polynomial methods;

  mps_boolean monic;

  int m;

  int degree;

  cplx_t * P;

  mpc_t * mP;

  mpq_t * mpqPr;

  mpq_t * mpqPi;

  int flags;
};
 
 
#endif /* _MPS_PRIVATE */
 
/* From here on you will find declaration of the public method
   available for the matrix polynomials */

mps_monomial_matrix_poly* mps_monomial_matrix_poly_new (mps_context * ctx,
                                                        int degree,
                                                        int m,
                                                        mps_boolean monic);
void mps_monomial_matrix_poly_free (mps_context * ctx,
                                    mps_polynomial * poly);

void mps_monomial_matrix_poly_add_flags (mps_context * ctx,
                                         mps_monomial_matrix_poly * mpoly,
                                         int flag);

void mps_monomial_matrix_poly_clear_flags (mps_context * ctx,
                                           mps_monomial_matrix_poly * mpoly,
                                           int flag);

void mps_monomial_matrix_poly_set_coefficient_d (mps_context * ctx,
                                                 mps_monomial_matrix_poly *mpoly,
                                                 int i,
                                                 cplx_t * matrix);

void mps_monomial_matrix_poly_set_coefficient_q (mps_context * ctx,
                                                 mps_monomial_matrix_poly *mpoly,
                                                 int i,
                                                 mpq_t * matrix_r,
                                                 mpq_t * matrix_i);

mps_boolean mps_monomial_matrix_poly_meval (mps_context * ctx,
                                            mps_polynomial * poly,
                                            mpc_t x,
                                            mpc_t value,
                                            rdpe_t error);

long int
mps_monomial_matrix_poly_raise_data (mps_context * ctx,
                                     mps_polynomial * p,
                                     long int wp);
 
#ifdef  __cplusplus
}
#endif
 
 
#endif /* __MPS_MONOMIAL_MATRIX_POLY_H */