 
 
 void compute_eigensystem_of_original_mapping_no_rotation(cublasHandle_t handle, user_map_vector Axb, void *user_struct, int k, int m, int N, char which[2], real *residualAV_real, real *residualAV_imag, real complex *eigenvaluesH_kk, real complex *eigenvectorsH_kk, real complex *eigenvectorsH_kk_sorted, cublasComplex *eigenvectorsA_d, real *V_d, real *V1_d, real *Q_Schur_d){
-
+    //on input eigenvectors are in V1_d
 
     cublasComplex *eigenvectorsH_d,  *eigenvectorsA_unsorted_d, *eigenvectorsA_sorted_d;
     real *Vcolres_d, *VRres_d;

        Arnoldi::to_host_from_device_real_cpy(V_imag_local, eigenvectors_imag_d, N, k, 1);

        print_matrix("EigVecA.dat", N, k, eigenvectorsA);
        print_matrix("V1.dat", N_rotational, k, V1_local);
        print_matrix("V_real.dat", N, k, V_real_local);//eigenvectors_real_d
        print_matrix("V_imag.dat", N, k, V_imag_local);//eigenvectors_imag_d
        print_matrix("V.dat", N_rotational, k, V_local);
        print_matrix("H.dat", m, m, H);
        print_matrix("H1.dat", m, m, H1);
        print_matrix("H2.dat", m, m, H2);
