00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057 #ifndef UVES_UTILS_POLYNOMIAL_H
00058 #define UVES_UTILS_POLYNOMIAL_H
00059
00060
00061
00062
00063
00064 #include <uves_propertylist.h>
00065 #include <cpl.h>
00066
00067
00068
00069
00070
00071 typedef struct _polynomial polynomial ;
00072
00073
00074
00075
00076
00077
00078 polynomial *uves_polynomial_new(const cpl_polynomial *pol);
00079 polynomial *uves_polynomial_new_zero(int dim);
00080 polynomial *uves_polynomial_duplicate(const polynomial *p);
00081
00082
00083 polynomial *uves_polynomial_convert_from_table(cpl_table *t);
00084 polynomial *uves_polynomial_collapse(const polynomial *p, int varno, double value);
00085 polynomial * uves_polynomial_fit_1d(const cpl_vector * x_pos,
00086 const cpl_vector * values,
00087 const cpl_vector * sigmas,
00088 int poly_deg,
00089 double * mse);
00090 polynomial *uves_polynomial_fit_2d(const cpl_bivector * xy_pos,
00091 const cpl_vector * values,
00092 const cpl_vector * sigmas,
00093 int poly_deg1,
00094 int poly_deg2,
00095 double * mse,
00096 double * red_chisq,
00097 polynomial ** variance);
00098
00099 polynomial *uves_polynomial_add_2d(const polynomial *p1, const polynomial *p2);
00100
00101 int uves_polynomial_get_degree(const polynomial *p);
00102
00103 void uves_polynomial_delete(polynomial **p);
00104 void uves_polynomial_delete_const(const polynomial **p);
00105
00106 cpl_table *uves_polynomial_convert_to_table(const polynomial *p);
00107 int uves_polynomial_get_dimension(const polynomial *p);
00108 void uves_polynomial_dump(const polynomial *p, FILE *stream);
00109 cpl_error_code uves_polynomial_shift(polynomial *p, int varno, double shift);
00110 cpl_error_code uves_polynomial_rescale(polynomial *p, int varno, double scale);
00111 double uves_polynomial_get_coeff_1d(const polynomial *p, int degree);
00112 double uves_polynomial_get_coeff_2d(const polynomial *p,
00113 int degree1, int degree2);
00114 double uves_polynomial_evaluate_1d(const polynomial *p, double x);
00115 double uves_polynomial_evaluate_2d(const polynomial *p, double x1, double x2);
00116 double uves_polynomial_solve_1d(const polynomial *p, double value,
00117 double guess, int multiplicity);
00118 double uves_polynomial_solve_2d(const polynomial *p, double value,
00119 double guess, int multiplicity,
00120 int varno, double x_value);
00121 double uves_polynomial_derivative_1d(const polynomial *p, double x);
00122 double uves_polynomial_derivative_2d(const polynomial *p, double x1,
00123 double x2, int varno);
00124 cpl_error_code uves_polynomial_derivative(polynomial *p, int varno);
00125
00126 #endif