     ellipse.projectiveTransform( mInvT, transformedEllipse );
 
     bool first_iteration = true;
+
+cerr << "TagIndex: " << tagIndex << " number of loops: " << iterations << endl;
+
     for( ; iterations>0; iterations-- ) {
         float neighSize = currentNeighbourSize * max( transformedEllipse.a(),
                                                       transformedEllipse.b() );

 */
__global__
void idComputeResult( const int         tagIndex,
                      NearbyPoint*      nPoint,
                      const CutStruct*  vCuts,
                      const CutSignals* allcut_signals,
                      const int         vCutsSize )
{
    int myPair = blockIdx.x * 32 + threadIdx.y;
    int j      = __float2int_rd( 1.0f + __fsqrt_rd(1.0f+8.0f*myPair) ) / 2;
    int i      = myPair - j*(j-1)/2;


#include "nearby_point.h"
#include "tag_cut.h"

#undef NO_ATOMIC

using namespace std;

namespace popart {
