lamino_ramp.cl 996 B

12345678910111213141516171819202122232425262728293031323334353637383940414243
  1. __kernel void lamino_ramp_create_filter(
  2. __global float * flt,
  3. const unsigned int width,
  4. const unsigned int fwidth,
  5. const unsigned int height,
  6. const float theta,
  7. const float tau)
  8. {
  9. const int idx = get_global_id(0);
  10. const int idy = get_global_id(1);
  11. const int index = idy * width + idx;
  12. const float scale = (float)width/(float)height;
  13. flt[index] = 0;
  14. if( idy > 0 ) return;
  15. const float sin2 = -sin(theta) / 2.;
  16. const float const_pi_tau= M_PI_F * M_PI_F * tau;
  17. const uint quatw = fwidth / 2;
  18. if(idx == 0)
  19. {
  20. flt[idx] =scale* -sin2 / ( 4 * tau);
  21. return;
  22. }
  23. if( (idx <= quatw) && ( (idx%2) != 0))
  24. {
  25. flt[idx] = scale*sin2/((float)idx*(float)idx*const_pi_tau);
  26. return;
  27. }
  28. int x = width - idx;
  29. if( (idx >= (width - quatw)) && ( (x%2) != 0))
  30. {
  31. flt[idx] = scale*sin2/((float)x*(float)x*const_pi_tau);
  32. return;
  33. }
  34. }