<div>Well, but as far asI have found, O(i) stands for &#39;simple momentum matrix elements&#39;.Then,</div>
<div>I take it for grant that (O(i))* O(j) the diagnal (if i = j) and off-diagnal (if i != j) components.</div>
<div>&nbsp;</div>
<div>So I Wonder eagerly that why it is necessary including the body below</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(1)=o2(1)+o1(1)*conjg(o1(1))&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(2)=o2(2)+o1(2)*conjg(o1(2))&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(3)=o2(3)+o1(3)*conjg(o1(3))&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (det(l).LT.0) then&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(4)=o2(4)+conjg(o1(2))*o1(1)&nbsp; 
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(5)=o2(5)+conjg(o1(3))*o1(1)&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(6)=o2(6)+conjg(o1(3))*o1(2)&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(4)=o2(4)+o1(2)*conjg(o1(1))&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(5)=o2(5)+o1(3)*conjg(o1(1))&nbsp; 
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; o2(6)=o2(6)+o1(3)*conjg(o1(2))&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; end if&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&nbsp;in the cycle signed by 772?In another word,What is&nbsp;its meaning from Physics?<br>Anyway, Thank you very much!<br><br>&nbsp;</div>
<div><span class="gmail_quote">2007/3/15, Laurence Marks &lt;<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>&gt;:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Many parts of the Wien2k code have been written for clarity rather<br>than speed; without knowing this part exactly this is probably the
<br>case because it is not a speed-critical part of the code.<br><br>--<br>Laurence Marks<br>Department of Materials Science and Engineering<br>MSE Rm 2036 Cook Hall<br>2220 N Campus Drive<br>Northwestern University<br>Evanston, IL 60208, USA
<br>Tel: (847) 491-3996 Fax: (847) 491-7820<br>email: L-marks at northwestern dot edu<br>Web: <a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>EMM2007 <a href="http://ns.crys.ras.ru/EMMM07/">http://ns.crys.ras.ru/EMMM07/
</a><br>_______________________________________________<br>Wien mailing list<br><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br></blockquote></div><br>