#include "global.h"
void propagate()
{
/* definition of e_i/** 6 2 5 * 3 0 1* 7 4 8*/
/* the only trick here is to have loop directions antiparallel to the lattice vectors */	register int x,y;	
  for(y=max_y-1;y>=0;y--)
   {++)	
       {			
        f[6][y][x]=f[6][y-1][x+1];	
       }	
     }		
   for(y=0;y<max_y-1;y++)	
     {
      for(x=max_x-1;x>0;x--)	
       {		
        f[8][y][x]=f[8][y+1][x-1];
       }		
      for(x=0;x<max_x-1;x++)		
       {			
        f[7][y][x]=f[7][y+1][x+1];