 				float3 origin_normal = input_normals ? fetch(input_normals, index) : make_float3(0.f, 0.f, 0.f); 				
 				origin_normal = normalized_safe(origin_normal); //normalize if non-zero
 
-				float3 rotation_axis = getRotationAxes(index, origin_normal);
+                float3 rotation_axis = AxesStrategy::getRotationAxes(index, origin_normal);
 				rotation_axis = normalized_safe(rotation_axis); //normalize if non-zero
 
 				const float eps = numeric_limits<float>::epsilon ();

	thrust::transform(beg, end, out, GtThan(min_neighb));
	cudaSafeCall( cudaGetLastError() );
	cudaSafeCall( cudaDeviceSynchronize() );	
}
		{	
			enum
			{
                CTA_SIZE = 192
			};
			
			int work_size;
