<div dir="ltr">Dear Profs. Blaha and Marks and wien users,<br><br>I have a question on how the new Mixer is implemented as it gives me<br>different results on different computers.  Different results mean that<br>the eigenvalues differ by 10^-6 ry, so the wavefunctions are even<br>worse. I am using wien2k_13 version, but there is no difference for<br>newer versions.<br><br><br>In FprojmemR1m.F, call DGESVD gets SVD of Vect, but SVD normally<br>produces an arbitrary sign for U and VT for the same Vect like<br>vect    1     1     1.7542057564<br>vect    1     2     1.7522638412<br>vect    2     1     1.7441856025<br>vect    2     2     1.7425017533<br><br>In one case, my U is<br><br><br>debug3c.outputm:MMT91 U VT    1     1    -0.7091042206   -0.7074733729<br>debug3c.outputm:MMT91 U VT    1     2    -0.7051036834   -0.7067399994<br>debug3c.outputm:MMT91 U VT    2     1    -0.7051036834   -0.7067399994<br>debug3c.outputm:MMT91 U VT    2     2     0.7091042206    0.7074733729<br><br>Another case has<br><br>debug11.outputm:MMT91 U VT    1     1     0.7091042206    0.7074733728<br>debug11.outputm:MMT91 U VT    1     2    -0.7051036834    0.7067399994<br>debug11.outputm:MMT91 U VT    2     1     0.7051036834   -0.7067399994<br>debug11.outputm:MMT91 U VT    2     2     0.7091042206    0.7074733728<br><br>I run 4 iterations to activate the mixer.  When I diff case.scf, these two<br>cases have no difference to the last digit,<br><br>But after the mixing, it differs a lot though the total energy is the<br>same (which is not my interest).<br><br>My questions are (1) whether the new mixer finally uses some products<br>like (U x VT) to generate the new charge density, and (2) whether you<br>have examples that I use to verify this arbitrary sign does not matter.<br><br>Thank you so much in advance!<br><br>Best wishes,<br><br>Guoping<br><br></div>