operaFit.h

Go to the documentation of this file.

#ifndef OPERAFIT_H
#define OPERAFIT_H

/*******************************************************************
 ****                LIBRARY FOR OPERA v1.0                     ****
 *******************************************************************
 Library name: operaFit
 Version: 1.0
 Description: This C library implements LM curve fitting for commonly 
 used functions such as polynomials, gaussians, etc.  
 Author(s): CFHT OPERA team
 Affiliation: Canada France Hawaii Telescope 
 Location: Hawaii USA
 Date: Aug/2011
 Contact: eder@cfht.hawaii.edu
 
 Copyright (C) 2011  Opera Pipeline team, Canada France Hawaii Telescope
 
 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, see:
 http://software.cfht.hawaii.edu/licenses
 -or-
 http://www.gnu.org/licenses/gpl-3.0.html
 ********************************************************************/

// $Date$
// $Id$
// $Revision$
// $Locker$
// $Log$


/*
 * The functions below uses operaLMFit library (LMFIT)
 */
void operaLMFitPolynomial(unsigned m_dat, double *x, double *y, int n_par, double par[], double *chi2, int verbose);
void operaLMFitGaussian(unsigned m_dat, double *x, double *y, double *a, double *x0, double *sig, double *chi2, int verbose);
void operaLMFit2DPolynomial(unsigned m_dat, double *x, double *y, double *fxy, int n_par, double par[], double *chi2, int verbose);
void operaLMFit2DGaussian(unsigned m_dat, double *x, double *y, double *fxy, double *a, double *x0, double *y0, double *sigx, double *sigy, double *chi2, int verbose);

/*
 * The functions below uses nr library, although all necessary functions to make them work are in operaFit.c
 */
// spline and 2D-spline interpolation and related functions
void operaFitSpline(unsigned nin, float *xin,float *yin, unsigned nout, float *xout, float *yout);
void operaFit2DSpline(unsigned nxin, float *xin, unsigned nyin, float *yin,double *fxyin, unsigned nxout, float *xout, unsigned nyout, float *yout, float *fxyout);
void spline(float x[], float y[], int n, float yp1, float ypn, float y2[]);
void splint(float xa[], float ya[], float y2a[], int n, float x, float *y);
void splin2(float x1a[], float x2a[], float **ya, float **y2a, int m, int n,float x1, float x2, float *y);
void splie2(float x1a[], float x2a[], float **ya, int m, int n, float **y2a);

// medfit and related functions
void operaMedfit(float xin[], float yin[], int ndata, float *a, float *b, float *abdev);
float rofunc(float b);
float select_nr(unsigned long k, unsigned long n, float arr[]);

// common functions to spline and medfit
float **matrix(long nrl, long nrh, long ncl, long nch);
void free_matrix(float **m, long nrl, long nrh, long ncl, long nch);
float *vector(long nl, long nh);
void free_vector(float *v, long nl, long nh);
void nrerror(char error_text[]);

#endif