mk_hex_sites

 creates the pattern of the hexagonal cells

  INPUT
   par         general parameter

  OUTPUT
   sites       center of the cells in complex format



 $Id: mk_hex_sites.m,v 1.1 2004/05/12 16:22:59 dalai Exp $

0001 % creates the pattern of the hexagonal cells
0002 %
0003 %  INPUT
0004 %   par         general parameter
0005 %
0006 %  OUTPUT
0007 %   sites       center of the cells in complex format
0008 %
0009 %
0010 %
0011 % $Id: mk_hex_sites.m,v 1.1 2004/05/12 16:22:59 dalai Exp $
0012 
0013 function sites = mk_hex_sites(km,lm,r)
0014 
0015 % the base
0016 base = exp(i*2/3*pi*(0:2));
0017 
0018 % rotation + redimension
0019 a = base*exp(i*pi/6)/sqrt(3);
0020 
0021 % a simplified version of RUNE stuff
0022 % km = par.km;
0023 % lm = par.lm;
0024 
0025 kv = 0:km-1;
0026 lv = 0:lm-1;
0027 
0028 [m1 m2] = meshgrid(kv*base(1),kv*base(2));
0029 kk_side = m1(:) + m2(:) + a(1);
0030 
0031 [m1 m2] = meshgrid(lv*base(1),lv*base(3));
0032 ll_side = m1(:) + m2(:) + a(3);
0033 
0034 [m1 m2] = meshgrid(kv*base(2),lv*base(3));
0035 lk_side = m1(:) + m2(:) + a(2);
0036 
0037 sites = [kk_side(:).' ll_side(:).' lk_side(:).'];
0038 sites = (sites - mean(sites)).*r;