<div dir="auto"><div>You are correct about what rint13.f does. It is a Simpson integrations using equal spacing (from the logarithm of the radial mesh). It is not written in the most clear code, you may want to compare it to charge*.f elsewhere which is cleaner.</div><div dir="auto"><br></div><div dir="auto">N.B., in parts of the code charge*.f has been changed to a numerically more stable integral, mainly lapw[0-2].<br><br><div data-smartmail="gmail_signature" dir="auto">---<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</a><br>Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent">www.cfw.org/100-percent</a><br>Co-Editor, Acta Cryst A<br>    </div><div class="gmail_extra" dir="auto"><br><div class="gmail_quote">On May 24, 2017 08:25, "Dara Goldar" <<a href="mailto:dara.goldar@gmail.com">dara.goldar@gmail.com</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div>
<div dir="ltr">Dear Wien2k-team
<div><br>
</div>
<div>I am currently implementing a program that will compute the overlap between various (L)APW(+lo) eigenfunctions extracted from Wien2k-simulations. As I understand, the subroutine RINT13 computes the radial integral over two-component radial functions on
 an exponential mesh, which is why I've spent some time now studying it. </div>
<div><br>
</div>
<div>There is one thing in particular about the code that I'm not able to make sense of. I've scissored out the code and attached it to this e-mail in case you have a newer version.</div>
<div><br>
</div>
<div>-------------------------- code from RINT13.F ----------------------------</div>
<div>
<div>59        GOTO 10</div>
<div>60   20 CONTINUE</div>
<div>61        P1 = RNOT*(C1*A(1)*X(1) + C2*B(1)*Y(1))</div>
<div>62        P2 = R1*(C1*A(J1)*X(J1) + C2*B(J1)*Y(J1))</div>
<div>63        S = 2*Z2 + 4*Z4 + R*(C1*A(J)*X(J) + C2*B(J)*Y(J)) + P2</div>
<div>64        S = (DX*S+P1)/3.0D+0</div>
<div>65        IF (J1 .GT. 1) S = S + 0.5D+0*DX*(P1+P2)</div>
<div>66!</div>
</div>
<div>67        RETURN</div>
<div>-------------------------- end code from RINT13.F -------------------------</div>
<div><br>
</div>
<div>On line 64, why is P1 not multiplied with DX? All other contributions to S are multiplied with DX, so I see no reason why P1 on line 64 shouldn't be. Apart from that, the code makes sense to me.</div>
<div><br>
</div>
<div>With regards</div>
<div style="font-size:12.800000190734863px">Dara Goldar, Master's student<br>
</div>
<div><span style="font-size:12.800000190734863px">NTNU/ Norwegian Defence Establishment</span>  </div>
</div>
</div>

</blockquote></div><br></div></div></div>