nr/nr.h

#ifndef _NR_H_
#define _NR_H_

#include "nr/binom.h"
#include "nr/gamma.h"
#include "nr/beta.h"

struct MarquadtMinimizer
{
        float **covar;
        float **alpha; 
        float *atry;
        float *beta;
        float *da;
        float **oneda;

        void (*funcs)(float, float, float [], float *, float [], int);
        
        float  alamda;
        int        mfit;
        float  ochisq; 
        float  chisq;
};

typedef struct MarquadtMinimizer MarquadtMinimizer;

enum {NR_MRQMIN_OK, NR_MRQMIN_SINGULAR_MATRIX};

void covsrt(float **covar, int ma, int ia[], int mfit);

int  gaussj(float **a, int n, float **b, int m);

void mrqcof(float x[], float y[], float sig[], int ndata, 
                        float a[], int ia[], int ma, 
                        float **alpha, float beta[], 
                        float *chisq,
                        void (*funcs)(float, float [], float *, float [], int));

void fgauss(float x, float a[], float *y, float dyda[], int na);

void mrqmin(float x[], float y[], float sig[], int ndata, 
                        float a[], int ia[], int ma, 
                        float **covar, float **alpha, 
                        float *chisq, 
                        void (*funcs)(float, float [], float *, float [], int), 
                        float *alamda);

int mrqmin3D(float x[], float y[], float** z, float** sig, int xdata, int ydata,
                          float a[], int ia[], int ma, 
                          float **covar, float **alpha, 
                          float *chisq,
                          void (*funcs)(float, float, float [], float *, float [], int), 
                          float *alamda);

int mrqmin3Dv2(float x[], float y[], float** z, float** sig, int xdata, int ydata,
                           float a[], int ia[] , int ma, 
                           MarquadtMinimizer *minimizer);

void MarquadtMinimizer_destroy(MarquadtMinimizer* minimizer, int ma);

int mrqcof3D(float x[], float y[], float** z, float** sig, int xdata, int ydata,  
                          float a[], int ia[], int ma, 
                          float **alpha, float beta[], 
                          float *chisq,
                          void (*funcs)(float, float, float[], float *, float [], int));

void four1 (float data[],  unsigned long nn, int isign); /* FFT algorithm */
void realft(float data[],  unsigned long n,  int isign);
void twofft(float data1[], float data2[],        float fft1[], float fft2[], unsigned long n);
void convlv(float data[],  unsigned long n,  float respns[], unsigned long m, int isign, float ans[]);

void savgol(float c[], int np, int nl, int nr, int ld, int m);
void lubksb(float **a, int n, int *indx, float b[]);
void ludcmp(float **a, int n, int *indx, float *d);

#endif

---