#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include "mucal.h"
#include "panel.h"

/* ********************************************** */

int facet_order (
   facet  FACET[],
   int    L_order[],
   int    distrib[],
   int    LEN,
   float  ring_width,
   int    layer) 
{
   int i=0, fx_min, fx, indx, num=0, num2=0, num3=0, iadd;

   float r, r_min;

   iadd = layer * N_CryPos;

   num = 0;
   for (fx=0; fx<LEN; fx++) {
      if (FACET[fx + iadd].nOption < 0) continue;
      num++;
   }
   printf("In facet_order: input facets: %4d, valid facets (nOption >=0): %4d\n", LEN, num);
 
again:
   r_min = 30000.0;
   for (fx=0; fx<LEN; fx++) {
      if (FACET[fx + iadd].nOption < 0) continue;
      r = FACET[fx + iadd].r;
      if (r_min > r) {r_min = r; fx_min = fx;}
   } 

   if (r_min < 30000.0) {
      L_order[i++] = fx_min;
      FACET[fx_min + iadd].nOption = -2;
      indx = (int)(FACET[fx_min + iadd].r);
      distrib[indx]++;
      num2++;
   } else {
      for (fx=0; fx<LEN; fx++) if (FACET[fx + iadd].nOption == -2) {num3++; FACET[fx + iadd].nOption = 0;}
      distrib[0] = 0;
      printf("Facets sorted %5d %5d %5d %5d, first: %5d, (orig #: %5d), r_min: %5.1f, fx_min: %4d, i: %4d\n", 
         LEN, num, num2, num3, 0, L_order[0], FACET[L_order[0] + iadd].r, fx_min, i);
      printf("Facets sorted %5d %5d %5d %5d,  last: %5d, (orig #: %5d), r_max: %5.1f\n", 
         LEN, num, num2, num3, fx_min, L_order[fx_min], FACET[fx_min + iadd].r);

      return (i);
   }
   if (i < N_CryPos -1) goto again; else {printf("In 'facet_order': index exceeds limit (%7d)\n", N_CryPos ); exit(0);}
}