    
 
     //check residual!
-    for(int i=0;i<k;i++){
-        Axb(user_struct, &V1_d[i*N], &Vcolres_d[i*N]);
-        Arnoldi::check_for_nans("IRA: Schur basis projeciton out in residual", N, &Vcolres_d[i*N]);
-    }   
+//    for(int i=0;i<k;i++){
+//        Axb(user_struct, &V1_d[i*N], &Vcolres_d[i*N]);
+//        Arnoldi::check_for_nans("IRA: Schur basis projeciton out in residual", N, &Vcolres_d[i*N]);
+//    }   
     Arnoldi::matrixMultMatrix_GPU(handle, N, k, k, V1_d, 1.0, Q_Schur_d, 0.0, VRres_d);
     for(int i=0;i<k;i++){
         Arnoldi::vectors_add_GPU(handle, N, -1.0, &Vcolres_d[i*N], &VRres_d[i*N]);

        print_matrix("Q.dat", m, m, Q); 
        print_matrix("EigVecH.dat", k, k, eigenvectorsH_kk_sorted);
        print_vector("EigH.dat", k, eigenvaluesH_kk);
        print_vector("f.dat", N, vec_f_local);  

        delete [] eigenvectorsA, vec_f_local, V_local, V1_local;
        delete [] V_real_local,V_imag_local;

    // Now store EigenvectorsH to GPU as cublasComplex.
    real_complex_to_cublas_complex(k*k, eigenvectorsH_kk_sorted,  eigenvectorsH_d);

    real_device_to_cublas_complex(N*k, V1_d, eigenvectorsA_unsorted_d); //check here V_d or V1_d!!!
    
    Arnoldi::to_device_from_host_int_cpy(sorted_list_d, sorted_list, k, 1, 1);

