PTEROSOR/QuantumPackage
10
0
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-12-23 21:03:49 +01:00
Code Releases Activity

complex orthonormalize mos

Browse Source
This commit is contained in:
Kevin Gasperich 2020-01-28 16:44:16 -06:00
parent 92c2a3961e
commit 1dc9c3ed0b
1 changed files with 14 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/mo_one_e_ints/orthonormalize.irp.f
View File

@ -1,11 +1,20 @@
subroutine orthonormalize_mos subroutine orthonormalize_mos
implicit none implicit none
integer :: m,p,s integer :: m,p,s
m = size(mo_coef,1) if (is_periodic) then
p = size(mo_overlap,1) m = size(mo_coef_complex,1)
call ortho_lowdin(mo_overlap,p,mo_num,mo_coef,m,ao_num) p = size(mo_overlap_complex,1)
mo_label = 'Orthonormalized' call ortho_lowdin_complex(mo_overlap_complex,p,mo_num,mo_coef_complex,m,ao_num)
SOFT_TOUCH mo_coef mo_label mo_label = 'Orthonormalized'
SOFT_TOUCH mo_coef_complex mo_label
!TODO: should we do anything with the separate real/imag parts of mo_coef_complex?
else
m = size(mo_coef,1)
p = size(mo_overlap,1)
call ortho_lowdin(mo_overlap,p,mo_num,mo_coef,m,ao_num)
mo_label = 'Orthonormalized'
SOFT_TOUCH mo_coef mo_label
endif
end end