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 #ifndef UVES_UTILS_POLYNOMIAL_H
00052 #define UVES_UTILS_POLYNOMIAL_H
00053
00054
00055
00056
00057
00058 #include <uves_propertylist.h>
00059 #include <cpl.h>
00060
00061
00062
00063
00064
00065 typedef struct _polynomial polynomial ;
00066
00067
00068
00069
00070
00071
00072 polynomial *uves_polynomial_new(const cpl_polynomial *pol);
00073 polynomial *uves_polynomial_new_zero(int dim);
00074 polynomial *uves_polynomial_duplicate(const polynomial *p);
00075 polynomial *uves_polynomial_convert_from_plist_midas(const uves_propertylist *plist,
00076 const char *regression_name);
00077 polynomial *uves_polynomial_convert_from_table(cpl_table *t);
00078 polynomial *uves_polynomial_collapse(const polynomial *p, int varno, double value);
00079 polynomial * uves_polynomial_fit_1d(const cpl_vector * x_pos,
00080 const cpl_vector * values,
00081 const cpl_vector * sigmas,
00082 int poly_deg,
00083 double * mse);
00084 polynomial *uves_polynomial_fit_2d(const cpl_bivector * xy_pos,
00085 const cpl_vector * values,
00086 const cpl_vector * sigmas,
00087 int poly_deg1,
00088 int poly_deg2,
00089 double * mse,
00090 double * red_chisq,
00091 polynomial ** variance);
00092
00093 polynomial *uves_polynomial_add_2d(const polynomial *p1, const polynomial *p2);
00094
00095 int uves_polynomial_get_degree(const polynomial *p);
00096
00097 void uves_polynomial_delete(polynomial **p);
00098 void uves_polynomial_delete_const(const polynomial **p);
00099
00100 cpl_table *uves_polynomial_convert_to_table(const polynomial *p);
00101 int uves_polynomial_get_dimension(const polynomial *p);
00102 void uves_polynomial_dump(const polynomial *p, FILE *stream);
00103 cpl_error_code uves_polynomial_shift(polynomial *p, int varno, double shift);
00104 cpl_error_code uves_polynomial_rescale(polynomial *p, int varno, double scale);
00105 double uves_polynomial_get_coeff_1d(const polynomial *p, int degree);
00106 double uves_polynomial_get_coeff_2d(const polynomial *p,
00107 int degree1, int degree2);
00108 double uves_polynomial_evaluate_1d(const polynomial *p, double x);
00109 double uves_polynomial_evaluate_2d(const polynomial *p, double x1, double x2);
00110 double uves_polynomial_solve_1d(const polynomial *p, double value,
00111 double guess, int multiplicity);
00112 double uves_polynomial_solve_2d(const polynomial *p, double value,
00113 double guess, int multiplicity,
00114 int varno, double x_value);
00115 double uves_polynomial_derivative_1d(const polynomial *p, double x);
00116 double uves_polynomial_derivative_2d(const polynomial *p, double x1,
00117 double x2, int varno);
00118 cpl_error_code uves_polynomial_derivative(polynomial *p, int varno);
00119
00120 #endif