         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, k, V1_local);
+        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, k, V_local);
+        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);

    // 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);

