Working but crashing on dealloc of command queue

qgis · Aug 8, 2018 · 6f9a72e · 6f9a72e
1 parent 51d1bb3
commit 6f9a72e
Showing 1 changed file with 111 additions and 0 deletions.
diff --git a/src/analysis/raster/slope.cl b/src/analysis/raster/slope.cl
@@ -0,0 +1,111 @@
+#pragma OPENCL EXTENSION cl_khr_fp64 : enable
+
+float calcFirstDer( float x11, float x21, float x31, float x12, float x22, float x32, float x13, float x23, float x33,
+                       float mInputNodataValue, float mOutputNodataValue, double mZFactor, double mCellSize )
+{
+ //the basic formula would be simple, but we need to test for nodata values...
+ //X: return (( (x31 - x11) + 2 * (x32 - x12) + (x33 - x13) ) / (8 * mCellSizeX));
+ //Y: return (((x11 - x13) + 2 * (x21 - x23) + (x31 - x33)) / ( 8 * mCellSizeY));
+
+ int weight = 0;
+ double sum = 0;
+
+
+ //first row
+ if ( x31 != mInputNodataValue && x11 != mInputNodataValue ) //the normal case
+ {
+   sum += ( x31 - x11 );
+   weight += 2;
+ }
+ else if ( x31 == mInputNodataValue && x11 != mInputNodataValue && x21 != mInputNodataValue ) //probably 3x3 window is at the border
+ {
+   sum += ( x21 - x11 );
+   weight += 1;
+ }
+ else if ( x11 == mInputNodataValue && x31 != mInputNodataValue && x21 != mInputNodataValue ) //probably 3x3 window is at the border
+ {
+   sum += ( x31 - x21 );
+   weight += 1;
+ }
+
+ //second row
+ if ( x32 != mInputNodataValue && x12 != mInputNodataValue ) //the normal case
+ {
+   sum += 2.0f * ( x32 - x12 );
+   weight += 4;
+ }
+ else if ( x32 == mInputNodataValue && x12 != mInputNodataValue && x22 != mInputNodataValue )
+ {
+   sum += 2.0f * ( x22 - x12 );
+   weight += 2;
+ }
+ else if ( x12 == mInputNodataValue && x32 != mInputNodataValue && x22 != mInputNodataValue )
+ {
+   sum += 2.0f * ( x32 - x22 );
+   weight += 2;
+ }
+
+ //third row
+ if ( x33 != mInputNodataValue && x13 != mInputNodataValue ) //the normal case
+ {
+   sum += ( x33 - x13 );
+   weight += 2;
+ }
+ else if ( x33 == mInputNodataValue && x13 != mInputNodataValue && x23 != mInputNodataValue )
+ {
+   sum += ( x23 - x13 );
+   weight += 1;
+ }
+ else if ( x13 == mInputNodataValue && x33 != mInputNodataValue && x23 != mInputNodataValue )
+ {
+   sum += ( x33 - x23 );
+   weight += 1;
+ }
+
+ if ( weight == 0 )
+ {
+   return mOutputNodataValue;
+ }
+
+ return sum / ( weight * mCellSize ) * mZFactor;
+}
+
+
+__kernel void processNineCellWindow( __global float *scanLine1,
+                                   __global float *scanLine2,
+                                   __global float *scanLine3,
+                                   __global float *resultLine,
+                                   __global float *mInputNodataValue,
+                                   __global float *mOutputNodataValue,
+                                   __global double *mZFactor,
+                                   __global double *mCellSizeX,
+                                   __global double *mCellSizeY
+                       ) {
+
+  // Get the index of the current element
+  int i = get_global_id(0);
+
+  // Do the operation
+  //return (( (x31 - x11) + 2 * (x32 - x12) + (x33 - x13) ) / (8 * mCellSizeX))
+  float derX = calcFirstDer(   scanLine1[i],   scanLine2[i],   scanLine3[i],
+                               scanLine1[i+1], scanLine2[i+1], scanLine3[i+1],
+                               scanLine1[i+2], scanLine2[i+2], scanLine3[i+2],
+                               *mInputNodataValue, *mOutputNodataValue, *mZFactor, *mCellSizeX
+                             );
+  //return (((x11 - x13) + 2 * (x21 - x23) + (x31 - x33)) / ( 8 * mCellSizeY));
+  float derY = calcFirstDer(   scanLine1[i+2], scanLine1[i+1], scanLine1[i],
+                               scanLine2[i+2], scanLine2[i+1], scanLine2[i],
+                               scanLine3[i+2], scanLine3[i+1], scanLine3[i],
+                               *mInputNodataValue, *mOutputNodataValue, *mZFactor, *mCellSizeY
+                             );
+  if ( derX == *mOutputNodataValue || derY == *mOutputNodataValue )
+  {
+    resultLine[i] = *mOutputNodataValue;
+  }
+  else
+  {
+    double res = sqrt( derX * derX + derY * derY );
+    res = atanpi( res );
+    resultLine[i] = (float)res * 180.0;
+  }
+}