//
//  Interp_Funct.cc
//  Brian Magill
//  10/9/2003
//
//  For interpolating an entire curve from a series of
//  curves.  Uses quadratic interpolation.
//
//----------------------------------------------------------------------
//

#include <gsl/gsl_multifit.h>
#include <gsl/gsl_fit.h>
#include <stdio.h>
#include <Interp_Funct.h>


int linear_fit(const double * x, const double * y, long n, 
                double * intercept, double * slope, double * cov00, double * cov01, 
                double * cov11, double * sumsq)
    {
  
  
    if( gsl_fit_linear (x, 1, y, 1, n, intercept, slope, 
                            cov00, cov01, cov11, sumsq) )
        {
        
        fprintf(stderr, "linear_fit: Problems with linear least squares fit!\n");        
        return -100;
        
        }

    return 0;
    }
    
int linear_fit_est(double x, double c0, double c1, double c00, double c01, 
                    double c11, double *y, double *y_err)
    {                    

    if( gsl_fit_linear_est (x, c0, c1, c00, c01, c11, y, y_err)  )
        {  
        
        fprintf(stderr, "linear_fit_est: Problems with linear least squares fit!\n");        
        return -100;
        
        }
    
    return 0;
    }