Ignore:
Timestamp:
09/14/10 10:44:50 (9 years ago)
Author:
rlentz
Message:

latest

File:
1 edited

Legend:

Unmodified
Added
Removed
  • OpenCLIterative3DDeconvolution/src/publication/sobel.cl

    r6961 r6962  
    66    int x = get_global_id(0); 
    77    int y = get_global_id(1); 
    8     float p[9]; 
     8    float p0, p1, p2, p3, p4, p5, p6, p7, p8 = 0; 
    99    int offset = y * width + x; 
    1010if( x < 1 || y < 1 || x > width - 2 || y > height - 2 ) 
     
    1414else 
    1515{ 
    16     p[0] = input[offset - width - 1] ; 
    17     p[1] = input[offset - width] ; 
    18     p[2] = input[offset - width + 1] ; 
    19     p[3] = input[offset - 1] ; 
    20     p[4] = input[offset] ; 
    21     p[5] = input[offset + 1] ; 
    22     p[6] = input[offset + width - 1] ; 
    23     p[7] = input[offset + width] ; 
    24     p[8] = input[offset + width + 1] ; 
     16    p0 = input[offset - width - 1] ; 
     17    p1 = input[offset - width] ; 
     18    p2 = input[offset - width + 1] ; 
     19    p3 = input[offset - 1] ; 
     20    p4 = input[offset] ; 
     21    p5 = input[offset + 1] ; 
     22    p6 = input[offset + width - 1] ; 
     23    p7 = input[offset + width] ; 
     24    p8 = input[offset + width + 1] ; 
    2525 
    26     float sum1 = p[0] + 2*p[1] + p[2] - p[6] - 2*p[7] - p[8];  //GY 
    27     float sum2 = p[0] + 2*p[3] + p[6] - p[2] - 2*p[5] - p[8];  //GX 
     26    float sum1 = p0 + 2*p1 + p2 - p6 - 2*p7 - p8;  //GY 
     27    float sum2 = p0 + 2*p3 + p6 - p2 - 2*p5 - p8;  //GX 
    2828    
    2929    output[offset] = sqrt(  sum1*sum1 + sum2*sum2 ); 
Note: See TracChangeset for help on using the changeset viewer.