sc_mpfit.h

Go to the documentation of this file.

/*
 *  This file is part of the SKYCORR software package.
 *  Copyright (C) 2009-2013 European Southern Observatory
 *
 *  This programme 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 2 of the License, or
 *  (at your option) any later version.
 *
 *  This programme 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 programme. If not, see <http://www.gnu.org/licenses/>.
 */
 
 
#ifdef __cplusplus
extern "C" {
#endif
 
/*****************************************************************************
 *                                INCLUDES                                   *
 ****************************************************************************/
 
/* Config header */
 
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
 
/* Sky correction headers */
 
#include <sc_basic.h>
#include <sc_modsky.h>
 
/*****************************************************************************
 *                                DEFINES                                    *
 ****************************************************************************/
 
#ifndef SC_MPFIT_H
#define SC_MPFIT_H
 
/* Definition of constants */
 
#define SC_DEFERRVAL 9.99
#define SC_CORRFAC_MIN 0.01
#define SC_CORRFAC_MAX 100.
#define SC_MINNLIN 3
#define SC_MAXRELFAC 2.
#define SC_MAXPARERR 0.5
 
/*****************************************************************************
 *                                GLOBALS                                    *
 ****************************************************************************/
 
/* Declaration of global variables */
 
extern int nfev;
extern cpl_boolean lastcall;
 
/*****************************************************************************
 *                                TYPEDEF                                    *
 ****************************************************************************/
 
/* Definition of structures */
 
typedef struct _scvars_ {
    cpl_table *scispec;
    cpl_table *skyspec;
    cpl_table *fitpar;
    cpl_vector *sinc;
    const cpl_parameterlist *parlist;
} scvars;
 
typedef struct _scpars_ {
    int n;
    double *p;
    mp_par *pars;
} scpars;
 
/*****************************************************************************
 *                               PROTOTYPES                                  *
 ****************************************************************************/
 
/* Declaration of functions */
 
cpl_error_code sc_mpfit(mp_result *result, cpl_table *scispec,
                        cpl_table *skyspec, cpl_table *fitpar,
                        const cpl_parameterlist *parlist);
int sc_mpfit_calcdev(int m, int n, double *p, double *dy, double **dvec,
                     void *vars);
cpl_error_code sc_mpfit_modinitpar(cpl_table *fitpar,
                                   cpl_table *scispec,
                                   cpl_table *skyspec,
                                   const cpl_parameterlist *parlist);
cpl_error_code sc_mpfit_setpar(scpars *fitpars,
                               const cpl_table *fitpar, const char fittype);
cpl_error_code sc_mpfit_allocmempar(scpars *fitpars, const int npar);
cpl_error_code sc_mpfit_freemempar(scpars *fitpars);
cpl_error_code sc_mpfit_allocmemresult(mp_result *result, const int m,
                                       const int n);
cpl_error_code sc_mpfit_initresult(mp_result *result, const int m,
                                   const int n);
cpl_error_code sc_mpfit_copyresult(mp_result *outresult,
                                   const mp_result *inresult);
cpl_error_code sc_mpfit_freememresult(mp_result *result);
cpl_error_code sc_mpfit_substbadfitpar(cpl_table *fitpar,
                                       const cpl_table *initfitpar,
                                       const cpl_parameterlist *parlist);
//cpl_error_code sc_mpfit_writeresults(const mp_result *result,
//                                     const cpl_table *scispec,
//                                     const cpl_table *fitpar,
//                                     const cpl_parameterlist *parlist,
//                                     const int calls, const double runtime);
 
#endif /* SC_MPFIT_H */
 
#ifdef __cplusplus
}
#endif