issue #1

include/gsl/gsl_deriv.h

@@ -0,0 +1,50 @@
+/* deriv/gsl_deriv.h
+ * 
+ * Copyright (C) 2000 David Morrison
+ * 
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 3 of the License, or (at
+ * your option) any later version.
+ * 
+ * This program is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+ * General Public License for more details.
+ * 
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
+ */
+
+#ifndef __GSL_DERIV_H__
+#define __GSL_DERIV_H__
+#include <gsl/gsl_math.h>
+
+#undef __BEGIN_DECLS
+#undef __END_DECLS
+#ifdef __cplusplus
+# define __BEGIN_DECLS extern "C" {
+# define __END_DECLS }
+#else
+# define __BEGIN_DECLS /* empty */
+# define __END_DECLS /* empty */
+#endif
+
+__BEGIN_DECLS
+
+int gsl_deriv_central (const gsl_function *f,
+                       double x, double h,
+                       double *result, double *abserr);
+
+int gsl_deriv_backward (const gsl_function *f,
+                        double x, double h,
+                        double *result, double *abserr);
+
+int gsl_deriv_forward (const gsl_function *f,
+                       double x, double h,
+                       double *result, double *abserr);
+
+__END_DECLS
+
+#endif /* __GSL_DERIV_H__ */

include/gsl/gsl_fit.h

@@ -0,0 +1,85 @@
+/* fit/gsl_fit.h
+ * 
+ * Copyright (C) 2000, 2007 Brian Gough
+ * 
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 3 of the License, or (at
+ * your option) any later version.
+ * 
+ * This program is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+ * General Public License for more details.
+ * 
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
+ */
+
+#ifndef __GSL_FIT_H__
+#define __GSL_FIT_H__
+
+#include <stdlib.h>
+#include <gsl/gsl_math.h>
+
+#undef __BEGIN_DECLS
+#undef __END_DECLS
+#ifdef __cplusplus
+# define __BEGIN_DECLS extern "C" {
+# define __END_DECLS }
+#else
+# define __BEGIN_DECLS /* empty */
+# define __END_DECLS /* empty */
+#endif
+
+__BEGIN_DECLS
+
+int gsl_fit_linear (const double * x, const size_t xstride,
+                    const double * y, const size_t ystride,
+                    const size_t n,
+                    double * c0, double * c1, 
+                    double * cov00, double * cov01, double * cov11, 
+                    double * sumsq);
+
+
+int gsl_fit_wlinear (const double * x, const size_t xstride,
+                     const double * w, const size_t wstride,
+                     const double * y, const size_t ystride,
+                     const size_t n,
+                     double * c0, double * c1, 
+                     double * cov00, double * cov01, double * cov11, 
+                     double * chisq);
+
+int
+gsl_fit_linear_est (const double x, 
+                    const double c0, const double c1, 
+                    const double cov00, const double cov01, const double cov11,
+                    double *y, double *y_err);
+
+
+int gsl_fit_mul (const double * x, const size_t xstride,
+                 const double * y, const size_t ystride,
+                 const size_t n,
+                 double * c1, 
+                 double * cov11, 
+                 double * sumsq);
+
+int gsl_fit_wmul (const double * x, const size_t xstride,
+                  const double * w, const size_t wstride,
+                  const double * y, const size_t ystride,
+                  const size_t n,
+                  double * c1, 
+                  double * cov11, 
+                  double * sumsq);
+
+
+int
+gsl_fit_mul_est (const double x, 
+                 const double c1, 
+                 const double cov11,
+                 double *y, double *y_err);
+
+__END_DECLS
+
+#endif /* __GSL_FIT_H__ */

include/gsl/gsl_machine.h

