forked from alexei-matveev/paragauss-gpl
-
Notifications
You must be signed in to change notification settings - Fork 0
/
calc_radial_potG.f90
73 lines (69 loc) · 2.72 KB
/
calc_radial_potG.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
!
! ParaGauss, a program package for high-performance computations of
! molecular systems
!
! Copyright (C) 2014 T. Belling, T. Grauschopf, S. Krüger,
! F. Nörtemann, M. Staufer, M. Mayer, V. A. Nasluzov, U. Birkenheuer,
! A. Hu, A. V. Matveev, A. V. Shor, M. S. K. Fuchs-Rohr, K. M. Neyman,
! D. I. Ganyushin, T. Kerdcharoen, A. Woiterski, A. B. Gordienko,
! S. Majumder, M. H. i Rotllant, R. Ramakrishnan, G. Dixit,
! A. Nikodem, T. Soini, M. Roderus, N. Rösch
!
! This program is free software; you can redistribute it and/or modify
! it under the terms of the GNU General Public License version 2 as
! published by the Free Software Foundation [1].
!
! This program is distributed in the hope that it will be useful, but
! WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
! General Public License for more details.
!
! [1] http://www.gnu.org/licenses/gpl-2.0.html
!
! Please see the accompanying LICENSE file for further information.
!
subroutine calc_radial_potG(L_a,L_b,L_c)
! used to cal PC potential grads
use calc_3center_module, &
j1m=>j1,j2m=>j2,j3m=>j3,i1m=>i1,i2m=>i2,m_upM=>m_up
implicit none
!------------ Declaration of formal parameters ---------------
integer(kind=i4_kind), intent(in) :: L_a,L_b,L_c
!** End of interface *****************************************
integer(kind=i4_kind) :: Lam,j1,j2,j3,i1,i2,i3,m,jj,m_up,llam,i3a,i3b
gamma_arg2(:)=( (gamma_arg(:,1)-xc(1))**2 + &
(gamma_arg(:,2)-xc(2))**2 + &
(gamma_arg(:,3)-xc(3))**2 ) * fact0
gamma_help(:,1:N_max+1)=gamma(N_max+1,gamma_arg2(:))*z
jj_fac=0.0_r8_kind
do LLam=llam_m,l_a+l_b+l_c
do m_up=0,mup_llam(LLam)
do m=0,min(m_up,LLam-llam_m)
jj_fac(:,m_up,LLam)=jj_fac(:,m_up,LLam) &
+bin_fac(m_up, m)*gamma_help(:,1+LLam-m)
enddo
enddo
enddo
jj=0
do j1=0,L_a
do j2=0,L_b
do j3=0,L_c
if(.NOT. even_triangle(j1,j2,j3) ) cycle
!if( 2*((j1+j2+j3)/2) /= (j1+j2+j3) ) cycle
Lam = (j1 + j2 + j3)/2
llam=L_a+L_b+L_c-Lam
do i1=0,L_a-j1
do i2=0,L_b-j2
do i3=0,L_c-j3
jj=jj+1
m_up=L_a-i1+i2+j2-Lam
i3a=L_c-i3+j1+i1-Lam
i3b=L_b-i2+j3+i3-Lam
radialNuc_mat(:,jj)=jj_fac(:,m_up,llam)*nucR_fac(:,m_up,i3a,i3b)
end do
end do
end do
end do
end do
end do
end subroutine calc_radial_potG