__global__ void kernel_Median3RegTri9( short *output, int i_dim, int j_dim)
{  
  
  // coordonnees absolues du point
  int j = __mul24(blockIdx.x,blockDim.x) + threadIdx.x ; 
  int i = __mul24(blockIdx.y,blockDim.y) + threadIdx.y ;
  
  /**************************************************************************
   *             tri(s)
   **************************************************************************/
  int a0, a1, a2, a3, a4, a5, a6, a7, a8 ;

  /********************************************************************************
   * affectation des valeurs
   ********************************************************************************/
  a0 = tex2D(tex_img_ins, j-1, i-1) ;
  a1 = tex2D(tex_img_ins, j  , i-1) ;
  a2 = tex2D(tex_img_ins, j+1, i-1) ;
  a3 = tex2D(tex_img_ins, j-1, i) ;
  a4 = tex2D(tex_img_ins, j  , i) ;
  a5 = tex2D(tex_img_ins, j+1, i) ;
  a6 = tex2D(tex_img_ins, j-1, i+1) ;
  a7 = tex2D(tex_img_ins, j  , i+1) ;
  a8 = tex2D(tex_img_ins, j+1, i+1) ;

  //tri selection
  bubTriReg9(&a0, &a1, &a2, &a3, &a4, &a5, &a6, &a7, &a8);
  
  //median au milieu !
  output[ __mul24(i, j_dim) +j ] = a4 ;
 
}