fast Interpolation!

bool InitLightTable()
{
  int X;


  // Allocate memory for the table
  LightTable = new float[16384]; // 2^16/sizeof(float)
  LightTable128 = new float[16384];
  if( !LightTable // !LightTable128 )
  {
    // Out of memory
    return false;
  }


  // Intensity = E/r^2*Sin(a) ,   Sin(a) = z/r
  //
  //                sqrt(X)
  // Lighted point *-------+
  //                 \ / a /
  //                   \   / z
  //                  r  \ / 
  //                       * Lightsource (Energy = E/constant)

	
  // Calculate intensity (E = 1, z = 1)
  for( X = 0; X < 16384; X++ )
  {
    LightTable[X] = (float)(1.0f/((float)X/100.0f+1.0f)/sqrt((float)X/100.0f+1.0f));
    LightTable128[X] = (float)(1.0f/((float)X*128.0f/100.0f+1.0f)/sqrt((float)X*128.0f/100.0f+1.0f));
  }


  // To get a desired light intesity at for example 
  // X = 10, E = 2 and z = 5, do either of:
  //
  // I = 2*LightTable[100*10/5^2]/5^2
  // I = 2*LightTable128[(100/128)*10/5^2]/5^2


  // Success
  return true;
}