 
     END_IL_TIMER(_time_knn);
     START_IL_TIMER();
+
+    if (opt.verbosity >= 1) {
+        std::cout << "done.\n Computing Pij matrix..." << std::flush;
+    }
 
     // Search Perplexity
     thrust::device_vector<float> pij_non_symmetric_device(num_points * num_neighbors);


    END_IL_TIMER(_time_init_low_dim);
    START_IL_TIMER();

    // FIT-TNSE Parameters
    int n_interpolation_points = 3;


    START_IL_TIMER();

    if (opt.verbosity >= 1) {
        std::cout << "Initializing cuda handles... " << std::flush;
    }

    // Construct the handles
    cublasHandle_t dense_handle;
    CublasSafeCall(cublasCreate(&dense_handle));
