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 16:05:15 +02:00
parent 9dfa4e6531
commit f3127f4300
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/utils/complex_orthogonalization.f90
View File

@ -7,17 +7,18 @@ subroutine complex_orthogonalize_matrix(N, vectors)
! Local variables ! Local variables
integer :: i, j integer :: i, j
complex*16,allocatable :: L(:,:),Linv(:,:) complex*16,allocatable :: L(:,:),Linv(:,:),tmp(:,:)
complex*16 :: proj complex*16 :: proj
complex*16 :: norm complex*16 :: norm
! Copy the input matrix to a temporary matrix ! Copy the input matrix to a temporary matrix
allocate(L(N,N),Linv(N,N)) allocate(L(N,N),Linv(N,N),tmp(N,N))
L = matmul(transpose(vectors) ,vectors) tmp = transpose(vectors)
L = matmul(tmp ,vectors)
call complex_cholesky_decomp(N,L) call complex_cholesky_decomp(N,L)
call complex_inverse_matrix(N,L,Linv) call complex_inverse_matrix(N,L,Linv)
vectors = matmul(vectors,transpose(Linv)) vectors = matmul(vectors,transpose(Linv))
deallocate(L,Linv) deallocate(L,Linv,tmp)
end subroutine end subroutine
subroutine complex_orthonormalize(N,vectors,A) subroutine complex_orthonormalize(N,vectors,A)