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

void high_derivs (int ix, int iy, struct fitpars *h, struct derivs *hdf, 
    float dfdr, float *drd, float dfdI, float *dIdfp)
{
    int n;

    for (n=1; n <= h->npars; n++) {
        if (h->ia[n] > 0) 
	    if (dIdfp != NULL)
                hdf->dpm[n][iy][ix] += (dfdr * drd[n] + dfdI * dIdfp[n]);
	    else
                hdf->dpm[n][iy][ix] += (dfdr * drd[n]);
    };
}