scemama/qp_plugins_scemama
1
0
Fork 0
mirror of https://gitlab.com/scemama/qp_plugins_scemama.git synced 2024-11-08 23:23:42 +01:00
Code Issues Releases Wiki Activity

Compare commits

base: scemama:46aadac71c2d715ecf5b25de8adba0968e73765b
Branches Tags
scemama:master
..
compare: scemama:568252c5068a8f14ac08daa27d16a71cb734d9bc
Branches Tags
scemama:master

2 Commits
46aadac71c ... 568252c506

Author SHA1 Message Date
Anthony Scemama
568252c506 Merge branch 'master' of gitlab.com:scemama/qp_plugins_scemama 2021-08-24 15:14:46 +02:00
Anthony Scemama
b7fcc7f9c3 Fixed normalization in trexio 2021-08-24 15:14:29 +02:00
1 changed files with 11 additions and 2 deletions
Show Stats Expand all files Collapse all files

13
devel/trexio/export_trexio.irp.f
View File

@ -155,11 +155,20 @@ program export_trexio
rc = trexio_write_ao_shell(f, ao_shell)
call check_success(rc)
integer :: i
integer :: i, pow0(3), powA(3), j, k, l, nz
double precision :: normA, norm0, C_A(3), overlap_x, overlap_z, overlap_y, c
nz=100
C_A(1) = 0.d0
C_A(2) = 0.d0
C_A(3) = 0.d0
allocate(factor(ao_num))
print *, ao_first_of_shell
if (ao_normalized) then
do i=1,ao_num
factor(i) = ao_coef_normalization_factor(i) / shell_normalization_factor( ao_shell(i) )
l = ao_first_of_shell(ao_shell(i))
factor(i) = (ao_coef_normalized(i,1)+tiny(1.d0))/(ao_coef_normalized(l,1)+tiny(1.d0))
enddo
else
factor(:) = 1.d0