fprintf(stderr, "Internal states (after %d iterations):\n", it);
fprintf(stderr, "A:\n%s\n", debug_eigen(A_).c_str());
fprintf(stderr, "maxA:\n%s\n", debug_eigen(maxA_).c_str());
fprintf(stderr, "C:\n%s\n", debug_eigen(C_).c_str());
fprintf(stderr, "C_shared:\n%s\n", debug_vector(C_shared_).c_str());
fprintf(stderr, "phi:\n%s\n", debug_eigen(phi_).c_str());
fprintf(stderr, "Internal states (after %d iterations):\n", it);
fprintf(stderr, "A:\n%s\n", debug_eigen(A_).c_str());
fprintf(stderr, "maxA:\n%s\n", debug_eigen(maxA_).c_str());
fprintf(stderr, "C:\n%s\n", debug_eigen(C_).c_str());
fprintf(stderr, "C_shared:\n%s\n", debug_vector(C_shared_).c_str());
fprintf(stderr, "phi:\n%s\n", debug_eigen(phi_).c_str());