5         result += 
sqrt( 
pow(centre.x - v_array[
i].x,2) + 
pow(centre.y - v_array[
i].y,2));
 __device__ double2 pow(double2 a, int b)
power operation for a double2 
 
double sepAvg(int2 *vArray, int2 centre, int length)
 
end if sqrt(sum([x(ii), y(ii)].^ 2))< radius %% ignore edges if(length(DT.vertexAttachments