@@ -0,0 +1,104 @@
+/* Author:  B. Gough and G. Jungman */
+#ifndef __GSL_MACHINE_H__
+#define __GSL_MACHINE_H__
+
+#include <limits.h>
+#include <float.h>
+
+/* magic constants; mostly for the benefit of the implementation */
+
+/* -*-MACHINE CONSTANTS-*-
+ *
+ * PLATFORM: Whiz-O-Matic 9000
+ * FP_PLATFORM: IEEE-Virtual
+ * HOSTNAME: nnn.lanl.gov
+ * DATE: Fri Nov 20 17:53:26 MST 1998
+ */
+#define GSL_DBL_EPSILON        2.2204460492503131e-16
+#define GSL_SQRT_DBL_EPSILON   1.4901161193847656e-08
+#define GSL_ROOT3_DBL_EPSILON  6.0554544523933429e-06
+#define GSL_ROOT4_DBL_EPSILON  1.2207031250000000e-04
+#define GSL_ROOT5_DBL_EPSILON  7.4009597974140505e-04
+#define GSL_ROOT6_DBL_EPSILON  2.4607833005759251e-03
+#define GSL_LOG_DBL_EPSILON   (-3.6043653389117154e+01)
+
+#define GSL_DBL_MIN        2.2250738585072014e-308
+#define GSL_SQRT_DBL_MIN   1.4916681462400413e-154
+#define GSL_ROOT3_DBL_MIN  2.8126442852362996e-103
+#define GSL_ROOT4_DBL_MIN  1.2213386697554620e-77
+#define GSL_ROOT5_DBL_MIN  2.9476022969691763e-62
+#define GSL_ROOT6_DBL_MIN  5.3034368905798218e-52
+#define GSL_LOG_DBL_MIN   (-7.0839641853226408e+02)
+
+#define GSL_DBL_MAX        1.7976931348623157e+308
+#define GSL_SQRT_DBL_MAX   1.3407807929942596e+154
+#define GSL_ROOT3_DBL_MAX  5.6438030941222897e+102
+#define GSL_ROOT4_DBL_MAX  1.1579208923731620e+77
+#define GSL_ROOT5_DBL_MAX  4.4765466227572707e+61
+#define GSL_ROOT6_DBL_MAX  2.3756689782295612e+51
+#define GSL_LOG_DBL_MAX    7.0978271289338397e+02
+
+#define GSL_FLT_EPSILON        1.1920928955078125e-07
+#define GSL_SQRT_FLT_EPSILON   3.4526698300124393e-04
+#define GSL_ROOT3_FLT_EPSILON  4.9215666011518501e-03
+#define GSL_ROOT4_FLT_EPSILON  1.8581361171917516e-02
+#define GSL_ROOT5_FLT_EPSILON  4.1234622211652937e-02
+#define GSL_ROOT6_FLT_EPSILON  7.0153878019335827e-02
+#define GSL_LOG_FLT_EPSILON   (-1.5942385152878742e+01)
+
+#define GSL_FLT_MIN        1.1754943508222875e-38
+#define GSL_SQRT_FLT_MIN   1.0842021724855044e-19
+#define GSL_ROOT3_FLT_MIN  2.2737367544323241e-13
+#define GSL_ROOT4_FLT_MIN  3.2927225399135965e-10
+#define GSL_ROOT5_FLT_MIN  2.5944428542140822e-08
+#define GSL_ROOT6_FLT_MIN  4.7683715820312542e-07
+#define GSL_LOG_FLT_MIN   (-8.7336544750553102e+01)
+
+#define GSL_FLT_MAX        3.4028234663852886e+38
+#define GSL_SQRT_FLT_MAX   1.8446743523953730e+19
+#define GSL_ROOT3_FLT_MAX  6.9814635196223242e+12
+#define GSL_ROOT4_FLT_MAX  4.2949672319999986e+09
+#define GSL_ROOT5_FLT_MAX  5.0859007855960041e+07
+#define GSL_ROOT6_FLT_MAX  2.6422459233807749e+06
+#define GSL_LOG_FLT_MAX    8.8722839052068352e+01
+
+#define GSL_SFLT_EPSILON        4.8828125000000000e-04
+#define GSL_SQRT_SFLT_EPSILON   2.2097086912079612e-02
+#define GSL_ROOT3_SFLT_EPSILON  7.8745065618429588e-02
+#define GSL_ROOT4_SFLT_EPSILON  1.4865088937534013e-01
+#define GSL_ROOT5_SFLT_EPSILON  2.1763764082403100e-01
+#define GSL_ROOT6_SFLT_EPSILON  2.8061551207734325e-01
+#define GSL_LOG_SFLT_EPSILON   (-7.6246189861593985e+00)
+
+/* !MACHINE CONSTANTS! */
+
+
+/* a little internal backwards compatibility */
+#define GSL_MACH_EPS  GSL_DBL_EPSILON
+
+
+
+/* Here are the constants related to or derived from
+ * machine constants. These are not to be confused with
+ * the constants that define various precision levels
+ * for the precision/error system.
+ *
+ * This information is determined at configure time
+ * and is platform dependent. Edit at your own risk.
+ *
+ * PLATFORM: WHIZ-O-MATIC
+ * CONFIG-DATE: Thu Nov 19 19:27:18 MST 1998
+ * CONFIG-HOST: nnn.lanl.gov
+ */
+
+/* machine precision constants */
+/* #define GSL_MACH_EPS         1.0e-15 */
+#define GSL_SQRT_MACH_EPS       3.2e-08
+#define GSL_ROOT3_MACH_EPS      1.0e-05
+#define GSL_ROOT4_MACH_EPS      0.000178
+#define GSL_ROOT5_MACH_EPS      0.00100
+#define GSL_ROOT6_MACH_EPS      0.00316
+#define GSL_LOG_MACH_EPS       (-34.54)
+
+
+#endif /* __GSL_MACHINE_H__ */