#include <stdio.h>
#include "structs.h"

void clearmodel (struct image *model, struct derivs *df, struct fitpars *fpar)
{
    void clear_fpar_images (struct derivs *df, struct fitpars *fpar);

    int i, j;
    struct fitpars *hptr;
    struct derivs *hdf;

    for (j=1; j <= model->naxes[2]; j++) {
        for (i=1; i <= model->naxes[1]; i++)
            model->z[j][i] = 0.;
    };

    while (fpar != NULL) {
	clear_fpar_images (df, fpar);

        df = df->next;
        fpar = fpar->next;
    };
}

void clear_fpar_images (struct derivs *df, struct fitpars *fpar)
{
    long i, j;
    int k;

    for (k=0; k<=fpar->npars; k++) {
        if (k==0 || fpar->ia[k] > 0) {
            for (j=1; j <= df->naxes[2]; j++) {
                for (i=1; i <= df->naxes[1]; i++)
		    if (df->dpm[k] != NULL)
                        df->dpm[k][j][i] = 0.;
            };
        };
    };

    if (fpar->high != NULL) 
        clear_fpar_images (df->high, fpar->high);
    
    if (fpar->next != NULL)
        clear_fpar_images (df->next, fpar->next);
}