LCPQ/qmcchem
10
1
Fork 1
mirror of https://gitlab.com/scemama/qmcchem.git synced 2024-12-22 04:13:31 +01:00
Code Releases Activity

Removed dead code

Browse Source
This commit is contained in:
Anthony Scemama 2016-01-14 17:07:59 +01:00
parent d878b5d33b
commit 55d3c7b68a
3 changed files with 1 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
ocaml/.gitignore vendored
View File

@ -1,6 +1,7 @@
*.cmi
*.cmx
*.o
*.cmo
ezfio.ml
Qptypes.ml
.ls_md5

2
src/SAMPLING/dmc_step.irp.f
View File

@ -135,7 +135,6 @@ for p in properties:
END_SHELL
! Brownian step
double precision :: p,q
real :: delta_x
logical :: accepted
@ -299,7 +298,6 @@ BEGIN_PROVIDER [ double precision, pop_weight, (0:dmc_projection+1) ]
! Population weight of DMC
END_DOC
pop_weight = 1.d0
pop_weight(dmc_projection) = 1.d0/dble(size(pop_weight))
END_PROVIDER
BEGIN_PROVIDER [ integer, trapped_walk, (walk_num_8) ]

49
src/SAMPLING/reconfigure.irp.f
View File

@ -1,52 +1,3 @@
subroutine reconfigure_old(ipos,w)
implicit none
integer, intent(inout) :: ipos(walk_num)
double precision, intent(in) :: w(walk_num)
integer :: kptab(walk_num), kmtab(walk_num)
double precision :: wp(walk_num), wm(walk_num)
double precision :: tmp
integer :: k, l
double precision :: qmc_ranf, rand
integer :: ipos_tmp(walk_num*2)
l=0
do k=1,walk_num
tmp = w(k)-1.d0
rand = qmc_ranf()
if (tmp >= 0.d0) then
l=l+1
ipos_tmp(l) = k
if (rand < tmp) then
l=l+1
ipos_tmp(l) = k
endif
else
if (rand > -tmp) then
l=l+1
ipos_tmp(l) = k
endif
endif
enddo
if (l>walk_num) then
do k=l,walk_num+1,-1
rand = qmc_ranf() * dble(walk_num)
ipos_tmp(int(rand)+1) = ipos_tmp(k)
enddo
else
do k=l+1,walk_num
rand = qmc_ranf() * dble(k-1)
ipos_tmp(k) = ipos_tmp(int(rand)+1)
enddo
endif
do k=1,walk_num
ipos(k) = ipos_tmp(k)
enddo
end
subroutine reconfigure(ipos,w)
implicit none
integer, intent(inout) :: ipos(*)