/* Last edited on 2008-08-21 11:16:30 by stolfi */
/***************************************************************************
 *   Copyright (C) 2008 by Douglas Castro   *
 *   douglas@localhost.localdomain   *
 *                                                                         *
 ***************************************************************************/
#include <bases.h>

int compare_doubles (const void *a, const void *b)
{
  double *ad = (double*)a;
  double *bd = (double*)b;
  return (int) (*ad - *bd);
}

void cria_base(VReg pac[], double xmin[], double xmax[], int profund, Preditor pred, double eps)
{
	int  k=5; // posteriormente k pode se tornar argumento da funcao
	//int ctr = (k*k)/2;
	int j;
	
// 	double vet[] = {43.3,2.0, 4.6, 3.2,5.2};
// 	int p = sizeof(vet)/sizeof(double);// da numero de casa do vetor vet

	if(profund>=28)
	{
		return;
	}
	// copiando o pacote
	double vet2[k*k];
	for(j=0;j<k*k;j++)
	{
		vet2[j]=pac[j].fv;
	}
	qsort (vet2, k*k, sizeof (double), compare_doubles);
	if (fabs(vet2[24]-vet2[0])<eps) 
	{ 
		return;
	}

	/* Precisa dividir pacote */

	/* Decide o eixo {ex} perpendicular ao corte: */
	int ex = profund % 2;
  
	int qual;

	for (qual = 0; qual <= 1; qual++)
	{ 
		interp_pacote(pac, xmin, xmax, ex, pred, qual);
		cria_base(pac, xmin, xmax, profund+1, pred, eps);
	}

}

void interp_pacote(VReg pac[],double xmin[], double xmax[], int ex,	Preditor pred,	int qual)
{
	int k=5;

	int i,j;
	if(ex == 0)
	{ // montando gabarito
		for(i=0; i<k; i++)
		{ 
			for(j=0; j<k; j++)
			{
				double fa, fb, fc;
				if(j==0)
				{
					fa = 0.0;
				}
				else
				{
					fa = pac[i*k+j-1].fv;
				}
				fb = pac[i*k+j].fv;
				if(j==4)
				{
					fc = 0.0;
				}
				else
				{
					fc = pac[i*k+j+1].fv;
				}

				pac[i*k+j].fv = pred(fa, fb, fc);
			}
		}
	}
	else
	{
    // montando gabarito
		int i,j;
		for(i=0;i<k;i++)
		{
			for(j=0;j<k;j++)
			{
				double fa, fb, fc;
				if(i==4)
				{
					fc = 0.0;
				}
				else
				{
					fc = pac[(i+1)*k+j].fv;
				}
				fb = pac[i*k+j].fv;
				if(i==0)
				{
					fa = 0.0;
				}
				else
				{
					fa = pac[(i-1)*k+j].fv;
				}
				pac[i*k+j].fv = pred(fa, fb, fc);
			}
		}
	}
}