PTEROSOR/QuAcK
10
1
Fork 0
mirror of https://github.com/pfloos/quack synced 2025-05-06 23:24:58 +02:00
Code Releases Activity

debugging

Browse Source
This commit is contained in:
Loris Burth 2025-04-02 15:59:20 +02:00
parent 8a3c738135
commit 9dfa4e6531
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/LR/complex_phRLR.f90
View File

@ -46,9 +46,11 @@ subroutine complex_phRLR(TDA,nS,Aph,Bph,EcRPA,Om,XpY,XmY)
RPA_matrix(1:nS,nS+1:2*nS) = Bph(:,:)
RPA_matrix(nS+1:2*nS,1:nS) = -Bph(:,:)
RPA_matrix(nS+1:2*nS,nS+1:2*nS) = -Aph(:,:)
call complex_diagonalize_matrix_without_sort(2*nS,RPA_matrix,OmOmminus)
call complex_sort_eigenvalues_RPA(2*nS,OmOmminus,RPA_matrix)
call complex_normalize_RPA(nS,RPA_matrix)
call complex_diagonalize_matrix(2*nS,RPA_matrix,OmOmminus)
call complex_orthogonalize_matrix(2*nS,RPA_matrix)
!call complex_diagonalize_matrix_without_sort(2*nS,RPA_matrix,OmOmminus)
!call complex_sort_eigenvalues_RPA(2*nS,OmOmminus,RPA_matrix)
!call complex_normalize_RPA(nS,RPA_matrix)
Om(:) = OmOmminus(1:nS)
if(maxval(abs(OmOmminus(1:nS)+OmOmminus(nS+1:2*nS))) > 1e-12) &
call print_warning('We dont find a Om and -Om structure as solution of the RPA. There might be a problem somewhere.')