#include <stdio.h>
#include <fcntl.h>
#include <errno.h>
#include <string.h>
#include <strings.h>
#include <math.h>
#include "hisview_view.h"
#include "hisview_macros.h"
#include "hisview_widgets.h"

#ifdef underscore
#define  hsv_to_rgb hsv_to_rgb_ 
#endif
#ifdef f2c
#define  hsv_to_rgb hsv_to_rgb__
#endif

#define MAXCOL 200
#define FIRSTCOL 50
#define MAXNHUE 10
#define MAXNSAT 4
#define MAXNVAL 5
#define INCRSAT 0.25
#define INCRVAL 0.2


static char dumstr[80];

static float hue[MAXCOL];
static float sat[MAXCOL];
static float val[MAXCOL];

float red[MAXCOL];
float green[MAXCOL];
float blue[MAXCOL];
float hue_to_red[MAXCOL];
float sat_to_green[MAXCOL];
float val_to_blue[MAXCOL];
float fullred[MAXCOL];
float fullgreen[MAXCOL];
float fullblue[MAXCOL];
float gray[MAXCOL];

unsigned long ncolors; 
int nhue, nsat, nval;
float valhue[MAXNHUE];
float valsat[MAXNSAT];
float valval[MAXNVAL];
unsigned char colmap[MAXCOL];

static unsigned long colors[MAXCOL];
static unsigned long firstcolors[FIRSTCOL];

static char *namec[MAXCOL];
static XColor exact_defs[MAXCOL];
static XColor first_defs[FIRSTCOL];

static char colnam[8];

void transfert_function()
{
  int k;

  for (k = 0; k < ncolors; k++)
    sprintf(namec[k], "#%02X%02X%02X",
	    (unsigned short) (255*red[k]),
	    (unsigned short) (255*green[k]),
	    (unsigned short) (255*blue[k]));



  for (k=0; k < ncolors; k++)
    if (!XParseColor(display,
		     cmap,
		     namec[k],
		     &exact_defs[k])){
      sprintf(dumstr,
	      "I can't parse col # %d, value %s.\n",
	      k, namec[k]);
      printf(dumstr);
    }

  XStoreColors(display, cmap, exact_defs, ncolors);

}


void init_luts()
{
  int i, j, k, n;
  float t, h;

  nhue = MAXNHUE;
  nsat = MAXNSAT;
  nval = MAXNVAL;

  for (i=0; i<nhue; i++) {
    valhue[i] = (270./nhue)*i;
    for (j=0; j < nsat*nval; j++)
      hue[nsat*nval*i + j] = valhue[i];
  }

  for (i=0; i<nval; i++)
    valval[i] = INCRVAL + INCRVAL*i;

  for (k=0; k<nhue; k++) {
    for (i=0; i<nval; i++) {
      for (j = 0; j < nsat; j++)
	val[nsat*nval*k + nsat*i + j] = valval[i];
    }
  }
  for (i=0; i<nsat; i++)
    valsat[i] = INCRSAT + INCRSAT*i;

  for (n=0; n<nhue; n++) {
    for (k = 0; k < nval; k++) {
      for (i=0; i<nsat; i++) {
	sat[nsat*nval*n + nsat*k + i] = valsat[i];
      }
    }
  }

  for (i=0; i<ncolors; i++) {
    hsv_to_rgb(&hue[i], &sat[i], &val[i], 
	       &hue_to_red[i], &sat_to_green[i], &val_to_blue[i]);
  
    red[i] = hue_to_red[i];
    green[i] = sat_to_green[i];
    blue[i] = val_to_blue[i];
 
  }

  t = 1.;
  for (i=0; i<ncolors; i++) {
    gray[i] = (float) i/(ncolors-1);
    h = gray[i]*270.;
    hsv_to_rgb(&h, &t, &t, 
	       &fullred[i], &fullgreen[i], &fullblue[i]);      


  }

  transfert_function();
}

void init_colors()
{
  int i;
  float t;

  ncolors = MAXCOL;



  if(!XAllocColorCells(display, cmap, NULL,NULL,NULL, firstcolors, FIRSTCOL)) {
    printf("Can't allocate colors\n"); fflush(stdout);
    exit(1);
  }

  for (i = 0; i < FIRSTCOL; i++)
    first_defs[i].pixel = firstcolors[i];

  XQueryColors(display, default_cmap, first_defs, FIRSTCOL);

  XStoreColors(display, cmap, first_defs, FIRSTCOL);



  if (!XAllocColorCells(display, cmap, NULL, NULL, NULL, colors, ncolors)) {
    printf("Can't allocate colors\n"); fflush(stdout);
    exit(1);
  }
  /*
  XInstallColormap(display, cmap);
  */

  for (i = 0; i < ncolors; i++) {
    exact_defs[i].pixel = colors[i];
    colmap[i] = (unsigned char) exact_defs[i].pixel;
    namec[i] = (char *) malloc(8);
  }
}