<font face="Default Sans Serif,Verdana,Arial,Helvetica,sans-serif" size="2">For WIEN2k_12, you have to look at lines 885-894 in hamilt.F.<br><br><font color="#990099">-----wien-bounces@zeus.theochem.tuwien.ac.at wrote: -----</font><div style="padding-left:5px;"><div style="padding-right:0px;padding-left:5px;border-left:solid black 2px;">To: A Mailing list for WIEN2k users &lt;wien@zeus.theochem.tuwien.ac.at&gt;<br>From: Peter Blaha <pblaha@theochem.tuwien.ac.at><br>Sent by: wien-bounces@zeus.theochem.tuwien.ac.at<br>Date: 05/03/2013 09:53AM<br>Subject: Re: [Wien] a possible bug in hamilt.F (local-local part)<br><br><div><font face="Courier New,Courier,monospace" size="2">Hi,<br><br>I'm not sure I understand your message correctly.<br><br>You are saying in WIEN2k_06 there was this error and the first local <br>orbital was not taken into account properly ???<br><br>What about WIEN2k_12 ?? Is the error still present or not ??<br><br>At least in _12 I do not have a do-loop 278 &nbsp;and lines 901-4 do not have <br>a do-loop either .... ????<br><br>Regards<br><br>On 10/22/2003 03:14 PM, Guo-ping Zhang wrote:<br>&gt; Dear Peter and Wien users,<br>&gt;<br>&gt; I noticed a potential bug &nbsp;in hamilt.F. The matrix misses one row.<br>&gt;<br>&gt; In _06 &nbsp;version, (same with _12 version (Line 901-904)<br>&gt;<br>&gt; &nbsp; &nbsp;DO 278 J = NV + NNLO - (MO1+L) - &amp;<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (2*L+1)*(JEQO-1) - (jlo-jlop)*(2*l+1)*mult(jneq), NV + NNLO &amp;<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -(jlo-jlop)*(MO1+L+1)-(jlo-jlop)*(JEQO-1)*(2*l+1)<br>&gt;<br>&gt; the low limit jumps incorrectly.<br>&gt;<br>&gt; Here is an example of case.output1 of _06 version.<br>&gt;<br>&gt; &nbsp; &nbsp; 642 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;-.1874358827D-03 &nbsp;0.0000000000D+00<br>&gt; J,I,IHELP,JNEQ &nbsp;sphere-local Line 891<br>&gt; &nbsp; &nbsp; 642 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;-.1874358827D-03 &nbsp;-.5997948246D-02<br>&gt; J,I,IHELP,JNEQ &nbsp;sphere-local Line 955<br>&gt; &nbsp; &nbsp; 644 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;0.0000000000D+00 &nbsp;0.0000000000D+00<br>&gt; J,I,IHELP,JNEQ, local-local 1<br>&gt; &nbsp; &nbsp; 644 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;0.2341672678D+00 &nbsp;0.0000000000D+00<br>&gt; J,I,IHELP,JNEQ, local-local 2<br>&gt; &nbsp; &nbsp; 644 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;0.2341672678D+00 &nbsp;0.0000000000D+00<br>&gt; J,I,IHELP,JNEQ, local-local 3<br>&gt;<br>&gt; See J index just runs from 642 to 644, without 643. In this case, from 1<br>&gt; to 642 is planewave, from 643 on is local orbital.<br>&gt;<br>&gt; This is the _12 results.<br>&gt;<br>&gt; Here is correct.<br>&gt;<br>&gt; &nbsp; &nbsp; 641 &nbsp; 643 &nbsp; &nbsp; 3 &nbsp; &nbsp; 1 &nbsp;-.1200259258D-01 &nbsp;0.0000000000D+00<br>&gt; J_g,I_g,IHELP,JNEQ, sphere-local<br>&gt; &nbsp; &nbsp; 642 &nbsp; 643 &nbsp; &nbsp; 3 &nbsp; &nbsp; 1 &nbsp;-.1200259258D-01 &nbsp;0.0000000000D+00<br>&gt; J_g,I_g,IHELP,JNEQ, sphere-local<br>&gt; &nbsp; &nbsp; 643 &nbsp; 643 &nbsp; &nbsp; 3 &nbsp; &nbsp; 1 &nbsp;0.2560228649D+00 &nbsp;0.0000000000D+00<br>&gt; J_g,I_g,IHELP,JNEQ, local-local 2<br>&gt;<br>&gt; Then, wrong ( 643 is missing ).<br>&gt;<br>&gt; &nbsp; &nbsp; 641 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;-.1867862609D-03 &nbsp;-.5977160350D-02<br>&gt; J_g,I_g,IHELP,JNEQ, sphere-local<br>&gt; &nbsp; &nbsp; 642 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;-.1867862609D-03 &nbsp;-.5977160350D-02<br>&gt; J_g,I_g,IHELP,JNEQ, sphere-local<br>&gt; &nbsp; &nbsp; 644 &nbsp; 644 &nbsp; &nbsp; 4 &nbsp; &nbsp; 1 &nbsp;0.2341672678D+00 &nbsp;-.1179069863D-17<br>&gt; J_g,I_g,IHELP,JNEQ, local-local 2<br>&gt;<br>&gt;<br>&gt; Would you please look into this?<br>&gt;<br>&gt; I can send you my debug code and case.struct file, if necessary.<br>&gt;<br>&gt; Thanks a lot!<br>&gt;<br>&gt; Best regards,<br>&gt; Guoping<br>&gt; _______________________________________________<br>&gt; Wien mailing list<br>&gt; Wien@zeus.theochem.tuwien.ac.at<br>&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt; SEARCH the MAILING-LIST at:<br>&gt; <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br><br>-- <br><br>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; P.Blaha<br>--------------------------------------------------------------------------<br>Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>Phone: +43-1-58801-165300 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; FAX: +43-1-58801-165982<br>Email: blaha@theochem.tuwien.ac.at &nbsp; &nbsp;WWW: <br><a href="http://info.tuwien.ac.at/theochem/">http://info.tuwien.ac.at/theochem/</a><br>--------------------------------------------------------------------------<br>_______________________________________________<br>Wien mailing list<br>Wien@zeus.theochem.tuwien.ac.at<br><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>SEARCH the MAILING-LIST at: &nbsp;<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br></font></div></pblaha@theochem.tuwien.ac.at></div></div><div></div></font>