#include <stdlib.h>
#include <string.h>
#include "structs.h"
#include "debug.h"

void sersic (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void exponential (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void gaussian (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void moffat (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void ferrer (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void king (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void nuker (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df);
void edgedisk (char objtype[], float a[], int ia[], int normtype, 
    struct fitpars *high, struct image *model, struct derivs *df); 

void psfunc (float a[], int ia[], struct image *model, struct derivs *df);
void sky (float a[], int ia[], struct image *model, struct derivs *df);

/*  void sinc (float a[], int ia[], struct fitpars *high,
    struct image *model, struct derivs *df);
void isophote (float a[], int ia[], int normtype, struct fitpars *high, 
    struct image *model, struct derivs *df);    */

void objselect (float a[], int ia[], int normtype, struct fitpars *high, char objtype[], 
    struct image *model, struct derivs *df)
{
    struct fitpars *wptr;

    if (strncmp (objtype,"devauc", 6)==0 || strncmp (objtype,"sersic", 6)==0) 
	sersic (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "expdisk", 7)==0 )
	exponential (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "gaussian", 8)==0 )
	gaussian (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "moffat", 6)==0 )
	moffat (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "ferrer", 6)==0 )
	ferrer (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "king", 4)==0 )
	king (objtype, a, ia, normtype, high, model, df);    

    if (strncmp (objtype, "psf", 3)==0 )
	psfunc (a, ia, model, df);

    if (strncmp (objtype, "sky", 3)==0)
	sky (a, ia, model, df);

    if (strncmp (objtype, "nuker", 5)==0 )
	nuker (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "edgedisk", 8)==0 )
	edgedisk (objtype, a, ia, normtype, high, model, df); 

/*
    if (strncmp (objtype, "isophote", 8)==0 )
	isophote (objtype, a, ia, normtype, high, model, df);

    if (strncmp (objtype, "sinc", 4)==0 )
	fsinc (a, ia, high, model, df); 			*/

}