2
c $Id: smd_util.F 19707 2010-10-29 17:59:36Z d3y133 $
5
SUBROUTINE smd_util_rebox(n,latt,rlatt,aaa)
10
double precision rlatt(3,3),latt(3,3)
11
double precision aaa(n,3)
14
double precision ssx,ssy,ssz,xss,yss,zss
27
> write(*,*) "rebox",aaa(i,1),aaa(i,2),aaa(i,3)
28
ssx=(rlatt(1,1)*aaa(i,1)+rlatt(1,2)*aaa(i,2)+rlatt(1,3)*aaa(i,3))
29
ssy=(rlatt(2,1)*aaa(i,1)+rlatt(2,2)*aaa(i,2)+rlatt(2,3)*aaa(i,3))
30
ssz=(rlatt(3,1)*aaa(i,1)+rlatt(3,2)*aaa(i,2)+rlatt(3,3)*aaa(i,3))
36
aaa(i,1)=(latt(1,1)*xss+latt(1,2)*yss+latt(1,3)*zss)
37
aaa(i,2)=(latt(2,1)*xss+latt(2,2)*yss+latt(2,3)*zss)
38
aaa(i,3)=(latt(3,1)*xss+latt(3,2)*yss+latt(3,3)*zss)
46
subroutine smd_util_print_force_array(un,na,
52
double precision fff(na,3)
56
write(un,*) i, fff(i,1),fff(i,2),fff(i,3)
62
subroutine smd_util_print_charge_array(un,na,lo,hi,
68
double precision fff(na)