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

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

/* Sequential numbering of individual facets in the lens. This version for three facet layers
 * The first three parameters in the call are 'ring_no' (the current hexagon ring number) layer number and facet sequential number.
 * The routine returns the initial facet number for the subsequent hexagonal facet ring.
 *
   Niels Lund April 2016 - December 2018

   panel_new[0] = Foc_len;
   panel_new[1] = lens_diam;  
   panel_new[2] = facet_diam;  
   panel_new[3] = x_focusU;
   panel_new[4] = y_focusU;
   panel_new[5] = x_focusD;
   panel_new[6] = y_focusD;
   panel_new[7] = cameraH;
   panel_new[8] = lensDiameter;
*/

int hexring(
   int   ring_no, 
   int   layer_no, 
   int   facet_no, 
   float panel_new[],       /* structure specifying the lens and facet geometry */
   facet *FACET)
{
   int i, j, k=0, iadd;
   static int last_facet_no=0;

   iadd = layer_no * N_CryPos;

   for (i=0; i<ring_no; i++) {		 /* 1 */
      FACET[facet_no + iadd].u = ring_no - i;
      FACET[facet_no + iadd].v = i;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    
   for (i=0; i<ring_no; i++) {		 /* 2 */
      FACET[facet_no + iadd].u = -i;
      FACET[facet_no + iadd].v = ring_no;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    
   for (i=0; i<ring_no; i++) {		 /* 3 */
      FACET[facet_no + iadd].u = -ring_no;
      FACET[facet_no + iadd].v = ring_no - i;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    
   for (i=0; i<ring_no; i++) {		 /* 4 */
      FACET[facet_no + iadd].u = -ring_no + i;
      FACET[facet_no + iadd].v = -i;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    
   for (i=0; i<ring_no; i++) {		 /* 5 */
      FACET[facet_no + iadd].u = i;
      FACET[facet_no + iadd].v = -ring_no;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    
   for (i=0; i<ring_no; i++) {		 /* 6 */
      FACET[facet_no + iadd].u = ring_no;
      FACET[facet_no + iadd].v = -ring_no + i;
      FACET[facet_no + iadd].ring_num = ring_no;
      facet_no = facet_geo(&(FACET[facet_no + iadd]), panel_new, layer_no, facet_no);
   }    

   last_facet_no = facet_no;
   return facet_no;
}

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