<div dir="ltr"><div><div><div><div>Dear All,<br><br></div>I'm calculating eigenvalues and eigenvectors of the matrix which has specific kind of symmetry.<br></div>Due to this symmetry I obtain the eigenvalues which are doubly degenerated. <br></div>So eg. eigeinvalue 'e1' has eigenvectors 'a1' and 'b1'. These eigenvectors are related to each other by the relation a1 = T b1, where T is a matrix (given for my problem).<br>So it is enough to calculate only one eigenvector for each eigenvalue (and the second one can be calculated by matvec operation). This situation has been described in <a href="http://dl.acm.org/citation.cfm?id=2494747">http://dl.acm.org/citation.cfm?id=2494747</a>.<br><br></div><div>How could I take advantage on this in EPSSolve in Jacobi-Davidson method? <br>Could I add two vectors to the subspace (the second one would be calculated by multiplying the first one by matrix T) in every iteration? Should I modify function "dvd_updateV_update_gen" in dvd_updatev.c ? <br><br></div><div>I would be very grateful for any suggestion.<br><br></div><div>Krzysztof<br></div><br><br><div><br>  <br></div></div>