<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>