I am new here. I try to understand the CASINO code related to ewald summation, but I confused about one question.
The demenstion of vecji in subroutine ewald_2d code and some called ewald_2d is different.:
SUBROUTINE
Code: Select all
SUBROUTINE ewald_2d(n,vecji,pote,field)
IMPLICIT NONE
INTEGER,INTENT(in) :: n
REAL(dp),INTENT(in) :: vecji(4,n)
....
do loop=1,n
d=vecji(1:2,loop)
z=vecji(3,loop)
.....
Code: Select all
if(ie>1)then
call ewald_2d(ie-1,vecji(1,1),pote1)
else
pote1=0.d0
How does it work?
Thanks!
Yecheng