// --Kernel-Funktion fuer Sobel-Operator

#define AT(x,y)		((y)*length+(x))

__global__ void sobel_operator(unsigned char* in, unsigned char* out, int length, int width) {

// --Wer bin ich? Bestimmung des Bildpunktes
// --Es wird jeweils um 1 inkrementiert, da auf dem Rand kein Sobel-Wert berechnet wird.
int x=blockDim.x*blockIdx.x+threadIdx.x+1;
int y=blockDim.y*blockIdx.y+threadIdx.y+1;

// --Anwendung des Operators (Achte auf Rand!)
if (x<length-2 && y<width-2) {

    // --Sobel-Operator in horizontaler Richtung
    int sh=1*in[AT(x-1,y+1)]-1*in[AT(x-1,y-1)]+
           2*in[AT(x,y+1)]-2*in[AT(x,y-1)]+
           1*in[AT(x+1,y+1)]-1*in[AT(x+1,y-1)];

    // --Sobel-Operator in vertikaler Richtung
    int sv=1*in[AT(x+1,y+1)]-1*in[AT(x-1,y+1)]+
           2*in[AT(x+1,y)]-2*in[AT(x-1,y)]+
           1*in[AT(x+1,y-1)]-1*in[AT(x-1,y-1)];

    //	--Zusaetzlich invertieren zur besseren Darstellung
    out[AT(x,y)]=255-(unsigned char)sqrt(float(sh)*float(sh)+float(sv)*float(sv));
    }
}