1/* fit/gsl_fit.h
2 *
3 * Copyright (C) 2000, 2007 Brian Gough
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 3 of the License, or (at
8 * your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful, but
11 * WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 * General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
18 */
19
20#ifndef __GSL_FIT_H__
21#define __GSL_FIT_H__
22
23#include <stdlib.h>
24#include <gsl/gsl_math.h>
25
26#undef __BEGIN_DECLS
27#undef __END_DECLS
28#ifdef __cplusplus
29# define __BEGIN_DECLS extern "C" {
30# define __END_DECLS }
31#else
32# define __BEGIN_DECLS /* empty */
33# define __END_DECLS /* empty */
34#endif
35
36__BEGIN_DECLS
37
38int gsl_fit_linear (const double * x, const size_t xstride,
39 const double * y, const size_t ystride,
40 const size_t n,
41 double * c0, double * c1,
42 double * cov00, double * cov01, double * cov11,
43 double * sumsq);
44
45
46int gsl_fit_wlinear (const double * x, const size_t xstride,
47 const double * w, const size_t wstride,
48 const double * y, const size_t ystride,
49 const size_t n,
50 double * c0, double * c1,
51 double * cov00, double * cov01, double * cov11,
52 double * chisq);
53
54int
55gsl_fit_linear_est (const double x,
56 const double c0, const double c1,
57 const double cov00, const double cov01, const double cov11,
58 double *y, double *y_err);
59
60
61int gsl_fit_mul (const double * x, const size_t xstride,
62 const double * y, const size_t ystride,
63 const size_t n,
64 double * c1,
65 double * cov11,
66 double * sumsq);
67
68int gsl_fit_wmul (const double * x, const size_t xstride,
69 const double * w, const size_t wstride,
70 const double * y, const size_t ystride,
71 const size_t n,
72 double * c1,
73 double * cov11,
74 double * sumsq);
75
76
77int
78gsl_fit_mul_est (const double x,
79 const double c1,
80 const double cov11,
81 double *y, double *y_err);
82
83__END_DECLS
84
85#endif /* __GSL_FIT_H__ */
86