1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
/*
* 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_CHEBYSHEV_H_
#define MPS_CHEBYSHEV_H_
MPS_BEGIN_DECLS
#define MPS_CHEBYSHEV_POLY_TYPE_NAME "mps_chebyshev_poly"
#define MPS_CHEBYSHEV_POLY(t) ((mps_chebyshev_poly*)t)
#define MPS_IS_CHEBYSHEV_POLY(t) mps_polynomial_check_type (t, "mps_chebyshev_poly")
typedef struct {
/**
* @brief Base implementation of a polynomial.
*/
mps_polynomial super;
/**
* @brief Floating point coefficients of the polynomial in the Chebyshev
* base
*/
cplx_t * fpc;
/**
* @brief DPE floating point coefficients of the polynomial in the Chebyshev
* base.
*/
cdpe_t * dpc;
/**
* @brief Multiprecision complex coefficients of the polynomial in the Chebyshev
* base.
*/
mpc_t * mfpc;
/**
* @brief Rational coefficients of the polynomial. These are the real parts of the
* coefficients.
*/
mpq_t * rational_real_coeffs;
/**
* @brief Ratinonal coefficients of the polynomial. These are the imaginary parts
* of the coefficients.
*/
mpq_t * rational_imag_coeffs;
/**
* @brief Leading coefficient of the polynomial.
*/
mpc_t lc;
/**
* @brief Internal mutex used to manage the change of precision.
*/
pthread_mutex_t precision_mutex;
} mps_chebyshev_poly;
/**
* @brief Create a new polynomial represented in the Chebyshev base with
* degree set to n.
*/
mps_chebyshev_poly * mps_chebyshev_poly_new (mps_context * ctx, int n, mps_structure structure);
/**
* @brief Set the coefficient relative to the i-th element of the Chebyshev
* base.
*
* This function takes a rational number as input, and is usable only if the
* Chebyshev polynomial is represented using rational coefficients.
*/
void mps_chebyshev_poly_set_coefficient_q (mps_context * ctx, mps_chebyshev_poly * poly, int i,
mpq_t real_part, mpq_t imag_part);
/**
* @brief Set the coefficient relative to the i-th element of the Chebyshev
* base.
*
* This function takes a multiprecision floating point number as input.
*/
void mps_chebyshev_poly_set_coefficient_f (mps_context * ctx, mps_chebyshev_poly * poly,
int i, mpc_t coeff);
/**
* @brief Set the coefficient of the i-th element of the Chebyshev base.
*
* This function takes an integer value as input.
*
* @param ctx The current mps_context
* @param poly The Chebyshev polynomial whose coefficient should be set.
* @param i The degree of the coefficient to set.
* @param real_coeff The real part of the new coefficient.
* @param imag_coeff The imaginary part of the new coefficient.
*/
void mps_chebyshev_poly_set_coefficient_i (mps_context * ctx, mps_chebyshev_poly * poly,
int i, long int real_coeff, long int imag_coeff);
mps_chebyshev_poly * mps_chebyshev_poly_read_from_stream (mps_context * ctx, mps_input_buffer * buffer,
mps_structure structure, mps_density density,
long int precision);
MPS_END_DECLS
#endif
|
